begin print("Warning: this example can be slow!") f = addfile("fh.0000_tl.press_gr.egrdsf.partial.nc","r") veg = f->VGTYP_96_SFC veg@lat2d = f->gridlat_96 veg@lon2d = f->gridlon_96 wks = gsn_open_wks("png","contour7f") cmap = read_colormap_file("BlAqGrYeOrRe") ; Read color map res = True ; Plot options desired. res@gsnMaximize = True ; Maximize plot in frame. res@cnFillOn = True ; Contour fill on res@cnLinesOn = False ; Turn off contour lines res@cnFillPalette = cmap(16:,:) ; Set color map res@cnFillMode = "RasterFill" ; Must use raster fill, ; or will run out of memory. res@cnLevelSelectionMode = "ExplicitLevels" ; Set contour levels res@cnLevels = ispan(0,14,1) ; explicitly. res@mpProjection = "LambertConformal" res@mpLambertMeridianF = -95 ; Central meridian res@mpLambertParallel1F = 33 ; Two parallels res@mpLambertParallel2F = 5 res@mpLimitMode = "LatLon" res@mpMinLatF = 24 ; Zoom in on res@mpMaxLatF = 50 ; a smaller res@mpMinLonF = -125 ; latitude and res@mpMaxLonF = -65 ; longitude area. res@pmTickMarkDisplayMode = "Always" ; Turn on map tickmarks res@mpOutlineBoundarySets = "GeophysicalAndUSStates" res@mpGeophysicalLineColor = "red" ; Color of continental outlines res@mpUSStateLineColor = "red" ; Color of US States map = gsn_csm_contour_map(wks,veg,res) ; Call the gsn_csm function for ; drawing contours over a map. end