[ncl-talk] Subregion for a polar plot

[Ipshita Majhi]

Dear NCL,

I want to extract a region from my polar plot and plot it. I read one
example

http://www.ncl.ucar.edu/Applications/Scripts/maponly_2.ncl

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.


;*******************************************
;This is to extract different regions of snow cover and the region

;10/02/2014
;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

;*******************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;*******************************************

a =addfile("~/Documents/Snow_cover_IMS/nhsce_v01r01_19661004_20140602.nc
","r")

sce=a->snow_cover_extent(:,:,:)
lat=a->latitude
lon=a->longitude


sce at lat2d=lat
sce at lon2d=lon

;Now trying to get the region of interests
;Western Europe(we)

we=sce(40:60,10:30,:)

;All of Europe

print(dimsizes(we))

;Plotting the region and creating a box there

;*;************************************************

  wks = gsn_open_wks("ps","Snow_IMS_Polar_Color")               ; open a ps
file

  res                             = True

  res at mpFillOn                    = False        ; turn off gray fill
  res at mpOutlineBoundarySets       = "National"   ; turn on country
boundaries
  res at mpGeophysicalLineColor      = "Navy"       ; color of cont. outlines
  res at mpGeophysicalLineThicknessF = 1.5          ; thickness of
outlines
  res at gsnPolar                    ="NH"
  res at mpMaxLatF                   = 40
  res at mpMinLatF                   = 60
  res at mpMaxLonF                   = 30
  res at mpMinLonF                   = 10


  plot = gsn_csm_contour_map_polar(wks,we,res)

;********************************************************

Best Regards
Ipshita
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141003/aaa19dc7/attachment.html