begin ; ; Note that "t" on this file has no attributes like "long_name" ; or "units". Add units attributes to the coordinate arrays and ; see what happens to plot. ; tf = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/Tstorm.cdf","r") ; Open netCDF file. T = tf->t(0,:,:) ; Read in temperature data. T&lon@units = "degrees_east" T&lat@units = "degrees_north" ; printVarSummary(T) wks = gsn_open_wks("png","contour1d") ; Open a png file called "contour1d.png" res = True ; Variable to hold plot options res@cnFillOn = True ; Contour fill on res@lbOrientation = "Vertical" ; Move labelbar to side res@cnFillPalette = "BlueYellowRed" ; Change the color map res@tiMainString = "'BlueYellowRed' color map" plot = gsn_csm_contour(wks,T,res) ; Call the gsn_csm function for ; drawing the contours. end