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

Land use geotiff file conversion to binary to run on geogrid

laurizio

New member
Hi!

I'm new to WRF. I created a new land use and land cover map based on the ESACCI data. I reclassified the data to IGBP-modified classes with 20 categories, so the values of each pixel are integers.
Convert GeoTiff file to binary via gdal:

> gdal_translate -of ENVI igbp-esacci.tif igbp-esacci.bin
> Input file size is 36000, 18000
> 0...10...20...30...40...50...60...70...80...90...100 – done.

Generating the files: igbp-esacci.bin, igbp_esacci.hdr and igbp_esacci.bin.aux.xml.

I renamed the binary file to 00001-36000. 00001-18000 and created the “index” file

“type=category
category_min=1
category_max=20
projection=polar_wgs84
dx=0.01
dy=-0.01
known_x=1.0
known_y=1.0
sab_lat=90
lon_known=180
word size = 1
tile_x=36000
tile_y=18000
tile_z=1
units="category"
missing_value=255
description="Noah-MP 20-category IGBP-MODIS based on ESACCI land use"
mminlu="MODIFIED_IGBP_MODIS_NOAH"
iswater=17
island=17
isice=15
isurbano=13”

Then I added to the GEOGRID.TBL file:

> landmask_water = esacci:17 # Calculate a land mask from this field
> interp_option = esacci:nearest_neighbor
> rel_path = esacci:esacci/

Added to namelist.input

> geog_data_res = “modis_30S+esacci”

When I run geogrid.exec the following error occurred when trying to open the file geo_em.d01.nc in ncview:
> Min and Max both 1e+20 for variable LANDUSEF
 
Top