problem reading/plotting data

Hello,
I am trying to read in a file to plot using gsn_csm_contour_map_ce,
and am having a problem with defining the dimensions of the file.
The script is below, and I've put in bold the lines where I am having
a problem. I also have included below the file attirbutes. I want
to plot the mean field, but am not sure how to define the dimensions
of the file (time, level001, latitude, longitude) when i read in the
data. I know this is basic, but I'm not sure which should be 0 as
opposed to :, but can't find information on this on the website.
Help would be very much appreciated!
Thank you,
Susan

;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
cdf_file = addfile("MeanTOMSAI19801991.nc","r")
;************************************************
; read in mean data
;************************************************
meanai = cdf_file->mean(:,0,:,:)
;************************************************
; plotting parameters
;************************************************
   wks = gsn_open_wks("ps","TOMSmean")
   res = True

gsn_define_colormap(wks,"ViBlGrWhYeOrRe") ; choose colormap

res_at_cnFillOn = True ; turn on color
res_at_gsnSpreadColors = True ; use full colormap
res_at_cnLinesOn = False ; turn off contour lines
res_at_lbLabelAutoStride = True ; automatic lb label
stride
;res_at_cnLevelSelectionMode = "ManualLevels" ; manually set cn levels
;res_at_cnMinLevelValF = 0 ; min level
;res_at_cnMaxLevelValF = 30 ; max level
;res_at_cnLevelSpacingF = 2 ; contour level spacing
SCRIPT TAKES MUCH LONGER TO RUN WITH INTERVAL AT .1

;res_at_mpMinLatF = -60 ; range to zoom in on
;res_at_mpMaxLatF = 60.
;res_at_mpMinLonF = -180.
;res_at_mpMaxLonF = 180.

plot = gsn_csm_contour_map_ce(wks,menai,res)

end

filename: MeanTOMSAI1980-1991
path: MeanTOMSAI1980-1991.nc
    file global attributes:
       title : ./bin/getmomentsfield TOMSData/toms2.ctl mean month 1
ave 12 begin 1980 end 1991 restrain 0 ./data/m497.ctl
       Conventions : COARDS
       GFDL_NCIR_VER_NO : 0.7
       missing_value : 3e+33
    dimensions:
       longitude = 288
       latitude = 120
       level001 = 1
       time = 1 // unlimited
    variables:
       float longitude ( longitude )
          long_name : Longitude values
          cartesian_axis : X
          units : degrees_E
          ipositive : 1

       float latitude ( latitude )
          long_name : Latitudes values
          cartesian_axis : Y
          units : degrees_N
          ipositive : 1

       float level001 ( level001 )
          long_name : Level values
          cartesian_axis : Z
          units : level
          ipositive : 1
          positive : up

       double time ( time )
          long_name : Time coordinate values
          cartesian_axis : T
          units : months since 1-01-15 00:00:00.0
          ipositive : 1
          calendar : Gregorian

       float mean ( time, level001, latitude, longitude )
          long_name : mean
          units : unknown
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float sd ( time, level001, latitude, longitude )
          long_name : standard deviation
          units : unknown
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float skew ( time, level001, latitude, longitude )
          long_name : skewness
          units : 1
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float kurt ( time, level001, latitude, longitude )
          long_name : kurtosis
          units : 1
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float adev ( time, level001, latitude, longitude )
          long_name : absolute deviation
          units : unknown
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float sdm ( time, level001, latitude, longitude )
          long_name : standard deviation / mean
          units : 1
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float min ( time, level001, latitude, longitude )
          long_name : minimum
          units : unknown
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

       float max ( time, level001, latitude, longitude )
          long_name : maximum
          units : unknown
          valid_range : <ARRAY>
          missing_value : 3e+33
          _FillValue : 3e+33

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Nov 16 2008 - 12:31:10 MST