i'm trying to run the following straightforward script. it runs
completely through with no errors - but i get an empty map with no SST
anomalies (i'm trying to plot the Reynolds SST anomaly data from last
week).
here is the script. am i doing something onviously wrong with respect to
setting the graphics control "res@"
;******************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
;************************************************
; read the file
;************************************************
in = addfile("/tigris/dt/hcullen/NCL/data/reynolds_ssta.cdf","r")
vars = getfilevarnames(in)
print(vars)
nvars = dimsizes(vars)
ssta = in->ssta
T= in->T
lat = in->Y
lon = in->X
printVarSummary(ssta)
;************************************************
; assign variables and reorder lat and lon
;************************************************
ssta!0="T"
ssta!1="lat"
ssta!2="lon"
;************************************************
; create plot
;************************************************
wks = gsn_open_wks("X11","hotspots")
gsn_define_colormap(wks,"gui_default")
res = True ; plot options desired
res@mpGridAndLimbOn = False ; turn on lat/lon lines
res@mpGridLatSpacingF = 30. ; spacing for lat lines
res@mpGridLonSpacingF = 30. ; spacing for lon lines
res@gsnSpreadColors = True
res@cnFillOn = True ; turn on color fill
res@cnInfoLabelOn = False ; turn off contour info
label
res@cnLinesOn = False ; turn off contour lines
res@cnLineLabelsOn = False ; turn off line labels
res@gsnAddCyclic = False ; data already has cyclic
point
res@mpOutlineOn = True ; turn on continental
outlines
res@mpOutlineBoundarySets = "National" ; add country boundaries
res@lbLabelBarOn = False ; No single label
bar
res@cnLevelSelectionMode = "ManualLevels" ; manual contour levels
res@cnMinLevelValF = -5.0 ; min level
res@cnMaxLevelValF = 5.0 ; max level
res@cnLevelSpacingF = 0.1 ; interval
printVarSummary(ssta)
plot = gsn_csm_contour_map_ce(wks,ssta(0,:,:), res)
end
-- Heidi M. Cullen Environmental & Societal Impacts Group National Center for Atmospheric Research P. O. Box 3000 Boulder, Colorado 80307-3000 USAFor Overnight, send to: 3450 Mitchell Lane Boulder, Colorado 80301 Phone:303.497.8132 - Fax:303.497.8125 - email:hcullen AT ucar.edu _______________________________________________ ncl-talk mailing list ncl-talk AT ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Mon Jul 29 2002 - 14:59:43 MDT