three maps in one plot

From: ousmane ndiaye <ondiaye70_at_nyahnyahspammersnyahnyah>
Date: Fri Mar 23 2012 - 10:29:54 MDT

Hi,
I'm drawing three maps on a plot but want each map to span the whole width
of the page.
Here is my plot and my script, what should I change
ousmane

==
; the model file
in = addfile(files(i),"r");
pr=in->tp_P11_L1_GLL0_acc(:,:,:,{12.5:17},{342:349}) ; extracting senegal
xm=dim_avg_n(dim_avg_n(pr,4),3) ;average XY over senegal

;************************************************
; plotting parameters
;************************************************
res=True
res@gsnFrame = False ; Do not advance page frame
res@gsnDraw = False ;so that individual plots are not drawn
res@trYMinF = 0 ; Set axes limits.
res@trYMaxF = 90;
res@trXMinF = 1
res@trXMaxF = 184

plot = new(3,graphic) ; create graphics array
d=new(3,graphic)

;=== start the three panels plots
do j=0,2
xmfcst= (xm(:,:,ix(j))-xm(:,:,ix(j)-1)) ;forecast lead
xmin = dim_min_n( xmfcst, 0 ) ;min value of all runs
xmax = dim_max_n( xmfcst, 0 ) ;max value of all runs
do k=0,183
        xp(k)=x(k)
        xp(2*184-1-k)=x(k)
        yp(k)=xmax(k)
        yp(2*184-1-k)=xmin(k)
end do

  res@tiMainString = tt(j) ; title
  res@tiYAxisString = "Rainfall in mm/day" ; y axis title
  res@tiXAxisString = "Date" ; x axis title
res@tmXBMode ="Explicit"
res@tmXBValues = (/1,32,62,93,124,154/) ; bottom tick mark
res@tmXBLabels = (/"1My","1Jn","1Jl","1Au","1Sp","1Oc"/) ; labels

res@xyMonoDashPattern = True ;By default, lines are dashed make them all
solid
res@xyLineColor = "Black"
res@xyLineThicknessF = 2 ; 3x as thick
plt = gsn_csm_y(wks,obs,res) ;plot obs
res@xyLineColor = "Blue"
plot(j) = gsn_csm_y(wks,dim_avg_n(xmfcst,0),res) ;plot ens_mean
overlay(plot(j),plt) ;overlay the two and put them in the first
; plot the shading
 gsres = True ; poly res
 gsres@tfPolyDrawOrder = "Predraw" ; draw this first
 gsres@gsFillColor = "SlateBlue" ; color chosen
 d(j) = gsn_add_polygon (wks,plot(j),xp,yp,gsres)
end do

;*******************************************
; draw panel plot with title
;*******************************************
  pres = True ; mod panel plot
  pres@gsnPanelBottom = 0.05 ; add some space at bottom
  pres@txString = model(i)
  pres@gsnMaximize = True ;Maximizing plots in a frame
  gsn_panel(wks,plot,(/3,1/),pres) ; create panel plot
frame(wks)
==

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

ECMWF_.png
Received on Fri Mar 23 10:30:15 2012

This archive was generated by hypermail 2.1.8 : Tue Mar 27 2012 - 08:50:24 MDT