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