Remote analysis with GrADSDODS

Diane Portis portis at ATMOS.UIUC.EDU
Mon Feb 26 16:11:49 EST 2007


Hello,

I have been trying out the remote analysis capabilities of gds.  I was
able to do it when the remote dataset was not in a subdirectory (e.g.
the ssta dataset that is located at http://cola8.iges.org:9090/dods),
however when the dataset was located in a subdirectory,  I received an
error message that the dataset did not exist or was not an netcdf file.

The first paste below shows that the file is readable with the sdfopen
command, the second paste lists my script to do remote analysis on that
dataset and the error message when I run the script. I suspect that it
is a syntax error.

Thanks for any insights,
Diane

---Opening the file with sdfopen -
okay-----------------------------------------------------------------------
> ga-> sdfopen
> http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete
> Scanning self-describing file:
> http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete
> SDF file
> http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete
> is open as file 1
> LON set to -220 -0.625
> LAT set to 0 89.625
> LEV set to 1000 1000
> Time values set: 1979:1:1:0 1979:1:1:0

--Trying to do analysis on server side-------------------------------------

script: climt1.gs
> baseurl = 'http://nomads.ncdc.noaa.gov:9091/dods/_expr_'
> datasets = 'NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete'
> expr='tloop(aave(tcdc,lon=-220,lon=0.625,lat=70,lat=89.625))'
> xlim='1:1'
> ylim='1:1'
> zlim='1:1'
> tlim='00Z01jan1979:00Z01may1979'
> 'sdfopen 'baseurl'{'datasets'}{'expr'}{'xlim','ylim','zlim','tlim'}'
> 'define var = result.1'
> 'd var'
> ------------------------------------------------------------------
> ga-> run climt1.gs
> nc_open failure:
> Unknown error 1001
> http://nomads.ncdc.noaa.gov:9091/dods/_expr_{NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete}
(broke this line up when pasting, so it could be read)
> {tloop(aave(tcdc,lon=-220,lon=0.625,lat=70,lat=89.625))}{1:1,1:1,1:1,00Z01jan1979:00Z01may1979}
> does not exist or is not a netCDF file.
>
> Couldn't ingest SDF metadata.
> If this was an HDF-SDS file, try gradshdf.
> DEFINE error:  no file open yet
> DISPLAY error:  no file open yet
> ga->



More information about the gradsusr mailing list