Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

(RESOLVED) Using next update GFS (2019, January): Error in ext_pkg_write_field

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

Meteo Consult

New member
Hi all,

I try to run WPS (version 3.5.1) with the next GFS (FV3-GFS) avalaible at http://para.nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/para/gfs.20181003/00/ and unfortunately I don't succeed to run metgrid.exe.

The error is :

2018-09-27 16:49:07.311 --- INFORM: GHT at level 200100.000000 already exists; leaving it alone.
2018-09-27 16:49:07.311 --- INFORM: LANDSEA at level 200100.000000 already exists; leaving it alone.
2018-09-27 16:49:07.311 --- INFORM: PSFC at level 200100.000000 already exists; leaving it alone.
2018-09-27 16:49:07.311 --- INFORM: RH at level 200100.000000 already exists; leaving it alone.
2018-09-27 16:49:07.458 --- Creating derived fields.
2018-09-27 16:49:07.458 --- INFORM: Going to create the field STC_WPS
2018-09-27 16:49:07.464 --- INFORM: STC_WPS at level 1.000000 already exists; leaving it alone.
2018-09-27 16:49:07.465 --- INFORM: STC_WPS at level 2.000000 already exists; leaving it alone.
2018-09-27 16:49:07.465 --- INFORM: STC_WPS at level 1.000000 already exists; leaving it alone.
2018-09-27 16:49:07.465 --- INFORM: STC_WPS at level 2.000000 already exists; leaving it alone.
2018-09-27 16:49:07.465 --- INFORM: STC_WPS at level 3.000000 already exists; leaving it alone.
2018-09-27 16:49:07.465 --- INFORM: STC_WPS at level 4.000000 already exists; leaving it alone.
2018-09-27 16:49:07.465 --- INFORM: Going to create the field SMC_WPS
2018-09-27 16:49:07.471 --- INFORM: SMC_WPS at level 1.000000 already exists; leaving it alone.
2018-09-27 16:49:07.471 --- INFORM: SMC_WPS at level 2.000000 already exists; leaving it alone.
2018-09-27 16:49:07.471 --- INFORM: SMC_WPS at level 1.000000 already exists; leaving it alone.
2018-09-27 16:49:07.472 --- INFORM: SMC_WPS at level 2.000000 already exists; leaving it alone.
2018-09-27 16:49:07.472 --- INFORM: SMC_WPS at level 3.000000 already exists; leaving it alone.
2018-09-27 16:49:07.472 --- INFORM: SMC_WPS at level 4.000000 already exists; leaving it alone.
2018-09-27 16:49:07.472 --- INFORM: Going to create the field SW
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SW000010 at level 200100.000000 to fill level 1.000000 of SW.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SW010040 at level 200100.000000 to fill level 2.000000 of SW.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SW040100 at level 200100.000000 to fill level 3.000000 of SW.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SW100200 at level 200100.000000 to fill level 4.000000 of SW.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SW000010 at level 200100.000000 to fill level 1.000000 of SW.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SW010200 at level 200100.000000 to fill level 2.000000 of SW.
2018-09-27 16:49:07.472 --- INFORM: Going to create the field SOILM
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILM000 at level 200100.000000 to fill level 1.000000 of SOILM.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILM005 at level 200100.000000 to fill level 2.000000 of SOILM.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILM020 at level 200100.000000 to fill level 3.000000 of SOILM.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILM040 at level 200100.000000 to fill level 4.000000 of SOILM.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILM160 at level 200100.000000 to fill level 5.000000 of SOILM.
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILM300 at level 200100.000000 to fill level 6.000000 of SOILM.
2018-09-27 16:49:07.472 --- INFORM: Going to create the field SOILT
2018-09-27 16:49:07.472 --- INFORM: Couldn't find SOILT000 at level 200100.000000 to fill level 1.000000 of SOILT.
2018-09-27 16:49:07.473 --- INFORM: Couldn't find SOILT005 at level 200100.000000 to fill level 2.000000 of SOILT.
2018-09-27 16:49:07.473 --- INFORM: Couldn't find SOILT020 at level 200100.000000 to fill level 3.000000 of SOILT.
2018-09-27 16:49:07.473 --- INFORM: Couldn't find SOILT040 at level 200100.000000 to fill level 4.000000 of SOILT.
2018-09-27 16:49:07.473 --- INFORM: Couldn't find SOILT160 at level 200100.000000 to fill level 5.000000 of SOILT.
2018-09-27 16:49:07.473 --- INFORM: Couldn't find SOILT300 at level 200100.000000 to fill level 6.000000 of SOILT.
2018-09-27 16:49:07.473 --- INFORM: Going to create the field PRES
2018-09-27 16:49:07.505 --- INFORM: PRES at level 200100.000000 already exists; leaving it alone.
2018-09-27 16:49:07.696 --- INFORM: Field LANDSEA.mask does not have a valid mask and will not be checked for missing values
2018-09-27 16:49:07.778 --- Initializing output module.
2018-09-27 16:49:07.781 --- ERROR: Error in ext_pkg_write_field


This is probably due to the addition of new fields in the gribs files.

I read that WPS 4.0 is supporting FV3-GFS but I would prefer not to have to complete all my WPS.

How should I change my Vtable and my Metgrid.TBL for the WPS to work (see attached files : View attachment Vtable.txt ; View attachment METGRID.TBL) ?

Tkanks for your help.

Regards,

Sebastien Poitevin
 
@Meteo Consult,
If you are interested in the code changes that allow V4.0 to work with FV3-GFS files, you can take a look at this code submission on GitHub:
https://github.com/wrf-model/WPS/commit/5f9c79a6b4f659cccace1bb960259b83ed7eae9f

The description is at the top, with code line modifications listed beneath. You could compare the modified files between V4.0 and V3.5.1 to see if it's possible to simply incorporate those changes into the 3.5.1 code. If not, another possibility may be to use 3.5.1 geo_em* output to run with V4.0 ungrib/metgrid. This way you will at least not have to rerun all of WPS.
 
Here are more details about processing FV3-GFS data :

Please add pmin to &ungrib to specify the minimum pressure level to process. Usually the value of pmin should be equal to or smaller than p_top_requested in the namelist.input. However, please ignore the option ec_rec_len = 26214508 , which is not applied to GFS data.

Sorry that somehow I didn't notice that you are trying to process FV3_GFS data. NCEP updated GFS output a few times in the past years, which makes the older version of WPS not work anymore. We suggest people to use the newest version of WPS.

Ming Chen
 
Thank you both for your help.

As explained above I want to avoid having to recompile the code because the model runs in an operational way.

On the other hand, by removing the fields at 0.4 mb FV3-GFS (by doing wgrib2 -not ": 0.4 mb" -grib_out ...), I managed to run my metgrid.exe without any problem.

I am testing my WRF 3.5.1 with. I had to change my namelist.input by replacing num_metgrid_levels = 32 by 34 but otherwise everything seems to work fine.

I just have to compare my new WRF forecasts from FV3-GFS with those from the operational GFS to make sure there is no anomaly.
 
Top