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

1/3 arc-second DEM national map data into WPS geogrid

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.

gjduine

New member
Hello,

I’m looking to incorporate 1/3 arc second DEM data into WRF from national map database (https://apps.nationalmap.gov/downloader/#/). Specifically, from this website I used 1/3 arc-second DEM under Elevation Products (3DEP). I used the .tif file, then used gdal commands to convert into binary format, and changing filenames to respect the geogrid rules considering 00001-10812.00001-10812 etc.

The gdal command was something in this style:
gdal –of ENVI foo.fit temp.dat

The projection works fine, but the values are in the order 10^9. See the 250 m domain (domain 3) in the attached figure. For reference, I also attached a figure that shows all desired domains, but using the standard 30s input domain. With some imagination I can see that using the 1/3 arc-second data the expected terrain elevation is present with a resolved valley. However, the values for terrain elevation are in the range of 10^9.

Using the .flt instead .tif and convert it to binary format gave similar results to what is posted here. I used some of the .hdr output to get some sense of the .tif file, and it shows that in the .tif data, the terrain elevation are in the expected value (0-1500 m). So, my conclusion was that either gdal conversion goes wrong, or my index file has wrong definitions. I’ve been looking into most options of the index file like endian, tile_bdr or word_size but that results in wrong projections, or even less normal values for terrain elevation.

I also looked into the convert_geotiff program that could produce the correct index file, but so far I have not been lucky to install this on the local machine.

Has anyone been successfull using this dataset? Or, perhaps if anyone has some kind of a hint to a solution that would be great.

Thanks for your help,

Gert-Jan
 

Attachments

  • sbareg_d01-d04_30s.png
    sbareg_d01-d04_30s.png
    246.7 KB · Views: 1,051
  • sbareg_d03_13s.png
    sbareg_d03_13s.png
    25.3 KB · Views: 1,051
  • index.txt
    376 bytes · Views: 35
Hi,
This may lie outside the scope of what we are able to support, but I've passed this along to a colleague who may know. If so, they will post a response. Otherwise, hopefully someone in the community will be able to assist.
 
Hi,

thanks for your reply. I just wanted to confirm that I managed to install and run convert_geotiff following procedures written up here https://wiki.openwfm.org/mediawiki/index.php/How_to_run_WRF-SFIRE_with_real_data and here https://wiki.openwfm.org/wiki/How_to_convert_data_for_Geogrid

so this uses the wrfxpy environment in python. Then, using the downloaded dem tif file described below, ran the code using ZSF for terrain elevation, which will give HGT_M input for geogrid. The program gives index file directly, which is perfect.

I'm also interested in using land cover at similar resolution (order of 10-30 meter). Nevertheless, there's somehow no conversion available for land use in the convert_geotiff (only fuel category for fire). I will follow up in a different post however as it considers a different problem.

Thanks,

Gert-Jan
 
Top