Opened 2 years ago
#323 new Defect
Wrong scaling of snow depth and convective precip when reading ECMWF input in GRIB2
Reported by: | pirmink | Owned by: | |
---|---|---|---|
Priority: | major | Milestone: | |
Component: | FP input data | Version: | FLEXPART 10.4 |
Keywords: | Cc: |
Description
Flexpart expects precipitation in mm/h and snow depth in m water equivalent. This is the same for GRIB1 and GRIB2.
For GRIB2 however, when the input field is SD or CP, the conversion_factor=1000. is set in readwind_ecmwf.f90. Later in this subroutine, the values of the SD and CP fields read from GRIB2 are divided by conversion_factor, i.e. by 1000.
This seems to be wrong, as flex_extract outputs precipitation in mm/h, which is what flexpart expects according to com_mod.f90 and get_wetscav.f90. SD is in m, as expected by flexpart according to getvdep.f90. In both cases no conversion is needed in flexpart.
Note: See
TracTickets for help on using
tickets.