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) Troubles trying to change from MODIS to USGS data

A

Arty

Member
Hello,

Although I read the User's Guide and did as explained, I can't make WPS to use USGS data instead of MODIS. I need to try and use USGS because of ndown/WRF.exe issue (using old coarse run under USGS as input via ndown.WRF.exe for a finer simulation, currently "segfaulting", maybe because of USGS vs. MODIS ?). See this thread for more info.

I tried different attributes for "geog_res_data" :

'usgs_lakes+default' --> going MODIS everytime, "don't find ..."
'usgs_lakes' --> always going to default, and finally fails when requiring first missing dataset
'usgs_lakes+30s' --> no better, cf. geogrid.log hereunder

There's two "source" directories for geog data (which I'm only allowed to read) on the computer I'm working on :

- /geog, which lists these directories :
albedo_ncep
clayfrac_5m
erod
greenfrac
greenfrac_fpar_modis
hangl
hanis
hasynw
hasys
hasysw
hasyw
hcnvx
hlennw
hlens
hlensw
hlenw
hslop
hstdv
hzmax
islope
lai_modis_10m
lai_modis_30s
lake_depth
landuse_10m
landuse_2m
landuse_30s
landuse_30s_with_lakes
landuse_5m
maxsnowalb
modis_landuse_20class_15s
modis_landuse_20class_30s
modis_landuse_21class_30s
nlcd2006_ll_30s
NUDAPT44_1km
orogwd_10m
orogwd_1deg
orogwd_20m
orogwd_2deg
orogwd_30m
README
sandfrac_5m
soiltemp_1deg
soiltype_bot_10m
soiltype_bot_2m
soiltype_bot_30s
soiltype_bot_5m
soiltype_top_10m
soiltype_top_2m
soiltype_top_30s
soiltype_top_5m
ssib_landuse_10m
ssib_landuse_5m
topo_10m
topo_2m
topo_30s
topo_5m
topo_gmted2010_30s
varsso
varsso_10m
varsso_2m
varsso_5m

- /geog4/WPS, which lists these directories :
albedo_modis
greenfrac_fpar_modis
lai_modis_10m
lai_modis_30s
maxsnowalb_modis
modis_landuse_20class_30s_with_lakes
orogwd_10m
orogwd_1deg
orogwd_20m
orogwd_2deg
orogwd_30m
soiltemp_1deg
soiltype_bot_30s
soiltype_top_30s
topo_gmted2010_30s
varsso
varsso_10m
varsso_2m
varsso_5m

