Error with precipitation fields while running FLEXPART
Hello, I have been attempting to run FLEXPART, but I am receiving the following error message after I run the program:
Welcome to FLEXPART Version 11.0 (2023-07-11)
Git: b6dd2bd Tue Jan 21 15:30:28 2025 +0100
FLEXPART is free software released under the GNU General Public License.
FLEXPART is running with ETA coordinates.
WARNING: SURF_ONLY in COMMAND will be deprecated and renamed SFC_ONLY
Continuing with SURF_ONLY...
WARNING: FILE COMMAND LRECOUTSTEP not provided,
value of LOUTSTEP will be used if RECEPTORS are
required.
WARNING: FILE COMMAND LRECOUTAVER not provided,
value of LOUTAVER will be used if RECEPTORS are
required.
WARNING: FILE COMMAND LRECOUTSTEP not provided,
value of LOUTSAMPLE will be used if RECEPTORS are
required.
----------------
INFORMATION: SUBGRIDSCALE TERRAIN EFFECT IS
NOT PARAMETERIZED DURING THIS SIMULATION.
----------------
Switch for LCM output LCMOUTPUT = F
*********** WARNING **********************************
* FLEXPART running in parallel mode *
* Number of uncertainty classes in *
* set to number of threads: 12 *
* All other computations are done with 12 threads. *
*******************************************************
ECMWF metdata detected
NXSHIFT is set to 359
grid dim: 721 361 138 138 137 138
ERROR STOP Conditions for precipitation interpolation not fulfilled! Please check number of precipitation fields per type in GRIB files and numpf parameter in par_mod.f90.
Error termination. Backtrace:
#0 0x14df80338960 in ???
#1 0x14df803394d9 in ???
#2 0x14df8033ab47 in ???
#3 0x5613e5747074 in __windfields_mod_MOD_gridcheck_ecmwf
at /home/kbloxam/flexpart/src/windfields_mod.f90:599
#4 0x5613e5756309 in read_options_and_initialise_flexpart_
at /home/kbloxam/flexpart/src/FLEXPART.f90:272
#5 0x5613e5756d01 in flexpart
at /home/kbloxam/flexpart/src/FLEXPART.f90:95
#6 0x5613e5609985 in main
at /home/kbloxam/flexpart/src/FLEXPART.f90:36This seems like it should be related to the precipitation field. I have tried going to the par_mod.f90 file and changed the numpf field from 1 to 3 and I still get the same message. I have even tried changing it to 4 to see if I get a different error message but I still get the same.
The data I am using was downloaded from ECMWF using the following CONTROL file input:
START_DATE 19950101
DTIME 1
TYPE AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN AN
TIME 00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23
STEP 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00
ACCTYPE FC
ACCTIME 06/18
ACCMAXSTEP 12
CLASS EA
STREAM OPER
GRID 0.5
LEFT -179.5
LOWER -90.
UPPER 90.
RIGHT 180.
LEVELIST 1/to/137
RESOL 359
ETA 1
CWC 1
PREFIX EAg
ECTRANS 1I should also point out that in building flex_extract I used ecCodes version 2.35. When I activate the environment I used to build flex_extract, start Python, and import eccodes, I get the message:
/home/kbloxam/flex_extract/pyflexextract/lib/python3.10/site-packages/gribapi/__init__.py:23: UserWarning: ecCodes 2.38.0 or higher is recommended. You are running version 2.35.0
However, outside of Python, if I type pip show eccodes, it says it is version 2.39.2, even when the environment has been activated. Perhaps there is a disconnect somewhere because when I have looked at the data that was downloaded using flex_extract there are some fields that aren't identified. For example,
"137[-] HYBL="Hybrid level"" "Specific humidity [kg/kg]" "[kg/kg]" "SPFH" "137-HYBL" '01-Jan-1995' '01-Jan-1995' '0 sec'
"1[-] HYBL="Hybrid level"" "(prodType 192, cat 201, subcat 31) [-]" "[-]" "unknown" "1-HYBL" '01-Jan-1995' '01-Jan-1995' '0 sec'
"2[-] HYBL="Hybrid level"" "(prodType 192, cat 201, subcat 31) [-]" "[-]" "unknown" "2-HYBL" '01-Jan-1995' '01-Jan-1995' '0 sec'
...and I am wondering if the ones that don't specify the name are related to the precipitation fields.
If you have any idea what the problem might be, I would greatly appreciate the advice.
Thank you,
Kevin Bloxam (He/Him) Postdoctoral Fellow Department of Physics | Faculty of Arts and Science | University of Toronto kevin.bloxam@utoronto.ca