<div dir="ltr">Dear NCL,<div><br></div><div>I figured out the lat and lon problem. I was not attaching the coordinates correctly. I have another problem. When I plot with significance the whole plot is hatched. I am attaching that plot with this email and the code as well. I would be grateful if you could guide me on this.</div><div><div>;*******************************************</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>;*******************************************</div><div><br></div><div>;============================================</div><div>; set up parameters</div><div>;============================================</div><div><br></div><div>;********************************************</div><div>;This is reading in Snow water equivalent which goes from 0 to 250 mm</div><div>; from 1982 to 2012</div><div>;********************************************</div><div><br></div><div>dira = "/Users/ipshita/Desktop/GlobsnowV2.0/version02/"</div><div>dirb = dira + "2011/"</div><div>;dira = "./"</div><div>;dirb = "./"</div><div>a = addfile(dira+"<a href="http://monthly_82_13.nc">monthly_82_13.nc</a>","r")</div><div>b = addfile(dirb+"<a href="http://GlobSnow_SWE_L3A_20110101_v2.0.nc">GlobSnow_SWE_L3A_20110101_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>;********************************************</div><div>;This is to extract 1982-2012 of swe data</div><div>;********************************************</div><div>swe_82_12=swe(0:371,:,:)</div><div><br></div><div>;********************************************</div><div><br></div><div>;********************************************</div><div>;This is to read in All India Rainfall</div><div>;For March , March , June</div><div>;********************************************</div><div><br></div><div>MJ_1982_2014=asciiread("~/Documents/NCL_files/SST/Monthly/MJ.txt",(/34,1/),"float")</div><div>MJ_1982_2012=MJ_1982_2014(0:30,0)</div><div><br></div><div>print(dimsizes(MJ_1982_2012))</div><div><br></div><div>;*******************************************</div><div>;Extracting monthly data for sce</div><div>;*******************************************</div><div><br></div><div>March_swe = new((/31,721,721/),float)</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>;=================================================</div><div>March_swe@_FillValue = -9.99e+08</div><div>;**********************************************</div><div><br></div><div>March_dt = dtrend_msg(March_swe&time,March_swe(lat|:,lon|:,time|:),True,False)</div><div><br></div><div>printMinMax(March_dt,False)</div><div><br></div><div>;****************************************</div><div>;Calculating Correlation</div><div>;***************************************</div><div>x1d=ndtooned(MJ_1982_2012)</div><div><br></div><div>corr_March=escorc(March_dt,x1d)</div><div>print(dimsizes(corr_March))</div><div><br></div><div>corr_March!0="lat" </div><div>corr_March!1="lon"</div><div> </div><div>corr_March&lat = swe&lat ; assign coordinate values and </div><div>corr_March&lon=swe&lon ; units attributes </div><div><br></div><div>copy_VarAtts(swe,corr_March) </div><div>copy_VarCoords(swe,corr_March)</div><div> </div><div>corr_March&lat@units = "degrees_north" </div><div>corr_March&lon@units = "degrees_east"</div><div>;**********************************************</div><div>;**********************************************</div><div>;Calculating significance</div><div>;**********************************************</div><div><br></div><div>prob_March=rtest(corr_March,31,0)</div><div>sig_March =100*(1-prob_March)</div><div>sig_March = mask(sig_March, sig_March.eq.100,False) </div><div>sig_March = mask(sig_March, sig_March.le.90,False)</div><div><br></div><div>sig_March !0="lat" </div><div>sig_March !1="lon"</div><div><br></div><div>sig_March&lat = swe&lat ; assign coordinate values and </div><div>sig_March&lon=swe&lon ; units attributes </div><div><br></div><div>copy_VarAtts(swe,sig_March) </div><div>copy_VarCoords(swe,sig_March)</div><div> </div><div> sig_March&lat@units = "degrees_north" </div><div>sig_March&lon@units = "degrees_east"</div><div><br></div><div>;**********************************************</div><div>;************************************************</div><div>; create plot</div><div>;************************************************</div><div> wks = gsn_open_wks("pdf","SWE_March_MJ_correlation_significance") ; </div><div> </div><div> res = True ; Plot modes desired.</div><div> res@gsnDraw = False ; Do not draw plot</div><div> res@gsnFrame = False ; Do not advance frome</div><div> res@gsnMaximize = True ; Maximize plot</div><div> res@cnFillOn = True ; color plot desired</div><div> res@cnFillPalette = "BrownBlue12"</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> </div><div> </div><div>;*******************************************</div><div>; georeferencing: plot on polar projection</div><div>;*******************************************</div><div> res@trGridType = "TriangularMesh" ; allow missing coordinates</div><div> res@gsnAddCyclic = False</div><div> res@cnFillMode = "RasterFill"</div><div><br></div><div> res@gsnPolar = "NH" ; specify the hemisphere</div><div> res@mpMinLatF = 35</div><div> </div><div> res@tiMMarchainString = "March SWE and MJ correlation"</div><div><br></div><div> plot = gsn_csm_contour_map_polar(wks,corr_March,res)</div><div> </div><div>;**********************************************</div><div>;Adding significance plot</div><div>;********************************************</div><div><br></div><div>;===============================================</div><div>;Significance</div><div>;================================================ </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@cnLinesOn = False ; no contour lines</div><div> res1@cnMonoFillColor = True</div><div> res1@cnMonoFillPattern = False</div><div> res1@lbLabelBarOn = False </div><div> res1@cnLevelSelectionMode = "ManualLevels" ; manually specify contour levels</div><div>res1@cnMinLevelValF = 90 ; min level</div><div>res1@cnMaxLevelValF = 100 ; max level</div><div> res1@cnLevelSpacingF = 1 ; contour interval</div><div> res1@gsnSpreadColors = False</div><div> res1@cnFillPatterns = (/17,17,17,17,17,17,17,17,17,17/)</div><div>res1@cnLineLabelsOn = False</div><div> </div><div>; scale = 1.0</div><div> ; dotsize = 0.005 </div><div> ; res@cnFillScaleF = scale</div><div> ; res@cnFillDotSizeF = dotsize</div><div><br></div><div><br></div><div> plot1 = gsn_csm_contour(wks,sig_March,res1)</div><div> </div><div> overlay (plot, plot1)</div><div><br></div><div> draw (plot)</div><div> frame(wks) </div><div> </div><div> </div></div><div>;************************************</div><div><br></div><div>I get the following warnings</div><div><br></div><div><div>warning:tiMMarchainString is not a valid resource in SWE_March_MJ_correlation_significance_contour at this time</div><div>warning:tiMMarchainString is not a valid resource in ezplot1p_map at this time</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>get_lon_values: Warning: The range your of longitude values is greater than 360.</div><div>warning:tmXBMode is not a valid resource in SWE_March_MJ_correlation_significance_contour.PlotManager at this time</div><div>warning:tmXBValues is not a valid resource in SWE_March_MJ_correlation_significance_contour.PlotManager at this time</div><div>warning:tmXBLabels is not a valid resource in SWE_March_MJ_correlation_significance_contour.PlotManager at this time</div><div>warning:tmXBMinorValues is not a valid resource in SWE_March_MJ_correlation_significance_contour.PlotManager at this time</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>get_lat_values: Warning: Your latitude values do not fall between -90 and 90 inclusive.</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>You will not get 'nice' latitude labels.</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 21, 2016 at 8:56 AM, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.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"><div><div><div><div>Hi,<br><br></div>What do you mean by ".... lat lon dimensions don't show up"? What do you get instead?<br><br></div>I think the reason for the warnings is that you are not plotting on top of a map, and thus those resources don't apply. Perhaps you want to plot using gsn_csm_contour_map() ?<br><br></div>HTH...<br></div>Rick<br><br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Sun, Nov 20, 2016 at 5:17 PM, Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Dear NCL,<div><br></div><div>I am trying to plot SWE data and the lat lon dimensions don't show up.</div><div>I would be grateful if you could guide me on this</div><div>******************************<wbr>***************************</div><div><div>;*****************************<wbr>***************</div><div>;This is reading in Snow water equivalent which goes from 0 to 250 mm</div><div>; from 1982 to 2012</div><div>;*****************************<wbr>***************</div><div><br></div><div>dira = "/Users/ipshita/Desktop/Globsn<wbr>owV2.0/version02/"</div><div>dirb = dira + "2011/"</div><div>;dira = "./"</div><div>;dirb = "./"</div><div>a = addfile(dira+"<a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a><wbr>","r")</div><div>b = addfile(dirb+"<a href="http://GlobSnow_SWE_L3A_20110101_v2.0.nc" target="_blank">GlobSnow_SWE_L3A<wbr>_20110101_v2.0.nc</a>","r")</div><div><br></div><div>swe = a->swe</div><div>time = a->time</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><br></div><div>;*****************************<wbr>***************</div><div>;This is to read in All India Rainfall</div><div>;For March , March , June</div><div>;*****************************<wbr>***************</div><div><br></div><div>MJ_1982_2014=asciiread("~/Docu<wbr>ments/NCL_files/SST/Monthly/<wbr>MJ.txt",(/34,1/),"float")</div><div>MJ_1982_2012=MJ_1982_2014(0:30<wbr>,0)</div><div><br></div><div>print(dimsizes(MJ_1982_2012))</div><div><br></div><div>;*****************************<wbr>**************</div><div>;Extracting monthly data for sce</div><div>;*****************************<wbr>**************</div><div><br></div><div>March_swe = new((/31,721,721/),float)</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,Marc<wbr>h_swe(lat|:,lon|:,time|:),True<wbr>,False)</div><div><br></div><div>printMinMax(March_dt,False)</div><div><br></div><div>;*****************************<wbr>***********</div><div>;Calculating Correlation</div><div>;*****************************<wbr>**********</div><div>x1d=ndtooned(MJ_1982_2012)</div><div><br></div><div>corr_March=escorc(March_dt,x1d<wbr>)</div><div>corr_March!0="lat" </div><div>corr_March!1="lon" </div><div>corr_March&lat=lat</div><div>corr_March&lon=lon</div><div><br></div><div>copy_VarCoords(swe,corr_March)<br></div><div><br></div><div><br></div><div>;*****************************<wbr>*****************</div><div><br></div><div>;*****************************<wbr>*******************</div><div>; create plot</div><div>;*****************************<wbr>*******************</div><div> wks = gsn_open_wks("x11","SWE_March_<wbr>MJ_correlation_significance") ; </div><div> </div><div> res = True ; Plot modes desired.</div><div> res@gsnMaximize = True ; Maximize plot</div><div> res@cnFillOn = True ; color plot desired</div><div> res@cnFillPalette = "BrownBlue12"</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> </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><br></div><div> res@gsnPolar = "NH" ; specify the hemisphere</div><div> res@mpMinLatF = 35</div><div> </div><div> res@tiMainString = "March SWE and MJ correlation"</div><div><br></div><div> plot = gsn_csm_contour(wks,corr_March<wbr>,res)</div><div> </div><div>;*****************************<wbr>**************************</div><div><br></div><div><div>warning:gsnPolar is not a valid resource in SWE_March_MJ_correlation_signi<wbr>ficance_contour at this time</div><div>warning:mpMinLatF is not a valid resource in SWE_March_MJ_correlation_signi<wbr>ficance_contour at this time</div><div><br></div></div><div><br></div><div class="m_-2299216546635779236m_-8777661468870570439gmail_signature"><div dir="ltr"><br></div></div>
</div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_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>
</div>