[gradsusr] gradsusr Digest, Vol 52, Issue 68
Yaqian He
yahe at mix.wvu.edu
Mon Jun 23 08:09:44 EDT 2014
Dear James,
Thanks for your reply.
I added the 'set lat 13 20.0' and 'set lon -15.0 20.0' after 'set dfile
2', however, it still shows that 'Error in gagchk: axis sizes are not the
same Operation error: Incompatable grids'. Could you give me more help?
Do I need to regrid my data?
Thanks very much!
Best wishes!
Yaqian
On Fri, Jun 20, 2014 at 2:51 PM, <gradsusr-request at gradsusr.org> wrote:
> Send gradsusr mailing list submissions to
> gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
> gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
> gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
> 1. Re: Fwd: How to calculate the difference between two NDVI
> datasets (James T. Potemra)
> 2. Re: Displaying a .nc file (David Stanfel)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 20 Jun 2014 08:39:36 -1000
> From: "James T. Potemra" <jimp at hawaii.edu>
> Subject: Re: [gradsusr] Fwd: How to calculate the difference between
> two NDVI datasets
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <53A47FE8.603 at hawaii.edu>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Yaqian,
>
> I think you need to specify all ranges for each file since the grids are
> different. You only define the lat/lon based on the first file
> (GIMMS). To be safe you should set the lat/lon/z/t for each data set.
> So, add the 'set lat 13 20' and 'set lon -15 20' after 'set dfile 2' to
> ensure your two grids are the same.
>
> Jim
>
> On 6/20/14 3:48 AM, Yaqian He wrote:
> > Hi,
> >
> > I want to calculate the average seasonal (Dec to Feb) mean
> > difference between GIMMS and MODIS NDVI, here is script I wrote:
> >
> > 'reinit'
> > *gimms
> > 'sdfopen /Users/geog493/ndvi/ndvi0.5deg.2006.nc
> > <http://ndvi0.5deg.2006.nc>'
> > 'sdfopen /Users/geog493/ndvi/ndvi0.5deg_modis_westafrica.nc
> > <http://ndvi0.5deg_modis_westafrica.nc>'
> > 'set lat 13 20.0'
> > 'set lon -15.0 20.0'
> > * DJF
> > 'define gy01=ave(ndvi.1,t=467,t=472)'
> > 'define gy02=ave(ndvi.1,t=491,t=496)'
> > 'define gy03=ave(ndvi.1,t=515,t=520)'
> > 'define gy04=ave(ndvi.1,t=539,t=544)'
> > 'define gy05=ave(ndvi.1,t=563,t=568)'
> > 'define gy06=ave(ndvi.1,t=587,t=592)'
> > 'define gyclim=(gy01+gy02+gy03+gy04+gy05+gy06)/6.'
> > *'modis ndvi'
> > 'set dfile 2'
> > 'set t 1 307'
> > * DJF
> > 'define my01=ave(ndvi.2,t=15,t=20)'
> > 'define my02=ave(ndvi.2,t=38,t=43)'
> > 'define my03=ave(ndvi.2,t=61,t=66)'
> > 'define my04=ave(ndvi.2,t=84,t=89)'
> > 'define my05=ave(ndvi.2,t=107,t=112)'
> > 'define my06=ave(ndvi.2,t=130,t=135)'
> > 'define myclim=(my01+my02+my03+my04+my05+my06)/6.'
> > 'define difyclim=gyclim-myclim'
> > 'd gyclim-myclim'
> >
> > However, it always shows the error "Error in gagchk: axis sizes are
> > not the same
> > Operation error: Incompatable grids"
> >
> > The two datasets have the same resolution (0.5X0.5), the original
> > information for the two datasets are:
> > modis
> > xdef 81 linear -20 0.5
> > ydef 81 linear 0 0.5
> > gimms
> > xdef 720 linear -179.75 0.5
> > ydef 304 linear -61.75 0.5
> > Are there something wrong here? Could you help me?
> >
> > Thanks in advance!
> >
> > Best wishes!
> > Yaqian
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20140620/f7419945/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Fri, 20 Jun 2014 14:51:20 -0400
> From: David Stanfel <dstanfel at ualberta.ca>
> Subject: Re: [gradsusr] Displaying a .nc file
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <
> CADML8btPBhW7oEkKpLJYLFpsjMESqDfPVf97hFgOdhfokD5iiw at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Yes, the axis stuff is pretty confusing, I've been looking around for what
> I should be using for it and still can't find it. I made the suggested
> changes to my control file:
>
> DSET ^atmos_month.nc.0000
> UNDEF -9.99e8_slp
> UNDEF -9.99e8_temp
> UNDEF -9.99e8_zsurf
> DTYPE netcdf
> TITLE 6 Months Sample Model Output
> *
> XDEF 144 LINEAR 0.0 2.5
> *
> YDEF 3 LINEAR 1 1
> *
> ZDEF 24 LINEAR 1 1
> *
> TDEF 1 LINEAR 01JAN0001 1dy
> *
> VARS 3
> slp 0 t,y,x sea level pressure (Pa)
> temp 24 t,z,y,x temperature (K)
> zsurf 0 y,z surface height (m)
> ENDVARS
>
> This was the result when it was called and when I tried to display the
> three variables that were specified:
>
> dstanfel:bin ninjawarrior343$ ./grads
>
>
> Grid Analysis and Display System (GrADS) Version 2.1.a2
>
> Copyright (c) 1988-2013 by the Institute for Global Environment and Society
> (IGES)
>
> GrADS comes with ABSOLUTELY NO WARRANTY
>
> See file COPYRIGHT for more information
>
>
> Config: v2.1.a2 little-endian readline grib2 netcdf hdf4-sds hdf5
> opendap-grids,stn geotiff shapefile cairo
>
> Issue 'q config' command for more detailed configuration information
>
> Landscape mode? ('n' for portrait): y
>
> GX Package Initialization: Size = 11 8.5
>
> ga-> open test.ctl
>
> Scanning description file: test.ctl
>
> Data file atmos_month.nc.0000 is open as file 1
>
> LON set to 0 360
>
> LAT set to 1 3
>
> LEV set to 1 1
>
> Time values set: 1:1:1:0 1:1:1:0
>
> E set to 1 1
>
> ga-> d temp
>
> *** WARNING! ***
>
> *** The I/O for this variable will be extremely slow
>
> *** because the chunks are too big to fit in the cache
>
> *** chunk size = 41472 bytes (144 * 3 * 24 * 1 * 4 bytes)
>
> *** cache size = -2017461983885971885 bytes
>
> *** There are two ways to control the cache size:
>
> *** add a CACHESIZE entry to the descriptor file
>
> *** or change the cache size scale factor with 'set cachesf'
>
> *** Please read http://iges.org/grads/gadoc/compression.html
>
> gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds dimension
> bound
>
> Data Request Error: Error for variable 'temp'
>
> Error ocurred at column 1
>
> DISPLAY error: Invalid expression
>
> Expression = temp
>
> ga-> zsurf
>
> Unknown command: zsurf
>
> ga-> d zsurf
>
> Open Error on Map Data Set: /usr/local/lib/grads/lowres
>
> Contouring: 0 to 2.7e+34 interval 3e+33
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> Error opening stroke character data set
>
> Data set names = /usr/local/lib/grads/font0.dat ; font0.dat
>
> ga-> d slp
>
> gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds dimension
> bound
>
> Data Request Error: Error for variable 'slp'
>
> Error ocurred at column 1
>
> DISPLAY error: Invalid expression
>
> Expression = slp
>
> ga->
>
> I tried changing the environment of where the font data is stored with:
>
> export GASCRP=$HOME/University/Research_job/supplibs
>
>
> but this did not change the error.
>
>
> Thanks,
>
>
> David
>
>
> On Fri, Jun 20, 2014 at 8:18 AM, Jennifer Adams <jma at cola.iges.org> wrote:
>
> > David,
> > You need to read the documentation about creating descriptors for netcdf
> > data.
> > http://iges.org/grads/gadoc/SDFdescriptorfile.html
> >
> > Your file has numerous coordinate variables, yet the ncdump output shows
> > them all to contain missing data, so I'm not sure how you can define the
> > proper grid. The 'lon' coordinate is promising, it has 144 values, so you
> > could try assuming it's a 2.5-degree grid and create an X axis like this:
> > xdef 144 linear 0 2.5
> > I don't know what the information in the attribute named
> > domain_decomposition means. Is it relevant?
> >
> > The lat axis has only 3 values ? I have no idea how to guess at that, so
> I
> > woud put an abstract axis definition, like this:
> > ydef 3 linear 1 1
> >
> > For the Z axis, most of the variables use pfull, but once again, there is
> > no information about what those values might be, so:
> > zdef 24 linear 1 1
> >
> > The time dimension has size 0, yet most of the variables in your file
> > appear to have a time dimension. I don't know how to interpret that,
> > perhaps just a default will do:
> > tdef 1 linear 01jan0001 1dy
> >
> > As for variables, if you extract all the variables that use the axes lon,
> > lat, or pfull, you get a long list. Here are three examples:
> >
> > float slp(time, lat, lon) ;
> > float temp(time, pfull, lat, lon) ;
> > float zsurf(lat, lon) ;
> >
> > Variable declarations for these three variables would look like this:
> >
> > vars 3
> > slp 0 t,y,x sea level pressure (Pa)
> > temp 24 t,z,y,x temperature (K)
> > zsurf 0 y,z surface height (m)
> > endvars
> >
> > You will also need
> > dtype netcdf
> > undef -9.99e8 missing_value
> >
> > The undef entry is important because the variables in your file have
> > different missing values, so you need to put the name of the attribute in
> > there so GrADS will get it right. The -9.99e8 can be any number -- the
> > attribute name is what's critical here.
> >
> > But I can't say for sure whether this will work because of your 0-size
> > time axis. And the mapping will be off until you figure out what the
> > latitude values are.
> > --Jennifer
> >
> > On Jun 19, 2014, at 1:50 PM, David Stanfel wrote:
> >
> > Here is the output of the ncdump:
> >
> > dstanfel at bugaboo:~/models/mom/work/ESM2M_pi-control_C2> ncdump -c
> > atmos_month.nc.0000
> > netcdf atmos_month.nc {
> > dimensions:
> > lon = 144 ;
> > lonb = 145 ;
> > lat = 3 ;
> > latb = 4 ;
> > time = UNLIMITED ; // (0 currently)
> > nv = 2 ;
> > pfull = 24 ;
> > phalf = 25 ;
> > variables:
> > double lon(lon) ;
> > lon:long_name = "longitude" ;
> > lon:units = "degrees_E" ;
> > lon:cartesian_axis = "X" ;
> > lon:domain_decomposition = 1, 144, 1, 144 ;
> > lon:edges = "lonb" ;
> > double lonb(lonb) ;
> > lonb:long_name = "longitude edges" ;
> > lonb:units = "degrees_E" ;
> > lonb:cartesian_axis = "X" ;
> > double lat(lat) ;
> > lat:long_name = "latitude" ;
> > lat:units = "degrees_N" ;
> > lat:cartesian_axis = "Y" ;
> > lat:domain_decomposition = 1, 90, 1, 3 ;
> > lat:edges = "latb" ;
> > double latb(latb) ;
> > latb:long_name = "latitude edges" ;
> > latb:units = "degrees_N" ;
> > latb:cartesian_axis = "Y" ;
> > latb:domain_decomposition = 1, 91, 1, 4 ;
> > double time(time) ;
> > time:long_name = "time" ;
> > time:units = "days since 0001-01-01 00:00:00" ;
> > time:cartesian_axis = "T" ;
> > time:calendar_type = "NOLEAP" ;
> > time:calendar = "NOLEAP" ;
> > time:bounds = "time_bounds" ;
> > double nv(nv) ;
> > nv:long_name = "vertex number" ;
> > nv:units = "none" ;
> > nv:cartesian_axis = "N" ;
> > double pfull(pfull) ;
> > pfull:long_name = "ref full pressure level" ;
> > pfull:units = "mb" ;
> > pfull:cartesian_axis = "Z" ;
> > pfull:positive = "down" ;
> > pfull:edges = "phalf" ;
> > double phalf(phalf) ;
> > phalf:long_name = "ref half pressure level" ;
> > phalf:units = "mb" ;
> > phalf:cartesian_axis = "Z" ;
> > phalf:positive = "down" ;
> > float fco2nat(time, lat, lon) ;
> > fco2nat:long_name = "flux of carbon dioxide" ;
> > fco2nat:units = "mol CO2/(m2 s)" ;
> > fco2nat:missing_value = 1.e+20f ;
> > fco2nat:_FillValue = 1.e+20f ;
> > fco2nat:cell_methods = "time: mean" ;
> > fco2nat:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float high_cld_amt(time, lat, lon) ;
> > high_cld_amt:long_name = "high cloud amount" ;
> > high_cld_amt:units = "percent" ;
> > high_cld_amt:missing_value = 1.e+20f ;
> > high_cld_amt:_FillValue = 1.e+20f ;
> > high_cld_amt:cell_methods = "time: mean" ;
> > high_cld_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float low_cld_amt(time, lat, lon) ;
> > low_cld_amt:long_name = "low cloud amount" ;
> > low_cld_amt:units = "percent" ;
> > low_cld_amt:missing_value = 1.e+20f ;
> > low_cld_amt:_FillValue = 1.e+20f ;
> > low_cld_amt:cell_methods = "time: mean" ;
> > low_cld_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float mid_cld_amt(time, lat, lon) ;
> > mid_cld_amt:long_name = "mid cloud amount" ;
> > mid_cld_amt:units = "percent" ;
> > mid_cld_amt:missing_value = 1.e+20f ;
> > mid_cld_amt:_FillValue = 1.e+20f ;
> > mid_cld_amt:cell_methods = "time: mean" ;
> > mid_cld_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tot_cld_amt(time, lat, lon) ;
> > tot_cld_amt:long_name = "total cloud amount" ;
> > tot_cld_amt:units = "percent" ;
> > tot_cld_amt:missing_value = 1.e+20f ;
> > tot_cld_amt:_FillValue = 1.e+20f ;
> > tot_cld_amt:cell_methods = "time: mean" ;
> > tot_cld_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float taubx(time, lat, lon) ;
> > taubx:long_name = "x base flux for grav wave drag" ;
> > taubx:units = "kg/m/s2" ;
> > taubx:missing_value = -999.f ;
> > taubx:_FillValue = -999.f ;
> > taubx:cell_methods = "time: mean" ;
> > taubx:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tauby(time, lat, lon) ;
> > tauby:long_name = "y base flux for grav wave drag" ;
> > tauby:units = "kg/m/s2" ;
> > tauby:missing_value = -999.f ;
> > tauby:_FillValue = -999.f ;
> > tauby:cell_methods = "time: mean" ;
> > tauby:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float udt_gwd(time, pfull, lat, lon) ;
> > udt_gwd:long_name = "u wind tendency for gravity wave drag" ;
> > udt_gwd:units = "m/s2" ;
> > udt_gwd:missing_value = -999.f ;
> > udt_gwd:_FillValue = -999.f ;
> > udt_gwd:cell_methods = "time: mean" ;
> > udt_gwd:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float vdt_gwd(time, pfull, lat, lon) ;
> > vdt_gwd:long_name = "v wind tendency for gravity wave drag" ;
> > vdt_gwd:units = "m/s2" ;
> > vdt_gwd:missing_value = -999.f ;
> > vdt_gwd:_FillValue = -999.f ;
> > vdt_gwd:cell_methods = "time: mean" ;
> > vdt_gwd:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float udt_rdamp(time, pfull, lat, lon) ;
> > udt_rdamp:long_name = "u wind tendency for Rayleigh damping" ;
> > udt_rdamp:units = "m/s2" ;
> > udt_rdamp:missing_value = -999.f ;
> > udt_rdamp:_FillValue = -999.f ;
> > udt_rdamp:cell_methods = "time: mean" ;
> > udt_rdamp:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float vdt_rdamp(time, pfull, lat, lon) ;
> > vdt_rdamp:long_name = "v wind tendency for Rayleigh damping" ;
> > vdt_rdamp:units = "m/s2" ;
> > vdt_rdamp:missing_value = -999.f ;
> > vdt_rdamp:_FillValue = -999.f ;
> > vdt_rdamp:cell_methods = "time: mean" ;
> > vdt_rdamp:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float bk(phalf) ;
> > bk:long_name = "vertical coordinate sigma value" ;
> > bk:units = "none" ;
> > bk:missing_value = 1.e+20f ;
> > bk:_FillValue = 1.e+20f ;
> > bk:cell_methods = "time: point" ;
> > float pk(phalf) ;
> > pk:long_name = "pressure part of the hybrid coordinate" ;
> > pk:units = "pascal" ;
> > pk:missing_value = 1.e+20f ;
> > pk:_FillValue = 1.e+20f ;
> > pk:cell_methods = "time: point" ;
> > float zsurf(lat, lon) ;
> > zsurf:long_name = "surface height" ;
> > zsurf:units = "m" ;
> > zsurf:missing_value = 1.e+20f ;
> > zsurf:_FillValue = 1.e+20f ;
> > zsurf:cell_methods = "time: point" ;
> > float cld_amt(time, pfull, lat, lon) ;
> > cld_amt:long_name = "cloud fraction" ;
> > cld_amt:units = "none" ;
> > cld_amt:missing_value = -1.e+10f ;
> > cld_amt:_FillValue = -1.e+10f ;
> > cld_amt:cell_methods = "time: mean" ;
> > cld_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float ice_wat(time, pfull, lat, lon) ;
> > ice_wat:long_name = "cloud ice water specific humidity" ;
> > ice_wat:units = "kg/kg" ;
> > ice_wat:missing_value = -1.e+10f ;
> > ice_wat:_FillValue = -1.e+10f ;
> > ice_wat:cell_methods = "time: mean" ;
> > ice_wat:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float liq_wat(time, pfull, lat, lon) ;
> > liq_wat:long_name = "cloud liquid specific humidity" ;
> > liq_wat:units = "kg/kg" ;
> > liq_wat:missing_value = -1.e+10f ;
> > liq_wat:_FillValue = -1.e+10f ;
> > liq_wat:cell_methods = "time: mean" ;
> > liq_wat:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float omega(time, pfull, lat, lon) ;
> > omega:long_name = "omega" ;
> > omega:units = "pa/sec" ;
> > omega:missing_value = -1.e+10f ;
> > omega:_FillValue = -1.e+10f ;
> > omega:cell_methods = "time: mean" ;
> > omega:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float ps(time, lat, lon) ;
> > ps:long_name = "surface pressure" ;
> > ps:units = "Pa" ;
> > ps:missing_value = -1.e+10f ;
> > ps:_FillValue = -1.e+10f ;
> > ps:cell_methods = "time: mean" ;
> > ps:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float sphum(time, pfull, lat, lon) ;
> > sphum:long_name = "specific humidity" ;
> > sphum:units = "kg/kg" ;
> > sphum:missing_value = -1.e+10f ;
> > sphum:_FillValue = -1.e+10f ;
> > sphum:cell_methods = "time: mean" ;
> > sphum:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float slp_dyn(time, lat, lon) ;
> > slp_dyn:long_name = "sea-level pressure" ;
> > slp_dyn:units = "mb" ;
> > slp_dyn:valid_range = 800.f, 1200.f ;
> > slp_dyn:missing_value = -1.e+10f ;
> > slp_dyn:_FillValue = -1.e+10f ;
> > slp_dyn:cell_methods = "time: mean" ;
> > slp_dyn:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float temp(time, pfull, lat, lon) ;
> > temp:long_name = "temperature" ;
> > temp:units = "deg_K" ;
> > temp:valid_range = 100.f, 350.f ;
> > temp:missing_value = -1.e+10f ;
> > temp:_FillValue = -1.e+10f ;
> > temp:cell_methods = "time: mean" ;
> > temp:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float ucomp(time, pfull, lat, lon) ;
> > ucomp:long_name = "zonal wind" ;
> > ucomp:units = "m/sec" ;
> > ucomp:valid_range = -330.f, 350.f ;
> > ucomp:missing_value = -1.e+10f ;
> > ucomp:_FillValue = -1.e+10f ;
> > ucomp:cell_methods = "time: mean" ;
> > ucomp:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float vcomp(time, pfull, lat, lon) ;
> > vcomp:long_name = "meridional wind" ;
> > vcomp:units = "m/sec" ;
> > vcomp:valid_range = -330.f, 350.f ;
> > vcomp:missing_value = -1.e+10f ;
> > vcomp:_FillValue = -1.e+10f ;
> > vcomp:cell_methods = "time: mean" ;
> > vcomp:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float drag_heat(time, lat, lon) ;
> > drag_heat:long_name = "drag coeff for heat" ;
> > drag_heat:units = "none" ;
> > drag_heat:missing_value = 1.e+20f ;
> > drag_heat:_FillValue = 1.e+20f ;
> > drag_heat:cell_methods = "time: mean" ;
> > drag_heat:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float drag_moist(time, lat, lon) ;
> > drag_moist:long_name = "drag coeff for moisture" ;
> > drag_moist:units = "none" ;
> > drag_moist:missing_value = 1.e+20f ;
> > drag_moist:_FillValue = 1.e+20f ;
> > drag_moist:cell_methods = "time: mean" ;
> > drag_moist:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float drag_mom(time, lat, lon) ;
> > drag_mom:long_name = "drag coeff for momentum" ;
> > drag_mom:units = "none" ;
> > drag_mom:missing_value = 1.e+20f ;
> > drag_mom:_FillValue = 1.e+20f ;
> > drag_mom:cell_methods = "time: mean" ;
> > drag_mom:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float evap(time, lat, lon) ;
> > evap:long_name = "evaporation rate" ;
> > evap:units = "kg/m2/s" ;
> > evap:missing_value = 1.e+20f ;
> > evap:_FillValue = 1.e+20f ;
> > evap:cell_methods = "time: mean" ;
> > evap:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float sphum_flux(time, lat, lon) ;
> > sphum_flux:long_name = "flux of specific humidity" ;
> > sphum_flux:units = "kg/kg kg air/(m2 s)" ;
> > sphum_flux:missing_value = 1.e+20f ;
> > sphum_flux:_FillValue = 1.e+20f ;
> > sphum_flux:cell_methods = "time: mean" ;
> > sphum_flux:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float ice_mask(time, lat, lon) ;
> > ice_mask:long_name = "fractional amount of sea ice" ;
> > ice_mask:units = "none" ;
> > ice_mask:valid_range = -0.01f, 1.01f ;
> > ice_mask:missing_value = 1.e+20f ;
> > ice_mask:_FillValue = 1.e+20f ;
> > ice_mask:cell_methods = "time: mean" ;
> > ice_mask:time_avg_info = "average_T1,average_T2,average_DT" ;
> > ice_mask:interp_method = "conserve_order1" ;
> > float land_mask(lat, lon) ;
> > land_mask:long_name = "fractional amount of land" ;
> > land_mask:units = "none" ;
> > land_mask:valid_range = -0.01f, 1.01f ;
> > land_mask:missing_value = 1.e+20f ;
> > land_mask:_FillValue = 1.e+20f ;
> > land_mask:cell_methods = "time: point" ;
> > land_mask:interp_method = "conserve_order1" ;
> > float lwflx(time, lat, lon) ;
> > lwflx:long_name = "net (down-up) longwave flux" ;
> > lwflx:units = "w/m2" ;
> > lwflx:missing_value = 1.e+20f ;
> > lwflx:_FillValue = 1.e+20f ;
> > lwflx:cell_methods = "time: mean" ;
> > lwflx:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float rh_ref(time, lat, lon) ;
> > rh_ref:long_name = "relative humidity at 2 m" ;
> > rh_ref:units = "percent" ;
> > rh_ref:missing_value = 1.e+20f ;
> > rh_ref:_FillValue = 1.e+20f ;
> > rh_ref:cell_methods = "time: mean" ;
> > rh_ref:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float shflx(time, lat, lon) ;
> > shflx:long_name = "sensible heat flux" ;
> > shflx:units = "w/m2" ;
> > shflx:missing_value = 1.e+20f ;
> > shflx:_FillValue = 1.e+20f ;
> > shflx:cell_methods = "time: mean" ;
> > shflx:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tau_x(time, lat, lon) ;
> > tau_x:long_name = "zonal wind stress" ;
> > tau_x:units = "pa" ;
> > tau_x:missing_value = 1.e+20f ;
> > tau_x:_FillValue = 1.e+20f ;
> > tau_x:cell_methods = "time: mean" ;
> > tau_x:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tau_y(time, lat, lon) ;
> > tau_y:long_name = "meridional wind stress" ;
> > tau_y:units = "pa" ;
> > tau_y:missing_value = 1.e+20f ;
> > tau_y:_FillValue = 1.e+20f ;
> > tau_y:cell_methods = "time: mean" ;
> > tau_y:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float t_ref(time, lat, lon) ;
> > t_ref:long_name = "temperature at 2 m" ;
> > t_ref:units = "deg_k" ;
> > t_ref:valid_range = 100.f, 400.f ;
> > t_ref:missing_value = 1.e+20f ;
> > t_ref:_FillValue = 1.e+20f ;
> > t_ref:cell_methods = "time: mean" ;
> > t_ref:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float t_ref_min(time, lat, lon) ;
> > t_ref_min:long_name = "temperature at 2 m" ;
> > t_ref_min:units = "deg_k" ;
> > t_ref_min:valid_range = 100.f, 400.f ;
> > t_ref_min:missing_value = 1.e+20f ;
> > t_ref_min:_FillValue = 1.e+20f ;
> > t_ref_min:cell_methods = "time: min" ;
> > t_ref_min:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float t_ref_max(time, lat, lon) ;
> > t_ref_max:long_name = "temperature at 2 m" ;
> > t_ref_max:units = "deg_k" ;
> > t_ref_max:valid_range = 100.f, 400.f ;
> > t_ref_max:missing_value = 1.e+20f ;
> > t_ref_max:_FillValue = 1.e+20f ;
> > t_ref_max:cell_methods = "time: max" ;
> > t_ref_max:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float t_surf(time, lat, lon) ;
> > t_surf:long_name = "surface temperature" ;
> > t_surf:units = "deg_k" ;
> > t_surf:valid_range = 100.f, 400.f ;
> > t_surf:missing_value = 1.e+20f ;
> > t_surf:_FillValue = 1.e+20f ;
> > t_surf:cell_methods = "time: mean" ;
> > t_surf:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float u_ref(time, lat, lon) ;
> > u_ref:long_name = "zonal wind component at 10 m" ;
> > u_ref:units = "m/s" ;
> > u_ref:valid_range = -400.f, 400.f ;
> > u_ref:missing_value = 1.e+20f ;
> > u_ref:_FillValue = 1.e+20f ;
> > u_ref:cell_methods = "time: mean" ;
> > u_ref:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float v_ref(time, lat, lon) ;
> > v_ref:long_name = "meridional wind component at 10 m" ;
> > v_ref:units = "m/s" ;
> > v_ref:valid_range = -400.f, 400.f ;
> > v_ref:missing_value = 1.e+20f ;
> > v_ref:_FillValue = 1.e+20f ;
> > v_ref:cell_methods = "time: mean" ;
> > v_ref:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float wind(time, lat, lon) ;
> > wind:long_name = "wind speed for flux calculations" ;
> > wind:units = "m/s" ;
> > wind:valid_range = 0.f, 400.f ;
> > wind:missing_value = 1.e+20f ;
> > wind:_FillValue = 1.e+20f ;
> > wind:cell_methods = "time: mean" ;
> > wind:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float IWP(time, lat, lon) ;
> > IWP:long_name = "Ice water path" ;
> > IWP:units = "kg/m2" ;
> > IWP:missing_value = 1.e+20f ;
> > IWP:_FillValue = 1.e+20f ;
> > IWP:cell_methods = "time: mean" ;
> > IWP:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float LWP(time, lat, lon) ;
> > LWP:long_name = "Liquid water path" ;
> > LWP:units = "kg/m2" ;
> > LWP:missing_value = 1.e+20f ;
> > LWP:_FillValue = 1.e+20f ;
> > LWP:cell_methods = "time: mean" ;
> > LWP:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float precip(time, lat, lon) ;
> > precip:long_name = "Total precipitation rate" ;
> > precip:units = "kg/m2/s" ;
> > precip:missing_value = 1.e+20f ;
> > precip:_FillValue = 1.e+20f ;
> > precip:cell_methods = "time: mean" ;
> > precip:time_avg_info = "average_T1,average_T2,average_DT" ;
> > precip:interp_method = "conserve_order1" ;
> > float prec_conv(time, lat, lon) ;
> > prec_conv:long_name = "Precipitation rate from convection" ;
> > prec_conv:units = "kg(h2o)/m2/s" ;
> > prec_conv:missing_value = 1.e+20f ;
> > prec_conv:_FillValue = 1.e+20f ;
> > prec_conv:cell_methods = "time: mean" ;
> > prec_conv:time_avg_info = "average_T1,average_T2,average_DT" ;
> > prec_conv:interp_method = "conserve_order1" ;
> > float prec_ls(time, lat, lon) ;
> > prec_ls:long_name = "Precipitation rate from strat cloud" ;
> > prec_ls:units = "kg/m2/s" ;
> > prec_ls:missing_value = 1.e+20f ;
> > prec_ls:_FillValue = 1.e+20f ;
> > prec_ls:cell_methods = "time: mean" ;
> > prec_ls:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float qdt_conv(time, pfull, lat, lon) ;
> > qdt_conv:long_name = "Spec humidity tendency from convection" ;
> > qdt_conv:units = "kg/kg/s" ;
> > qdt_conv:missing_value = -999.f ;
> > qdt_conv:_FillValue = -999.f ;
> > qdt_conv:cell_methods = "time: mean" ;
> > qdt_conv:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float qdt_ls(time, pfull, lat, lon) ;
> > qdt_ls:long_name = "Spec humidity tendency from strat cloud" ;
> > qdt_ls:units = "kg/kg/s" ;
> > qdt_ls:missing_value = -999.f ;
> > qdt_ls:_FillValue = -999.f ;
> > qdt_ls:cell_methods = "time: mean" ;
> > qdt_ls:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float rh(time, pfull, lat, lon) ;
> > rh:long_name = "relative humidity" ;
> > rh:units = "percent" ;
> > rh:missing_value = -999.f ;
> > rh:_FillValue = -999.f ;
> > rh:cell_methods = "time: mean" ;
> > rh:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float snow_conv(time, lat, lon) ;
> > snow_conv:long_name = "Frozen precip rate from convection" ;
> > snow_conv:units = "kg(h2o)/m2/s" ;
> > snow_conv:missing_value = 1.e+20f ;
> > snow_conv:_FillValue = 1.e+20f ;
> > snow_conv:cell_methods = "time: mean" ;
> > snow_conv:time_avg_info = "average_T1,average_T2,average_DT" ;
> > snow_conv:interp_method = "conserve_order1" ;
> > float snow_ls(time, lat, lon) ;
> > snow_ls:long_name = "Frozen precip rate from strat cloud" ;
> > snow_ls:units = "kg/m2/s" ;
> > snow_ls:missing_value = 1.e+20f ;
> > snow_ls:_FillValue = 1.e+20f ;
> > snow_ls:cell_methods = "time: mean" ;
> > snow_ls:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_conv(time, pfull, lat, lon) ;
> > tdt_conv:long_name = "Temperature tendency from convection" ;
> > tdt_conv:units = "deg_K/s" ;
> > tdt_conv:missing_value = -999.f ;
> > tdt_conv:_FillValue = -999.f ;
> > tdt_conv:cell_methods = "time: mean" ;
> > tdt_conv:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_ls(time, pfull, lat, lon) ;
> > tdt_ls:long_name = "Temperature tendency from strat cloud" ;
> > tdt_ls:units = "deg_K/s" ;
> > tdt_ls:missing_value = -999.f ;
> > tdt_ls:_FillValue = -999.f ;
> > tdt_ls:cell_methods = "time: mean" ;
> > tdt_ls:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float WVP(time, lat, lon) ;
> > WVP:long_name = "Column integrated water vapor" ;
> > WVP:units = "kg/m2" ;
> > WVP:missing_value = 1.e+20f ;
> > WVP:_FillValue = 1.e+20f ;
> > WVP:cell_methods = "time: mean" ;
> > WVP:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float alb_sfc(time, lat, lon) ;
> > alb_sfc:long_name = "surface albedo" ;
> > alb_sfc:units = "percent" ;
> > alb_sfc:missing_value = -999.f ;
> > alb_sfc:_FillValue = -999.f ;
> > alb_sfc:cell_methods = "time: mean" ;
> > alb_sfc:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lwdn_sfc(time, lat, lon) ;
> > lwdn_sfc:long_name = "LW flux down at surface" ;
> > lwdn_sfc:units = "watts/m2" ;
> > lwdn_sfc:missing_value = -999.f ;
> > lwdn_sfc:_FillValue = -999.f ;
> > lwdn_sfc:cell_methods = "time: mean" ;
> > lwdn_sfc:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lwdn_sfc_clr(time, lat, lon) ;
> > lwdn_sfc_clr:long_name = "clear skyLW flux down at surface" ;
> > lwdn_sfc_clr:units = "watts/m2" ;
> > lwdn_sfc_clr:missing_value = -999.f ;
> > lwdn_sfc_clr:_FillValue = -999.f ;
> > lwdn_sfc_clr:cell_methods = "time: mean" ;
> > lwdn_sfc_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lwup_sfc(time, lat, lon) ;
> > lwup_sfc:long_name = "LW flux up at surface" ;
> > lwup_sfc:units = "watts/m2" ;
> > lwup_sfc:missing_value = -999.f ;
> > lwup_sfc:_FillValue = -999.f ;
> > lwup_sfc:cell_methods = "time: mean" ;
> > lwup_sfc:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lwup_sfc_clr(time, lat, lon) ;
> > lwup_sfc_clr:long_name = "clear skyLW flux up at surface" ;
> > lwup_sfc_clr:units = "watts/m2" ;
> > lwup_sfc_clr:missing_value = -999.f ;
> > lwup_sfc_clr:_FillValue = -999.f ;
> > lwup_sfc_clr:cell_methods = "time: mean" ;
> > lwup_sfc_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float netrad_toa(time, lat, lon) ;
> > netrad_toa:long_name = "net radiation (lw + sw) at toa" ;
> > netrad_toa:units = "watts/m2" ;
> > netrad_toa:missing_value = -999.f ;
> > netrad_toa:_FillValue = -999.f ;
> > netrad_toa:cell_methods = "time: mean" ;
> > netrad_toa:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float netrad_toa_clr(time, lat, lon) ;
> > netrad_toa_clr:long_name = "clear skynet radiation (lw + sw) at toa" ;
> > netrad_toa_clr:units = "watts/m2" ;
> > netrad_toa_clr:missing_value = -999.f ;
> > netrad_toa_clr:_FillValue = -999.f ;
> > netrad_toa_clr:cell_methods = "time: mean" ;
> > netrad_toa_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float olr(time, lat, lon) ;
> > olr:long_name = "outgoing longwave radiation" ;
> > olr:units = "watts/m2" ;
> > olr:missing_value = -999.f ;
> > olr:_FillValue = -999.f ;
> > olr:cell_methods = "time: mean" ;
> > olr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float olr_clr(time, lat, lon) ;
> > olr_clr:long_name = "clear skyoutgoing longwave radiation" ;
> > olr_clr:units = "watts/m2" ;
> > olr_clr:missing_value = -999.f ;
> > olr_clr:_FillValue = -999.f ;
> > olr_clr:cell_methods = "time: mean" ;
> > olr_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float qo3(time, pfull, lat, lon) ;
> > qo3:long_name = "ozone mixing ratio" ;
> > qo3:units = "kg/kg" ;
> > qo3:missing_value = -999.f ;
> > qo3:_FillValue = -999.f ;
> > qo3:cell_methods = "time: mean" ;
> > qo3:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float qo3_col(time, lat, lon) ;
> > qo3_col:long_name = "ozone column" ;
> > qo3_col:units = "DU" ;
> > qo3_col:missing_value = -999.f ;
> > qo3_col:_FillValue = -999.f ;
> > qo3_col:cell_methods = "time: mean" ;
> > qo3_col:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_sfc(time, lat, lon) ;
> > swdn_sfc:long_name = "SW flux down at surface" ;
> > swdn_sfc:units = "watts/m2" ;
> > swdn_sfc:missing_value = -999.f ;
> > swdn_sfc:_FillValue = -999.f ;
> > swdn_sfc:cell_methods = "time: mean" ;
> > swdn_sfc:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_sfc_clr(time, lat, lon) ;
> > swdn_sfc_clr:long_name = "clear skySW flux down at surface" ;
> > swdn_sfc_clr:units = "watts/m2" ;
> > swdn_sfc_clr:missing_value = -999.f ;
> > swdn_sfc_clr:_FillValue = -999.f ;
> > swdn_sfc_clr:cell_methods = "time: mean" ;
> > swdn_sfc_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_sfc(time, lat, lon) ;
> > swup_sfc:long_name = "SW flux up at surface" ;
> > swup_sfc:units = "watts/m2" ;
> > swup_sfc:missing_value = -999.f ;
> > swup_sfc:_FillValue = -999.f ;
> > swup_sfc:cell_methods = "time: mean" ;
> > swup_sfc:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_sfc_clr(time, lat, lon) ;
> > swup_sfc_clr:long_name = "clear skySW flux up at surface" ;
> > swup_sfc_clr:units = "watts/m2" ;
> > swup_sfc_clr:missing_value = -999.f ;
> > swup_sfc_clr:_FillValue = -999.f ;
> > swup_sfc_clr:cell_methods = "time: mean" ;
> > swup_sfc_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_toa(time, lat, lon) ;
> > swdn_toa:long_name = "SW flux down at TOA" ;
> > swdn_toa:units = "watts/m2" ;
> > swdn_toa:missing_value = -999.f ;
> > swdn_toa:_FillValue = -999.f ;
> > swdn_toa:cell_methods = "time: mean" ;
> > swdn_toa:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_toa_clr(time, lat, lon) ;
> > swdn_toa_clr:long_name = "clear skySW flux down at TOA" ;
> > swdn_toa_clr:units = "watts/m2" ;
> > swdn_toa_clr:missing_value = -999.f ;
> > swdn_toa_clr:_FillValue = -999.f ;
> > swdn_toa_clr:cell_methods = "time: mean" ;
> > swdn_toa_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_toa(time, lat, lon) ;
> > swup_toa:long_name = "SW flux up at TOA" ;
> > swup_toa:units = "watts/m2" ;
> > swup_toa:missing_value = -999.f ;
> > swup_toa:_FillValue = -999.f ;
> > swup_toa:cell_methods = "time: mean" ;
> > swup_toa:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_toa_clr(time, lat, lon) ;
> > swup_toa_clr:long_name = "clear skySW flux up at TOA" ;
> > swup_toa_clr:units = "watts/m2" ;
> > swup_toa_clr:missing_value = -999.f ;
> > swup_toa_clr:_FillValue = -999.f ;
> > swup_toa_clr:cell_methods = "time: mean" ;
> > swup_toa_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_lw(time, pfull, lat, lon) ;
> > tdt_lw:long_name = "temperature tendency for LW radiation" ;
> > tdt_lw:units = "deg_K/sec" ;
> > tdt_lw:missing_value = -999.f ;
> > tdt_lw:_FillValue = -999.f ;
> > tdt_lw:cell_methods = "time: mean" ;
> > tdt_lw:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_lw_clr(time, pfull, lat, lon) ;
> > tdt_lw_clr:long_name = "clear skytemperature tendency for LW radiation" ;
> > tdt_lw_clr:units = "deg_K/sec" ;
> > tdt_lw_clr:missing_value = -999.f ;
> > tdt_lw_clr:_FillValue = -999.f ;
> > tdt_lw_clr:cell_methods = "time: mean" ;
> > tdt_lw_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_sw(time, pfull, lat, lon) ;
> > tdt_sw:long_name = "temperature tendency for SW radiation" ;
> > tdt_sw:units = "deg_K/sec" ;
> > tdt_sw:missing_value = -999.f ;
> > tdt_sw:_FillValue = -999.f ;
> > tdt_sw:cell_methods = "time: mean" ;
> > tdt_sw:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_sw_clr(time, pfull, lat, lon) ;
> > tdt_sw_clr:long_name = "clear skytemperature tendency for SW radiation" ;
> > tdt_sw_clr:units = "deg_K/sec" ;
> > tdt_sw_clr:missing_value = -999.f ;
> > tdt_sw_clr:_FillValue = -999.f ;
> > tdt_sw_clr:cell_methods = "time: mean" ;
> > tdt_sw_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float netlw_lin_trop(time, lat, lon) ;
> > netlw_lin_trop:long_name = "net LW flux at linear tropopause" ;
> > netlw_lin_trop:units = "watts/m2" ;
> > netlw_lin_trop:missing_value = -999.f ;
> > netlw_lin_trop:_FillValue = -999.f ;
> > netlw_lin_trop:cell_methods = "time: mean" ;
> > netlw_lin_trop:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float netlw_lin_trop_clr(time, lat, lon) ;
> > netlw_lin_trop_clr:long_name = "clear skynet LW flux at linear
> tropopause"
> > ;
> > netlw_lin_trop_clr:units = "watts/m2" ;
> > netlw_lin_trop_clr:missing_value = -999.f ;
> > netlw_lin_trop_clr:_FillValue = -999.f ;
> > netlw_lin_trop_clr:cell_methods = "time: mean" ;
> > netlw_lin_trop_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float netlw_200hPa(time, lat, lon) ;
> > netlw_200hPa:long_name = "net LW flux at 200 hPa" ;
> > netlw_200hPa:units = "watts/m2" ;
> > netlw_200hPa:missing_value = -999.f ;
> > netlw_200hPa:_FillValue = -999.f ;
> > netlw_200hPa:cell_methods = "time: mean" ;
> > netlw_200hPa:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float netlw_200hPa_clr(time, lat, lon) ;
> > netlw_200hPa_clr:long_name = "clear skynet LW flux at 200 hPa" ;
> > netlw_200hPa_clr:units = "watts/m2" ;
> > netlw_200hPa_clr:missing_value = -999.f ;
> > netlw_200hPa_clr:_FillValue = -999.f ;
> > netlw_200hPa_clr:cell_methods = "time: mean" ;
> > netlw_200hPa_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_lin_trop(time, lat, lon) ;
> > swdn_lin_trop:long_name = "SW flux down at linear tropopause" ;
> > swdn_lin_trop:units = "watts/m2" ;
> > swdn_lin_trop:missing_value = -999.f ;
> > swdn_lin_trop:_FillValue = -999.f ;
> > swdn_lin_trop:cell_methods = "time: mean" ;
> > swdn_lin_trop:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_lin_trop_clr(time, lat, lon) ;
> > swdn_lin_trop_clr:long_name = "clear skySW flux down at linear
> tropopause"
> > ;
> > swdn_lin_trop_clr:units = "watts/m2" ;
> > swdn_lin_trop_clr:missing_value = -999.f ;
> > swdn_lin_trop_clr:_FillValue = -999.f ;
> > swdn_lin_trop_clr:cell_methods = "time: mean" ;
> > swdn_lin_trop_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_200hPa(time, lat, lon) ;
> > swdn_200hPa:long_name = "SW flux down at 200 hPa" ;
> > swdn_200hPa:units = "watts/m2" ;
> > swdn_200hPa:missing_value = -999.f ;
> > swdn_200hPa:_FillValue = -999.f ;
> > swdn_200hPa:cell_methods = "time: mean" ;
> > swdn_200hPa:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_200hPa_clr(time, lat, lon) ;
> > swdn_200hPa_clr:long_name = "clear skySW flux down at 200 hPa" ;
> > swdn_200hPa_clr:units = "watts/m2" ;
> > swdn_200hPa_clr:missing_value = -999.f ;
> > swdn_200hPa_clr:_FillValue = -999.f ;
> > swdn_200hPa_clr:cell_methods = "time: mean" ;
> > swdn_200hPa_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_lin_trop(time, lat, lon) ;
> > swup_lin_trop:long_name = "SW flux up at linear tropopause" ;
> > swup_lin_trop:units = "watts/m2" ;
> > swup_lin_trop:missing_value = -999.f ;
> > swup_lin_trop:_FillValue = -999.f ;
> > swup_lin_trop:cell_methods = "time: mean" ;
> > swup_lin_trop:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_lin_trop_clr(time, lat, lon) ;
> > swup_lin_trop_clr:long_name = "clear skySW flux up at linear tropopause"
> ;
> > swup_lin_trop_clr:units = "watts/m2" ;
> > swup_lin_trop_clr:missing_value = -999.f ;
> > swup_lin_trop_clr:_FillValue = -999.f ;
> > swup_lin_trop_clr:cell_methods = "time: mean" ;
> > swup_lin_trop_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_200hPa(time, lat, lon) ;
> > swup_200hPa:long_name = "SW flux up at 200 hPa" ;
> > swup_200hPa:units = "watts/m2" ;
> > swup_200hPa:missing_value = -999.f ;
> > swup_200hPa:_FillValue = -999.f ;
> > swup_200hPa:cell_methods = "time: mean" ;
> > swup_200hPa:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swup_200hPa_clr(time, lat, lon) ;
> > swup_200hPa_clr:long_name = "clear skySW flux up at 200 hPa" ;
> > swup_200hPa_clr:units = "watts/m2" ;
> > swup_200hPa_clr:missing_value = -999.f ;
> > swup_200hPa_clr:_FillValue = -999.f ;
> > swup_200hPa_clr:cell_methods = "time: mean" ;
> > swup_200hPa_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float mc(time, phalf, lat, lon) ;
> > mc:long_name = "Cumulus Mass Flux from RAS" ;
> > mc:units = "kg/m2/s" ;
> > mc:missing_value = -999.f ;
> > mc:_FillValue = -999.f ;
> > mc:cell_methods = "time: mean" ;
> > mc:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float qdt_vdif(time, pfull, lat, lon) ;
> > qdt_vdif:long_name = "Spec humidity tendency from vert diff" ;
> > qdt_vdif:units = "kg/kg/s" ;
> > qdt_vdif:missing_value = -999.f ;
> > qdt_vdif:_FillValue = -999.f ;
> > qdt_vdif:cell_methods = "time: mean" ;
> > qdt_vdif:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tdt_vdif(time, pfull, lat, lon) ;
> > tdt_vdif:long_name = "Temperature tendency from vert diff" ;
> > tdt_vdif:units = "deg_K/s" ;
> > tdt_vdif:missing_value = -999.f ;
> > tdt_vdif:_FillValue = -999.f ;
> > tdt_vdif:cell_methods = "time: mean" ;
> > tdt_vdif:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float udt_vdif(time, pfull, lat, lon) ;
> > udt_vdif:long_name = "Zonal wind tendency from vert diff" ;
> > udt_vdif:units = "m/s2" ;
> > udt_vdif:missing_value = -999.f ;
> > udt_vdif:_FillValue = -999.f ;
> > udt_vdif:cell_methods = "time: mean" ;
> > udt_vdif:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float vdt_vdif(time, pfull, lat, lon) ;
> > vdt_vdif:long_name = "Meridional wind tendency from vert diff" ;
> > vdt_vdif:units = "m/s2" ;
> > vdt_vdif:missing_value = -999.f ;
> > vdt_vdif:_FillValue = -999.f ;
> > vdt_vdif:cell_methods = "time: mean" ;
> > vdt_vdif:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float co2(time, pfull, lat, lon) ;
> > co2:long_name = "carbon dioxide" ;
> > co2:units = "kg/kg" ;
> > co2:missing_value = -1.e+10f ;
> > co2:_FillValue = -1.e+10f ;
> > co2:cell_methods = "time: mean" ;
> > co2:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float co2_dvmr(time, pfull, lat, lon) ;
> > co2_dvmr:long_name = "carbon dioxide (dry vmr)" ;
> > co2_dvmr:units = "mol/mol" ;
> > co2_dvmr:missing_value = -1.e+10f ;
> > co2_dvmr:_FillValue = -1.e+10f ;
> > co2_dvmr:cell_methods = "time: mean" ;
> > co2_dvmr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float slp(time, lat, lon) ;
> > slp:long_name = "sea level pressure" ;
> > slp:units = "pa" ;
> > slp:missing_value = 1.e+20f ;
> > slp:_FillValue = 1.e+20f ;
> > slp:cell_methods = "time: mean" ;
> > slp:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float wind_ref(time, lat, lon) ;
> > wind_ref:long_name = "absolute value of wind at 10 m" ;
> > wind_ref:units = "m/s" ;
> > wind_ref:valid_range = -400.f, 400.f ;
> > wind_ref:missing_value = 1.e+20f ;
> > wind_ref:_FillValue = 1.e+20f ;
> > wind_ref:cell_methods = "time: mean" ;
> > wind_ref:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float rh_ref_cmip(time, lat, lon) ;
> > rh_ref_cmip:long_name = "relative humidity at 2 m" ;
> > rh_ref_cmip:units = "percent" ;
> > rh_ref_cmip:missing_value = 1.e+20f ;
> > rh_ref_cmip:_FillValue = 1.e+20f ;
> > rh_ref_cmip:cell_methods = "time: mean" ;
> > rh_ref_cmip:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float q_ref(time, lat, lon) ;
> > q_ref:long_name = "specific humidity at 2 m" ;
> > q_ref:units = "kg/kg" ;
> > q_ref:missing_value = -1.f ;
> > q_ref:_FillValue = -1.f ;
> > q_ref:cell_methods = "time: mean" ;
> > q_ref:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float snow_tot(time, lat, lon) ;
> > snow_tot:long_name = "Frozen precip rate from all sources" ;
> > snow_tot:units = "kg(h2o)/m2/s" ;
> > snow_tot:missing_value = 1.e+20f ;
> > snow_tot:_FillValue = 1.e+20f ;
> > snow_tot:cell_methods = "time: mean" ;
> > snow_tot:time_avg_info = "average_T1,average_T2,average_DT" ;
> > snow_tot:interp_method = "conserve_order1" ;
> > float cld_amt_2d(time, lat, lon) ;
> > cld_amt_2d:long_name = "total cloud amount" ;
> > cld_amt_2d:units = "percent" ;
> > cld_amt_2d:missing_value = 1.e+20f ;
> > cld_amt_2d:_FillValue = 1.e+20f ;
> > cld_amt_2d:cell_methods = "time: mean" ;
> > cld_amt_2d:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float WP_all_clouds(time, lat, lon) ;
> > WP_all_clouds:long_name = "Total water path -- all clouds" ;
> > WP_all_clouds:units = "kg/m2" ;
> > WP_all_clouds:missing_value = 1.e+20f ;
> > WP_all_clouds:_FillValue = 1.e+20f ;
> > WP_all_clouds:cell_methods = "time: mean" ;
> > WP_all_clouds:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float IWP_all_clouds(time, lat, lon) ;
> > IWP_all_clouds:long_name = "Ice water path -- all clouds" ;
> > IWP_all_clouds:units = "kg/m2" ;
> > IWP_all_clouds:missing_value = 1.e+20f ;
> > IWP_all_clouds:_FillValue = 1.e+20f ;
> > IWP_all_clouds:cell_methods = "time: mean" ;
> > IWP_all_clouds:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float conv_cld_base(time, lat, lon) ;
> > conv_cld_base:long_name = "pressure at convective cloud base" ;
> > conv_cld_base:units = "Pa" ;
> > conv_cld_base:missing_value = -999.f ;
> > conv_cld_base:_FillValue = -999.f ;
> > conv_cld_base:cell_methods = "time: mean" ;
> > conv_cld_base:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float conv_cld_top(time, lat, lon) ;
> > conv_cld_top:long_name = "pressure at convective cloud top" ;
> > conv_cld_top:units = "Pa" ;
> > conv_cld_top:missing_value = -999.f ;
> > conv_cld_top:_FillValue = -999.f ;
> > conv_cld_top:cell_methods = "time: mean" ;
> > conv_cld_top:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float conv_freq(time, lat, lon) ;
> > conv_freq:long_name = "frequency of convection" ;
> > conv_freq:units = "number" ;
> > conv_freq:missing_value = -999.f ;
> > conv_freq:_FillValue = -999.f ;
> > conv_freq:cell_methods = "time: mean" ;
> > conv_freq:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tot_cloud_area(time, pfull, lat, lon) ;
> > tot_cloud_area:long_name = "Cloud area -- all clouds" ;
> > tot_cloud_area:units = "percent" ;
> > tot_cloud_area:missing_value = -999.f ;
> > tot_cloud_area:_FillValue = -999.f ;
> > tot_cloud_area:cell_methods = "time: mean" ;
> > tot_cloud_area:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tot_liq_amt(time, pfull, lat, lon) ;
> > tot_liq_amt:long_name = "Liquid amount -- all clouds" ;
> > tot_liq_amt:units = "kg/kg" ;
> > tot_liq_amt:missing_value = -999.f ;
> > tot_liq_amt:_FillValue = -999.f ;
> > tot_liq_amt:cell_methods = "time: mean" ;
> > tot_liq_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float tot_ice_amt(time, pfull, lat, lon) ;
> > tot_ice_amt:long_name = "Ice amount -- all clouds" ;
> > tot_ice_amt:units = "kg/kg" ;
> > tot_ice_amt:missing_value = -999.f ;
> > tot_ice_amt:_FillValue = -999.f ;
> > tot_ice_amt:cell_methods = "time: mean" ;
> > tot_ice_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float rh_cmip(time, pfull, lat, lon) ;
> > rh_cmip:long_name = "relative humidity" ;
> > rh_cmip:units = "percent" ;
> > rh_cmip:missing_value = -999.f ;
> > rh_cmip:_FillValue = -999.f ;
> > rh_cmip:cell_methods = "time: mean" ;
> > rh_cmip:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float qo3v(time, pfull, lat, lon) ;
> > qo3v:long_name = "ozone mole fraction" ;
> > qo3v:units = "1.e-9" ;
> > qo3v:missing_value = -999.f ;
> > qo3v:_FillValue = -999.f ;
> > qo3v:cell_methods = "time: mean" ;
> > qo3v:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float q_star(time, lat, lon) ;
> > q_star:long_name = "moisture scale" ;
> > q_star:units = "kg water/kg air" ;
> > q_star:missing_value = 1.e+20f ;
> > q_star:_FillValue = 1.e+20f ;
> > q_star:cell_methods = "time: mean" ;
> > q_star:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float rh_min(time, pfull, lat, lon) ;
> > rh_min:long_name = "relative humidity" ;
> > rh_min:units = "percent" ;
> > rh_min:missing_value = -999.f ;
> > rh_min:_FillValue = -999.f ;
> > rh_min:cell_methods = "time: min" ;
> > rh_min:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float rh_max(time, pfull, lat, lon) ;
> > rh_max:long_name = "relative humidity" ;
> > rh_max:units = "percent" ;
> > rh_max:missing_value = -999.f ;
> > rh_max:_FillValue = -999.f ;
> > rh_max:cell_methods = "time: max" ;
> > rh_max:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float droplets(time, pfull, lat, lon) ;
> > droplets:long_name = "Droplet number concentration" ;
> > droplets:units = "/cm3" ;
> > droplets:missing_value = -999.f ;
> > droplets:_FillValue = -999.f ;
> > droplets:cell_methods = "time: mean" ;
> > droplets:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float enth_ls_col(time, lat, lon) ;
> > enth_ls_col:long_name = "Column enthalpy tendency from strat cloud" ;
> > enth_ls_col:units = "W/m2" ;
> > enth_ls_col:missing_value = 1.e+20f ;
> > enth_ls_col:_FillValue = 1.e+20f ;
> > enth_ls_col:cell_methods = "time: mean" ;
> > enth_ls_col:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float wat_ls_col(time, lat, lon) ;
> > wat_ls_col:long_name = "Column total water tendency from strat cloud" ;
> > wat_ls_col:units = "kg/m2/s" ;
> > wat_ls_col:missing_value = 1.e+20f ;
> > wat_ls_col:_FillValue = 1.e+20f ;
> > wat_ls_col:cell_methods = "time: mean" ;
> > wat_ls_col:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float enth_conv_col(time, lat, lon) ;
> > enth_conv_col:long_name = "Column enthalpy tendency from convection" ;
> > enth_conv_col:units = "W/m2" ;
> > enth_conv_col:missing_value = 1.e+20f ;
> > enth_conv_col:_FillValue = 1.e+20f ;
> > enth_conv_col:cell_methods = "time: mean" ;
> > enth_conv_col:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float wat_conv_col(time, lat, lon) ;
> > wat_conv_col:long_name = "Column total water tendency from convection" ;
> > wat_conv_col:units = "kg(h2o)/m2/s" ;
> > wat_conv_col:missing_value = 1.e+20f ;
> > wat_conv_col:_FillValue = 1.e+20f ;
> > wat_conv_col:cell_methods = "time: mean" ;
> > wat_conv_col:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_vis_dir(time, lat, lon) ;
> > swdn_vis_dir:long_name = "downward direct visible sfc sw flux" ;
> > swdn_vis_dir:units = "watts/m2" ;
> > swdn_vis_dir:missing_value = -999.f ;
> > swdn_vis_dir:_FillValue = -999.f ;
> > swdn_vis_dir:cell_methods = "time: mean" ;
> > swdn_vis_dir:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_vis_dif(time, lat, lon) ;
> > swdn_vis_dif:long_name = "downward diffuse visible sfc sw flux" ;
> > swdn_vis_dif:units = "watts/m2" ;
> > swdn_vis_dif:missing_value = -999.f ;
> > swdn_vis_dif:_FillValue = -999.f ;
> > swdn_vis_dif:cell_methods = "time: mean" ;
> > swdn_vis_dif:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_vis_clr(time, lat, lon) ;
> > swdn_vis_clr:long_name = "downward visible sfc sw flux clear sky" ;
> > swdn_vis_clr:units = "watts/m2" ;
> > swdn_vis_clr:missing_value = -999.f ;
> > swdn_vis_clr:_FillValue = -999.f ;
> > swdn_vis_clr:cell_methods = "time: mean" ;
> > swdn_vis_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_tot_dir(time, lat, lon) ;
> > swdn_tot_dir:long_name = "downward direct total sfc sw flux" ;
> > swdn_tot_dir:units = "watts/m2" ;
> > swdn_tot_dir:missing_value = -999.f ;
> > swdn_tot_dir:_FillValue = -999.f ;
> > swdn_tot_dir:cell_methods = "time: mean" ;
> > swdn_tot_dir:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float swdn_tot_dif(time, lat, lon) ;
> > swdn_tot_dif:long_name = "downward diffuse total sfc sw flux" ;
> > swdn_tot_dif:units = "watts/m2" ;
> > swdn_tot_dif:missing_value = -999.f ;
> > swdn_tot_dif:_FillValue = -999.f ;
> > swdn_tot_dif:cell_methods = "time: mean" ;
> > swdn_tot_dif:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float flux_sw_vis(time, lat, lon) ;
> > flux_sw_vis:long_name = "net visible sfc sw flux" ;
> > flux_sw_vis:units = "watts/m2" ;
> > flux_sw_vis:missing_value = -999.f ;
> > flux_sw_vis:_FillValue = -999.f ;
> > flux_sw_vis:cell_methods = "time: mean" ;
> > flux_sw_vis:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float sw_dn_t_dir_clr(time, lat, lon) ;
> > sw_dn_t_dir_clr:long_name = "downward clearsky direct total sfc sw flux"
> ;
> > sw_dn_t_dir_clr:units = "watts/m2" ;
> > sw_dn_t_dir_clr:missing_value = -999.f ;
> > sw_dn_t_dir_clr:_FillValue = -999.f ;
> > sw_dn_t_dir_clr:cell_methods = "time: mean" ;
> > sw_dn_t_dir_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float sw_dn_t_dif_clr(time, lat, lon) ;
> > sw_dn_t_dif_clr:long_name = "downward clearsky diffuse total sfc sw
> flux" ;
> > sw_dn_t_dif_clr:units = "watts/m2" ;
> > sw_dn_t_dif_clr:missing_value = -999.f ;
> > sw_dn_t_dif_clr:_FillValue = -999.f ;
> > sw_dn_t_dif_clr:cell_methods = "time: mean" ;
> > sw_dn_t_dif_clr:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float LWP2(time, lat, lon) ;
> > LWP2:long_name = "In-cloud liquid water path of stratiform clouds" ;
> > LWP2:units = "kg/m2" ;
> > LWP2:missing_value = -999.f ;
> > LWP2:_FillValue = -999.f ;
> > LWP2:cell_methods = "time: mean" ;
> > LWP2:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float z_pbl(time, lat, lon) ;
> > z_pbl:long_name = "depth of planetary boundary layer" ;
> > z_pbl:units = "m" ;
> > z_pbl:missing_value = 1.e+20f ;
> > z_pbl:_FillValue = 1.e+20f ;
> > z_pbl:cell_methods = "time: mean" ;
> > z_pbl:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float conv_cloud_area(time, pfull, lat, lon) ;
> > conv_cloud_area:long_name = "Convective cloud area" ;
> > conv_cloud_area:units = "percent" ;
> > conv_cloud_area:missing_value = -999.f ;
> > conv_cloud_area:_FillValue = -999.f ;
> > conv_cloud_area:cell_methods = "time: mean" ;
> > conv_cloud_area:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float conv_liq_amt(time, pfull, lat, lon) ;
> > conv_liq_amt:long_name = "Convective cloud liquid amount" ;
> > conv_liq_amt:units = "kg/kg" ;
> > conv_liq_amt:missing_value = -999.f ;
> > conv_liq_amt:_FillValue = -999.f ;
> > conv_liq_amt:cell_methods = "time: mean" ;
> > conv_liq_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float conv_ice_amt(time, pfull, lat, lon) ;
> > conv_ice_amt:long_name = "Convective cloud ice amount" ;
> > conv_ice_amt:units = "kg/kg" ;
> > conv_ice_amt:missing_value = -999.f ;
> > conv_ice_amt:_FillValue = -999.f ;
> > conv_ice_amt:cell_methods = "time: mean" ;
> > conv_ice_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lsc_cloud_area(time, pfull, lat, lon) ;
> > lsc_cloud_area:long_name = "Large-scale cloud area" ;
> > lsc_cloud_area:units = "percent" ;
> > lsc_cloud_area:missing_value = -999.f ;
> > lsc_cloud_area:_FillValue = -999.f ;
> > lsc_cloud_area:cell_methods = "time: mean" ;
> > lsc_cloud_area:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lsc_liq_amt(time, pfull, lat, lon) ;
> > lsc_liq_amt:long_name = "Large-scale cloud liquid amount" ;
> > lsc_liq_amt:units = "kg/kg" ;
> > lsc_liq_amt:missing_value = -999.f ;
> > lsc_liq_amt:_FillValue = -999.f ;
> > lsc_liq_amt:cell_methods = "time: mean" ;
> > lsc_liq_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > float lsc_ice_amt(time, pfull, lat, lon) ;
> > lsc_ice_amt:long_name = "Large-scale cloud ice amount" ;
> > lsc_ice_amt:units = "kg/kg" ;
> > lsc_ice_amt:missing_value = -999.f ;
> > lsc_ice_amt:_FillValue = -999.f ;
> > lsc_ice_amt:cell_methods = "time: mean" ;
> > lsc_ice_amt:time_avg_info = "average_T1,average_T2,average_DT" ;
> > double average_T1(time) ;
> > average_T1:long_name = "Start time for average period" ;
> > average_T1:units = "days since 0001-01-01 00:00:00" ;
> > average_T1:missing_value = 1.e+20 ;
> > average_T1:_FillValue = 1.e+20 ;
> > double average_T2(time) ;
> > average_T2:long_name = "End time for average period" ;
> > average_T2:units = "days since 0001-01-01 00:00:00" ;
> > average_T2:missing_value = 1.e+20 ;
> > average_T2:_FillValue = 1.e+20 ;
> > double average_DT(time) ;
> > average_DT:long_name = "Length of average period" ;
> > average_DT:units = "days" ;
> > average_DT:missing_value = 1.e+20 ;
> > average_DT:_FillValue = 1.e+20 ;
> > double time_bounds(time, nv) ;
> > time_bounds:long_name = "time axis boundaries" ;
> > time_bounds:units = "days" ;
> > time_bounds:missing_value = 1.e+20 ;
> > time_bounds:_FillValue = 1.e+20 ;
> >
> > // global attributes:
> > :filename = "atmos_month.nc.0000" ;
> > :NumFilesInSet = 30 ;
> > :title = "ESM2M_pi-control_C2" ;
> > :grid_type = "regular" ;
> > :grid_tile = "N/A" ;
> > data:
> >
> > lon = _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _, _ ;
> >
> > lonb = _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> > _,
> > _, _ ;
> >
> > lat = _, _, _ ;
> >
> > latb = _, _, _, _ ;
> >
> > nv = _, _ ;
> >
> > pfull = _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _,
> > _, _ ;
> >
> > phalf = _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _,
> > _,
> > _, _ ;
> > }
> >
> > Thanks,
> >
> > David
> >
> >
> > On Thu, Jun 19, 2014 at 1:41 PM, Jennifer Adams <jma at cola.iges.org>
> wrote:
> >
> >> David,
> >> Can you post the output from ncdump -c on your netcdf file?
> >> --Jennifer
> >>
> >> On Jun 19, 2014, at 1:14 PM, David Stanfel wrote:
> >>
> >> I made the changes that were recommended with the following results. My
> >> X11 window does pop up when I run grads but it stays black. This is what
> >> happens when I try to display a variable:
> >>
> >> ga-> open cont.ctl
> >> Scanning description file: cont.ctl
> >> Data file atmos_month.nc.0000 is open as file 1
> >> LON set to 0 360
> >> LAT set to -88.29 88.29
> >> LEV set to 1000 1000
> >> Time values set: 0:1:1:0 0:1:1:0
> >> E set to 1 1
> >> ga-> q file
> >> File 1 : 6 Months Sample Model Output
> >> Descriptor: cont.ctl
> >> Binary: atmos_month.nc.0000
> >> Type = Gridded
> >> Xsize = 96 Ysize = 80 Zsize = 18 Tsize = 804 Esize = 1
> >> Number of Variables = 7
> >> u 18 t,z,y,x u
> >> v 18 t,z,y,x v
> >> w 18 t,z,y,x v
> >> t 18 t,z,y,x v
> >> sh 18 t,z,y,x v
> >> rh 18 t,z,y,x v
> >> z 18 t,z,y,x v
> >> ga-> d u
> >> Error: Variable u not in netcdf file
> >> Data Request Error: Error for variable 'u'
> >> Error ocurred at column 1
> >> DISPLAY error: Invalid expression
> >> Expression = u
> >> ga-> d z
> >> Error: Variable z not in netcdf file
> >> Data Request Error: Error for variable 'z'
> >> Error ocurred at column 1
> >> DISPLAY error: Invalid expression
> >> Expression = z
> >>
> >> It still does not recognize that that variable is in the netcdf file for
> >> some reason.
> >>
> >> Thanks,
> >>
> >> David
> >>
> >>
> >> On Thu, Jun 19, 2014 at 10:22 AM, Hosmay Lopez <hlopez at rsmas.miami.edu>
> >> wrote:
> >>
> >>> Try this CTL file:
> >>>
> >>> The only difference is how you defined the variables' arguments.
> >>>
> >>>
> >>> DSET ^atmos_month.nc.0000
> >>> UNDEF 1.e-10
> >>> DTYPE netcdf
> >>> TITLE 6 Months Sample Model Output
> >>> *
> >>> XDEF 96 LINEAR 0.0 3.75
> >>> *
> >>> YDEF 80 GAUSR30 1
> >>> *
> >>> ZDEF 18 LEVELS 1000 950 900 850 800 750 700 650 600 550 500 450 400
> 350
> >>> 300
> >>> 250 200 150
> >>> *
> >>> TDEF 804 LINEAR JAN0 1MO
> >>> *
> >>> VARS 7
> >>> u 18 t,z,y,x u
> >>> v 18 t,z,y,x v
> >>> w 18 t,z,y,x v
> >>> t 18 t,z,y,x v
> >>> sh 18 t,z,y,x v
> >>> rh 18 t,z,y,x v
> >>> z 18 t,z,y,x v
> >>> ENDVARS
> >>>
> >>> Good luck!
> >>>
> >>> Hosmay Lopez
> >>> UM/CIMAS/RSMAS NOAA/AOML/PHOD
> >>> 4301 Rickenbacker Causeway,
> >>> Miami, FL 33149, USA
> >>> Phone: (305) 361-4329
> >>>
> >>> ------------------------------
> >>> *From:* gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org]
> >>> on behalf of David Stanfel [dstanfel at ualberta.ca]
> >>> *Sent:* Thursday, June 19, 2014 12:07 AM
> >>> *To:* GrADS Users Forum
> >>> *Subject:* Re: [gradsusr] Displaying a .nc file
> >>>
> >>> It did prompt a list of variables including u, v, w, and z. It still
> >>> doesn't want to display them.
> >>>
> >>> ga-> q file
> >>> File 1 : 6 Months Sample Model Output
> >>> Descriptor: cont.ctl
> >>> Binary: atmos_month.nc.0000
> >>> Type = Gridded
> >>> Xsize = 96 Ysize = 80 Zsize = 18 Tsize = 804 Esize = 1
> >>> Number of Variables = 7
> >>> u 18 0 u
> >>> v 18 0 v
> >>> w 18 0 v
> >>> t 18 0 v
> >>> sh 18 0 v
> >>> rh 18 0 v
> >>> z 18 0 v
> >>> ga-> d v
> >>> Data Request Error: Error for variable 'v'
> >>> Error ocurred at column 1
> >>> DISPLAY error: Invalid expression
> >>> Expression = v
> >>> ga-> d sh
> >>> Error: Variable sh not in netcdf file
> >>> Data Request Error: Error for variable 'sh'
> >>> Error ocurred at column 1
> >>> DISPLAY error: Invalid expression
> >>> Expression = sh
> >>>
> >>> Thanks,
> >>>
> >>> David
> >>>
> >>>
> >>> On Wed, Jun 18, 2014 at 11:58 PM, James T. Potemra <jimp at hawaii.edu>
> >>> wrote:
> >>>
> >>>> If you enter "q file" at the GrADS prompt it should list the
> variables.
> >>>>
> >>>>
> >>>> On 6/18/14 5:52 PM, David Stanfel wrote:
> >>>>
> >>>> How do I know what variable name will work? I've tried u, v, w, and z
> >>>> and still only get:
> >>>>
> >>>> ga-> d v
> >>>> Error: Variable v not in netcdf file
> >>>> Data Request Error: Error for variable 'v'
> >>>> Error ocurred at column 1
> >>>> DISPLAY error: Invalid expression
> >>>> Expression = v
> >>>>
> >>>> Thanks,
> >>>>
> >>>> David
> >>>>
> >>>>
> >>>> On Wed, Jun 18, 2014 at 11:13 PM, Hosmay Lopez <
> hlopez at rsmas.miami.edu>
> >>>> wrote:
> >>>>
> >>>>> David,
> >>>>>
> >>>>> You should let GrADS know that the file is a netCDF. You do that by
> >>>>> adding the following text line in your control file:
> >>>>>
> >>>>> DTYPE netcdf
> >>>>>
> >>>>> Insert such after the UNDEF line.
> >>>>>
> >>>>> Good luck!
> >>>>>
> >>>>> Hosmay Lopez
> >>>>> UM/CIMAS/RSMAS NOAA/AOML/PHOD
> >>>>> 4301 Rickenbacker Causeway,
> >>>>> Miami, FL 33149, USA
> >>>>> Phone: (305) 361-4329
> >>>>>
> >>>>> ------------------------------
> >>>>> *From:* gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org
> ]
> >>>>> on behalf of David Stanfel [dstanfel at ualberta.ca]
> >>>>> *Sent:* Wednesday, June 18, 2014 10:58 PM
> >>>>> *To:* GrADS Users Forum
> >>>>> *Subject:* [gradsusr] Displaying a .nc file
> >>>>>
> >>>>> So I figured out how to use a descriptor file but I am still having
> >>>>> troubles opening the file. My descriptor file is as follows:
> >>>>>
> >>>>> DSET ^atmos_month.nc.0000
> >>>>> UNDEF 1.e-10
> >>>>> TITLE 6 Months Sample Model Output
> >>>>> *
> >>>>> XDEF 96 LINEAR 0.0 3.75
> >>>>> *
> >>>>> YDEF 80 GAUSR30 1
> >>>>> *
> >>>>> ZDEF 18 LEVELS 1000 950 900 850 800 750 700 650 600 550 500 450 400
> >>>>> 350 300
> >>>>> 250 200 150
> >>>>> *
> >>>>> TDEF 804 LINEAR JAN0 1MO
> >>>>> *
> >>>>> VARS 7
> >>>>> u 18 0 u
> >>>>> v 18 0 v
> >>>>> w 18 0 v
> >>>>> t 18 0 v
> >>>>> sh 18 0 v
> >>>>> rh 18 0 v
> >>>>> z 18 0 v
> >>>>> ENDVARS
> >>>>>
> >>>>> It was able to open it without issue:
> >>>>>
> >>>>> ga-> open cont.ctl
> >>>>> Scanning description file: cont.ctl
> >>>>> Data file atmos_month.nc.0000 is open as file 1
> >>>>> LON set to 0 360
> >>>>> LAT set to -88.29 88.29
> >>>>> LEV set to 1000 1000
> >>>>> Time values set: 0:1:1:0 0:1:1:0
> >>>>> E set to 1 1
> >>>>>
> >>>>> However, I am having trouble displaying the data. I tried to use the
> >>>>> display command (d) with this result:
> >>>>>
> >>>>> ga-> d atmos_month.nc.0000
> >>>>> Syntax error: Bad file number for variable atmos_month
> >>>>> Error ocurred at column 1
> >>>>> DISPLAY error: Invalid expression
> >>>>> Expression = atmos_month.nc.0000
> >>>>>
> >>>>> Is my syntax off?
> >>>>>
> >>>>> Thanks,
> >>>>>
> >>>>> David
> >>>>>
> >>>>> _______________________________________________
> >>>>> gradsusr mailing list
> >>>>> gradsusr at gradsusr.org
> >>>>> http://gradsusr.org/mailman/listinfo/gradsusr
> >>>>>
> >>>>>
> >>>>
> >>>>
> >>>> _______________________________________________
> >>>> gradsusr mailing listgradsusr at gradsusr.orghttp://
> gradsusr.org/mailman/listinfo/gradsusr
> >>>>
> >>>>
> >>>>
> >>>> _______________________________________________
> >>>> gradsusr mailing list
> >>>> gradsusr at gradsusr.org
> >>>> http://gradsusr.org/mailman/listinfo/gradsusr
> >>>>
> >>>>
> >>>
> >>> _______________________________________________
> >>> gradsusr mailing list
> >>> gradsusr at gradsusr.org
> >>> http://gradsusr.org/mailman/listinfo/gradsusr
> >>>
> >>>
> >> _______________________________________________
> >> gradsusr mailing list
> >> gradsusr at gradsusr.org
> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >>
> >>
> >> --
> >> Jennifer M. Adams
> >> Center for Ocean-Land-Atmosphere Studies (COLA)
> >> 111 Research Hall, Mail Stop 2B3
> >> George Mason University
> >> 4400 University Drive
> >> Fairfax, VA 22030
> >>
> >>
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> gradsusr mailing list
> >> gradsusr at gradsusr.org
> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >>
> >>
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > --
> > Jennifer M. Adams
> > Center for Ocean-Land-Atmosphere Studies (COLA)
> > 111 Research Hall, Mail Stop 2B3
> > George Mason University
> > 4400 University Drive
> > Fairfax, VA 22030
> >
> >
> >
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20140620/d5dc051d/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 52, Issue 68
> ****************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140623/4c52e9c8/attachment-0001.html
More information about the gradsusr
mailing list