[SciPy-User] Can create, but not access, netCDF variables with no dimension

[Oleksandr Huziy]

Do you have nco installed?
You could select onlyt p0 into a separate file, then I could look into it?

ncks -v p0 <your complete file> only_p0.nc

I think it is also possible to do with cdo.

In any case it is better to upgrade scipy, though I am not sure if netcdf
has changed since 1.10.

Cheers



2014-04-14 12:04 GMT-04:00 Mannucci, Anthony J (335G) <
anthony.j.mannucci at jpl.nasa.gov>:

>   Thanks for your help using getValue. That should help with the example
> I provided. However, the variable does not show up in the keys(). Here is
> an ncdump excerpt.
>
>  More information: I'm using scipy 0.10. Perhaps this issue has been
> resolved in a later version?
>
>  -Tony
>
>  netcdf apr2011storm_s008D110429_0100_110501_0000 {
> dimensions:
>         time = UNLIMITED ; // (48 currently)
>         lon = 72 ;
>         lat = 36 ;
>         lev = 29 ;
>         ilev = 29 ;
>         mlon = 81 ;
>         mlat = 97 ;
>         mlev = 32 ;
>         imlev = 32 ;
>         mtimedim = 3 ;
>         latlon = 2 ;
>         dtidedim = 2 ;
>         sdtidedim = 10 ;
>         datelen = 24 ;
>         filelen = 1024 ;
> variables:
>         double time(time) ;
>                 time:long_name = "time" ;
>                 time:units = "minutes since 2011-3-21 0:0:0" ;
>                 time:initial_year = 2011 ;
>                 time:initial_day = 80 ;
>                 time:initial_mtime = 80, 0, 0 ;
>         double lon(lon) ;
>                 lon:long_name = "geographic longitude (-west, +east)" ;
>                 lon:units = "degrees_east" ;
>         double lat(lat) ;
>                 lat:long_name = "geographic latitude (-south, +north)" ;
>                 lat:units = "degrees_north" ;
>         double lev(lev) ;
>                 lev:long_name = "midpoint levels" ;
>                 lev:short_name = "ln(p0/p)" ;
>                 lev:units = "" ;
>                 lev:positive = "up" ;
>                 lev:standard_name = "atmosphere_ln_pressure_coordinate" ;
>                 lev:formula_terms = "p0: p0 lev: lev" ;
>                 lev:formula = "p(k) = p0 * exp(-lev(k))" ;
>         double ilev(ilev) ;
>                 ilev:long_name = "interface levels" ;
>                 ilev:short_name = "ln(p0/p)" ;
>                 ilev:units = "" ;
>                 ilev:positive = "up" ;
>                 ilev:standard_name = "atmosphere_ln_pressure_coordinate" ;
>                 ilev:formula_terms = "p0: p0 lev: ilev" ;
>                 ilev:formula = "p(k) = p0 * exp(-ilev(k))" ;
>
>  <<< More stuff. The following don't show up>>>
>          double p0 ;
>                 p0:long_name = "Reference pressure" ;
>                 p0:units = "millibars" ;
>         double p0_model ;
>                 p0_model:long_name = "Reference pressure (as used by the
> model)" ;
>                 p0_model:units = "microbars" ;
>
>
>   --
>  Tony Mannucci
> Supervisor, Ionospheric and Atmospheric Remote Sensing Group
>  Mail-Stop 138-308,                     Tel > (818) 354-1699
>  Jet Propulsion Laboratory,              Fax > (818) 393-5115
>  California Institute of Technology,     Email >
> Tony.Mannucci at jpl.nasa.gov
>  4800 Oak Grove Drive,
> http://scienceandtechnology.jpl.nasa.gov/people/a_mannucci/
>  Pasadena, CA 91109
>
>


-- 
Sasha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20140414/08db4db6/attachment.html>