Finally the geogrid.exe fails, cf. geogrid.log hereunder :
Parsed 29 entries in GEOGRID.TBL
Processing domain 1 of 1
ERROR: Could not open /home3/datawork/pmauger/DATA/geog/albedo_modis/index
application called MPI_Abort(MPI_COMM_WORLD, 33279808) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=33279808
:
system msg for write_line failure : Bad file descriptor
2023-02-25 21:03:43.531 --- START_MONTH = 0
2023-02-25 21:03:43.531 --- START_DAY = 0
2023-02-25 21:03:43.531 --- START_HOUR = 0
2023-02-25 21:03:43.531 --- START_MINUTE = 0
2023-02-25 21:03:43.531 --- START_SECOND = 0
2023-02-25 21:03:43.531 --- END_YEAR = 0
2023-02-25 21:03:43.531 --- END_MONTH = 0
2023-02-25 21:03:43.531 --- END_DAY = 0
2023-02-25 21:03:43.531 --- END_HOUR = 0
2023-02-25 21:03:43.531 --- END_MINUTE = 0
2023-02-25 21:03:43.532 --- END_SECOND = 0
2023-02-25 21:03:43.532 --- START_DATE = 2005-02-02_00:00:00
2023-02-25 21:03:43.532 --- END_DATE = 2005-02-07_00:00:00
2023-02-25 21:03:43.532 --- INTERVAL_SECONDS = 21600
2023-02-25 21:03:43.532 --- IO_FORM_GEOGRID = 2
2023-02-25 21:03:43.532 --- OPT_OUTPUT_FROM_GEOGRID_PATH = /home3/datawork/pmauger/CONFIGS/THT5/WRF_FILES/WPS_DATA
2023-02-25 21:03:43.532 --- SUBGRID_RATIO_X = 1
2023-02-25 21:03:43.532 --- SUBGRID_RATIO_Y = 1
2023-02-25 21:03:43.532 --- DEBUG_LEVEL = 0
2023-02-25 21:03:43.532 --- ACTIVE_GRID = .TRUE.
2023-02-25 21:03:43.532 --- NOCOLONS = .FALSE.
2023-02-25 21:03:43.532 --- /
2023-02-25 21:03:43.532 --- &GEOGRID
2023-02-25 21:03:43.532 --- PARENT_ID = 1
2023-02-25 21:03:43.532 --- PARENT_GRID_RATIO = 1
2023-02-25 21:03:43.532 --- I_PARENT_START = 1
2023-02-25 21:03:43.532 --- J_PARENT_START = 1
2023-02-25 21:03:43.532 --- S_WE = 1
2023-02-25 21:03:43.532 --- E_WE = 51
2023-02-25 21:03:43.532 --- S_SN = 1
2023-02-25 21:03:43.532 --- E_SN = 51
2023-02-25 21:03:43.532 --- GEOG_DATA_RES = usgs_lakes+30s
2023-02-25 21:03:43.532 --- DX = 7500.000000
2023-02-25 21:03:43.532 --- DY = 7500.000000
2023-02-25 21:03:43.532 --- MAP_PROJ = mercator
2023-02-25 21:03:43.532 --- POLE_LAT = 90.000000
2023-02-25 21:03:43.532 --- POLE_LON = 0.000000
2023-02-25 21:03:43.532 --- REF_LAT = -17.500000
2023-02-25 21:03:43.532 --- REF_LON = -149.500000
2023-02-25 21:03:43.532 --- REF_X = 100000002004087734272.000000
2023-02-25 21:03:43.532 --- REF_Y = 100000002004087734272.000000
2023-02-25 21:03:43.532 --- TRUELAT1 = -17.500000
2023-02-25 21:03:43.532 --- TRUELAT2 = 0.000000
2023-02-25 21:03:43.532 --- STAND_LON = -149.500000
2023-02-25 21:03:43.532 --- GEOG_DATA_PATH = /home3/datawork/pmauger/DATA/geog
2023-02-25 21:03:43.532 --- OPT_GEOGRID_TBL_PATH = /home3/datahome/pmauger/WRF/../WPS/geogrid
2023-02-25 21:03:43.532 --- /
2023-02-25 21:03:43.535 --- Processing domain 1 of 1
2023-02-25 21:03:43.535 --- INFORM: For HGT_M, couldn't find interpolator sequence for resolution usgs_lakes.
2023-02-25 21:03:43.535 --- INFORM: For HGT_M, couldn't find interpolator sequence for resolution 30s.
2023-02-25 21:03:43.535 --- INFORM: Using default interpolator sequence for HGT_M.
2023-02-25 21:03:43.535 --- INFORM: For HGT_M, couldn't find usgs_lakes data source.
2023-02-25 21:03:43.536 --- INFORM: For HGT_M, couldn't find 30s data source.
2023-02-25 21:03:43.536 --- INFORM: Using default data source for HGT_M.
2023-02-25 21:03:43.536 --- INFORM: Using usgs_lakes interpolator sequence for LANDUSEF.
2023-02-25 21:03:43.536 --- INFORM: Using usgs_lakes data source for LANDUSEF.
2023-02-25 21:03:43.536 --- INFORM: For SOILTEMP, couldn't find interpolator sequence for resolution usgs_lakes.
2023-02-25 21:03:43.536 --- INFORM: For SOILTEMP, couldn't find interpolator sequence for resolution 30s.
2023-02-25 21:03:43.537 --- INFORM: Using default interpolator sequence for SOILTEMP.
2023-02-25 21:03:43.537 --- INFORM: For SOILTEMP, couldn't find usgs_lakes data source.
2023-02-25 21:03:43.537 --- INFORM: For SOILTEMP, couldn't find 30s data source.
2023-02-25 21:03:43.538 --- INFORM: Using default data source for SOILTEMP.
2023-02-25 21:03:43.538 --- INFORM: For SOILCTOP, couldn't find interpolator sequence for resolution usgs_lakes.
2023-02-25 21:03:43.538 --- INFORM: Using 30s interpolator sequence for SOILCTOP.
2023-02-25 21:03:43.538 --- INFORM: For SOILCTOP, couldn't find usgs_lakes data source.
2023-02-25 21:03:43.538 --- INFORM: Using 30s data source for SOILCTOP.
2023-02-25 21:03:43.538 --- INFORM: For SOILCBOT, couldn't find interpolator sequence for resolution usgs_lakes.
2023-02-25 21:03:43.538 --- INFORM: Using 30s interpolator sequence for SOILCBOT.
2023-02-25 21:03:43.538 --- INFORM: For SOILCBOT, couldn't find usgs_lakes data source.
2023-02-25 21:03:43.538 --- INFORM: Using 30s data source for SOILCBOT.
2023-02-25 21:03:43.538 --- INFORM: For ALBEDO12M, couldn't find interpolator sequence for resolution usgs_lakes.
2023-02-25 21:03:43.538 --- INFORM: For ALBEDO12M, couldn't find interpolator sequence for resolution 30s.
2023-02-25 21:03:43.538 --- INFORM: Using default interpolator sequence for ALBEDO12M.
2023-02-25 21:03:43.538 --- INFORM: For ALBEDO12M, couldn't find usgs_lakes data source.
2023-02-25 21:03:43.538 --- INFORM: For ALBEDO12M, couldn't find 30s data source.
2023-02-25 21:03:43.538 --- INFORM: Using default data source for ALBEDO12M.
2023-02-25 21:03:43.538 --- ERROR: Could not open /home3/datawork/pmauger/DATA/geog/albedo_modis/index

