LANDUSE issues

[jcoutinho]

Hello,

I am trying to update the land-use dataset in WPS; however, the values in LU_INDEX and LANDUSEF are not being updated correctly.

So far, I have:

• Converted my land-use data to the required binary format;
• Updated the namelist.wps file;
• Modified the GEOGRID.TBL.ARW file.

Despite these changes, the output fields do not appear to reflect the new land-use data.

Am I missing any additional steps in the process? Any guidance would be greatly appreciated.

Thank you for your help.

&geogrid
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 1201,
e_sn = 1261,
geog_data_res = 'mylandusedata+default',
dx = 4000,
dy = 4000,
map_proj = 'mercator',
ref_lat = -13.50,
ref_lon = -54,
truelat1 = -13.5,
truelat2 = 0.0,
stand_lon = -54.0,
geog_data_path = '/scratch/jcoutinh/Build_WRF/WPS_GEOG'
/


name=LANDUSEF
priority=1
dest_type=categorical
z_dim_name=land_cat
dominant = LU_INDEX
landmask_water = mylandusedata:17
landmask_water = nlcd2006_9s:17 # Calculate a landmask from this field
landmask_water = nlcd2006_30s:17 # Calculate a landmask from this field
landmask_water = nlcd2011_9s:17 # Calculate a landmask from this field
landmask_water = nlcd2006:17 # Calculate a landmask from this field
landmask_water = ssib_10m:16 # Calculate a landmask from this field
landmask_water = ssib_5m:16 # Calculate a landmask from this field
landmask_water = modis_15s:17 # Calculate a landmask from this field
landmask_water = modis_15s_lake:17,21 # Calculate a landmask from this field
landmask_water = modis_30s:17 # Calculate a landmask from this field
landmask_water = modis_30s_lake:17,21 # Calculate a landmask from this field
landmask_water = usgs_30s:16 # Calculate a landmask from this field
landmask_water = usgs_lakes:16,28 # Calculate a landmask from this field
landmask_water = modis_lakes:17,21 # Calculate a landmask from this field
landmask_water = modis_2010:17 # Calculate a landmask from this field
landmask_water = modis_2010_lake:17,21 # Calculate a landmask from this field
landmask_water = modis_2050:17 # Calculate a landmask from this field
landmask_water = modis_2050_lake:17,21 # Calculate a landmask from this field
landmask_water = usgs_2m:16 # Calculate a landmask from this field
landmask_water = usgs_5m:16 # Calculate a landmask from this field
landmask_water = usgs_10m:16 # Calculate a landmask from this field
landmask_water = lowres:17,21 # Calculate a landmask from this field
landmask_water = default:17,21 # Calculate a landmask from this field
interp_option = mylandusedata:nearest_neighbor
interp_option = nlcd2006_9s:average_gcell(0.0)
interp_option = nlcd2006_30s:average_gcell(0.0)
interp_option = nlcd2011_9s:average_gcell(0.0)
interp_option = nlcd2006:nearest_neighbor
interp_option = ssib_10m:four_pt
interp_option = ssib_5m:four_pt
interp_option = modis_15s:nearest_neighbor
interp_option = modis_15s_lake:nearest_neighbor
interp_option = modis_30s:nearest_neighbor
interp_option = modis_30s_lake:nearest_neighbor
interp_option = usgs_30s:nearest_neighbor
interp_option = usgs_lakes:nearest_neighbor
interp_option = modis_lakes:nearest_neighbor
interp_option = modis_2010:nearest_neighbor
interp_option = modis_2010_lake:nearest_neighbor
interp_option = modis_2050:nearest_neighbor
interp_option = modis_2050_lake:nearest_neighbor
interp_option = usgs_2m:four_pt
interp_option = usgs_5m:four_pt
interp_option = usgs_10m:four_pt
interp_option = lowres:average_gcell(4.0)+four_pt
interp_option = default:nearest_neighbor
rel_path = mylandusedata:mylandusedata/
rel_path = nlcd2006_9s:nlcd2006_ll_9s/
rel_path = nlcd2006_30s:nlcd2006_ll_30s/
rel_path = nlcd2011_9s:nlcd2011_ll_9s/
rel_path = nlcd2006:nlcd2006_ll_30s/
rel_path = ssib_10m:ssib_landuse_10m/
rel_path = ssib_5m:ssib_landuse_5m/
rel_path = modis_15s:modis_landuse_20class_15s/
rel_path = modis_15s_lake:modis_landuse_20class_15s_with_lakes/
rel_path = modis_30s:modis_landuse_20class_30s/
rel_path = modis_30s_lake:modis_landuse_20class_30s_with_lakes/
rel_path = usgs_30s:landuse_30s/
rel_path = usgs_lakes:landuse_30s_with_lakes/
rel_path = modis_lakes:modis_landuse_21class_30s/
rel_path = modis_2010:slucm/modis_landuse_20class_15s_2010s/
rel_path = modis_2010_lake:slucm/modis_landuse_20class_15s_2010s_with_lakes/
rel_path = modis_2050:slucm/modis_landuse_20class_15s_2050s/
rel_path = modis_2050_lake:slucm/modis_landuse_20class_15s_2050s_with_lakes/
rel_path = usgs_2m:landuse_2m/
rel_path = usgs_5m:landuse_5m/
rel_path = usgs_10m:landuse_10m/
rel_path = lowres:modis_landuse_20class_5m_with_lakes/
rel_path = default:modis_landuse_20class_30s_with_lakes/

My index file:
type = categorical
signed = yes
projection = regular_ll
dx = 0.036
dy = 0.036
known_x = 1.0
known_y = 1.0
known_lat = -34.0
known_lon = -74.0
row_order = top_bottom
tile_x = 1111
tile_y = 1095
tile_z = 1
wordsize = 2
missing_value = 0
description = mylanduse data 4km

 

[kwerner]

Hi, See this post that discusses creating a separate section in GEOGRID.TBL for the landuse you wish to take priority - giving it a priority of 2. See if those changes make any difference in your results.

 

[huyang]

Hi, See this post that discusses creating a separate section in GEOGRID.TBL for the landuse you wish to take priority - giving it a priority of 2. See if those changes make any difference in your results.

Hi ,

I followed the method you shared, but I keep getting this error message:

ERROR: In GEOGRID.TBL, dominant category at entry 3 cannot have the same name as dominant category at entry 4.

Could you tell me how to resolve this issue? Looking forward to your reply.

 

[jcoutinho]

Hi, See this post that discusses creating a separate section in GEOGRID.TBL for the landuse you wish to take priority - giving it a priority of 2. See if those changes make any difference in your results.

Hi!
Thank you for your reply. The previous errors have been addressed. However it, looks like that the model is still not reading my new land use map.
Please find attached a compressed file containing the namelist.wps, index, binary land use map, and GEOGRID.TBL.ARW.
Could you please take a look and let me know if you have any suggestions on what might be causing this issue?