begin ; ; This grid came to us via Dave Randall, Todd Ringler, and Ross Heikes ; of CSU. The data for this mesh were originally downloaded from: ; http://kiwi.atmos.colostate.edu/BUGS/geodesic/interpolate.html ; ; It will be contoured using NCL's triangular mesh capability. ; f = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/hswm_d000000p000.g2.nc","r") ke = f->kinetic_energy(2,:) wks = gsn_open_wks("png","contour5c") ; Open "contour5c.png" png file. res = True ; Plot mods desired res@cnFillOn = True ; Turn on contour fill res@cnFillPalette = "ncview_default" ; Change the color map res@mpProjection = "Orthographic" ; Change map projection res@mpCenterLonF = -60 ; Change center longitude res@mpPerimOn = False ; Turn off map perimeter res@tiMainString = "~F26~Contouring 1-dimensional x,y,z data" ;---These resources are what set the triangular mesh. r2d = 57.2957795 ; radians to degrees res@sfXArray = f->grid_center_lon * r2d res@sfYArray = f->grid_center_lat * r2d res@sfXCellBounds = f->grid_corner_lon * r2d res@sfYCellBounds = f->grid_corner_lat * r2d map = gsn_csm_contour_map(wks,ke,res) ; Call the gsn_csm function for ; drawing contours over a map. end