I understand geogrid.exe can't find " albedo_modis " as I changed from /geog/ directory to /geog4/WPS/ directory to access USGS data. Nevertheless, it seems that from the beginning geogrid.exe won't find what it needs to use USGS data.

Should I copy all /geog/* & /geog4/WPS/* directories to another unique directory and update my path ? Do I miss something in the existing directories ?
I read other threads (as this recent one) but couldn't find a workaround.

Thanks for your help.
 
Last edited:
Hi,
To use USGS data, the correct setting in the namelist.wps file is
Code: 
geog_data_res = 'usgs_lakes+default','usgs_lakes+default',

and to use that data, plus the default, you need to make sure you have the static fields listed in the first table (with the blue outline) on this page, in addition to the data for "landuse_30s_with_lakes." all in the same directory. However, an alternative is to modify your WPS/geogrid/GEOGRID.TBL file so that the setting for usgs_lakes points to an absolute path. For e.g., in the GEOGRID.TBL, change this line

Code: 
rel_path =     usgs_lakes:landuse_30s_with_lakes/

to

Code: 
abs_path =     /path-to-directory-containing-landuse_30s_with_lakes-data/usgs_lakes:landuse_30s_with_lakes/

and then for the "geog_data_path" setting, you can keep it set to the directory where you have all the other fields.
 
Hello,

I succeded putting all sub-directories (which were separated in two directories) in one directory ; this was indeed a path problem. Nevertheless, I just tried to run ndown with the new WPS/REAL ouputs and I still get PSFC of 10^22 magnitude ; so we can exclude the MODIS/USGS landuse base as root cause for this problem.

Thanks for your help.
 
kwerner said:
Hi,
To use USGS data, the correct setting in the namelist.wps file is
Code: 
geog_data_res = 'usgs_lakes+default','usgs_lakes+default',

and to use that data, plus the default, you need to make sure you have the static fields listed in the first table (with the blue outline) on this page, in addition to the data for "landuse_30s_with_lakes." all in the same directory. However, an alternative is to modify your WPS/geogrid/GEOGRID.TBL file so that the setting for usgs_lakes points to an absolute path. For e.g., in the GEOGRID.TBL, change this line

Code: 
rel_path =     usgs_lakes:landuse_30s_with_lakes/

to

Code: 
abs_path =     /path-to-directory-containing-landuse_30s_with_lakes-data/usgs_lakes:landuse_30s_with_lakes/

and then for the "geog_data_path" setting, you can keep it set to the directory where you have all the other fields.
Click to expand...
Hi
To use USGS 24 data instead of MODIS 21,the correct setting in the namelist.wps file is
Code: 
geog_data_res='default','default',

and to use that default, correspondingly you may change the file setting in LANDUSEF category of the GEOGRID.TBL in the geogrid directory as

from
Code: 
 rel_path =  default:modis_landuse_20class_30s_with_lakes/

to
Code: 
abs_path =  default:{/path_to_usgs_lakes_data_place_directory}/
 
Last edited:
You must log in or register to reply here.
Top