Opened 2 years ago

Last modified 2 years ago

#231 assigned Defect

RELEASE time variation problem

Reported by: ahilboll Owned by: ignacio
Priority: major Milestone: FLEXPART 10
Component: FP coding/compilation Version: FLEXPART 10.3
Keywords: Cc:


I'm using FLEXPART 10 from the download

When running a simulation with a RELEASE that has a duration longer than ltimesync / 2, and when I don't specify PPOINT_HOUR and PPOINT_DOW in the SPECIES file, then no particles ever get released.

I introduced some debugging statements in the code in order to find out what's the problem, and the problem seems to be that average_timecorrect in file releaseparticles.f90 has some really low value ~4E-09, which leads to rfraction being ~1E-04, which in combination with npart(i)=1000000 leads to numrel=0.

When I set the PPOINT_HOUR and PPOINT_DOW to 1.0 (multiple times, i.e., for each of the 24 hours / 7 days) in the SPECIES_NNN file, then everything works as expected.

Change History (2)

comment:1 Changed 2 years ago by ahilboll

comment:2 Changed 2 years ago by pesei

  • Component changed from FP other to FP coding/compilation
  • Milestone set to FLEXPART 10
  • Owner set to ignacio
  • Status changed from new to assigned

See also #115 - no idea whether this has been fixed or not, as the v9 code review was superseded by v10 and not completed. I think that readspecies.f90 would in general deserve some code refactoring.

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