Opened 5 months ago

Last modified 3 months ago

#209 accepted Support

Out of bounds error

Reported by: mitchekc Owned by: pesei
Priority: major Milestone: FLEXPART_WRF_3.4_FPbase_9
Component: FP coding/compilation Version: FLEXPART-WRF
Keywords: Cc:

Description (last modified by pesei)

When I run flexpart-wrf using the command ./flexwrf_intel_omp I get out of bounds errors when my WRF geographical grid is compared to the grid in flexwrf.input. Where does FLEXPART-WRF look to find the WRF geographical grid to compare to? I have tried multiple settings in flexwrf.input and I can't determine which specific fields in the flexwrf.input file are causing the issue.

Attachments (6)

output.rtf (46.7 KB) - added by mitchekc 5 months ago.
slurm-8579762.out (4.7 KB) - added by mitchekc 5 months ago.
output information
flexwrf.input.docx (17.3 KB) - added by mitchekc 5 months ago.
flexwrf.input file
slurm-8708506.out (2.3 KB) - added by mitchekc 5 months ago.
error file following compile with backtrace
namelist.input (4.2 KB) - added by mitchekc 5 months ago.
namelist.wps (1.4 KB) - added by mitchekc 5 months ago.

Download all attachments as: .zip

Change History (15)

Changed 5 months ago by mitchekc

comment:1 Changed 5 months ago by pesei

  • Component changed from FP other to FP coding/compilation
  • Description modified (diff)
  • Milestone set to FLEXPART_WRF_3.4_FPbase_9
  • Owner set to pesei
  • Status changed from new to accepted

comment:2 follow-up: Changed 5 months ago by pesei

Can you please post the error message, including backtrace with line numbers?

Changed 5 months ago by mitchekc

output information

comment:3 in reply to: ↑ 2 Changed 5 months ago by mitchekc

Replying to pesei:

Can you please post the error message, including backtrace with line numbers?

I have attached the output file. My interpretation of the error message is that my WRF output file data and my FLEXPART-WRF geographical domain are different. The thing I can't figure out is: where is FLEXPART-WRF looking to compare with the WRF data? Am I missing a file? Should my namelist.wps or namelist.input files from WRF be placed in a certain folder?

I have also attached my flexpart.input file for reference.

Changed 5 months ago by mitchekc

flexwrf.input file

comment:4 Changed 5 months ago by pesei

What's the time interval in your meteorological input?

And please recompile FLEXPART-WRF with -fbacktrace (or corresponding option in other compiler if not gfortran) so that we can see at which line the problems are encountered (repost the output).

Changed 5 months ago by mitchekc

error file following compile with backtrace

comment:5 Changed 5 months ago by mitchekc

Hi there
I recompiled with backtrace. The issue seems to be that my FLEXPART domain is different from my WRF domain. I'm not sure how to fix the problem since I don't know the source of the issue. I have attached my namelist.wps and namelist.input files. These files are where the domain for WRF are defined.

Changed 5 months ago by mitchekc

Changed 5 months ago by mitchekc

comment:6 Changed 5 months ago by mitchekc

Specific questions I have about the FLEXPART-WRF output file include:

  1. I get this notice in my file and I don't know what it means. It seems to say that only dispersion in the air is calculated. Is that correct?
INFORMATION: SUBGRIDSCALE TERRAIN EFFECT IS NOT PARAMETERIZED DURING THIS SIMULATION.


  1. # of vertical levels in WRF data
    n_bottom_top & "true" nuvz:     34     34
    nwz &     "augmented" nuvz:     35     35

    nwzmax, nuvzmax, nzmax    :     61     61     61
    nconvlevmax, nconvlev     :     60     33
    nx, ny, nxmax, nymax      :    147    121    250    190

Where does FLEXPART-WRF get this information from? It seems like it comes from one of the WRF namelist files. Does FLEXPART-WRF need the namelist files to work correctly or does FLEXPART-WRF get this information from the netCDF input data files?

  1. readoutgrid diagnostics
xtmp1, ytmp1 (in) -7.8316902161E+01  4.5500301361E+01
xtmp2, ytmp2 (in)  2.1383683594E+04  1.4807500000E+04
 -78.5038452148     outgrid_swlon
  45.6768264771     outgrid_swlat
 -78.1782531738     outgrid_nelon
  45.7616157532     outgrid_nelat
out_xm0, out_ym0  -7.8316902161E+01  4.5500301361E+01
dxout,   dyout     1.4600000000E+02  1.2100000000E+02
  #### FLEXPART MODEL ERROR! PART OF OUTPUT    ####
  #### GRID IS OUTSIDE MODEL DOMAIN. CHANGE    ####
  #### OUTGRID IN INPUT FILE                   ####

This error message makes it clear that I need to change the OUTGRID section of the input file. I don't know how to change it. WRF defines domain location using the latitude and longitude at the centre of the grid. FLEXPART-WRF is looking for the Lower Left point and grid distance or upper right point. I am trying to define a FLEXPART-WRF grid that starts at a lower left point of (-78.3169,45.5003), goes 147 km to the East and 121 km to the North. Is this what I have defined in the flexwrf.input file?

 xmet0 (WRF), out_xm0 (FLX) :   0.0000000E+00 <  -78.31690     ?
 ymet0 (WRF), out_ym0 (FLX) :   0.0000000E+00 <   45.50030     ?
 xr1 (WRF),   xr (FLX)      :    147000.0     >   21383.68     ?
 yr1 (WRF),   yr (FLX)      :    121000.0     >   14807.50     ?

What is this trying to communicate? It seems like it is trying to say that the specified grid size is greater than a grid size that FLEXPART-WRF has generated. What input data or files does FLEXPART-WRF use to determine xr (FLX) and yr (FLX)?

Thank you for your help

Last edited 3 months ago by pesei (previous) (diff)

comment:7 Changed 5 months ago by mitchekc

The target longitude/latitude I have given above should be in eastern Ontario, Canada.

comment:8 follow-up: Changed 5 months ago by mitchekc

I think the problem may be due to the variable xmet0 and ymet0 not being initialized. I can't find where these variables are assigned a value.

Last edited 3 months ago by pesei (previous) (diff)

comment:9 in reply to: ↑ 8 Changed 3 months ago by pesei

Replying to mitchekc:

I think the problem may be due to the variable xmet0 and ymet0 not being initialized. I can't find where these variables are assigned a value.

The makefile contains an FFLAG -finit-local-zero, and the gfortran manual says about it:
The -finit-local-zero option instructs the compiler to initialize local INTEGER, REAL, and COMPLEX variables to zero, LOGICAL variables to false, and CHARACTER variables to a string of null bytes.

Did you remove this flag?

Note: See TracTickets for help on using tickets.
hosted by ZAMG