<div dir="ltr"><div>Dear All,<br><br>I am trying to plot a polygon on a polar plot but it is not showing up<br><br>;************************************************<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot; <br>;************************************************<br>begin<br>;************************************************<br>; read in netCDF file and access data<br>a =addfile(&quot;~/Documents/Snow_cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc">nhsce_v01r01_19661004_20140602.nc</a>&quot;,&quot;r&quot;);<br>  sce=a-&gt;snow_cover_extent(0,:,:)<br> lat=a-&gt;latitude<br> lon=a-&gt;longitude<br><br><br> sce@lat2d=lat<br> sce@lon2d=lon<br>                                 ; <br>;************************************************<br>; create polygon (remember the earth is a sphere)<br>;************************************************<br>  xpt = new(20,float)<br>  ypt = new(20,float)<br><br>  xpt(0:9)     = fspan(120,-120,10)  ; lon around edge<br>  xpt(10:14)   = -180                ; up 240<br>  xpt(15:19)   = 120                 ; down 120<br><br>  ypt(0:9)     = 40                  ; lat around edge<br>  ypt(10:14)   = fspan(40.,90.,5)    ; lat up to pole<br>  ypt(15:19)   = fspan(90.,40.,5)    ; lat down from pole<br>;************************************************<br>; create plots<br>;************************************************<br>  wks = gsn_open_wks(&quot;ps&quot; ,&quot;sce_polar_polygon&quot;)              ; create ps plot     <br>   res            = True                          ; plot mods desired<br>  res@gsnPolar   = &quot;NH&quot;                          ; specify the hemisphere<br>  plot = gsn_csm_contour_map_polar(wks,sce,res)    ; create the plot<br><br><br> dum1 = gsn_add_polygon(wks,plot,xpt,ypt,False) ; add polygon<br>end<br>;***********************************************<br></div>Best Regards<br>Ipshita<br></div>