ARC GIS Process GEOGRID error

[Tamara Pletzer]

Hi, 

I was following the ARC GIS youtube tutorial for the Croton Test Case and got an error for the Process GEOGRID File step. 

Here is the error log:

Begining processing on Thu Mar 04 09:06:30 2021

64-bit: False

Input parameters:

    Parameter: Input GEOGRID File: C:\Users\pleta922\Desktop\GIS_Preprocessor_Test_Data\Croton_Lambert\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\pleta922\Desktop\GIS_Preprocessor_Test_Data\Croton_Lambert\NED_30m_02b_Croton.tif

    Parameter: Regridding (nest) Factor: 4

    Parameter: Number of routing grid cells to define stream: 32

    Parameter: OVROUGHRTFAC Value: 1

    Parameter: RETDEPRTFAC Value: 1

    Parameter: Channel Initiation Points Feature Class: None

    Parameter: Output ZIP File: C:\Users\pleta922\Desktop\GIS_Preprocessor_Test_Data\Croton_Lambert\NED30m_rf4_th32.zip

None

Step 1: NetCDF Conversion initiated...

    Map Projection: Lambert Conformal Conic

    Using MOAD_CEN_LAT for latitude of origin.

    Using Standard Parallel 2 in Lambert Conformal Conic map projection.

    Georeferencing step completed without error.

    New grid spacing: dx=250.0, dy=250.0

    New dimensions: rows=64, cols=60

    The GEOGRID File resolution is 1000.0sm

    Proj4: +proj=lcc +units=m +a=6370000.0 +b=6370000.0 +lat_1=30.0 +lat_2=60.0 +lat_0=41.471 +lon_0=-97.0 +x_0=0 +y_0=0 +k_0=1.0 +nadgrids=@null +wktext  +no_defs 

    GeoTransform: 1841999.59153 1000.0 0 278496.978218 0 -1000.0

    Created projection definition from input NetCDF GEOGRID file.

Creating CF-netCDF File.

    Esri PE String: PROJCS['Sphere_Lambert_Conformal_Conic',GEOGCS['GCS_Sphere',DATUM['D_Sphere',SPHEROID['Sphere',6370000.0,0.0]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Lambert_Conformal_Conic'],PARAMETER['false_easting',0.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',-97.0],PARAMETER['standard_parallel_1',30.0],PARAMETER['standard_parallel_2',60.0],PARAMETER['latitude_of_origin',41.471],UNIT['Meter',1.0]];-35691800 -29233700 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision

    Map Projection of input raster : lambert_conformal_conic

    Dimensions created after  0.01 seconds.

    Starting Process: Building to XMap/YMap

    Conversion of input raster to XMap/YMap completed without error.

    Coordinate variables and variable attributes set after  0.01 seconds.

    netCDF global attributes set after  0.02 seconds.

Step 2 initiated...

    The GEOGRID File resolution is 1000.0sm

    The High-resolution dataset will be 250.0m

    WKT for CRS 1:

      PROJCS['Sphere_Lambert_Conformal_Conic',GEOGCS['GCS_Sphere',DATUM['D_Sphere',SPHEROID['Sphere',6370000.0,0.0]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Lambert_Conformal_Conic'],PARAMETER['false_easting',0.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',-97.0],PARAMETER['standard_parallel_1',30.0],PARAMETER['standard_parallel_2',60.0],PARAMETER['latitude_of_origin',41.471],UNIT['Meter',1.0]];-35691800 -29233700 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision

    WKT for CRS 2:

      PROJCS['NAD_1983_Albers',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['central_meridian',-96.0],PARAMETER['Standard_Parallel_1',29.5],PARAMETER['Standard_Parallel_2',45.5],PARAMETER['latitude_of_origin',23.0],UNIT['Meter',1.0]],VERTCS['Unknown VCS',VDATUM['Unknown'],PARAMETER['Vertical_Shift',0.0],PARAMETER['Direction',1.0],UNIT['User_Defined_Unit',0.01]];-16901100 -6972200 266467840.990852;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision

    CRS differs between datasets.

    Custom geotransformation will be necessary.

    ArcGIS version 10.6 detected. Transformation: GeoTransform_Null_

I've attached the ArcCatalog Results messages as well. I have the newest wrf hydro preprocessing tools from github and spatial analyst extension box is ticked.

Any help is appreciated thanks!

 

[Kevin]

Tamara,

If you have access to ArcGIS Pro, I would try that, as we have found that platform to be more stable than ArcGIS Desktop. In any case, have you tried to run the Flow Accumulation tool (Spatial Analyst Tools > Hydrology > Flow Accumulation) to see if the tool is working properly? According to the screenshot you provided, the problem appears to be in the specification of the input parameter "flow_direction_type", although the ArcGIS Desktop 10.6 documentation for that tool shows it is a valid optional parameter name:

     desktop.arcgis.com/en/arcmap/10.6/tools/spatial-analyst-toolbox/flow-accumulation.htm#

If the Flow Accumulation geoprocessing tool works, then I would try the same thing on the python console (within ArcGIS), specifying the "flow_direction_type" manually using syntax like the following:

    arcpy.CheckOutExtension("Spatial")

    from arcpy.sa import *

    outFlowAccumulation = FlowAccumulation( "flowdir"  , '#', data_type='FLOAT', flow_direction_type="D8")

To make the above code snippet work, you'll need a raster layer 'flowdir', which is a flow direction raster. You can create that using the Flow Direction tool or see if you can scavenge the flow direction raster from your failed tool run, in the 'scratchdir' directory.

Let me know what you find out,

Kevin

[Tamara Pletzer]

Hi Kevin,

Thank you for your speedy response. I managed to get it to work with ArcGIS Pro, thank you for the tip!