Opened 5 years ago
Last modified 5 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: |
Description
I'm using FLEXPART 10 from the download https://www.flexpart.eu/downloads/63
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 5 years ago by ahilboll
comment:2 Changed 5 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.
Just checked, and it seems this got addressed by https://www.flexpart.eu/changeset/0ca497657ed0c132e49e911adc17490636fdc780/flexpart.git ...