Opened 14 months ago

Last modified 12 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 14 months ago.
slurm-8579762.out (4.7 KB) - added by mitchekc 14 months ago.
output information
flexwrf.input.docx (17.3 KB) - added by mitchekc 14 months ago.
flexwrf.input file
slurm-8708506.out (2.3 KB) - added by mitchekc 14 months ago.
error file following compile with backtrace
namelist.input (4.2 KB) - added by mitchekc 14 months ago.
namelist.wps (1.4 KB) - added by mitchekc 14 months ago.

Download all attachments as: .zip

Change History (15)

Changed 14 months ago by mitchekc

comment:1 Changed 14 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 14 months ago by pesei

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

Changed 14 months ago by mitchekc

output information

comment:3 in reply to: ↑ 2 Changed 14 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 14 months ago by mitchekc

flexwrf.input file

comment:4 Changed 14 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 14 months ago by mitchekc

error file following compile with backtrace

comment:5 Changed 14 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 14 months ago by mitchekc

Changed 14 months ago by mitchekc

comment:6 Changed 14 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 12 months ago by pesei (previous) (diff)

comment:7 Changed 13 months ago by mitchekc

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

comment:8 follow-up: Changed 13 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 12 months ago by pesei (previous) (diff)

comment:9 in reply to: ↑ 8 Changed 12 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