request
Don Hooper
hoop at COLORADO.EDU
Tue Jun 14 13:02:15 EDT 2005
Anshu Prakash Mishra,
I'm forwarding your query to the GRADSUSR mailing list, where all such
queries should go in the first place. If you aren't subscribed, check
the GrADS web page for instructions.
> BUt when i am truying to plot it in grdas,
> it plots for only on emonth time while file is made for 12 months climatology,
There's a good reason for this, and the clue comes in the output from the
sdfopen command:
> SDF file has no discernable time coordinate.
The units attribute on the time variable in this NetCDF file is not a
udunits-based time unit. Thus, it doesn't conform to the COARDS metadata
convention, thus sdfopen cannot handle this time axis. So, it makes up
one that has length one.
If you wish to use the time axis of this file, you will have to either write
a DDF (Data Descriptor File) for the xdfopen command, something like:
DSET ^Trop_Atl_Drifter_Monthly_Climatology.nc
TDEF time 12 LINEAR 0Z01JAN0001 1MO
If you put those two lines into an ASCII text file, and then give the
name of that ASCII text file to the xdfopen command at the GrADS prompt,
that should work. Note that the "^" on the DSET indicates that the .nc file
is in the same directory as the ASCII text file. If this is not the case,
you will have to put the full path of the .nc file on the DSET instead of
the "^" character.
It would also be possible to use this file with a "DTYPE NETCDF" DDF;
I'll leave it for others to give the URLs referring the gadoc pages
that are relevant to such an exercise.
Note that you are using a bad build of GrADS; these output lines:
> gxstart=0, gxstop=85, gystart=0, gystop=41.
> gzstart=0, gzstop=1, gtstart=0, gtstop=1.
indicate a build of GrADS that should *never* have been released to the
public. I strongly urge you to obtain a newer build from the COLA website:
http://grads.iges.org/grads/
-Hoop
> From anshu_786 at rediffmail.com Tue Jun 14 10:18:54 2005
> Date: 14 Jun 2005 16:20:24 -0000
> From: "anshu prakash mishra" <anshu_786 at rediffmail.com>
> Reply-To: "anshu prakash mishra" <anshu_786 at rediffmail.com>
> To: hoop at COLORADO.EDU
> Subject: request
>
> Dear sir,
> i have a netcdf file of monthly climatology given below.
> the ncdump -h shows the details of the file given below:
>
> SGI_2100 36% ls *.nc
> Trop_Atl_Drifter_Annual_Climatology.nc
> Trop_Atl_Drifter_Monthly_Climatology.nc
> SGI_2100 37% ncdump -h Trop_Atl_Drifter_Monthly_Climatology.nc
> netcdf Trop_Atl_Drifter_Monthly_Climatology {
> dimensions:
> time = 12 ;
> latitude = 41 ;
> longitude = 85 ;
> variables:
> double time(time) ;
> time:name = "time" ;
> time:long_name = "Climatological month" ;
> time:generic_name = "time" ;
> time:type = "EVEN" ;
> time:units = "Climatological month" ;
> time:data_min = 1. ;
> time:data_max = 12. ;
> time:F_format = "f9.0" ;
> time:C_format = "%9.0f" ;
> double latitude(latitude) ;
> latitude:name = "lat" ;
> latitude:long_name = "latitude of U/V grid" ;
> latitude:type = "UNEVEN" ;
> latitude:units = "degrees_north" ;
> latitude:C_format = "%8.4f" ;
> latitude:F_format = "f8.4" ;
> latitude:epic_code = 500. ;
> double longitude(longitude) ;
> longitude:name = "lon" ;
> longitude:long_name = "latitude of U/V grid" ;
> longitude:type = "UNEVEN" ;
> longitude:units = "degrees_east" ;
> longitude:C_format = "%8.4f" ;
> longitude:F_format = "f8.4" ;
> longitude:epic_code = 500. ;
> double u(latitude, longitude, time) ;
> u:long_name = "Zonal Current" ;
> u:units = "m/s" ;
> u:data_min = -999. ;
> u:data_max = 999. ;
> u:C_format = "%10.4f" ;
> u:F_format = "f10.4" ;
> u:_FillValue = -9.99e+33 ;
> u:epic_code = 318. ;
> double v(latitude, longitude, time) ;
> v:long_name = "Meridional Current" ;
> v:units = "m/s" ;
> v:data_min = -999. ;
> v:data_max = 999. ;
> v:C_format = "%10.4f" ;
> v:F_format = "f10.4" ;
> v:_FillValue = -9.99e+33 ;
> v:epic_code = 319. ;
> double sst(latitude, longitude, time) ;
> sst:long_name = "Sea Surface Temperature" ;
> sst:units = "degreeC" ;
> sst:data_min = -999. ;
> sst:data_max = 999. ;
> sst:C_format = "%10.4f" ;
> sst:F_format = "f10.4" ;
> sst:_FillValue = -9.99e+33 ;
> sst:epic_code = 25. ;
> double ug(latitude, longitude, time) ;
> ug:long_name = "Geostrophic Zonal Current" ;
> ug:units = "m/s" ;
> ug:data_min = -999. ;
> ug:data_max = 999. ;
> ug:C_format = "%10.4f" ;
> ug:F_format = "f10.4" ;
> ug:_FillValue = -9.99e+33 ;
> ug:epic_code = 340. ;
> double vg(latitude, longitude, time) ;
> vg:long_name = "Geostrophic Meridional Current" ;
> vg:units = "m/s" ;
> vg:data_min = -999. ;
> vg:data_max = 999. ;
> vg:C_format = "%10.4f" ;
> vg:F_format = "f10.4" ;
> vg:_FillValue = -9.99e+33 ;
> vg:epic_code = 341. ;
> double count(latitude, longitude) ;
> count:long_name = "Sea Surface Temperature" ;
> count:units = "drifter days per square deg" ;
> count:data_min = -999. ;
> count:data_max = 999. ;
> count:C_format = "%10.4f" ;
> count:F_format = "f10.4" ;
> count:_FillValue = -9.99e+33 ;
> double uE(latitude, longitude) ;
> uE:long_name = "Zonal Current Standard Error" ;
> uE:units = "m/s" ;
> uE:data_min = -999. ;
> uE:data_max = 999. ;
> uE:C_format = "%10.4f" ;
> uE:F_format = "f10.4" ;
> uE:_FillValue = -9.99e+33 ;
> double vE(latitude, longitude) ;
> vE:long_name = "Meridional Current Standard Error" ;
> vE:units = "m/s" ;
> vE:data_min = -999. ;
> vE:data_max = 999. ;
> vE:C_format = "%10.4f" ;
> vE:F_format = "f10.4" ;
> vE:_FillValue = -9.99e+33 ;
> double sstE(latitude, longitude) ;
> sstE:long_name = "Sea Surface Temperature Standard Error" ;
> sstE:units = "degreeC" ;
> sstE:data_min = -999. ;
> sstE:data_max = 999. ;
> sstE:C_format = "%10.4f" ;
> sstE:F_format = "f10.4" ;
> sstE:_FillValue = -9.99e+33 ;
> double ugE(latitude, longitude) ;
> ugE:long_name = "Geostrophic Zonal Current Standard Error" ;
> ugE:units = "m/s" ;
> ugE:data_min = -999. ;
> ugE:data_max = 999. ;
> ugE:C_format = "%10.4f" ;
> ugE:F_format = "f10.4" ;
> ugE:_FillValue = -9.99e+33 ;
> double vgE(latitude, longitude) ;
> vgE:long_name = "Geostrophic Meridional Current Standard Error" ;
> vgE:units = "m/s" ;
> vgE:data_min = -999. ;
> vgE:data_max = 999. ;
> vgE:C_format = "%10.4f" ;
> vgE:F_format = "f10.4" ;
> vgE:_FillValue = -9.99e+33 ;
>
> // global attributes:
> :HISTORY = "Initial release. 01Nov2003" ;
> :CONVENTIONS = "COARDS" ;
> :DATA_ORIGIN = "http://www.aoml.noaa.gov.edu/phod/~lumpkin/drifter_climatology.html" ;
> :DATA_AUTHOR = "Dr Rick Lumpkin" ;
> :FILE_AUTHOR = "Derrick Snowden" ;
> :FILE_CREATION_DATE = "14:24 18-Jan-2005" ;
> :FILE_SOURCE = "Atlantic Oceanographic and Meteorological Laboratory/PhOD" ;
> :ACKNOWLEDGEMENT01 = "This data is distributed in good faith." ;
> :ACKNOWLEDGEMENT02 = "It is an estimate and may not be considered exact." ;
> :ACKNOWLEDGEMENT03 = "If you choose to use this data, the source can be " ;
> :ACKNOWLEDGEMENT04 = "acknowledged by citing:" ;
> :ACKNOWLEDGEMENT05 = "Lumkin, Rick (2003). Decomposition of surface" ;
> :ACKNOWLEDGEMENT06 = "drifter observations in the Atlantic Ocean." ;
> :ACKNOWLEDGEMENT07 = "Geophysical Research Letters, V30, N14 Pp1753-1756." ;
> :missing_value = -9.99e+33 ;
> :_FillValue = -9.99e+33 ;
> }
>
> BUt when i am truying to plot it in grdas,
> it plots for only on emonth time while file is made for 12 months climatology,
> sdfopen open give following details:
> ga-> !ls *.nc
> Trop_Atl_Drifter_Annual_Climatology.nc
> Trop_Atl_Drifter_Monthly_Climatology.nc
> ga-> sdfopen Trop_Atl_Drifter_Monthly_Climatology.nc
> Scanning self-describing file: Trop_Atl_Drifter_Monthly_Climatology.nc
> SDF file has no discernable time coordinate.
> Found 11 data variables in SDF file.
> Found displayable variable u with 0 levels in SDF file.
> SDF: The double precision values of the variable u
> will be automatically converted to single precision.
> Found displayable variable v with 0 levels in SDF file.
> SDF: The double precision values of the variable v
> will be automatically converted to single precision.
> Found displayable variable sst with 0 levels in SDF file.
> SDF: The double precision values of the variable sst
> will be automatically converted to single precision.
> Found displayable variable ug with 0 levels in SDF file.
> SDF: The double precision values of the variable ug
> will be automatically converted to single precision.
> Found displayable variable vg with 0 levels in SDF file.
> SDF: The double precision values of the variable vg
> will be automatically converted to single precision.
> Found displayable variable count with 0 levels in SDF file.
> SDF: The double precision values of the variable count
> will be automatically converted to single precision.
> Found displayable variable ue with 0 levels in SDF file.
> SDF: The double precision values of the variable ue
> will be automatically converted to single precision.
> Found displayable variable ve with 0 levels in SDF file.
> SDF: The double precision values of the variable ve
> will be automatically converted to single precision.
> Found displayable variable sste with 0 levels in SDF file.
> SDF: The double precision values of the variable sste
> will be automatically converted to single precision.
> Found displayable variable uge with 0 levels in SDF file.
> SDF: The double precision values of the variable uge
> will be automatically converted to single precision.
> Found displayable variable vge with 0 levels in SDF file.
> SDF: The double precision values of the variable vge
> will be automatically converted to single precision.
> SDF file Trop_Atl_Drifter_Monthly_Climatology.nc is open as file 1
> LON set to -70 14
> LAT set to -20 20
> LEV set to 0 0
> Time values set: 1:1:1:0 1:1:1:0
> ga->ga-> set t 1
> Time values set: 1:1:1:0 1:1:1:0
> ga-> d sst
> gxstart=0, gxstop=85, gystart=0, gystop=41.
> gzstart=0, gzstop=1, gtstart=0, gtstop=1.
> Contouring: 17 to 29 interval 1
> ga-> set t 2
> Time values set: 1:1:1:0 1:1:1:0
> ga-> d sst
> Data Request Warning: Request beyond file limits
> Entire grid contents are set to missing data
> Dimension ranges are: X = 1 85 Y = 1 41 Z = 1 1 T = 2 1
> Warning issued for variable = sst
> Cannot contour grid - all undefined values
> ga->
>
> plhelp me where is problem an dhow i can plot it for all the 12 months in
> Grads. i thought but could not resolve it.
> pl help me,
> with best wishes,
> Anshu
More information about the gradsusr
mailing list