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

geogrid error (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.

jokervtv

New member
Dear Sir or Madam,

After I successfully generated all the execution files (WPS and WRF), I received the error message "ERROR: Error in ext_pkg_write_field" after running geogrid.exe. Attach the log at the end of the article.

If you can help, thank you very much.

--------
Code:
2020-06-04 22:17:12.291 ---  *** Starting program geogrid.exe ***
2020-06-04 22:17:12.291 --- Using the following namelist variables:
2020-06-04 22:17:12.291 --- &SHARE
2020-06-04 22:17:12.291 ---   WRF_CORE         = ARW
2020-06-04 22:17:12.291 ---   MAX_DOM          = 1
2020-06-04 22:17:12.291 ---   START_YEAR       = 0
2020-06-04 22:17:12.291 ---   START_MONTH      = 0
2020-06-04 22:17:12.291 ---   START_DAY        = 0
2020-06-04 22:17:12.291 ---   START_HOUR       = 0
2020-06-04 22:17:12.291 ---   START_MINUTE     = 0
2020-06-04 22:17:12.291 ---   START_SECOND     = 0
2020-06-04 22:17:12.291 ---   END_YEAR         = 0
2020-06-04 22:17:12.291 ---   END_MONTH        = 0
2020-06-04 22:17:12.291 ---   END_DAY          = 0
2020-06-04 22:17:12.291 ---   END_HOUR         = 0
2020-06-04 22:17:12.291 ---   END_MINUTE       = 0
2020-06-04 22:17:12.292 ---   END_SECOND       = 0
2020-06-04 22:17:12.292 ---   START_DATE       = 2006-08-16_12:00:00
2020-06-04 22:17:12.292 ---   END_DATE         = 2006-08-16_18:00:00
2020-06-04 22:17:12.292 ---   INTERVAL_SECONDS = 21600
2020-06-04 22:17:12.292 ---   IO_FORM_GEOGRID  = 2
2020-06-04 22:17:12.292 ---   OPT_OUTPUT_FROM_GEOGRID_PATH = ./
2020-06-04 22:17:12.292 ---   SUBGRID_RATIO_X  = 1
2020-06-04 22:17:12.292 ---   SUBGRID_RATIO_Y  = 1
2020-06-04 22:17:12.292 ---   DEBUG_LEVEL      = 0
2020-06-04 22:17:12.292 ---   ACTIVE_GRID      = .TRUE.
2020-06-04 22:17:12.292 ---   NOCOLONS         = .FALSE.
2020-06-04 22:17:12.292 --- /
2020-06-04 22:17:12.292 --- &GEOGRID
2020-06-04 22:17:12.292 ---   PARENT_ID         = 1
2020-06-04 22:17:12.292 ---   PARENT_GRID_RATIO = 1
2020-06-04 22:17:12.292 ---   I_PARENT_START    = 1
2020-06-04 22:17:12.292 ---   J_PARENT_START    = 1
2020-06-04 22:17:12.292 ---   S_WE              = 1
2020-06-04 22:17:12.292 ---   E_WE              = 74
2020-06-04 22:17:12.292 ---   S_SN              = 1
2020-06-04 22:17:12.292 ---   E_SN              = 61
2020-06-04 22:17:12.292 ---   GEOG_DATA_RES     = default
2020-06-04 22:17:12.292 ---   DX                = 30000.000000
2020-06-04 22:17:12.292 ---   DY                = 30000.000000
2020-06-04 22:17:12.292 ---   MAP_PROJ          = lambert
2020-06-04 22:17:12.292 ---   POLE_LAT          = 90.000000
2020-06-04 22:17:12.292 ---   POLE_LON          = 0.000000
2020-06-04 22:17:12.292 ---   REF_LAT           = 34.830002
2020-06-04 22:17:12.292 ---   REF_LON           = 81.029999
2020-06-04 22:17:12.292 ---   REF_X             = 100000002004087734272.000000
2020-06-04 22:17:12.292 ---   REF_Y             = 100000002004087734272.000000
2020-06-04 22:17:12.292 ---   TRUELAT1          = 30.000000
2020-06-04 22:17:12.292 ---   TRUELAT2          = 60.000000
2020-06-04 22:17:12.292 ---   STAND_LON         = 45.000000
2020-06-04 22:17:12.292 ---   GEOG_DATA_PATH    = /home/jschao/WPS_GEOG/
2020-06-04 22:17:12.292 ---   OPT_GEOGRID_TBL_PATH = geogrid/
2020-06-04 22:17:12.293 --- /
2020-06-04 22:17:12.299 --- Processing domain 1 of 1
2020-06-04 22:17:12.299 --- INFORM: Using default interpolator sequence for HGT_M.
2020-06-04 22:17:12.299 --- INFORM: Using default data source for HGT_M.
2020-06-04 22:17:12.300 --- INFORM: Using default interpolator sequence for LANDUSEF.
2020-06-04 22:17:12.300 --- INFORM: Using default data source for LANDUSEF.
2020-06-04 22:17:12.300 --- INFORM: Using default interpolator sequence for SOILTEMP.
2020-06-04 22:17:12.300 --- INFORM: Using default data source for SOILTEMP.
2020-06-04 22:17:12.300 --- INFORM: Using default interpolator sequence for SOILCTOP.
2020-06-04 22:17:12.300 --- INFORM: Using default data source for SOILCTOP.
2020-06-04 22:17:12.300 --- INFORM: Using default interpolator sequence for SOILCBOT.
2020-06-04 22:17:12.300 --- INFORM: Using default data source for SOILCBOT.
2020-06-04 22:17:12.301 --- INFORM: Using default interpolator sequence for ALBEDO12M.
2020-06-04 22:17:12.301 --- INFORM: Using default data source for ALBEDO12M.
2020-06-04 22:17:12.301 --- INFORM: Using default interpolator sequence for GREENFRAC.
2020-06-04 22:17:12.301 --- INFORM: Using default data source for GREENFRAC.
2020-06-04 22:17:12.301 --- INFORM: Using default interpolator sequence for LAI12M.
2020-06-04 22:17:12.301 --- INFORM: Using default data source for LAI12M.
2020-06-04 22:17:12.301 --- INFORM: Using default interpolator sequence for SNOALB.
2020-06-04 22:17:12.301 --- INFORM: Using default data source for SNOALB.
2020-06-04 22:17:12.302 --- INFORM: Using default interpolator sequence for CON.
2020-06-04 22:17:12.302 --- INFORM: Using default data source for CON.
2020-06-04 22:17:12.302 --- INFORM: Using default interpolator sequence for VAR.
2020-06-04 22:17:12.302 --- INFORM: Using default data source for VAR.
2020-06-04 22:17:12.302 --- INFORM: Using default interpolator sequence for OA1.
2020-06-04 22:17:12.302 --- INFORM: Using default data source for OA1.
2020-06-04 22:17:12.302 --- INFORM: Using default interpolator sequence for OA2.
2020-06-04 22:17:12.302 --- INFORM: Using default data source for OA2.
2020-06-04 22:17:12.303 --- INFORM: Using default interpolator sequence for OA3.
2020-06-04 22:17:12.303 --- INFORM: Using default data source for OA3.
2020-06-04 22:17:12.303 --- INFORM: Using default interpolator sequence for OA4.
2020-06-04 22:17:12.303 --- INFORM: Using default data source for OA4.
2020-06-04 22:17:12.303 --- INFORM: Using default interpolator sequence for OL1.
2020-06-04 22:17:12.303 --- INFORM: Using default data source for OL1.
2020-06-04 22:17:12.303 --- INFORM: Using default interpolator sequence for OL2.
2020-06-04 22:17:12.303 --- INFORM: Using default data source for OL2.
2020-06-04 22:17:12.304 --- INFORM: Using default interpolator sequence for OL3.
2020-06-04 22:17:12.304 --- INFORM: Using default data source for OL3.
2020-06-04 22:17:12.304 --- INFORM: Using default interpolator sequence for OL4.
2020-06-04 22:17:12.304 --- INFORM: Using default data source for OL4.
2020-06-04 22:17:12.304 --- INFORM: Using default interpolator sequence for VAR_SSO.
2020-06-04 22:17:12.304 --- INFORM: Using default data source for VAR_SSO.
2020-06-04 22:17:12.304 --- INFORM: Using default interpolator sequence for LAKE_DEPTH.
2020-06-04 22:17:12.304 --- INFORM: Using default data source for LAKE_DEPTH.
2020-06-04 22:17:12.305 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/lake_depth/index for field LAKE_DEPTH
2020-06-04 22:17:12.305 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.305 --- INFORM: Using default interpolator sequence for URB_PARAM.
2020-06-04 22:17:12.305 --- INFORM: Using default data source for URB_PARAM.
2020-06-04 22:17:12.305 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/NUDAPT44_1km/index for field URB_PARAM
2020-06-04 22:17:12.305 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.305 --- INFORM: Using default interpolator sequence for FRC_URB2D.
2020-06-04 22:17:12.305 --- INFORM: Using default data source for FRC_URB2D.
2020-06-04 22:17:12.305 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/urbfrac_nlcd2011/index for field FRC_URB2D
2020-06-04 22:17:12.305 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.305 --- INFORM: Using default interpolator sequence for IMPERV.
2020-06-04 22:17:12.305 --- INFORM: Using default data source for IMPERV.
2020-06-04 22:17:12.305 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/nlcd2011_imp_ll_9s/index for field IMPERV
2020-06-04 22:17:12.305 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.305 --- INFORM: Using default interpolator sequence for CANFRA.
2020-06-04 22:17:12.305 --- INFORM: Using default data source for CANFRA.
2020-06-04 22:17:12.305 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/nlcd2011_can_ll_9s/index for field CANFRA
2020-06-04 22:17:12.305 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.305 --- INFORM: Using default interpolator sequence for EROD.
2020-06-04 22:17:12.306 --- INFORM: Using default data source for EROD.
2020-06-04 22:17:12.306 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/erod/index for field EROD
2020-06-04 22:17:12.306 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.306 --- INFORM: Using default interpolator sequence for CLAYFRAC.
2020-06-04 22:17:12.306 --- INFORM: Using default data source for CLAYFRAC.
2020-06-04 22:17:12.306 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/clayfrac_5m/index for field CLAYFRAC
2020-06-04 22:17:12.306 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.306 --- INFORM: Using default interpolator sequence for SANDFRAC.
2020-06-04 22:17:12.306 --- INFORM: Using default data source for SANDFRAC.
2020-06-04 22:17:12.306 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/sandfrac_5m/index for field SANDFRAC
2020-06-04 22:17:12.306 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.306 --- INFORM: Using default interpolator sequence for IRRIGATION.
2020-06-04 22:17:12.306 --- INFORM: Using default data source for IRRIGATION.
2020-06-04 22:17:12.306 --- INFORM: Could not read 'index' file /home/jschao/WPS_GEOG/irrigation/fao/index for field IRRIGATION
2020-06-04 22:17:12.306 --- INFORM: This field is optional and will not be processed.
2020-06-04 22:17:12.325 --- ERROR: Error in ext_pkg_write_field
------------
 
kwerner said:
Hi,
Can you attach your namelist.wps file so that I can take a look?

Thanks for your reply. And the namelist.wps is here:

Code:
&share
 wrf_core = 'ARW',
 max_dom = 1,
 start_date = '2006-08-16_12:00:00','2006-08-16_12:00:00',
 end_date   = '2006-08-16_18:00:00','2006-08-16_12:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,   1,
 parent_grid_ratio =   1,   3,
 i_parent_start    =   1,  31,
 j_parent_start    =   1,  17,
 e_we              =  74, 112,
 e_sn              =  61,  97,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! The default datasets used to produce the MAXSNOALB and ALBEDO12M
 ! fields have changed in WPS v4.0. These fields are now interpolated
 ! from MODIS-based datasets.
 !
 ! To match the output given by the default namelist.wps in WPS v3.9.1,
 ! the following setting for geog_data_res may be used:
 !
 ! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default', 
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = 'default','default',
 dx = 30000,
 dy = 30000,
 map_proj = 'lambert',
 ref_lat   =  34.83,
 ref_lon   = 81.03,
 truelat1  =  30.0,
 truelat2  =  60.0,
 stand_lon = 45.0,
 geog_data_path = '/home/jschao/WPS_GEOG/'
/

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE'
 io_form_metgrid = 2, 
/
 
@jokervtv Which compiler and compiler version are you using? Which versions of the netCDF C and Fortran libraries do you have? I'll see if I can reproduce the issue locally.
 
Thank you for your reminder.

I checked the version of HDF and found that geogrid.exe links two versions of HDF. This is my question. One is compiled by myself and the other is provided by conda. Finally, I installed WRF using docker and it worked well.

Thanks again.
 
Top