<div dir="ltr"><div>Dear All</div><div>I am trying to make 6plots BUT only two are being made. If i set </div><div><br></div><div>plot(0) = gsn_csm_contour_map_ce(wks,precSeaClm(ns,:,:), res)  ; create plot<br>    plot(1) = gsn_csm_contour_map_ce(wks,prec1SeaClm(ns,:,:), res)  ; create plot</div><div><br></div><div>I get only 3plots</div><div>BUT IF I SET </div><div>plot(ns) = gsn_csm_contour_map_ce(wks,precSeaClm(ns,:,:), res)  ; create plot<br>    plot(ns) = gsn_csm_contour_map_ce(wks,prec1SeaClm(ns,:,:), res)  ; create plot </div><div>I only get 2plots..</div><div><br></div><div>kindly help</div><div><br></div><div>Below is my script full script and also attached...</div><div><br></div><div>;***********************************************<br>; climo_3.ncl<br>;<br>; Concepts illustrated:<br>;   - Calculating monthly climatologies<br>;   - Calculating seasonal totals for 36yrs<br>;   - Creating a color map using named colors<br>;   - Paneling four plots on a page<br>;   - Adding a common labelbar to paneled plots<br>;   - Explicitly setting contour levels<br>;   - Changing the center longitude for a cylindrical equidistant projection<br>;   - Turning off the individual contour labelbar<br>;   - Adding a common title to paneled plots<br>;   - Turning off map fill<br>;   - Creating a center subtitle<br>;<br>;***********************************************<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl&quot;</div><div>begin<br>;*******************************************************<br>; Read the file for GPCC<br>;*******************************************************<br>  ymStrt = 198101 <br>  ymLast = 199012<br>  a    = addfile(&quot;GPCC.nc&quot;,&quot;r&quot;)<br> <br>  TIME = a-&gt;time<br>  YYYYMM = cd_calendar(TIME, -1)<br>  iStrt  = ind(YYYYMM.eq.ymStrt)<br>  iLast  = ind(YYYYMM.eq.ymLast)</div><div>  prec = a-&gt;precip(iStrt:iLast,:,:)              ; (time,lat,lon)<br>  delete(prec@units)<br>  delete(prec@long_name)<br>  printMinMax(prec,True)<br>  <br>  printVarSummary(prec)<br>   time = a-&gt;time                                ; (time) ==&gt; YYMM<br>  ntim = dimsizes(time)<br>  precClm = clmMonTLL(prec)<br>   printVarSummary(precClm)</div><div>;*******************************************************<br>; Read the file for GPCP<br>;*******************************************************<br>  ymStrt = 198101 <br>  ymLast = 199012<br>  b    = addfile(&quot;GPCP.nc&quot;,&quot;r&quot;)<br>  ;TIME = b-&gt;time<br>  time1 =b-&gt;time<br>  delete(TIME)<br>  yyyymm = cd_calendar(time1, -1)<br>  iStrt  = ind(yyyymm.eq.ymStrt)<br>  iLast  = ind(yyyymm.eq.ymLast)</div><div>  prec1 = b-&gt;precip(iStrt:iLast,:,:)              ; (time,lat,lon)<br>  delete(prec1@units)<br>  delete(prec1@long_name)<br>  printMinMax(prec1,True)<br>  <br>  printVarSummary(prec1)<br>   time1 = b-&gt;time                             ; (time) ==&gt; YYMM<br>  ntim1 = dimsizes(time1)<br>  prec1Clm = clmMonTLL(prec1)<br>   printVarSummary(prec1Clm)</div><div>;************************************************<br>; Compute the monthly climatology using a function in contributed.ncl<br>;************************************************<br>  precMonClm = clmMonTLL(prec)   ; (12,:,:)<br>  precSeaClm = runave_n_Wrap(precMonClm, 3, -1, 0)  ; unweighted seasonal average<br>  printVarSummary(precSeaClm)<br>  <br>   prec1MonClm = clmMonTLL(prec1)   ; (12,:,:)<br>  prec1SeaClm = runave_n_Wrap(prec1MonClm, 3, -1, 0)  ; unweighted seasonal average<br>  printVarSummary(prec1SeaClm)<br>  <br>   season_index = (/3,6,10/)  ; MAM, JJA, OND <br>  season_label = (/ &quot;MAM&quot;, &quot;JJA&quot;,&quot;OND&quot;/)<br>  <br>;************************************************<br>; create colors<br>;*************************************************<br>  wks = gsn_open_wks(&quot;x11&quot;, &quot;clim_gpcp&quot;)          ; open  a work station<br>                                          <br>   <br> colors = (/ (/255,255,255/),(/255,255,255/), (/244,255,244/), \<br>  (/217,255,217/), (/163,255,163/), (/106,255,106/), \<br>  (/43,255,106/), (/0,224,0/), (/0,134,0/),(/255,255,0/),\<br>  (/255,127,0/) /) * 1.0 ; we multiply by 1 to make colors float<br>  <br>    colors = colors/255.               ; normalize (required by NCL)</div><div> ;************************************************<br>; create panel plots<br>;*************************************************<br> plot   = new ( 6, graphic)                    ; create graphical array<br> <br>    res                      = True               ; plot options desired<br>    <br>  ;res@cnFillDrawOrder      = &quot;Predraw&quot;    <br>  res@cnFillOn             = True               ; turn on color fill<br>  res@cnInfoLabelOn        = False              ; turn off contour info label<br>  res@cnLinesOn            = False              ; turn off contour lines<br>  res@cnLineLabelsOn       = False              ; turn off line labels<br>  res@cnFillPalette        =colors<br>  res@cnLevelSelectionMode = &quot;ExplicitLevels&quot;   ; set explicit contour levels<br> res@cnLevels        =(/0.1,0.2,0.4,0.8,1.6,3.2,6.4,12.8,25,30,40,45,60,75,83,100/)</div><div>                             <br>  res@mpFillOn             = False              ; turn off gray continents<br>  res@mpCenterLonF         = 180                ; Centers the plot at 180 <br> <br>  res@lbLabelBarOn         = False              ; No single label bar</div><div>  res@gsnDraw              = False<br>  res@gsnFrame             = False</div><div>  resP                     = True               ; panel options<br>  resP@txString            = &quot;Seasons&quot;   ; common title<br>  resP@gsnMaximize         = True               ; maximize image<br>  resP@gsnPanelLabelBar    = True               ; Add common label bar <br>  <br>  ;****************ZOOM INTO YOUR LOCATION*********************<br>  <br>  res@mpLimitMode           = &quot;LatLon&quot;<br> ;res@mpOutlineOn    =True<br> ;res@mpOutlineSpecifiers     =True<br>  res@mpMinLonF            = 21.               ; select a subregion<br>  res@mpMaxLonF            = 54.<br>  res@mpMinLatF            = -12. <br>  res@mpMaxLatF            = 22.5<br>  <br>  <br>  ;res@mpLandFillColor      = &quot;background&quot;       ; color of land<br>  <br>  ;*********************Masking*********************************<br> </div><div>      res@mpFillDrawOrder = &quot;Postdraw&quot;  <br>    <br>  <br>      resMap= True <br>      resMap@mpPerimOn = False <br>     resMap@mpGridAndLimbOn = False <br>    resMap@gsnTickMarksOn = False <br>    resMap@mpDataSetName = &quot;Earth..3&quot; <br>     resMap@mpOutlineOn = True <br>   <br>  <br>  res@mpFillOn                    = False        ; turn off gray fill<br>  res@mpOutlineBoundarySets       = &quot;National&quot;   ; turn on country boundaries<br>   <br>  resMap@mpFillOn = True </div><div>  resMap@mpOceanFillColor = &quot;transparent&quot; <br>  resMap@mpLandFillColor = &quot;transparent&quot; <br>  resMap@mpInlandWaterFillColor = &quot;transparent&quot; <br>   <br>   res@mpDataBaseVersion = &quot;MediumRes&quot;            ; choose higher resolution<br>  </div><div>   do ns=0,dimsizes(season_index)-1<br>    <br>    res@gsnCenterString   = season_label(ns)+&quot;:&quot;+(ymStrt/100) +&quot;-&quot;+(ymLast/100)<br> <br>    ;res@gsnCenterString   = &quot;&quot;<br>    plot(0) = gsn_csm_contour_map_ce(wks,precSeaClm(ns,:,:), res)  ; create plot<br>    plot(1) = gsn_csm_contour_map_ce(wks,prec1SeaClm(ns,:,:), res)  ; create plot<br>    <br>    <br>   end do</div><div>  gsn_panel(wks,plot,(/3,3/),resP)<br> <br> end</div><div><br></div></div>