Implementing Urban LCZs as landuse data source

Dataset specific topics involving WPS.

Implementing Urban LCZs as landuse data source

Postby amandanmc » Sat Feb 24, 2018 12:36 am

Dear all,
I've been trying for a while now to create a landuse dataset based on the URBAN Local Climate Zones (LCZ) following instructions from the paper "Urbanized WRF modelling using WUDAPT", but I'm not having any success. Procedure says I should extract my classification (both urban and natural classes) from SAGA in .txt and use the script rd_wr_binary.F90 to project/interpolate the data on a regular grid as needed by WPS. Program is also responsible for zeroing the natural classes of the LCZ, leaving the urban classes only. The final result should be an urban area classified according to the LCZs and the natural classes from USGS.

After running rd_wr_binary.exe, I got a file which I renamed as 00001-00874.00001-00828 (00001-ncols.00001-nrows). The index file related to the new landuse has the following info:

type=categorical
category_min=31
category_max=40
projection=regular_ll
dx=0.0010987375
dy=0.0010986719
known_x=1.0
known_y=1.0
known_lat=-20.39015
known_lon=-44.51095
wordsize=1
tile_x=874 #ncols
tile_y=828 #nrows
tile_z=1
units="category"
description="10-category LCZ"

Paper says I should also add a line the line “missing_value=0.”, but when I do, geogrid.exe doesn't read my landuse file.

I also added a few lines on GEOGRID.TBL to point to my new landuse dataset as below:
===============================
name=LANDUSEF
priority=2
dest_type=categorical
z_dim_name=land_cat
interp_option = wudapt:nearest_neighbor
abs_path= wudapt:/share/dados/amandanmc/modelagens/lcz_120m/
===============================
The variable geog_data_res in my namelist.wps is set as = 'wudapt+usgs_30s'

While running geogrid.exe I get a lots of warning like these:

WARNING: In source tile /share/dados/amandanmc/modelagens/lcz_120m/00001-00874.00001-00828, point (1, 1) has an invalid category of 0
WARNING: Attempted to assign an invalid category 0 to grid point (1, 1)

And when I check Geogrid.log, I see that, apparently, my dataset is not being used:
2018-02-23 19:55:11.773 --- INFORM: Using wudapt interpolator sequence for LANDUSEF.
2018-02-23 19:55:11.773 --- INFORM: Using wudapt data source for LANDUSEF.
2018-02-23 19:55:11.773 --- INFORM: For LANDUSEF, couldn't find interpolator sequence for resolution wudapt.
2018-02-23 19:55:11.773 --- INFORM: Using usgs_30s interpolator sequence for LANDUSEF.
2018-02-23 19:55:11.773 --- INFORM: For LANDUSEF, couldn't find wudapt data source.
2018-02-23 19:55:11.773 --- INFORM: Using usgs_30s data source for LANDUSEF.

Checking LU_INDEX from geo_em on VERDI I see that my urban classes are there (despite the warnings saying wudapt data source weren’t found) but the natural classes from USGS are not being used. I understand my dataset has priority over USGS and it,probably assigning strange values to the remain cells, but as geogrid.exe is not running when I add the line "missing_value=0." to the index, how can I solve the problem?

If anybody had any similar issue trying to implement the urban LCZs on WRF or have any idea on why to proceed, the help would be very much appreciated!! Thank you in advance.
amandanmc
 
Posts: 1
Joined: Fri Feb 23, 2018 10:06 pm

Return to Working with Various Datasets

Who is online

Users browsing this forum: Bing [Bot] and 6 guests