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

Reading in land use + impervious fraction data

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.

dashie

New member
Hello, I'm trying to get geogrid to read two data: one impervious fraction, and one land use (only urban categories 31-33; the idea is to merge it with other data later), both from ICLUS. Geogrid is not doing that: it outputs zero for the impervious fraction, and it outputs 30 for the land use (only because I set the category_min to 31).

Things I've tried:
- Recompiling
- Running on a different system (Which is what worked for me last time)
- Running a different version of WPS

I'm sure:
- The datum isn't all zero
- My domain is in the bounds of these data's domain.

I've attached the two data. Any help you could lend would be much appreciated!

index files:

type=continuous
projection=regular_ll
dx=0.0010230381016821
dy=0.0010230381016821
known_x=1.0
known_y=1.0
known_lat=32.545985734331
known_lon=-113.20778097593
wordsize=1
tile_x=2193
tile_y=1821
tile_z=1
scale_factor=0.01
units="fraction"
description="ICLUS2 impervious fraction"

---------------------

type=categorical
category_min=31
category_max=33
projection=regular_ll
dx=0.000883023796
dy=0.000883023796
known_x=1.0
known_y=1.0
known_lon=-113.20776968534
known_lat=32.546872618955
wordsize=1
missing_value=0.
tile_x=2540
tile_y=2109
tile_z=1
units="category"
description="3-category urban ICLUS2 urban (90m)"

-----------------
relevant namelist.wps section:

&geogrid
parent_id = 1,1,2,
parent_grid_ratio = 1,3,3,
i_parent_start = 1,50,64,
j_parent_start = 1,50,64,
e_we = 151,151,100,
e_sn = 151,151,100,
geog_data_res = 'iclus+nlcd2011_9s','iclus+nlcd2011_9s','iclus+nlcd2011_9s'
dx = 9000,
dy = 9000,
map_proj = 'lambert',
ref_lat = 33.43417,
ref_lon = -112.05111,
truelat1 = 33.43417,
truelat2 = 33.43417,
stand_lon = -112.05111,

----------------

Geogrid.tbl entries:

name=LANDUSEF
priority=1
dest_type=categorical
z_dim_name=land_cat
dominant = LU_INDEX
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 = 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 = 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 = 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 = iclus:four_pt
interp_option = default:nearest_neighbor
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 = 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 = iclus:iclus/lu/
rel_path = default:modis_landuse_20class_30s_with_lakes/
===============================
name=FRC_URB2D
priority=1
optional=yes
dest_type=continuous
fill_missing = 0.
interp_option=default: four_pt
rel_path=default:iclus/imp/
flag_in_output=FLAG_FRC_URB2D
===============================
 

Attachments

  • 00001-02540.00001-02109.gz
    567.5 KB · Views: 1
  • 00001-02540.00001-02109.gz
    196 KB · Views: 1
Top