<div dir="ltr">Hi Ipshita,<div>I think you want to set</div><div><span style="font-size:12.8px"> opt@gsnShadeHigh = 17</span></div><div><span style="font-size:12.8px">as opposed to </span></div><div><span style="font-size:12.8px"> opt@gsnShadeLow = 17 </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">as you want to shade above your given high contour of 90. </span><span style="font-size:12.8px">You might want to set res2@cnLinesOn = False as well. </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">If that does not fix it let ncl-talk know.</span></div><div><span style="font-size:12.8px">Adam</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 13, 2017 at 8:39 PM, 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">Dear NCL,<div><br></div><div>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.</div><div><br></div><div><br></div><div><div>;*****************************<wbr>**************</div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"</div><div>;*****************************<wbr>**************</div><div><br></div><div>;*****************************<wbr>**************</div><div>MJ_1982_2014=asciiread("/<wbr>Volumes/Ipshita_extra/<wbr>Documents/NCL_files/SST/<wbr>Monthly/MJ.txt",(/33,1/),"<wbr>float")</div><div>MJ_1982_2012=MJ_1982_2014(0:<wbr>31,0)</div><div><br></div><div>print(dimsizes(MJ_1982_2012))</div><div><br></div><div>;*****************************<wbr>********************</div><div>;Reading in SWE data</div><div>;*****************************<wbr>********************</div><div>a = addfile("<a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a>","r"<wbr>)</div><div>b = addfile("<a href="http://GlobSnow_SWE_L3A_20130507_v2.0.nc" target="_blank">GlobSnow_SWE_L3A_<wbr>20130507_v2.0.nc</a>","r")</div><div><br></div><div>swe = a->swe</div><div>time = a->time</div><div><br></div><div>lat = b->lat</div><div>lon = b->lon</div><div>swe@lat2d = lat</div><div>swe@lon2d = lon</div><div><br></div><div>print(dimsizes(lat))</div><div>print(dimsizes(lon))</div><div><br></div><div>swe_FillValue=-9.99e+08</div><div><br></div><div>utc_date = cd_calendar(time, 0)</div><div>year=utc_date(:,0)</div><div>month=utc_date(:,1)</div><div><br></div><div>;*****************************<wbr>***************</div><div>;This is to extract 1982-2012 of swe data</div><div>;*****************************<wbr>***************</div><div>swe_82_12=swe(0:371,:,:)</div><div><br></div><div>;*****************************<wbr>***************</div><div>;*****************************<wbr>**************</div><div>;Extracting monthly data for sce</div><div>;*****************************<wbr>**************</div><div><br></div><div>March_swe = swe(0:31,:,:)</div><div><br></div><div>March_swe@_FillValue = -9.99e+08</div><div><br></div><div>do nyr=0,371,12</div><div><br></div><div>March_swe(nyr/12,:,:) =swe_82_12(nyr+2,:,:)</div><div><br></div><div>end do</div><div>;=============================<wbr>====================</div><div>March_swe@_FillValue = -9.99e+08</div><div>;*****************************<wbr>*****************</div><div><br></div><div>March_dt = dtrend_msg(March_swe&time,<wbr>March_swe(lat|:,lon|:,time|:),<wbr>True,False)</div><div><br></div><div>printMinMax(March_dt,False)</div><div>MJ_dt=dtrend(MJ_1982_2012,<wbr>False)</div><div>;*****************************<wbr>***********</div><div>;Calculating Correlation</div><div>;*****************************<wbr>**********</div><div><br></div><div>corr_March=escorc(March_dt,MJ_<wbr>dt)</div><div><br></div><div>copy_VarCoords(swe(0,:,:),<wbr>corr_March) </div><div>copy_VarAtts(swe(0,:,:),corr_<wbr>March)</div><div><br></div><div>;*****************************<wbr>*****************</div><div>;Calculating significance</div><div>;*****************************<wbr>*****************</div><div><br></div><div>prob_March=rtest(corr_March,<wbr>31,0)</div><div><br></div><div>sig_March =100*(1-prob_March)</div><div><br></div><div> copy_VarCoords(swe(0,:,:),<wbr>sig_March) </div><div> copy_VarAtts(swe(0,:,:),sig_<wbr>March)</div><div><br></div><div> sig_March = mask(sig_March, sig_March.ge.100,False) </div><div> sig_March = mask(sig_March, sig_March.lt.90,False)</div><div><br></div><div><br></div><div>;*****************************<wbr>*******************</div><div>; create plot</div><div>;*****************************<wbr>*******************</div><div> </div><div> wks = gsn_open_wks("x11","swe_march_<wbr>mj_cor_sig") ; </div><div><br></div><div> gsn_define_colormap(wks, "BlueYellowRed") </div><div> </div><div> </div><div> ;****************************<wbr>********************</div><div> res = True ; Plot modes desired.</div><div> res@gsnMaximize = True ; Maximize plot</div><div> res@gsnDraw = False ; don't draw</div><div> res@gsnFrame = False ; don't advance frame </div><div> res@cnFillOn = True ; color plot desired</div><div> res@cnLinesOn = False ; turn off contour lines</div><div> res@pmLabelBarWidthF = 0.9 ; make wider</div><div> res@pmLabelBarHeightF = 0.1 ; default is taller</div><div> res@lbLabelFontHeightF = .018 ; default is HUGE</div><div> res@lbLabelBarOn =True ; turn off individual cb's</div><div><br></div><div> </div><div>;*****************************<wbr>**************</div><div>; georeferencing: plot on polar projection</div><div>;*****************************<wbr>**************</div><div> res@trGridType = "TriangularMesh" ; allow missing coordinates</div><div> res@gsnAddCyclic = False</div><div> res@cnFillMode = "RasterFill"</div><div> res@cnInfoLabelOn = False ; turn off info label</div><div> res@gsnPolar = "NH" ; specify the hemisphere</div><div> res@mpMinLatF = 40</div><div> res@mpMaxLatF = 70</div><div> res@mpMinLonF = -30</div><div> res@mpMinLonF = -180</div><div> plot = gsn_csm_contour_map_polar(wks,<wbr>corr_March, res) </div><div> ;****************************<wbr>*****************</div><div> </div><div> res2=True</div><div> res2@gsnDraw = False ; Do not draw plot</div><div> res2@gsnFrame = False ; Do not advance frome </div><div> res2@cnFillOn = True</div><div> res2@cnMonoFillColor = True</div><div> res2@cnMonoFillPattern = False</div><div> res2@lbLabelBarOn = False ; turn off label bar</div><div> res2@cnLevelSelectionMode = "ManualLevels" ; manually specify contour levels</div><div> res2@cnMinLevelValF = 90 ; min level</div><div> res2@cnMaxLevelValF = 100 ; max level</div><div> res2@cnLevelSpacingF = 1 ; contour interval</div><div> res2@gsnSpreadColors = False</div><div> </div><div><br></div><div> plot2=gsn_csm_contour(wks,sig_<wbr>March, res2) ; contours are at 90,100 </div><div><br></div><div><br></div><div> opt = True</div><div> opt@gsnShadeFillType = "pattern"</div><div> opt@gsnShadeLow = 17 </div><div><br></div><div>plot2=gsn_contour_shade(plot2,<wbr>-999.,90.,opt) ; </div><div> </div><div> overlay (plot, plot2)</div><div><br></div><div> draw (plot)</div><div> frame(wks) </div><div><br></div><div> </div><div><br></div><div>;=============================<wbr>====================</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>;</div><div><br></div><div><br></div>-- <br><div class="m_-665176293034648781gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><div><div><div><br></div><div><br></div><div>Ipshita Majhi<br></div>PhD Candidate<br></div>University of Alaska , Fairbanks<br></div>Atmospheric Science Department<br></div><a href="tel:(907)%20978-4220" value="+19079784220" target="_blank">(907)978-4220</a> <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div><div dir="ltr"><br></div><div dir="ltr">******************************<wbr>******************************<wbr>******************************<wbr>******************</div><div dir="ltr"><span style="color:rgb(29,33,41);font-family:Helvetica,Arial,sans-serif;font-size:14px">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!</span><br style="color:rgb(29,33,41);font-family:Helvetica,Arial,sans-serif;font-size:14px"><span style="color:rgb(29,33,41);font-family:Helvetica,Arial,sans-serif;font-size:14px">-Sri Sri</span><br></div><div dir="ltr"><span style="color:rgb(29,33,41);font-family:Helvetica,Arial,sans-serif;font-size:14px">******************************<wbr>******************************<wbr>******************************<wbr>**********</span></div></div></div></div></div></div></div></div>
</div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>