<div dir="ltr"><div>Dear NCL,<br><br>I want to extract a region from my polar plot and plot it. I read one example<br><br><a href="http://www.ncl.ucar.edu/Applications/Scripts/maponly_2.ncl">http://www.ncl.ucar.edu/Applications/Scripts/maponly_2.ncl</a><br><br></div>I applied it to polar plot but it is not working. There is approach of polygon and other things for polar plot is it possible to have a section extracted using min and max latitude. I would be grateful is somebody could assist me.<br><br><br><div><div><div>;*******************************************<br>;This is to extract different regions of snow cover and the region<br><br>;10/02/2014<br>;Western Eurasia 10W-30E , All of Eurasia, 0-190 E, 20-90 N, Southern Eurasia 0-190 E, 20-50 N, Himalayas 60-105 E, 30-45 N<br><br>;*******************************************<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>;*******************************************<br><br>a =addfile("~/Documents/Snow_cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc">nhsce_v01r01_19661004_20140602.nc</a>","r")<br><br>sce=a->snow_cover_extent(:,:,:)<br>lat=a->latitude<br>lon=a->longitude<br><br><br>sce@lat2d=lat<br>sce@lon2d=lon<br><br>;Now trying to get the region of interests<br>;Western Europe(we)<br><br>we=sce(40:60,10:30,:)<br><br>;All of Europe<br><br>print(dimsizes(we))<br><br>;Plotting the region and creating a box there<br><br>;*;************************************************<br> <br> wks = gsn_open_wks("ps","Snow_IMS_Polar_Color")              ; open a ps file<br> <br> res                            = True <br> <br> res@mpFillOn                   = False       ; turn off gray fill<br> res@mpOutlineBoundarySets      = "National"  ; turn on country boundaries<br> res@mpGeophysicalLineColor     = "Navy"      ; color of cont. outlines<br> res@mpGeophysicalLineThicknessF = 1.5         ; thickness of outlines            <br> res@gsnPolar                   ="NH"              <br> res@mpMaxLatF                  = 40        <br> res@mpMinLatF                  = 60<br> res@mpMaxLonF                  = 30<br> res@mpMinLonF                  = 10<br> <br> <br> plot = gsn_csm_contour_map_polar(wks,we,res)<br><br>;********************************************************<br><br></div><div>Best Regards<br>Ipshita<br></div><div><br><br></div></div></div></div>