[ncl-talk] In significance plot stippling dot not showing up

Ipshita Majhi ipmajhi at alaska.edu
Thu Sep 14 12:14:45 MDT 2017


Thank you , I worked and I am grateful for it.

Best Regards
Ipshita

On Thu, Sep 14, 2017 at 7:27 AM, Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Ipshita,
> I think you want to set
>  opt at gsnShadeHigh = 17
> as opposed to
>   opt at gsnShadeLow = 17
>
> as you want to shade above your given high contour of 90. You might want
> to set res2 at cnLinesOn = False as well.
>
> If that does not fix it let ncl-talk know.
> Adam
>
>
>
> On Wed, Sep 13, 2017 at 8:39 PM, Ipshita Majhi <ipmajhi at alaska.edu> wrote:
>
>> Dear NCL,
>>
>> I am plotting correlation and significance in the same plot. I am not
>> getting the dots , I am getting a contour instead and I am not sure where I
>> am making a mistake. If someone could guide me then I will be
>> grateful.Listed below is my code.
>>
>>
>> ;*******************************************
>> 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"
>> ;*******************************************
>>
>> ;*******************************************
>> MJ_1982_2014=asciiread("/Volumes/Ipshita_extra/Documents/
>> NCL_files/SST/Monthly/MJ.txt",(/33,1/),"float")
>> MJ_1982_2012=MJ_1982_2014(0:31,0)
>>
>> print(dimsizes(MJ_1982_2012))
>>
>> ;*************************************************
>> ;Reading in SWE data
>> ;*************************************************
>> a = addfile("monthly_82_13.nc","r")
>> b = addfile("GlobSnow_SWE_L3A_20130507_v2.0.nc","r")
>>
>> swe  = a->swe
>> time = a->time
>>
>> lat  = b->lat
>> lon  = b->lon
>> swe at lat2d = lat
>> swe at lon2d = lon
>>
>> print(dimsizes(lat))
>> print(dimsizes(lon))
>>
>> swe_FillValue=-9.99e+08
>>
>> utc_date = cd_calendar(time, 0)
>> year=utc_date(:,0)
>> month=utc_date(:,1)
>>
>> ;********************************************
>> ;This is to extract 1982-2012 of swe data
>> ;********************************************
>> swe_82_12=swe(0:371,:,:)
>>
>> ;********************************************
>> ;*******************************************
>> ;Extracting monthly data for sce
>> ;*******************************************
>>
>> March_swe = swe(0:31,:,:)
>>
>> March_swe at _FillValue = -9.99e+08
>>
>> do nyr=0,371,12
>>
>> March_swe(nyr/12,:,:) =swe_82_12(nyr+2,:,:)
>>
>> end do
>> ;=================================================
>> March_swe at _FillValue = -9.99e+08
>> ;**********************************************
>>
>> March_dt = dtrend_msg(March_swe&time,March_swe(lat|:,lon|:,time|:),True
>> ,False)
>>
>> printMinMax(March_dt,False)
>> MJ_dt=dtrend(MJ_1982_2012,False)
>> ;****************************************
>> ;Calculating Correlation
>> ;***************************************
>>
>> corr_March=escorc(March_dt,MJ_dt)
>>
>> copy_VarCoords(swe(0,:,:),corr_March)
>> copy_VarAtts(swe(0,:,:),corr_March)
>>
>> ;**********************************************
>> ;Calculating significance
>> ;**********************************************
>>
>> prob_March=rtest(corr_March,31,0)
>>
>> sig_March =100*(1-prob_March)
>>
>>  copy_VarCoords(swe(0,:,:),sig_March)
>>  copy_VarAtts(swe(0,:,:),sig_March)
>>
>>  sig_March = mask(sig_March, sig_March.ge.100,False)
>>  sig_March = mask(sig_March, sig_March.lt.90,False)
>>
>>
>> ;************************************************
>> ; create plot
>> ;************************************************
>>
>>  wks = gsn_open_wks("x11","swe_march_mj_cor_sig")               ;
>>
>>  gsn_define_colormap(wks, "BlueYellowRed")
>>
>>
>>  ;************************************************
>>     res                     = True              ; Plot modes desired.
>>     res at gsnMaximize         = True              ; Maximize plot
>>     res at gsnDraw             = False            ; don't draw
>>     res at gsnFrame            = False            ; don't advance frame
>>     res at cnFillOn            = True              ; color plot desired
>>     res at cnLinesOn           = False             ; turn off contour lines
>>     res at pmLabelBarWidthF    = 0.9             ; make wider
>>     res at pmLabelBarHeightF   = 0.1               ; default is taller
>>     res at lbLabelFontHeightF  = .018              ; default is HUGE
>>     res at lbLabelBarOn        =True           ; turn off individual cb's
>>
>>
>> ;*******************************************
>> ; georeferencing: plot on polar projection
>> ;*******************************************
>>     res at trGridType   = "TriangularMesh"        ; allow missing
>> coordinates
>>     res at gsnAddCyclic = False
>>     res at cnFillMode   = "RasterFill"
>>     res at cnInfoLabelOn       = False       ; turn off info label
>>     res at gsnPolar   = "NH"                          ; specify the
>> hemisphere
>>     res at mpMinLatF  = 40
>>     res at mpMaxLatF  = 70
>>     res at mpMinLonF  = -30
>>     res at mpMinLonF  = -180
>>     plot = gsn_csm_contour_map_polar(wks,corr_March, res)
>>  ;*********************************************
>>
>>   res2=True
>>   res2 at gsnDraw              = False           ; Do not draw plot
>>   res2 at gsnFrame             = False           ; Do not advance frome
>>   res2 at cnFillOn             = True
>>   res2 at cnMonoFillColor      = True
>>   res2 at cnMonoFillPattern    = False
>>   res2 at lbLabelBarOn        = False                            ; turn off
>> label bar
>>   res2 at cnLevelSelectionMode = "ManualLevels"   ; manually specify
>> contour levels
>>   res2 at cnMinLevelValF       = 90             ; min level
>>   res2 at cnMaxLevelValF       = 100             ; max level
>>   res2 at cnLevelSpacingF      =  1              ; contour interval
>>   res2 at gsnSpreadColors      = False
>>
>>
>>   plot2=gsn_csm_contour(wks,sig_March, res2) ; contours are at 90,100
>>
>>
>>   opt     = True
>>   opt at gsnShadeFillType = "pattern"
>>   opt at gsnShadeLow = 17
>>
>> plot2=gsn_contour_shade(plot2,-999.,90.,opt) ;
>>
>>  overlay (plot, plot2)
>>
>>  draw (plot)
>>  frame(wks)
>>
>>
>>
>> ;=================================================
>>
>>
>>
>>
>>
>> ;
>>
>>
>> --
>>
>>
>> Ipshita Majhi
>> PhD Candidate
>> University of Alaska , Fairbanks
>> Atmospheric Science Department
>> (907)978-4220 <(907)%20978-4220> ipmajhi at alaska.edu
>>
>> ************************************************************
>> ************************************************
>> Happiness comes when we stop complaining about the troubles we have and
>> offer thanks for all the troubles we don't have. Life is a gift!
>> -Sri Sri
>> ************************************************************
>> ****************************************
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726 <(303)%20497-1726>
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>



-- 


Ipshita Majhi
PhD Candidate
University of Alaska , Fairbanks
Atmospheric Science Department
(907)978-4220 ipmajhi at alaska.edu

************************************************************************************************************
Happiness comes when we stop complaining about the troubles we have and
offer thanks for all the troubles we don't have. Life is a gift!
-Sri Sri
****************************************************************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170914/ea41c5c5/attachment.html>


More information about the ncl-talk mailing list