Create preliminary Hydrologic Routing Inputs for WRF Hydro single catchment

[LukeParky]

[Martin20494]

• Necessary data:

  • geo_em.d01.nc: domain information
  • gauge_station.csv: coordinates of gauge stations
  • DEM: DEM of the area of interest

• Necessary GIS Preprocessor files: (Source: https://github.com/NCAR/wrf_hydro_gis_preprocessor/tree/master/wrfhydro_gis)

  • Build_Routing_Stack.py: To create Routing network
  • Create_Domain_Boundary_Shapefile.py: To create domain shapefile
  • Forecast_Point_Tools.py: To create gauge station points
  • Examine_Outputs_of_GIS_Preprocessor.py: To create some of outputs of GIS Preprocessor for checking

• Required Parameters: (Source: https://github.com/NCAR/wrf_hydro_training/tree/main/lessons/training)
   -i -WRF/WPS GEOGRID file (geo_em.d0*.nc)
   -d -High-resolution Elevation raster file (Esri GRID, GeoTIFF, VRT, etc.)
   -R -Regridding Factor – nesting relationship of routing:land surface model grid cells
   -t -Minimum basin area threshold (in routing grid cells)
   -o -Output ZIP File containing all script outputs

• Optional Parameters: (Source: https://github.com/NCAR/wrf_hydro_training/tree/main/lessons/training)
   --CSV -Station Locations location file (.csv)
   -b -Option to mask channel grids not contributing to provided station locations
   -r -Reach based (Muskingum / Muskingum-Cunge) routing option
   -l -Lake Polygons (polygon feature class or .shp)
   -O -OVROUGHRTFAC – Multiplier on Manning's roughness for overland flow. default=1.0
   -T -RETDEPRTFAC – Multiplier on maximum retention depth before flow is routed as overland flow. default=1.0
   -LKSATFAC – (script global variable) Multiplier on saturated hydraulic conductivity in lateral flow direction. default=1000.0
   --starts -Path to point shapefile or feature class containing channel initiation locations (overrides -t parameter)
   --gw -Path to polygon shapefile or feature class containing prescribed groundwater basin locations

---

• Example for creating routing stack for Mataura river, run this on Anaconda Powershell Prompt with activated environment (here is wrfh_gis_env_002):
  python Build_Routing_Stack.py -i C:\Users\mng42\wrf_wps\example_nz_014\geo_em.d01.nc --CSV C:\Users\mng42\wrf_wps\example_nz_014\river_source.csv -b False -r True -d C:\Users\mng42\wrf_wps\example_nz_014\dem_mataura_filled.tif -R 10 -t 50 -o C:\Users\mng42\wrf_wps\example_nz_014\output\mataura_test_001.zip

• For checking domain shape, run this on Anaconda Powershell Prompt with activated environment (here is wrfh_gis_env_002):
  python Create_Domain_Boundary_Shapefile.py -i C:\Users\mng42\wrf_wps\example_nz_014\geo_em.d01.nc -o C:\Users\mng42\wrf_wps\example_nz_014\output

• For checking gauge stations, run this on Anaconda Powershell Prompt with activated environment (here is wrfh_gis_env_002):
  python Forecast_Point_Tools.py -i C:\Users\mng42\wrf_wps\example_nz_013\river_source.csv -s C:\Users\mng42\wrf_wps\example_nz_013\output\points.shp -o C:\Users\mng42\wrf_wps\example_nz_013\output

• For checking some of other outputs from GIS Preprocessor functions that used python, run this on Anaconda Powershell Prompt with activated environment (here is wrfh_gis_env_002):
  python Examine_Outputs_of_GIS_Preprocessor.py -i C:\Users\mng42\wrf_wps\example_nz_012\output\mataura_test_001.zip -o C:\Users\mng42\wrf_wps\example_nz_012\output\raster_outputs

[Martin20494]

Initial script to create routing stack:

import os


def create_routing_stack(
        domain_nc,
        dem_tif,
        resolution,
        network_dense,
        output_zip=None
):
    command = (
        f'python Build_Routing_Stack.py '
        f'-i "{domain_nc}" '
        f'-b False '
        f'-r True '
        f'-d "{dem_tif}" '
        f'-R {resolution} '
        f'-t {network_dense} '
        f'-o "{output_zip}"'
    )
    os.system(command)

[Martin20494]

The "-b False" should be removed and at the moment the highest resolution for channel grid can be used is 50 m. Higher than this should be considered with machine memory. Furthermore, the locations of gauge stations should be found too to use the forecast points csv (depends on the amount of work, this task can be created as an issue, but at the moment, it is not necessary).