Re: trouble with opendap reading CMIP3 through the IRI Data Library

From: Wei Huang <huangwei_at_nyahnyahspammersnyahnyah>
Date: Fri Jan 03 2014 - 13:32:26 MST

Alessandra,

Try this script:

 f=addfile("http://iridl.ldeo.columbia.edu/expert/SOURCES/.WCRP/.CMIP3/%28ipcc4/20c3m%29%40%40/%28ipcc4/20c3m/miroc3_2_medres%29%40%40/.pcmdi.ipcc4.miroc3_2_medres.20c3m.run2.atm.mo.xml/zg/9.8/mul/ta/1005.7/mul/add/hus/2501000/mul/add/plev/30000/10000/RANGEEDGES/plev/integralgrid/plev/differences//name//deltap/def%5Bplev%5Dweighted-average/time/1/monthlyAverage/time/%28Dec%201900%29%28Jan%201999%29RANGE/time/3/runningAverage/time/12/splitstreamgrid/time/%2Fmon/renameGRID/time2/%2FT/renameGRID/-999/setmissing_value/dods", "r")
 print(f)

It should produce:
Variable: f
Type: file
filename: dods
path: http://iridl.ldeo.columbia.edu/expert/SOURCES/.WCRP/.CMIP3/%28ipcc4/20c3m%29%40%40/%28ipcc4/20c3m/miroc3_2_medres%29%40%40/.pcmdi.ipcc4.miroc3_2_medres.20c3m.run2.atm.mo.xml/zg/9.8/mul/ta/1005.7/mul/add/hus/2501000/mul/add/plev/30000/10000/RANGEEDGES/plev/integralgrid/plev/differences//name//deltap/def%5Bplev%5Dweighted-average/time/1/monthlyAverage/time/%28Dec%201900%29%28Jan%201999%29RANGE/time/3/runningAverage/time/12/splitstreamgrid/time/%2Fmon/renameGRID/time2/%2FT/renameGRID/-999/setmissing_value/dods
   file global attributes:
      Conventions : IRIDL
   dimensions:
      T = 98
      lat = 64
      lon = 128
      mon = 12
   variables:
      double lon ( lon )
         standard_name : longitude
         modulus : 360
         bounds : lon_bnds
         uniquename : lon
         units : degree_east
         axis : X
         gridtype : 1

      double lat ( lat )
         standard_name : latitude
         bounds : lat_bnds
         uniquename : lat
         units : degree_north
         axis : Y
         gridtype : 0

      float mon ( mon )
         pointwidth : 3
         calendar : 360
         modulo : 12
         modulus : 12
         gridtype : 1
         units : months since 1960-01-01

      float T ( T )
         calendar : 360
         pointwidth : 12
         gridtype : 0
         units : months since 1960-01-01

      float asum ( T, mon, lat, lon )
         original_name : geopotential height
         long_name : Geopotential Height
         standard_name : geopotential_height
         dodsname : zg
         cell_methods : time: mean (interval: 20 minutes comment: approx. value of variable model time-step)
         file_missing_value : 1e+20
         history : Boxes with less than 0.0% dropped
Boxes with less than 100.0% dropped
Averaged in time with overlapping interval 3
Boxes with less than 0.0% dropped
         units : unitless
         missing_value : -999
         _FillValue : -999

Otherwise, your NCL version is not opendap-enabled.

Wei

huangwei@ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Jan 3, 2014, at 1:10 PM, Alessandra Giannini <alesall@iri.columbia.edu> wrote:

>
> Hello everyone,
>
> I am having trouble reading CMIP3 model output through the IRI Data Library, using opendap, but am at a total loss as to why.
>
> This is the error message I get:
>
> ncvarget: ncid 327680; varid 4: NetCDF: I/O failure
> fatal:NetCdf: An error occurred while attempting to read variable (asum) from file (http://iridl.ldeo.columbia.edu/expert/SOURCES/.WCRP/.CMIP3/%28ipcc4/20c3m%29%40%40/%28ipcc4/20c3m/miroc3_2_medres%29%40%40/.pcmdi.ipcc4.miroc3_2_medres.20c3m.run2.atm.mo.xml/zg/9.8/mul/ta/1005.7/mul/add/hus/2501000/mul/add/plev/30000/10000/RANGEEDGES/plev/integralgrid/plev/differences//name//deltap/def%5Bplev%5Dweighted-average/time/1/monthlyAverage/time/%28Dec%201900%29%28Jan%201999%29RANGE/time/3/runningAverage/time/12/splitstreamgrid/time/%2Fmon/renameGRID/time2/%2FT/renameGRID/-999/setmissing_value/dods)
>
> I am pretty sure the data is fine. If you go to the link:
> <http://iridl.ldeo.columbia.edu/expert/SOURCES/.WCRP/.CMIP3/%28ipcc4/20c3m%29%40%40/%28ipcc4/20c3m/miroc3_2_medres%29%40%40/.pcmdi.ipcc4.miroc3_2_medres.20c3m.run2.atm.mo.xml/zg/9.8/mul/ta/1005.7/mul/add/hus/2501000/mul/add/plev/30000/10000/RANGEEDGES/plev/integralgrid/plev/differences//name//deltap/def%5Bplev%5Dweighted-average/time/1/monthlyAverage/time/%28Dec%201900%29%28Jan%201999%29RANGE/time/3/runningAverage/time/12/splitstreamgrid/time/%2Fmon/renameGRID/time2/%2FT/renameGRID/>
> you can visualize the data by clicking any of the plot icons.
>
> You will note that I am not simply reading data. Through the Data Library I can manipulate the data, and in this case I perform a vertical integral of moist static energy, i.e.
> gz + cT +Lq, from 300 to 100 hPa.
>
> I do not have an issue when I simply read precipitation in a similar way, from this url:
> <http://iridl.ldeo.columbia.edu/expert/SOURCES/.WCRP/.CMIP3/%28ipcc4/20c3m%29%40%40/%28ipcc4/20c3m/gfdl_cm2_0%29%40%40/.pcmdi.ipcc4.gfdl_cm2_0.20c3m.run3.atm.mo.xml/.pr/time/1/monthlyAverage/time/%28Dec%201900%29%28Jan%201999%29RANGE/time/3/runningAverage/time/12/splitstreamgrid/time//mon/renameGRID/time2//T/renameGRID/>
>
> All this is making me think that either the connection times out, or the memory is insufficient. But I can't get my hands around the cause, or solution for that matter...
> I would appreciate it if anyone could interpret this line specifically:
>
> ncvarget: ncid 327680; varid 4: NetCDF: I/O failure
>
> thanks much, and happy new year!
>
> best, alessandra
>
>
>
>
>
>
>
>
>
>
>
> --
> Alessandra Giannini
> IRI for Climate and Society - The Earth Institute at Columbia University
> P.O. Box 1000, Palisades NY 10964-8000
> phone/fax: +1 845 680-4473/4864 - email: alesall@iri.columbia.edu
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jan 03 13:32:41 2014

This archive was generated by hypermail 2.1.8 : Mon Jan 06 2014 - 13:09:53 MST