<div dir="ltr"><div>Hi James</div><div>I am on the right now as you can see on the attached plot. <br></div><div><br></div><div>I am also sending script i used.</div><div><br></div><div>Thank you</div><div><br></div><div><br></div><div><br>  <br>   fname = "<a href="http://tas_2006.nc">tas_2006.nc</a>"<br>    f     = addfile(fname,"r")<br>        t       = f->tas(:, {-15:15}, {5:35})   <br>    <br>    lon = t&lon<br>    lat = t&lat<br>    <br>  <br>  <br>   pltDir = "./"                ; where output will be written<br>   pltName= "panel29"                                   <br>   pltType= "png"               ; ps, eps, png, pdf, x11 .....<br><br>;****************************************************************************<br>; Miscellaneous <br>;****************************************************************************<br>  nrow = 3<br>  ncol = 4                           <br><br><br>  <br>;****************************************************************************<br>; Graphic resources<br>;****************************************************************************<br>  pltPath = pltDir+pltName<br>  wks     = gsn_open_wks(pltType ,pltPath) <br><br>  res                       = True     ; plot mods desired<br>  res@gsnDraw               = False<br>  res@gsnFrame              = False<br><br>  res@cnFillOn              = True     ; turn on color fill<br>  res@cnLinesOn             = False    ; turn of contour lines<br>  res@cnLineLabelsOn        = False    ; turn of contour line labels<br>  ;res@cnLevelSelectionMode  = "ManualLevels"    ; set manual contour levels<br><br>  res@lbLabelBarOn          = False<br>  <br>  <br><br> res@mpProjection                = "CylindricalEquidistant"       ; projection<br> res@mpLimitMode                 = "Corners"        ; method to zoom<br> res@mpLeftCornerLatF            = min(lat)<br> res@mpLeftCornerLonF            = min(lon)<br> res@mpRightCornerLatF           = max(lat)       <br> res@mpRightCornerLonF           = max(lon)<br> res@tfDoNDCOverlay              = True             ; do no<br> res@gsnAddCyclic                = False            ; regional data, don't add pt<br> <br><br>  res@tmYLLabelsOn          = False    ; turn off lat labels<br>  res@tmXBLabelsOn          = False    ; turn off lon labels<br><br>  res@mpDataBaseVersion     = "MediumRes"<br>  res@mpFillOn              = True    <br>  res@mpOutlineOn           = True<br>  ;res@mpFillColors          = (/-1,0,-1,-1/)       ; mask out all non-land areas<br>  res@mpFillOn                    = False        ; turn off gray fill<br><br>  res@mpFillDrawOrder       = "PostDraw"<br>  res@mpProjection          = "CylindricalEquidistant"<br>  <br>  <br>  ;Ressource column 1<br>   res1 = True<br>   res1 = res<br>   res1@cnFillPalette        = "WhiteYellowOrangeRed"<br>   res1@lbOrientation      = "Horizontal"        ; orientation<br>   ;res1@pmLabelBarWidthF   = 0.8<br>   ;res1@pmLabelBarHeightF  = 0.03<br>  <br>   res1@lbLabelBarOn       = False<br>   res1@lbBoxEndCapStyle      = "TriangleBothEnds"<br>   <br>   <br>   ;Ressource columns 2 & 3<br>     res2 = True<br>     res2 = res<br>     res2@cnFillPalette        = "MPL_YlGn"<br>     res2@lbOrientation      = "Horizontal"        ; orientation<br>     res2@pmLabelBarWidthF   = 1.1<br>     ;res2@pmLabelBarOrthogonalPosF  = -0.02<br>     res2@pmLabelBarParallelPosF  = -0.02<br><br>     res2@lbBoxEndCapStyle      = "TriangleBothEnds"<br>   <br><br>    ;Ressource column 4  <br>     res3 = True<br>     res3 = res<br>     res3@cnFillPalette        = "MPL_RdPu"<br>     res3@lbOrientation      = "Horizontal"        ; orientation<br>     ;res3@pmLabelBarWidthF   = 0.8<br>     <br>     res3@lbBoxEndCapStyle      = "TriangleBothEnds"<br>     <br><br>;****************************************************************************<br>; Generate panel entries<br>;****************************************************************************<br>                               ; only one time per file<br>  plot = new(nrow*ncol , "graphic")<br>  <br>  ;res@cnFillPalette        = "WhiteYellowOrangeRed"<br>    ;res1@tmYLLabelsOn          = True    ; turn off lat labels<br>   plot(0) = gsn_csm_contour_map(wks, t(0, :, :), res1)<br>   <br>     res2@tmYLLabelsOn          = False    ; turn off lat labels<br>   plot(1) = gsn_csm_contour_map(wks, t(1, :, :), res2)<br>   plot(2) = gsn_csm_contour_map(wks, t(2, :, :), res2)<br>   plot(3) = gsn_csm_contour_map(wks, t(3, :, :), res3)<br>   <br>   <br>    ;res1@tmYLLabelsOn          = True   <br>    plot(4) = gsn_csm_contour_map(wks, t(4, :, :), res1)<br>    res2@tmYLLabelsOn          = False    ; turn off lat labels<br>    plot(5) = gsn_csm_contour_map(wks, t(5, :, :),  res2)<br>    plot(6) = gsn_csm_contour_map(wks, t(6, :, :), res2)<br>    plot(7) = gsn_csm_contour_map(wks, t(7, :, :),  res3)<br>   <br>   <br>    res1@lbLabelBarOn          = True<br>    ;res2@tmYLLabelsOn          = True <br>    ;res1@tmXBLabelsOn          = True    ; turn off lon labels<br>    ;res2@tmXBLabelsOn          = True    ; turn off lon labels<br><br>    plot(8) = gsn_csm_contour_map(wks, t(8, :, :), res1)<br>    res2@tmYLLabelsOn          = False    ; turn off lat labels<br>    plot(9) = gsn_csm_contour_map(wks, t(9, :, :),  res2)<br>    <br>      res2@lbLabelBarOn       = True<br>    plot(10) = gsn_csm_contour_map(wks, t(10, :, :),res2)<br>    <br>    res3@lbLabelBarOn       = True<br>   plot(11) = gsn_csm_contour_map(wks, t(11, :, :), res3)<br> <br> <br> <br><br>  <br>;****************************************************************************<br>; Draw the panel and label bars<br>;****************************************************************************<br>  resP = True<br>  ;resP@gsnPanelMainString = "Land Only: Tmax: 2 Colorbars"<br>  resP@gsnMaximize = True<br>  resP@gsnFrame    = False<br> ;resP@gsnPanelRight  = 0.98<br>  ;resP@gsnPanelBottom = 0.045 <br>  <br>   resP@gsnPanelFigureStrings= (/"1","1","1","1","2","2","2","2","3","3","3","3"/) ; add strings to panel<br>    resP@amJust   = "BottomLeft"<br><br>  gsn_panel(wks,plot,(/nrow,ncol/),resP)  ; plots<br>  <br>  <br><br>  frame(wks)<br><br></div><div><br></div><div><br></div></div>