ncar / wrf_hydro_arcgis_preprocessor

WRF-Hydro ArcGIS Pre-processor

ARC-GIS TOOLS FOR PREPARING WRF-HYDRO ROUTING GRIDS

These tools are designed to help WRF-Hydro users create surface input data for WRF-Hydro. Presently these tools consist of scripts for use with the ArcGIS Geographical Information System. This is a stand-alone set of scripts which ArcGIS users can install and run locally on their own systems. ArcGIS is a licensed and proprietary piece of software and we recognize that not all WRF-Hydro users will have access to it.

The purpose of the WRF Hydro GIS Pre-Processing Tool is to create the data layers for terrestrial overland flow, subsurface flow and channel routing processes required by WRF Hydro. The outputs from these tools are geospatial and tabular data layers. This processing workflow for creating WRF-Hydro routing grids is available to users as an ArcGIS Python Toolbox. Instructions on how to install and use the stand-alone tool are provided in the .pdf file contained within the release.

wrf_hydro_arcgis_preprocessor's Issues

LKSATFAC mention + description missing from documentation

Just noting that this variable is missing from the documentation (noted by a collaborator). Appears in the main Technical Description though so at least there's info out there.

Utilities > Add reach-based routing option does not update RouteLink file with gages

The utility for adding reach-based routing parameter files to a routing stack does not add gages defined in the frxst_pts variable in Fulldom_hires.nc to the RouteLink.nc file's 'gages' variable.

wrfhydro Preprocessor: opening netcdf's in QGIS

Hi,

I see that the Fulldom_hires.nc follows CF-compliancy, but I can't seem to get them in the correct location in QGIS.. I see this as crs.spatial_ref attribute:

PROJCS["Sphere_Mercator",GEOGCS["GCS_Sphere",DATUM["D_Sphere",SPHEROID["Sphere",6370000.0,0.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",8.0],PARAMETER["Standard_Parallel_1",53.5],UNIT["Meter",1.0]];-11903700 -17965200 10000;-100000 10000;-100000 10000;0,001;0,001;0,001;IsHighPrecision

Which CRS do I need to indicate in QGIS to get the data at the right location?

thanks!

Kobe

HydroSheds dem data into preprocessor output wrong direction and accumulation model files

I'm using latest hydro preprocessing tool V5.0 with ArcGIS 10.6. When I put WRF
generated geogrid file and HydroSHEDS DEM into this tool, the river channel is wrong and the flow direction, flow accumulation data are wrong according to hydrosheds flow direction and accumulation data.

So ,I have two questions. First, where did the errors come from? how can I avoid them. Second, if I want to replace the result data with the data from hydrosheds, is it possible? As I can see, the result from tool is not in latlon projection while hydrosheds one is.

Thanks!

Error on argis 10.5.1 when

Attached run settings and log
![run_settings](https://user-images.githubusercontent.com/24871303/37336007-bfd42870-26b8-11e8-8b36-f2e16
20180310-1.zip.log
130bcb2.png)

wrf_hydro_arcgis_preprocessor v5.0

Hi,
Please i want to download the version 5.0, can you send me the link.
Thank you.

Error when I add lakes using Add Lake Parameter Arcmap 10.8

Hi all,
I use Arcmap 10.8 and I created GIS_preprocessing output data without lake data. But when I insert the lake then I conclude the error in the GIS_preprocessing module.
For example my input and output data are attached. Could you please test Add Lake Parameter with my lake data and output data? I don't know what's the problem. When I insert output data and lakes data into Add Lake Parameter option, I've got the following error:

DEM file link: https://www.dropbox.com/sh/hmpwobbz9qixxpe/AACdhcfa9RMIWaBZelRsUMN-a/HydroSHEDS_DEM/DEM_30s_GRID?dl=0&preview=eu_dem_30s_grid.zip&subfolder_nav_tracking=1

output.zip
lakes.zip
input.zip

max_elevs

I am not sure if only my computer has this issue. But when running wrf_hydro_functions.py at line 2306: "elevRange = {key:max_elevs[key]-val for key, val in min_elevs.items()} # Get lake depths", I get an error reported as: NameError: global name 'max_elevs' is not defined. It doesn't make much sense to me, since I know it has been defined previously. But if I convert this line into a complete for loop, the error is gone and the program can resume. Here it is an example workaround that works for me:

elevRange = {}
for key,val in min_elevs.items(): elevRange.update({key:max_elevs[key]-val})