<div dir="ltr"><br><div class="gmail_quote"><div dir="ltr"><br><div class="gmail_quote"><br><br><div dir="ltr"><div>Hi,<br></div><div>I don't know how to get rid of white spaces around the boundary. Could someone help me to fix this? I have attached a output figure for your review.<br></div><div><br></div><div>Thanks<br></div><div> <br></div><div><span class=""><br>latS  = 40.<br> latN  = 49.<br> lonL  = -116.<br></span> lonR  = -90.<br><br><br></div><div><div class="h5">wks = gsn_open_wks("pdf","eof")<br> gsn_define_colormap(wks,"BlWhRe")      ; choose colormap<br> plot = new(neof,graphic)               ; create graphic array<br>                                         ; only needed if paneling<br>; EOF patterns<br><br> res                     = True<br> res@gsnDraw             = False       ; don't draw yet<br> res@gsnFrame            = False       ; don't advance frame yet<br><br>;---This resource not needed in V6.1.0<br> res@gsnSpreadColors     = True        ; spread out color table<br><br> res@gsnAddCyclic        = False       ; plotted dataa are not cyclic<br><br> res@mpFillOn            = False       ; turn off map fill<br>res@mpLimitMode          = "LatLon"<br> res@mpMinLatF           = latS        ; zoom in on map<br> res@mpMaxLatF           = latN<br> res@mpMinLonF           = lonL<br> res@mpMaxLonF           = lonR<br> res@cnFillOn            = True        ; turn on color fill<br> res@cnLinesOn           = False       ; True is default<br>;res@cnLineLabelsOn      = False       ; True is default<br> res@lbLabelBarOn        = False       ; turn off individual lb's<br> res@mpPerimOn             = True                   ; draw box around map<br>res@mpGeophysicalLineThicknessF = 3.0<br>res@mpGeophysicalLineColor = "Black"; (/22/)<br>;res@mpNationalLineThicknessF = 3.0<br>res@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; add state boundaries<br>res@mpNationalLineColor = res@mpGeophysicalLineColor<br>res@mpUSStateLineThicknessF = 3.0<br>res@mpUSStateLineColor = res@mpGeophysicalLineColor<br><br><br>                                         ; set symmetric plot min/max<br> symMinMaxPlt(eof, 16, False, res)      ; contributed.ncl<br>; panel plot only resources<br> resP                    = True        ; modify the panel plot<br> resP@gsnMaximize        = False        ; large format<br> resP@gsnPanelLabelBar   = True        ; add common colorbar<br> resP@lbLabelAutoStride  = True        ; auto stride on labels<br><br> yStrt                   = yyyymm(0)/100<br> yLast                   = yyyymm(nyrs-1)/100<br>; resP@txString           = "SLP: "+season+": "+yStrt+"-"+yLast<br><br>;*******************************************<br>; first plot<br>;*******************************************<br> do n=0,neof-1<br>    res@gsnLeftString = "EOF "+(n+1)<br>    res@gsnRightString = sprintf("%5.1f", eof@pcvar(n)) +"%"<br>    plot(n)=gsn_csm_contour_map_ce(wks,eof(n,:,:),res)<br> end do<br>;*******************************************<br>; second plot<br>;*******************************************<br>; EOF time series [bar form]<br><br> rts          = True<br> rts@gsnDraw  = False      ; don't draw yet<br> rts@gsnFrame = False      ; don't advance frame yet<br> rts@gsnScale = True       ; force text scaling<br><br>; these four rtsources allow the user to stretch the plot size, and<br>; decide exactly where on the page to draw it.<br><br> rts@vpHeightF = 0.40       ; Changes the aspect ratio<br> rts@vpWidthF = 0.85<br> rts@vpXF     = 0.10       ; change start locations<br> rts@vpYF     = 0.75       ; the plot<br><br><br> rts@tiYAxisString = "mm/day"    rts@gsnYRefLine          = 0.             ; reference line<br> rts@gsnXYBarChart        = True           ; create bar chart<br> rts@gsnAboveYRefLineColor = "red"          ; above ref line fill red<br> rts@gsnBelowYRefLineColor = "blue"         ; below ref line fill blue<br><br>; panel plot only resources<br> rtsP                     = True           ; modify the panel plot<br> rtsP@gsnMaximize         = True           ; large format<br> rtsP@txString            = "Prec: "+season+": "+yStrt+"-"+yLast<br><br> year = yyyymm/100<br><br>; create individual plots<br> do n=0,neof-1<br>    rts@gsnLeftString = "EOF "+(n+1)<br>    rts@gsnRightString = sprintf("%5.1f", eof@pcvar(n)) +"%"<br>    plot(n) = gsn_csm_xy (wks,year,eof_ts(n,:),rts)<br> end do<br> gsn_panel(wks,plot,(/neof,1/),rtsP)    ; now draw as one plot<br>end<br>    <br><br></div></div></div>
</div><br></div>
</div><br></div>