[ncl-talk] Zoom plot
Tabish Ansari
tabishumaransari at gmail.com
Tue Nov 10 14:53:13 MST 2015
Hi
I have a 3d global data for ozone (time, lat, lon) and I want to plot only
over south asia. I am able to get a nice plot when plotting the full data
(with the world map in the background), but couldn't achieve the same for
south asia.
I am currently using the gsn_csm_contour_map function.
My script is below:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
a=addfile("/data2/ACCMIP/processed/acchist/o3tropcol_2000Trop/UM-CAM/r1i1p1/v2/o3tropcol_ACCMIP-monthly_UM-CAM_acchist_r1i1p1_200101-200112.nc","r")
wks = gsn_open_wks("x11","ozone-jan")
ozone = a->o3tropcol(1,10:30,20:35)
res=True
res at cnFillOn= True ; turn on color
gsn_define_colormap(wks,"gui_default") ; choose a colormap
res at cnLineLabelsOn = False ; Turn off contour line labels
res at cnLinesOn = False ; Turn off contour line labels
res at cnLevelSelectionMode = "ManualLevels" ; manually set the contour
levels with the following 3 resources
res at cnMinLevelValF = 0 ; set the minimum contour level
res at cnMaxLevelValF = 52 ; set the maximum contour
level
res at cnLevelSpacingF = 2 ; set the interval between contours
plot= gsn_csm_contour_
map
(wks,ozone,res)
end
Thanks in advance.
Tabish
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151110/28519de3/attachment.html
More information about the ncl-talk
mailing list