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