<div dir="ltr">Hi all,<div><br></div><div>I have three plots and want to create one plot with three panels how can I do that. I've attached the figures here and also the script for them.</div><div><br></div><div><br></div><div><div><br></div><div> wks = gsn_open_wks("eps","atm_sfc_fine")</div><div>;gsn_define_colormap(wks,"BlAqGrYeOrRe")</div><div>;gsn_define_colormap(wks,"BlWhRe")</div><div> plot = new(2,graphic)   </div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>  res                     = True               ; plot mods desired</div><div>  res@cnFillOn            = True              ; turn on color   </div><div>  res@gsnSpreadColors     = True               ; use full color map</div><div>  res@cnLinesOn           = False              ; no contour lines</div><div>  res@cnLineLabelsOn      = False              ; no line labels</div><div><br></div><div><br></div><div> ; res@gsnAddCyclic        = False </div><div>  res@gsnDraw             = False           ; don't draw</div><div>  res@gsnFrame            = False           ; don't advance frame</div><div>  res@cnInfoLabelOn       = False           ; turn off cn info label</div><div>  res@cnFillOn            = True            ; turn on color</div><div>  res@gsnSpreadColors     = True            ; spread out color table</div><div>  res@gsnSpreadColorStart = 2               ; color to startt</div><div><br></div><div>  res@lbLabelBarOn        = False           ; turn off individual cb's</div><div><br></div><div>;  res@cnLevelSelectionMode =  "AutomaticLevels"   ; AutomaticLevels</div><div><br></div><div><br></div><div>  res@cnLevelSelectionMode =  "ManualLevels"   ; AutomaticLevels</div><div>  res@cnMinLevelValF       = -40.              ; min level</div><div>  res@cnMaxLevelValF       =  40.               ; max level</div><div>  res@cnLevelSpacingF      =  4               ; interval</div><div><br></div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>res@gsnLeftString       = "Solar radiation absorption in the atmosphere"         ; change left string</div><div>plot(0) = gsn_csm_contour_map_ce(wks,Absp,res)</div><div>res@gsnLeftString       = "Surface net flux"         ; change left string</div><div>plot(1) = gsn_csm_contour_map_ce(wks,sfc,res)</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>;************************************************</div><div>; create panel</div><div>;************************************************</div><div>  resP                  = True                    ; modify the panel plot</div><div><br></div><div>  resP@gsnPanelLabelBar = True                    ; add common colorbar</div><div>  resP@gsnPanelFigureStrings = (/"a)","b)"/)      ; add strings to panel</div><div>;resP@txString         = "Annual SDE changes at the surface"     </div><div>resP@lbTitleString    = "W/m^2"                ; title string</div><div>resP@lbTitlePosition  = "Right"              ; title position</div><div>resP@lbTitleFontHeightF= .02                ; make title smaller</div><div>resP@lbTitleDirection = "Across"             ; title direction</div><div>gsn_panel(wks,plot,(/1,2/),resP)  </div><div><br></div><div><br></div><div>delete(plot)</div><div>delete(res)</div><div>delete(resP)</div><div>delete(wks)</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><div><br></div><div><br></div><div>;************************************************</div><div>; plotting parameters</div><div>;************************************************</div><div> wks   = gsn_open_wks ("eps","sfc_lat_summer")                  ; send graphics to PNG file</div><div><br></div><div> res                   = True                       ; plot mods desired</div><div>;res@tiMainString      = "Vertical profile of atmospheric heating rate"             ; add title</div><div>res@tiYAxisString = "Surface net flux (W/m^2)"  ; xaxis string</div><div>res@tiXAxisString = "Latitude";Pressure levels (hPa)"</div><div><br></div><div>res@trYMinF     = -10;0<span style="white-space:pre">          </span>; set minimum Y-axis value</div><div>res@trYMaxF     = 1;12250;1200<span style="white-space:pre">           </span>; set maximum Y-axis value</div><div>;res@trXMinF     = -0.5      ; set minimum X-axis value</div><div>;res@trXMaxF     = 1.2      ; set maximum X-axis value</div><div>; res@trYReverse = True                     ; reverse Y-axis</div><div>res@xyDashPattern  = 0                         ; Make curves all solid</div><div>res@xyLineColors      = "blue"          ; change line color</div><div>res@xyLineThicknesses = 4.0</div><div> plot  = gsn_csm_xy (wks,lat,sfc_summer2,res) ; create plot</div><div>;************************************************</div><div>delete(wks)</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div>Thanks,</div><div>Ana</div><div><br></div></div>