[ncl-talk] help
george luke otieno
gotieno2000 at gmail.com
Tue Jul 28 03:48:51 MDT 2015
Dear All; Trying to plot 3 seasons MAM,OND,JJA together. My data has
3dimension ordered like (time,lat,lon)
This is the complain from this line ....*res at gsnCenterString =
season(ns)+":"+time(0)/100 +"-"+ time(ntim-1)/100*
*fatal:Number of subscripts do not match number of dimensions of
variable,(1) Subscripts used, (2) Subscripts
expectedfatal:["Execute.c":8578]:Execute: Error occurred at or near line
160 in file trial*
*Below is the script..*
*;***********************************************; climo_3.ncl;; Concepts
illustrated:; - Calculating monthly climatologies; - Calculating
interannual variability; - Calculating standard deviations of monthly
means; - Creating a color map using named colors; - Paneling four plots
on a page; - Adding a common labelbar to paneled plots; - Explicitly
setting contour levels; - Changing the center longitude for a cylindrical
equidistant projection; - Turning off the individual contour labelbar;
- Adding a common title to paneled plots; - Turning off map fill; -
Creating a center
subtitle;;***********************************************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"*
*begin;************************************************; Read the
file;************************************************ a =
addfile("GPCP.nc","r") prec =
a->precip(0:431,:,:) ; (time,lat,lon)
printVarSummary(prec) time = a->time ;
(time) ==> YYMM ntim = dimsizes(time) ymStrt =
198001 ; start yyyymm ymLast =
199912 ; last yrStrt = ymStrt/100 yrLast
= ymLast/100 ; start yyyymm wksType =
"x11" wksName = "test" ; ": "+yrStrt+"_"+yrLast
;************************************************; plots ; Compute
3-month
climatologies;************************************************
; (2,6) season = (/ (/ 2, 3, 4/) \ ; March-May[MAM] ,
(/ 5, 6, 7/) \ ; June-August [JJA] , (/ 9, 10, 11/) /)
; Oct-Dec [OND] i_season = season - 1 ; NCL indices
season_label = (/ "March-May", "June-August","October-December"/) plot =
new ( 3, "graphic") wks = gsn_open_wks(wksType, wksName)
gsn_define_colormap(wks,"BlueYellowRed")
;************************************************; Compute the MAM,JJA,OND
climatology using a function in
contributed.ncl;************************************************; precClm
= month_to_season(prec) ; monthly climatology
;************************************************; create
colors;************************************************* wks =
gsn_open_wks("x11","test") ; open a workstation
colors = (/ (/255,255,255/),(/0,0,0/),(/255,255,255/), (/244,255,244/), \
(/217,255,217/), (/163,255,163/), (/106,255,106/), \ (/43,255,106/),
(/0,224,0/), (/0,134,0/),(/255,255,0/),\ (/255,127,0/) /) * 1.0 ; we
multiply by 1 to make colors float colors = colors/255. ;
normalize (required by NCL)*
* ;colors = (/"white","Balck","White","RoyalBlue"\ ; back/fore ground
; ,"lightSkyblue","lightseagreen","palegreen"/) ; choose
colors for color map gsn_define_colormap(wks,
colors) ; generate new color
map;************************************************; create panel
plots;************************************************* plot = new ( 3,
graphic) ; create graphical array*
res = True ; plot options desired
res at cnFillOn = True ; turn on color fill
res at cnInfoLabelOn = False ; turn off contour info
label
res at cnLinesOn = False ; turn off contour lines
res at cnLineLabelsOn = False ; turn off line labels
res at cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour
levels
res at cnLevels = (/0.0,0.5,1.0,5.0,10,15,20,25\ ; set unequal
contour levels
,30,35/)
res at mpFillOn = False ; turn off gray continents
res at mpCenterLonF = 180 ; Centers the plot at 180
res at lbLabelBarOn = False ; No single label bar
res at gsnDraw = False
res at gsnFrame = False
resP = True ; panel options
resP at txString = "GPCC Observation" ; common title
resP at gsnMaximize = True ; maximize image
resP at gsnPanelLabelBar = True ; Add common label bar
;****************ZOOM INTO YOUR LOCATION*********************
res at mpLimitMode = "LatLon"
mpOutlineOn =True
mpOutlineSpecifiers =True
res at mpMinLonF = 21. ; select a subregion
res at mpMaxLonF = 54.
res at mpMinLatF = -12.
res at mpMaxLatF = 22.5
res at mpLandFillColor = "background" ; color of land
res at cnFillDrawOrder = "Predraw"
;*********************Masking*********************************
res at mpAreaMaskingOn = True
res at mpMaskAreaSpecifiers = "Land"
res at mpOceanFillColor = 0
res at mpFillDrawOrder = "Postdraw"
resMap= True
resMap at mpPerimOn = False
resMap at mpGridAndLimbOn = False
resMap at gsnTickMarksOn = False
resMap at mpDataSetName = "Earth..3"
resMap at mpDataBaseVersion = "MediumRes"
resMap at mpOutlineOn = True
res = True
res at mpFillOn = False ; turn off gray fill
res at mpOutlineBoundarySets = "National" ; turn on country
boundaries
;res at mpGeophysicalLineColor = "Navy" ; color of cont. outlines
resMap at mpFillOn = True
resMap at mpOceanFillColor = "transparent"
resMap at mpLandFillColor = "transparent"
resMap at mpInlandWaterFillColor = "transparent"
res at mpDataBaseVersion = "MediumRes" ; choose higher resolution
;delete(precClm at long_name) ; do not want
do ns =0,1
do ns=0,1 ; loop over the months
*res at gsnCenterString = season(ns)+":"+time(0)/100 +"-"+
time(ntim-1)/100*
plot(i) = gsn_csm_contour_map_ce(wks,precClm(ns,:,:), res) ; create
plot
end do
gsn_panel(wks,plot,(/2,2/),resP)
;************************************************
; create second panel plot
;*************************************************
delete(precStd at long_name)
; Std Deviations
do ns=0,1 ; loop over the months
res at gsnCenterString = season(ns)+":"+time(0)/100 +"-"+
time(ntim-1)/100
plot(i) = gsn_csm_contour_map_ce(wks,precStd(ns,:,:), res) ; create
plot
end do
gsn_panel(wks,plot,(/2,2/),resP)
end do
end
Best Regards
George Otieno
Intergrated Climate System Modelling Lab.
Department of Environmental Atmospheric Sciences
608-737, Room 4302, Chungmu-Building
Pukyong National University,
45,Yongso-ro,Nam-gu, Busan, South Korea
Tel. +82-51-629-6643
Fax+82-51-629-7991
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150728/853e359f/attachment.html
More information about the ncl-talk
mailing list