Error for ArcGIS Tools for Preparing WRF-Hydro Routing Grids

38 views
Skip to first unread message

강구협(일반대학원 건설환경공학과)

unread,
Dec 28, 2023, 2:38:32 PM12/28/23
to wrf-hydro_users
Hi, I am a newbie to WRF-Hydro and I am struggling with ArcGIS Tools for Preparing WRF-Hydro. My ArcGIS is pro version 3.2 and I keep getting an error in Process GEOGRID File. Below is the full error and I have shared the geo_em.d01.nc I used. And I shared a capture of the hydroshed that I used. Any help would be greatly appreciated.
Best Regards,
Guhyup

Error messages
Begining processing on Fri Dec 29 04:21:18 2023
64-bit: True
Input parameters:
    Parameter: Input GEOGRID File: C:\Users\chi\Desktop\Guhyup\WRF-Hydro\geo_em.d01.nc
    Parameter: Forecast Points (CSV): None
    Parameter: Mask CHANNELGRID variable to forecast basins?: false
    Parameter: Create reach-based routing (RouteLink) files?: false
    Parameter: Create lake parameter (LAKEPARM) file?: false
    Parameter: Reservoirs Shapefile or Feature Class: None
    Parameter: Input Elevation Raster: C:\Users\chi\Desktop\Guhyup\WRF-Hydro\hyd_as_dem_15s.tif
    Parameter: Regridding (nest) Factor: 10
    Parameter: Number of routing grid cells to define stream: 100
    Parameter: OVROUGHRTFAC Value: 1
    Parameter: RETDEPRTFAC Value: 1
    Parameter: Channel Initiation Points Feature Class: None
    Parameter: Output ZIP File: C:\Users\chi\Documents\ArcGIS\Projects\MyProject2\test.zip
None
Step 1: NetCDF Conversion initiated...
    Map Projection: Mercator
    Using MOAD_CEN_LAT for latitude of origin.
    Georeferencing step completed without error.
    The GEOGRID File resolution is 1500.0sm
    Proj4: +proj=merc +units=m +a=6370000.0 +b=6370000.0 +lon_0=126.65 +lat_ts=35.055
    GeoTransform: -69749.52107774622 1500.0 0 3511571.0158267864 0 -1500.0
    Created projection definition from input NetCDF GEOGRID file.
Creating CF-netCDF File.
    Esri PE String: 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",126.65],PARAMETER["Standard_Parallel_1",35.055],UNIT["Meter",1.0]];-16381900 -24723900 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision
    Map Projection of input raster : mercator
    Dimensions created after  0.00 seconds.
    Coordinate variables and variable attributes set after  0.01 seconds.
    netCDF global attributes set after  0.01 seconds.
Step 2 initiated...
    The GEOGRID File resolution is 1500.0sm
    The High-resolution dataset will be 150.0m
    WKT for CRS 1:
      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",126.65],PARAMETER["Standard_Parallel_1",35.055],UNIT["Meter",1.0]];-16381900 -24723900 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision
    WKT for CRS 2:
      GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]];-400 -400 11258999068426.2;-100000 10000;-100000 10000;8.98315284119521E-09;0.001;0.001;IsHighPrecision
    CRS differs between datasets.
    Custom geotransformation will be necessary.
    Transformation: GeoTransform_Null_WRFHydro
    MakeRasterLayer process completed without error.
    The coarse grid has 135 rows and 93 columns.
    The input elevation grid (before reprojection) has 438 rows and 368 columns.
    Projecting input elevation data to WRF coordinate system.
    ArcGIS version 3.2 found. Using Custom Geotransformation (GeoTransform_Null_WRFHydro)
    Finished projecting input elevation data to WRF coordinate system.
    The fine grid (before ExtractByMask) has 1362 rows and 941 columns.
    Fine Grid has 1350 rows and 930 columns.
    Step 2 completed without error.
    Deriving geocentric coordinates on routing grid from bilinear interpolation of geogrid coordinates.
Creating CF-netCDF File.
    Esri PE String: 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",126.65],PARAMETER["Standard_Parallel_1",35.055],UNIT["Meter",1.0]];-16381900 -24723900 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision
    Map Projection of input raster : mercator
    Dimensions created after  0.00 seconds.
    Coordinate variables and variable attributes set after  0.03 seconds.
    Proceeding to add LATITUDE and LONGITUDE variables after  0.03 seconds.
    Variables populated after  0.03 seconds.
    Process completed without error.
    LATITUDE and LONGITUDE variables and variable attributes set after  0.04 seconds.
    netCDF global attributes set after  0.04 seconds.
    Process: landuse written to output netCDF.
Step 4 initiated...
    Channelgrid will not be masked to basins.
    Reach-based routing files will not be created.
    Process: FLOWDIRECTION written to output netCDF.
    Flow accumulation will be thresholded to build channel pixels.
    Process: FLOWACC written to output netCDF.
    Process: TOPOGRAPHY written to output netCDF.
    Process: RETDEPRTFAC written to output netCDF.
    Process: OVROUGHRTFAC written to output netCDF.
    Process: LKSATFAC written to output netCDF.
    Process: STREAMORDER written to output netCDF.
    Process: frxst_pts was empty. Constant value raster created.
    Process: gw_basns was empty. Constant value raster created.
    Process: LAKEGRID written to output netCDF.
    Process: CHANNELGRID written to output netCDF.
    Step 4 completed without error.
Beginning to build 2D groundwater basin inputs
  Building groundwater inputs using FullDom LINKID local basins
    Generating LINKID grid from CHANNELGRID and FLOWDIRECTION
    Creating layer from netCDF variable CHANNELGRID
Finished building groundwater basin grids in 26.82 seconds
geo_em.d01.nc
tif.jpg

Kevin

unread,
Dec 28, 2023, 11:33:25 PM12/28/23
to wrf-hydro_users, 강구협(일반대학원 건설환경공학과)
It looks like the script got most of the way through the GIS Pre-processing but may have failed when building the groundwater inputs. Do you have a screenshot of the traceback information so we can see what the exact error was?

Thanks,

Kevin

Kevin

unread,
Dec 29, 2023, 12:12:48 PM12/29/23
to wrf-hydro_users, Kevin, 강구협(일반대학원 건설환경공학과)
Also, the script may have generated most of the inputs you need for WRF-Hydro. Check the output directory for the outputs, possibly in a `/scratchdir` temporary directory. My guess is that the groundwater inputs are missing, because this is where there error occurred. You can generate these inputs separately using a different script in the GIS Pre-processing tools, called "Build_Groundwater_Inputs.py". 

Kevin

Reply all
Reply to author
Forward
0 new messages