Opened 5 years ago

Last modified 4 years ago

#214 accepted Defect

global OUTGRID is rejected

Reported by: pesei Owned by: pesei
Priority: major Milestone: FLEXPART 10
Component: FP coding/compilation Version: FLEXPART 9.0.2
Keywords: Cc:

Description (last modified by pesei)

readoutgrid.f90 checks whether the outgrid falls within the boundaries of the meteo domain. However, a global meteo domain is a special case which is not considered here.

For example, with 0.75 deg meteo resolution (ERA-Interim), the first grid point has a longitude of -179.25 deg.
For a global OUTGRID, outlon0 has to be set to -180. This is rejected by

  if ((outlon0+eps.lt.xlon0).or.(outlat0+eps.lt.ylat0) &
       .or.(xr.gt.xr1+eps).or.(yr.gt.yr1+eps)) then

Change History (6)

comment:1 Changed 5 years ago by pesei

  • Owner changed from petra to pesei
  • Status changed from new to accepted

The most simple workaround is to comment out the whole if block of this check. But we should handle this properly!

comment:2 Changed 5 years ago by pesei

  • Description modified (diff)

comment:3 Changed 5 years ago by pesei

  • Description modified (diff)

comment:4 Changed 5 years ago by pesei

  • Description modified (diff)

comment:5 Changed 5 years ago by pesei

  • Description modified (diff)

comment:6 Changed 4 years ago by pesei

A similar issue was reported in July 2020 on the mailing list under the subjecgt OUTGRID including dateline. It should be considered when working on that ticket.

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