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" load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/ut_string.ncl" ;; tickmark.ncl begin ;; Generate a nylat x nxlon array of randaom data, with a minimum of -20 and maximum of 50 ny = 121 nx = 240 test_var = generate_2d_array(15,15,-20.0,50.0,0,(/ny,nx/)) ;; Generate 1D lat/lon coordinate arrays and attach to the data, including units lat = fspan(-90.0,90.0,ny) lon = fspan(0.0,358.5,nx) test_var!0 = "lat" test_var!1 = "lon" test_var&lat = lat test_var&lat@units = "degrees_north" test_var&lon = lon test_var&lon@units = "degrees_east" ;; Open a workstation and plot a region on a cylindrical equidistant (default projection) map wks = gsn_open_wks("eps","tickmark_test") gsn_define_colormap(wks,"BlueRed") res = True res@gsnSpreadColors = True res@gsnDraw = False res@gsnFrame = False res@cnLinesOn = False res@cnFillOn = True res@cnInfoLabelOn = False res@cnLevelSelectionMode = "ManualLevels" res@cnMaxLevelValF = 50.0 res@cnMinLevelValF = -20.0 res@cnLevelSpacingF = 5.0 res@cnFillDrawOrder = "Predraw" res@cnLineLabelsOn = False res@tmXBOn = True res@tmXTOn = False res@tmYLOn = True res@tmYROn = False res@tmXBBorderOn = True res@tmXTBorderOn = True res@tmYLBorderOn = True res@tmYRBorderOn = True res@gsnCenterString = "Tickmark test" res@gsnCenterStringFontHeightF = 0.015 res@gsnLeftString = "" res@gsnRightString = "" res@mpMaxLatF = 60.0 ; select subregion res@mpMinLatF = 0.0 res@mpMinLonF = 40. res@mpMaxLonF = 200. res@mpCenterLonF = 120 res@mpCenterLatF = 30 res@mpOutlineDrawOrder = "PostDraw" res@mpFillOn = False res@gsnAddCyclic = False res@mpGridAndLimbOn = True ; print map gridlines res@mpGridLineDashPattern = 2 res@mpGridLatSpacingF = 20.0 res@mpGridLonSpacingF = 30.0 plot = gsn_csm_contour_map(wks,test_var(:,:),res) draw(plot) frame(wks) end