;************************************************* ; conLab_3.ncl ;************************************************ 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 ;************************************************ a = addfile("uv300.nc","r") ;************************************************ ; read in zonal winds ;************************************************ u = a->U(1,:,:) ;************************************************ ; create default plot ;************************************************ wks = gsn_open_wks("ps","conLab") ; open a ps file gsn_define_colormap(wks,"gui_default") ; choose colormap res = True res@cnFillOn = True ; turn on color res@gsnSpreadColors = True ; use full range of colormap res@mpMinLatF = -30 ; zoom map res@mpMaxLatF = 30 res@mpMinLonF = 0 res@mpMaxLonF = 70 res@cnLineLabelsOn = True ; turn on line labels res@cnLineLabelBackgroundColor = "white" ; white bckgrnd around label plot = gsn_csm_contour_map_ce(wks,u,res) ; create a cylindrical equidistant plot end