Adding Digital Model Terrain 25 meters

[jmcabrera]

Hi everyone!

I'm trying to add a topography height with a resolution of 25 meters.

I'm using the library GitHub - openwfm/convert_geotiff: A commandline utility for converting GeoTIFF files for use in WRF to covert the .tif.

The code I'm using is:

./convert_geotiff -w 4 -t 1500 -u "meters MSL" - "25m topography" 136_MDT25_GC.tif

After executing it I get some warnings:

proj_create_operation_factory_context: /home/administrator/miniconda3/share/proj/proj.db contains DATABASE.LAYOUT.VERSION.MINOR = 2 whereas a number >= 3 is expected. It comes from another PROJ installation.
pj_obj_create: /home/administrator/miniconda3/share/proj/proj.db contains DATABASE.LAYOUT.VERSION.MINOR = 2 whereas a number >= 3 is expected. It comes from another PROJ installation.
proj_create_operation_factory_context: /home/administrator/miniconda3/share/proj/proj.db contains DATABASE.LAYOUT.VERSION.MINOR = 2 whereas a number >= 3 is expected. It comes from another PROJ installation.
pj_obj_create: /home/administrator/miniconda3/share/proj/proj.db contains DATABASE.LAYOUT.VERSION.MINOR = 2 whereas a number >= 3 is expected. It comes from another PROJ installation.

However the files are created.

After that, I modify the GEOGRID.TBL.ARW, add the geodata "mtd_25+default" and run geogrid.

The running complete "correctly" however when I see the files geo_em.nc with the ncview and plot the variable HGT_M I get:

min and max both 0 for variable HGT_M

Thanks!

 

[jmcabrera]

Hi,
I could solve the problem changing the projectio from mercator to regular_ll.

Now seem that work for domain3 (240m) before with default and after withe the new DMT


and the domain4 (48m)


However, in the domain2, I'm getting some strange values, as you see in the range the values came from -3500 to 4536.73, I'm don't know why ir making this "border" around the island.


This is the .tif file which values came from 0 to 1910:

 

[jmcabrera]

I fix the problem modifying the GEOGRID.TBL

I add a water mask, however, I'm not sure if this will couse problems with lakes.

name = HGT_M
priority = 1
masked = water
dest_type = continuous
smooth_option = smth-desmth_special; smooth_passes=1
fill_missing=0.

interp_option = mtd_25:four_pt+average_4pt
interp_option = gmted2010_30s:four_pt+average_4pt
interp_option = gtopo_30s:average_gcell(4.0)+four_pt+average_4pt
interp_option = gtopo_2m:four_pt
interp_option = gtopo_5m:four_pt
interp_option = gtopo_10m:four_pt
interp_option = lowres:average_gcell(4.0)+four_pt
interp_option = default:average_gcell(4.0)+four_pt+average_4pt
interp_option = water:nearest_neighbor

rel_path = mtd_25:mtd_25/
rel_path = gmted2010_30s:topo_gmted2010_30s/
rel_path = gtopo_30s:topo_30s/
rel_path = gtopo_2m:topo_2m/
rel_path = gtopo_5m:topo_5m/
rel_path = gtopo_10m:topo_10m/
rel_path = lowres:topo_gmted2010_5m/
rel_path = default:topo_gmted2010_30s/
rel_path = water:water_mask/

And in geog_data_res i write:

geog_data_res = 'default', 'mtd_25+default', 'mtd_25+default', 'mtd_25+default'