gridcheck_ecmwf.f90: fix of small bug and improvements

[anphi]

1. check on field dimensions uses undefined variable nx

2. If ECMWF grib files contain only GRIB2 formatted fields, FLEXPART stops with an error message due to missing parameters from GRIB1 messages.

***ERROR: input file needs to contain GRiB1 formatted messages

Update 2022-02-09: In v10.4 the problem is not that a GRIB1 field is required, but that this messages is emitted if no T2 field is found (which is used to extract certain parameters). Rewriting the error message is probably good enough as a solution.

[cluekenwinkels]

Hello, I recently started using Flexpart and recieved the same ERROR in my text case. I followed the same steps as described in ​ticket 303. Is there a solution yet? Could you also give information about the flex_extact version you are using?

[pesei]

We'll try to offer at least a quick fix soon.

[pesei]

Closely related: #199 gridcheck: checking array dimensions needs to be improved

[anphi]

Dear cluekenwinkels,

without any more detail to your specific case I am unable to reproduce the error.
Just from checking the code, I have the suspicion that your EN files are missing the 2m temperature field. I am not going into much detail yet. Could you please provide me with more information?
I take it, you are using Flexpart 10.4? What flex_extract version do you use and do you still have the log file from the extraction of the EN files you are currently using? Please provide it to me (attaching it to this ticket would be good). And could you please send me the grib_ls output from at least one of your EN files please? I need to know what is actually in the files. Maybe you can attach one EN file to this ticket. That would be enough.

Thank you.
Anne

[pesei]

If your EN file is big, you can upload it here: ​https://bokubox.boku.ac.at/#f6dfa649844127b8da86fa8856618ce4 (max size 15 GB)
(NB: only I will be able to download)

[cluekenwinkels]

I followed the documentation for flex_extract 7.1.2, so this should be my version: ​https://www.flexpart.eu/flex_extract/installation.html#ref-download
The Flexpart version is 10.4.
The error only occurs using ERA5 data as I checked runs with ERA-Interim data which worked fine. The only flex_extract output I have left are the files with names EA########, all other files were deleted by flex_extract. I uploaded one of these files and the error and output files of the corresponding flex_extract run. Where can I find the EN files?
(I tried it with both settings: FORMAT=GRIB1 and FORMAT=GRIB2, and the final error remains the same.)

[pesei]

Thank you for the files. the EA-File is what we need (for many years we called the FLEXPART input files EN*, more recentely identiers specific to the data set are used, such as EA for ERA5, just old habit to say EN).

There are two uploads without file name - is this a stopped upload? They are not downloadable, so I deleted them.

[cluekenwinkels]

Ok perfect. There should be 3 files: the EA*, one .o* and one .e* file (the output and error file respectively). I did not try to upload any other files.

[anphi]

I need to tell you that the grib file looks completely wrong in the surface section. There are a couple of fields which are too much and others are completely missing. I do not remember seeing this specific combination of surface fields before and I cannot reproduce the full problem right now.

So, generally, there is a problem with flex_extract or your setup but not with FLEXPART.

Please upload the CONTROL file you used for the extraction and maybe you can also upload the run_local.sh file.

And I need to make sure that we are using the same version of flex_extract.
Did you use the tar file from the University of Vienna Phaidra website or did you clone the master branch of the git repository? And can you please confirm that you use 7.1.2 and not 7.1.1?

[cluekenwinkels]

I used this download link: git clone --single-branch --branch master ​https://www.flexpart.eu/gitmob/flex_extract
I could not find the version number within the files, but it should yield the most recent one, right?
I cannot upload data anymore to the bokubox link, so here are the first lines of my CONTROL file:

START_DATE 20180809
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.28125  
LEFT -24.
LOWER 9.875
UPPER 74.
RIGHT 60.375
LEVELIST 1/to/137
RESOL 799
ETA 1
FORMAT GRIB1
PREFIX EA
CWC 1
RRINT 1
ECTRANS 1

This should be the first lines of the run_local.sh I used (I changed it since then but the parameters should be right):

QUEUE=' '
START_DATE='20170102'
DATE_CHUNK=None
JOB_CHUNK=None
BASETIME=None
STEP=None
LEVELIST=None
AREA=None
INPUTDIR='./Workspace/path_to_output'
OUTPUTDIR=None
PP_ID=None
JOB_TEMPLATE=' '
CONTROLFILE='/work/bb1170/RUN/b381737/software/flexpart_v10.4_3d7eebf/preprocess/flex_extract/Run/Control/CONTROL_EA5_grib1'
DEBUG=0
REQUEST=0
PUBLIC=1

[anphi]

Thank you! I need some time to figure out what the problem is. I am on it. I'll let you know as soon as it's solved and how you should proceed.

And by the way, you used RRINT=1. This is not yet supported in Flexpart10.4. I currently have a beta version under test so you need to either use RRINT=0 or wait for the release of my beta version.

Best,
Anne