[ncl-talk] polygon now showing up in the polar plot

[Ipshita Majhi]

Thank you

On Fri, Oct 3, 2014 at 3:22 PM, Xin Xi <xin.xi30 at gmail.com> wrote:

> this is a draw order issue.
> set at the beginning:
> res at gsnDraw = False
> then
> plot =gsn_csm_contour_map(...)
> draw(plot)
> dum1=gsn_add_polygon(...)
>
>
> On Friday, October 3, 2014, Ipshita Majhi <ipmajhi at alaska.edu> wrote:
>
>> Dear All,
>>
>> I am trying to plot a polygon on a polar plot but it is not showing up
>>
>> ;************************************************
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> ;************************************************
>> begin
>> ;************************************************
>> ; read in netCDF file and access data
>> a =addfile("~/Documents/Snow_cover_IMS/nhsce_v01r01_19661004_20140602.nc
>> ","r");
>>   sce=a->snow_cover_extent(0,:,:)
>>  lat=a->latitude
>>  lon=a->longitude
>>
>>
>>  sce at lat2d=lat
>>  sce at lon2d=lon
>>                                  ;
>> ;************************************************
>> ; create polygon (remember the earth is a sphere)
>> ;************************************************
>>   xpt = new(20,float)
>>   ypt = new(20,float)
>>
>>   xpt(0:9)     = fspan(120,-120,10)  ; lon around edge
>>   xpt(10:14)   = -180                ; up 240
>>   xpt(15:19)   = 120                 ; down 120
>>
>>   ypt(0:9)     = 40                  ; lat around edge
>>   ypt(10:14)   = fspan(40.,90.,5)    ; lat up to pole
>>   ypt(15:19)   = fspan(90.,40.,5)    ; lat down from pole
>> ;************************************************
>> ; create plots
>> ;************************************************
>>   wks = gsn_open_wks("ps" ,"sce_polar_polygon")              ; create ps
>> plot
>>    res            = True                          ; plot mods desired
>>   res at gsnPolar   = "NH"                          ; specify the hemisphere
>>   plot = gsn_csm_contour_map_polar(wks,sce,res)    ; create the plot
>>
>>
>>  dum1 = gsn_add_polygon(wks,plot,xpt,ypt,False) ; add polygon
>> end
>> ;***********************************************
>> Best Regards
>> Ipshita
>>
>
>
> --
> Sent from Gmail Mobile
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141003/7dffe027/attachment.html