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

landuse one value only

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.

kaelel18

Member
Hi. I am trying to change the land cover using MODIS land cover data. Previously, I was able to do this without any problem but right now I am experiencin the same problem again even if I do what I previously did. Here are things I did.

1.. I downloaded the landcover data (see uploaded 2015modis.tar.xz) then transformed it into a binary file with its corresponding index.
Code:
[code][quote]
type = categorical
category_min=1
category_max=17
signed = yes
signed = yes
projection = regular_ll
dx = 0.004166666666306337
dy = 0.00416666666625
known_x = 1.0
known_y = 1.0
known_lat = 29.999999997
known_lon = 100.004166658
wordsize = 2
endian = big
tile_x = 12026
tile_y = 7200
tile_z = 1
missing_value = 255
row_order = bottom_top
units = "category"
description = "MODIS 2015"
mminlu="MODIFIED_IGBP_MODIS_NOAH""[/quote]
[/code]

2. I added a new entry in my geogrid.tbl

Code:
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 =        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 =        default:modis_landuse_20class_15s/
        interp_option =        ES:nearest_neighbor
	abs_path =	       ES:/home/gil/WRF/WPS_GEOG/2015_modis/
/

3. I added this to my namelist.input

Code:
&share
 wrf_core = 'ARW',
 max_dom = 3,
 start_date = '2018-04-21_00:00:00', '2018-04-21_00:00:00', '2018-04-21_00:00:00', 
 end_date   = '2018-04-30_00:00:00', '2018-04-30_00:00:00', '2018-04-30_00:00:00', 
 interval_seconds = 21600,
 io_form_geogrid = 2,
 debug_level = 0,
/

&geogrid
 parent_id         = 1,1,2,
 parent_grid_ratio = 1,3,3,
 i_parent_start    = 1,29,47,
 j_parent_start    = 1,22,40,
 e_we          = 100,130,112,
 e_sn          = 105,124,118,
 geog_data_res = 'ASTER+ES','ASTER+ES','ASTER+ES',
 dx =  9000,
 dy = 9000,
 map_proj =  'mercator',
 ref_lat   = 15.474,
 ref_lon   = 121.125,
 truelat1  = 15.474,
 truelat2  = 0,
 stand_lon = 121.125,
 geog_data_path = '/home/gil/WRF/WPS_GEOG/',
 opt_geogrid_tbl_path = '/home/gil/WRF/WPS/geogrid/',
 ref_x = 50.0,
 ref_y = 52.5,
/

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

&metgrid
 fg_name = 'FILE',
 io_form_metgrid = 2,
 opt_metgrid_tbl_path = '/home/gil/WRF/Test/SLUCM/',
/

&mod_levs
 press_pa = 201300 , 200100 , 100000 ,
             95000 ,  90000 ,
             85000 ,  80000 ,
             75000 ,  70000 ,
             65000 ,  60000 ,
             55000 ,  50000 ,
             45000 ,  40000 ,
             35000 ,  30000 ,
             25000 ,  20000 ,
             15000 ,  10000 ,
              5000 ,   1000
 /


But looking at LU_INDEX I get,

min max both 17 for LU_INDEX

Why does it only show one value? Thanks for the answer
 

Attachments

  • modis2015.tar
    165.2 MB · Views: 70
Can you verify that your land use dataset covers your WRF domain geographically? Based on the 'index' file, it looks like the south-west corner of the dataset begins at ~30 N latitude, while the center latitude of your domain is at 15.474 N latitude.
 
Hi mgduda. Thank you for the response. I generated the domain grids and overlayed it over the land use which I used to convert to binary. The domain is located inside the land cover.


Capture.PNG
 
Your 'index' file specifies that the south-west pixel of the dataset is at (29.999999997 lat, 100.004166658 lon); is this location actually that of the north-west pixel? If so, you could try setting
Code:
known_y = 7200.0
dy = -0.00416666666625
in your index file.
 
Hi mgduda. Apologies for the late response. Your suggestion worked. I am now able to run it right. Thanks!

Hi Ming Chen. No, im not using it. I'm using it for a previous simulation. I was not able to edit it out. :D
 
Top