<div dir="ltr"><div style="font-size:14px">Hi Ipshita,</div><div style="font-size:14px"><br></div><div style="font-size:14px">I think the sample size of your data for doing the significant test is small (sX = 7 and sY = 4). </div><div style="font-size:14px">If you set the minimum contour intervals for showing the significant test to be 95, probably none of the data pass this confidence level.</div><div style="font-size:14px"><br></div><div style="font-size:14px"><div style="font-size:14px"> res1@cnLevelSelectionMode = "ManualLevels" ; manually specify contour levels</div><div style="font-size:14px"> res1@cnMinLevelValF = 95 ; min level</div><div style="font-size:14px"> res1@cnMaxLevelValF = 100 ; max level</div><div style="font-size:14px"> res1@cnLevelSpacingF = 1 ; contour interval</div><div style="font-size:14px"><br></div><div style="font-size:14px">You may use printMinMax(<span style="font-size:14px">alpha,True) to see the maximum value of alpha.</span></div><div style="font-size:14px"><span style="font-size:14px"><br></span></div><div style="font-size:14px"><span style="font-size:14px">Cheers,</span></div><div style="font-size:14px">Hoffman</div><div style="font-size:14px"><br></div><div style="font-size:14px"><span style="font-size:14px"><br></span></div><div style="font-size:14px"><span style="font-size:14px"><br></span></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-11-23 2:38 GMT+01:00 Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear NCl,</div><div><br></div><div>I am trying to overlay two plots. But the significance plot does not show up, and there are no errors as such. I would be grateful if you could guide me on this.</div><div><br></div><div><br></div><div><br></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>;Reading in Storm track data by Peter Beinek 2016</div><div>;*****************************<wbr>********************</div><div> a=addfile("~/Documents/NCL_<wbr>files/Storm_tracks/<a href="http://cyclcnt_monthly_80-14.nc" target="_blank">cyclcnt_<wbr>monthly_80-14.nc</a>","r")</div><div><br></div><div> vNames = getfilevarnames(a) ; get variable names of all groups on file</div><div> </div><div> ;print (vNames) ; print all variable names on file</div><div> </div><div> dens=a->dens</div><div> time=a->time</div><div> lat=a->lat</div><div> lon=a->lon</div><div>;*****************************<wbr>*********************</div><div>;Doing composite analysis for the years</div><div>;with high and low monsoon</div><div>;High rainfall years: 1990,1994,1999,2000,2001,2008,<wbr>2013</div><div>;Low rainfall years: 1987,1991,2009,2012</div><div>;*****************************<wbr>**********************</div><div><br></div><div>;*****************************<wbr>**********************</div><div>;Now extracting the years with high rainfall</div><div>;*****************************<wbr>**********************</div><div><br></div><div>;********************</div><div>;Extracting for 1990</div><div>;********************</div><div>st1990=dens({119:130},:,:)</div><div><br></div><div>;*******************</div><div>;Extracting for 1994</div><div>;*******************</div><div>st1994=dens({167:178},:,:)</div><div>;********************</div><div>;1999</div><div>;*******************</div><div>st1999=dens({227:238},:,:)</div><div>;********************</div><div>;2000</div><div>;*********************</div><div>st2000=dens({239:250},:,:)</div><div>;*********************</div><div>;2001</div><div>;********************</div><div>st2001=dens({251:262},:,:)</div><div>;********************</div><div>;2008</div><div>;********************</div><div>st2008=dens({335:346},:,:)</div><div>;********************</div><div>;2013</div><div>;*********************</div><div>st2013=dens({395:406},:,:)</div><div>;*********************</div><div><br></div><div>;*****************************<wbr>************</div><div>;This averaging all the high years</div><div>;*****************************<wbr>************</div><div><br></div><div>st_highyrs= [/st1990,st1994,st1999,st2000,<wbr>st2001,st2008,st2013/] </div><div><br></div><div> ListSetType(st_highyrs, "join") </div><div> jab = st_highyrs[:] </div><div>; print(dimsizes(jab)) </div><div> zAvg = dim_avg_n_Wrap(jab,0) ; ==> zAvg(12,nlat,nlon)</div><div>; print(dimsizes(zAvg))</div><div><br></div><div>zAvg!1="lat" ; name dimensions </div><div>zAvg!2="lon" </div><div>zAvg&lat = dens&lat ; assign coordinate values and </div><div>zAvg&lon =dens&lon ; units attributes </div><div> copy_VarAtts(dens,zAvg) </div><div> copy_VarCoords_1(dens,zAvg)</div><div> zAvg&lat@units = "degrees_north" </div><div> zAvg&lon@units = "degrees_east" </div><div><br></div><div>;*****************************<wbr>****************</div><div>;Now reading in the years with low rainfall data </div><div>;For storm track 1987,1991,2009,2012</div><div>;*****************************<wbr>*****************</div><div><br></div><div>;********************</div><div>;1987</div><div>;*********************</div><div>st1987=dens({83:94},:,:)</div><div>;*********************</div><div>;********************</div><div>;1991</div><div>;*********************</div><div>st1991=dens({131:142},:,:)</div><div>;*********************</div><div>;********************</div><div>;2009</div><div>;*********************</div><div>st2009=dens({335:346},:,:)</div><div>;*********************</div><div>;*********************</div><div>;2012</div><div>;*********************</div><div>st2012=dens({383:394},:,:)</div><div>;*********************</div><div>;*****************************<wbr>************</div><div>;This averaging all the low years</div><div>;*****************************<wbr>************</div><div><br></div><div>st_lowyrs= [/st1987,st1991,st2009,st2012/<wbr>] </div><div><br></div><div>ListSetType(st_lowyrs, "join") </div><div>jab1 = st_lowyrs[:] </div><div>;print(dimsizes(jab1))</div><div><br></div><div>;*****************************<wbr>*************</div><div>;This is to calculate the std</div><div>;*****************************<wbr>*************</div><div> zStd1=dim_stddev_n_Wrap( jab1, 0 ) </div><div><br></div><div>zAvg1 = dim_avg_n(jab1,0) ; ==> zAvg(12,nlat,nlon)</div><div><br></div><div>zAvg1!1="lat" ; name dimensions </div><div>zAvg1!2="lon" </div><div>zAvg1&lat = dens&lat ; assign coordinate values and </div><div>zAvg1&lon =dens&lon ; units attributes </div><div> copy_VarAtts(dens,zAvg1) </div><div> copy_VarCoords_1(dens,zAvg1)</div><div> zAvg1&lat@units = "degrees_north" </div><div> zAvg1&lon@units = "degrees_east" </div><div><br></div><div><br></div><div>;*****************************<wbr>*************</div><div>;Difference between high and low years</div><div>;*****************************<wbr>*************</div><div><br></div><div>zDiff=zAvg-zAvg1</div><div><br></div><div>zDiff!1="lat" ; name dimensions </div><div>zDiff!2="lon" </div><div>zDiff&lat = dens&lat ; assign coordinate values and </div><div>zDiff&lon =dens&lon ; units attributes </div><div> copy_VarAtts(dens,zDiff) </div><div> copy_VarCoords_1(dens,zDiff)</div><div> zDiff&lat@units = "degrees_north" </div><div> zDiff&lon@units = "degrees_east" </div><div> zDiff= smth9(zDiff, 0.5, 0.25, True)</div><div><br></div><div> </div><div> ;****************************<wbr>*************</div><div> ;Calculating significance test</div><div> siglvl = 0.05</div><div> aveX =dim_avg_n_Wrap(jab,0) </div><div> aveY = dim_avg_n_Wrap(jab1,0) </div><div> varX = dim_variance_n_Wrap( jab, 0 ) </div><div> varY = dim_variance_n_Wrap( jab1, 0 ) </div><div> sX = 7 </div><div> sY = 4</div><div> </div><div> alpha = 100.*(1. - ttest(aveX,varX,sX, aveY,varY,sY, True, False)) </div><div> printMinMax(alpha,False)</div><div> </div><div>; sig_ alpha = mask( alpha, alpha.eq.100,False) </div><div>;sig_ alpha = mask( alpha, alpha.le.90,False)</div><div> </div><div> </div><div>alpha!0="lat" ; name dimensions </div><div>alpha!1="lon" </div><div>alpha&lat = dens&lat ; assign coordinate values and </div><div>alpha&lon =dens&lon ; units attributes </div><div> copy_VarAtts(dens,alpha) </div><div> copy_VarCoords_1(dens,alpha)</div><div>alpha&lat@units = "degrees_north" </div><div> alpha&lon@units = "degrees_east" </div><div> </div><div>;*****************************<wbr>*********************</div><div>; Now plotting the storm track</div><div>;*****************************<wbr>*********************</div><div>;*****************************<wbr>******************** </div><div> plot = new(2,graphic) </div><div> </div><div> ;*****************************<wbr>*******************</div><div>; create plot</div><div>;*****************************<wbr>*******************</div><div> wks = gsn_open_wks("X11","Composite_<wbr>high_low_storm_tracks_sig_dec_<wbr>april") ; </div><div><br></div><div> gsn_define_colormap(wks, "BlueWhiteOrangeRed") </div><div> </div><div>;*****************************<wbr>*******************</div><div>;Composite</div><div>;*****************************<wbr>*******************</div><div> res = True</div><div> res@gsnDraw = False ; don't draw</div><div> res@gsnFrame = False ; don't advance frame</div><div> res@cnFillOn = True ; turn on color</div><div> res@cnLinesOn = False ; turn off contour lines</div><div> res@cnLevelSpacingF = 0.1 ; contour interval</div><div> res@cnLevelSelectionMode = "ManualLevels";</div><div> res@cnMinLevelValF = -0.5 ; min level</div><div> res@cnMaxLevelValF =0.5</div><div> res@mpMinLatF = 30 ; choose a subregion</div><div> res@mpMaxLatF = 90.</div><div> res@mpMinLonF = 0 ; choose a subregion</div><div> res@mpMaxLonF = 180.</div><div> res@lbLabelBarOn = False ; turn off individual cb's</div><div> </div><div>;*****************************<wbr>************************</div><div>;Significance</div><div>;*****************************<wbr>************************</div><div> res1 = True</div><div> res1@gsnDraw = False ; Do not draw plot</div><div> res1@gsnFrame = False ; Do not advance frome </div><div> res1@cnFillOn = True</div><div> res1@cnMonoFillColor = True</div><div> res1@cnMonoFillPattern = False</div><div> res1@lbLabelBarOn = False ; turn off label bar</div><div> res1@cnLevelSelectionMode = "ManualLevels" ; manually specify contour levels</div><div> res1@cnMinLevelValF = 95 ; min level</div><div> res1@cnMaxLevelValF = 100 ; max level</div><div> res1@cnLevelSpacingF = 1 ; contour interval</div><div> res1@gsnSpreadColors = False</div><div> res1@cnInfoLabelOn = False ; turn off info label</div><div> res@cnLineLabelsOn = False</div><div> res1@cnFillPatterns = (/17,17,17,17,17,17,17,17,17,<wbr>17/)</div><div><br></div><div> </div><div> scale = 1.0</div><div> dotsize = 0.003 </div><div> res1@cnFillScaleF = scale</div><div> res1@cnFillDotSizeF = dotsize</div><div> ;*****************************<wbr>***********************</div><div> </div><div> </div><div> </div><div> plotA= gsn_csm_contour_map(wks,zDiff(<wbr>3,:,:), res) ;</div><div> plotB = gsn_csm_contour(wks,alpha(3,:<wbr>,:),res1) ; create plot</div><div> overlay(plotA,plotB) ; result will be plotA</div><div> plot(0) = plotA </div><div> </div><div> plotC = gsn_csm_contour_map(wks,zDiff(<wbr>11,:,:), res) ;</div><div> plotD = gsn_csm_contour(wks,alpha(11,<wbr>:,:),res1) ; create plot</div><div> </div><div> overlay (plotC, plotD) </div><div> plot(1) = plotC </div><div> </div><div><br></div><div> </div><div> </div><div> resP = True ; modify the panel plot</div><div> resP@gsnFrame = False ; don't advance panel plot</div><div> </div><div> </div><div>resP@gsnPanelFigureStrings= (/ "April" ,"December"/) ;</div><div>resP@<wbr>gsnPanelFigureStringsFontHeigh<wbr>tF=0.01</div><div>resP@amJust = "TopLeft"</div><div>resP@<wbr>gsnPanelFigureStringsBackgroun<wbr>dFillColor = (/-1/)</div><div>resP@txString = "Storm tracks difference between high and low MJ years" </div><div>resP@gsnPanelLabelBar = True ; add common colorbar</div><div>resP@lbLabelFontHeightF = 0.01 ; make labels smaller</div><div>resP@lbOrientation = "vertical" ; vertical label bar</div><div> ;</div><div> gsn_panel(wks,plot,(/2,1/),<wbr>resP) </div><div><br></div><div><br></div><div> frame(wks) </div><div>;=============================<wbr>==============</div><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_6610889102505543268gmail_signature"><div dir="ltr"><div><div><div><div>Ipshita Majhi<br></div>PhD Candidate<br></div>University of Alaska , Fairbanks<br></div>Atmospheric Science Department<br></div>(907)978-4220 <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</font></span></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></div>