NetCDF Analysis of Temperature & Precipitation Data + ModelBuilder

4925
16
06-07-2017 03:37 PM
ChrisMcCloud
New Contributor

I'm working on a project examining potential environmental changes in the Chesapeake Bay watershed and the resulting potential impacts on local jurisdictions and planning.  As part of this, I have a number of NetCDF files containing temperature and precipitation data projections across 93 years.  How do I, first, separate out certain years as raster files with attribute tables containing the temperature and precipitation data, and second, automate the process using ModelBuilder so that it is easier to duplicate the process across different jurisdictions and regions? I'm using ArcGIS Advanced 10.4.1.   netcdf data and model builder netcdf

0 Kudos
16 Replies
ChrisMcCloud
New Contributor

That error appeared both when I ran the Copy Raster tool directly on the output of Make NetCDF Raster and also when I attempted to run a Clip after attempting to project it into NAD 83/MD State Plane.  I've wondered about the table name & table not found portions of the error myself; there's no obvious cause for it, with no attribute table involved. I've also wondered about the table not found error.  

0 Kudos
ChrisMcCloud
New Contributor

I'll try and post a screenshot of the model as well, presuming that's what you mean by "job stream," as I'm not entirely sure if that's the case, keeping in mind that I have tried several different orders of the tools as well as several different tools, in some cases, and continued to get the same result (aside from the situation described below, which removed this problem but seems to have created another in the process).  

0 Kudos
ChrisRingo
Occasional Contributor

You could export the model to python and attach the code.  Even more helpful, you could run that code in the python window and provide the results.

Also, did you download the data from a site to which you could provide a link?  

At this point unless I can see exactly what you're doing and what data you're doing it to, I really can't speculate on what might be going on.  Maybe someone else has an idea from what you've provided so far, but I would need to see exactly what you're working with.

0 Kudos
ChrisMcCloud
New Contributor

I'll look into doing the code export as soon as I get my external drive with my data on it hooked up again, probably within an hour or so.  

The FTP link to the download for the data I am currently using is ftp://gdo-dcp.ucllnl.org/pub/dcp/archive/cmip5/bcsd/BCSD/bcc-csm1-1/rcp26/mon/r1i1p1/pr/.

0 Kudos
ChrisMcCloud
New Contributor

The output from the ModelBuilder process is as follows:

Executing (Make NetCDF Raster Layer (25)): MakeNetCDFRasterLayer "D:\Todd Research\GIS Datasets\Climate_Projections\CMIP_5_Data2\Monthly\RCP_4.5\Precipitation\BCSD_0.125deg_pr_Amon_bcc-csm1-1_rcp45_r1i1p1_200601-209912.nc" pr longitude latitude Precip_bcc-csm1-1_RCP45_Dec2025 # "time '12/16/2025 12:00:00 PM'" BY_VALUE
Start Time: Wed Jun 21 14:52:43 2017
Succeeded at Wed Jun 21 14:52:43 2017 (Elapsed Time: 0.12 seconds)
Executing (Project): Project "D:\Todd Research\Processing\Chesapeake_Bay.gdb\Somerset_County_MD_2016" "D:\Todd Research\Processing\Chesapeake_Bay.gdb\Somerset_County_MD_2016_WGS84" PROJCS['WGS_1984_World_Mercator',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Mercator'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',0.0],PARAMETER['Standard_Parallel_1',0.0],UNIT['Meter',1.0]] WGS_1984_(ITRF00)_To_NAD_1983 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]] NO_PRESERVE_SHAPE # NO_VERTICAL
Start Time: Wed Jun 21 14:52:43 2017
Succeeded at Wed Jun 21 14:52:44 2017 (Elapsed Time: 1.66 seconds)
Executing (Clip (159)): Clip Precip_bcc-csm1-1_RCP45_Dec2025 "283.763542 37.886605 284.460367 38.298944" "D:\Todd Research\Processing\Tangier_Smith_Islands.gdb\Precip_bcc-csm1-1_RCP45_Dec2025_Somerset_MD_WGS84" "D:\Todd Research\Processing\Chesapeake_Bay.gdb\Somerset_County_MD_2016_WGS84" 1.000000e+020 NONE NO_MAINTAIN_EXTENT
Start Time: Wed Jun 21 14:52:45 2017
ERROR 999999: Error executing function.
The table name is invalid.
No spatial reference exists.
FDO error: -2147220948
The table was not found. [Precip_bcc-csm1-1_RCP45_Dec2025_Somerset_MD_WGS84]
Failed to execute (Clip (159)).
Failed at Wed Jun 21 14:52:45 2017 (Elapsed Time: 0.14 seconds)

0 Kudos
IanMurray
Frequent Contributor

As for your most recent error, your output file(Precip_bcc-csm1-1_RCP45_Dec2025_Somerset_MD_WGS84)  is an ESRI Grid which can only be 13 characters or less, you need to add a file extension to allow for a longer name(I'd keep going with .tif).

I guess that also begs the question why your initial raster is even created since it is also apparently a Grid ( I don't see a file extension) that is too long of a file name.  (Precip_bcc-csm1-1_RCP45_Dec2025)

0 Kudos
ChrisMcCloud
New Contributor

I tried running Project Raster on the output of the Make NetCDF Raster Tool, projecting it into the already-applied coordinate system, and then running the time series extract script from there. That seems to resolve the spatial error, as I did not get that.  However, it also seems to remove the multidimensional component of the NetCDF raster in the process of running the project at that point, so only one raster is generated by the script, as opposed to one per NetCDF time series entry.  So that appears to solve one problem and create another. 

0 Kudos