<div dir="ltr">Hi Ipshita,<div><br></div><div>The routine gsn_csm_contour_map_polar is a specialized routine that draws a circular stereographic projection with the center of the plot at either the north or the south pole. The only lat/lon limit resource that makes sense in this context is mpMinLatF, which can limit the plot from its maximal possible extent which I think is the equator.</div><div><br></div><div>If you want to limit the plot to a specific region you should use a different routine, probably the more general gsn_csm_contour_map. And instead of limiting using the lat and lon resources (which for non-rectangular projections will probably include a larger area than you want), you might want to try</div><div>using the "Corners" limit mode. It might not be what you really want, but here is the same projection using the general map routine:</div><div><br></div>res@mpCenterLatF = 90<br>res@mpProjection = "stereographic"<br>plot = gsn_csm_contour_map(wks,we,res)<br><br><br>Add these resources to restrict the plot to a region:<br><br> res@mpLimitMode = "corners"<br> res@mpLeftCornerLatF = 40<br> res@mpLeftCornerLonF = 10<br> res@mpRightCornerLatF = 60<br> res@mpRightCornerLonF = 30<br><div><br></div><div>You will probably need to play with these to get what you are really looking for.</div><div> -dave</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 3, 2014 at 11:00 AM, Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">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" target="_blank">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<span class="HOEnZb"><font color="#888888"><br>Ipshita<br></font></span></div><div><br><br></div></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>