<div dir="ltr">Hi Mary,<div><br></div><div>I transfered the files . They are named <a href="http://GlobSnow_SWE_L3A_20110101_v2.0.nc">GlobSnow_SWE_L3A_20110101_v2.0.nc</a></div><div>and <a href="http://monthly_82_13.nc">monthly_82_13.nc</a>.</div><div><br></div><div>Thank you for your help.</div><div>Best Regards</div><div>Ipshita</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 2, 2016 at 1:49 PM, Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;</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 class="gmail_default" style="font-size:small">The only thing I can think of at this point is perhaps there&#39;s an issue with the longitude values crossing from 0 to 360 abruptly.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Can you provide your data file on our ftp site?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP" target="_blank">http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP</a><span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888"><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 2, 2016 at 3:30 PM, Ipshita Majhi <span dir="ltr">&lt;<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Mary,<div><br></div><div>I noticed that my plot has straight lines that run through it. I removed part of those lines when I commented out res@trGridType = &quot;TriangularMesh.</div><div>But some of the lines still remain.</div><div>How do I remove those lines?</div><div>I am attaching my code and the picture of the map with this email.</div><div>**************************************************************************</div><div><span><div>;*******************************************</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;</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>a=addfile(&quot;/Users/ipshita/Desktop/GlobsnowV2.0/version02/<a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a>&quot;,&quot;r&quot;)</div></span><div>b=addfile(&quot;/Users/ipshita/Desktop/GlobsnowV2.0/version02/2011/<a href="http://GlobSnow_SWE_L3A_20110101_v2.0.nc" target="_blank">GlobSnow_SWE_L3A_20110101_v2.0.nc</a>&quot;,&quot;r&quot;)</div><div>swe=a-&gt;swe</div><div><br></div><div>time=a-&gt;time</div><div><br></div><div>lat=b-&gt;lat</div><div>lon=b-&gt;lon</div><div><br></div><div>swe@lat2d=lat</div><div>swe@lon2d=lon</div><div><br></div><div>printMinMax(lat,False)</div><div>printMinMax(lon,False)</div><span><div>swe_FillValue=-9.99e+08</div><div><br></div><div>printVarSummary(swe)</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><br></div><div>;*******************************************</div><div>; Create plots</div><div>;*******************************************</div><div>    </div></span><div>    wks = gsn_open_wks(&quot;png&quot;,&quot;ease_map&quot;)  </div><span><div>    gsn_define_colormap(wks,&quot;amwg&quot;)             ; Change color map</div><div><br></div><div>    res                     = True              ; Plot modes desired.</div><div>    res@gsnMaximize         = True              ; Maximize plot</div><div><br></div><div>;---This resource not needed in V6.1.0</div><div><br></div><div>    res@gsnSpreadColors     = True              ; use full range of colormap</div><div><br></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><br></div><div>;    res@cnLevelSelectionMode= &quot;ExplicitLevels&quot; ; set explicit contour levels</div><div>;    res@cnLevels            = (/0,1,5,10,25,100,200,300,400/)</div><div>    </div><div>;************************************************</div><div>;No georeferencing: simple contour. Draws faster</div><div>;************************************************</div></span><div>    res@tiMainString = &quot;Polar Plot&quot;</div><div>   ; plot = gsn_csm_contour(wks,swe(0,:,:),res)    ; contour, no map</div><span><div><br></div><div>;*******************************************</div><div>; georeferencing: plot on polar projection</div><div>;*******************************************</div></span><div>  ;  res@trGridType = &quot;TriangularMesh&quot;        ; allow missing coordinates</div><div>  res@gsnAddCyclic = False</div><span><div>    res@gsnPolar   = &quot;NH&quot;                          ; specify the hemisphere</div><div>    res@mpMinLatF  = 35</div><div>    </div><div>    res@tiMainString = &quot;gsn_csm_contour_map_polar&quot;</div><div><br></div><div>    plot = gsn_csm_contour_map_polar(wks,swe(0,:,:),res)   </div></span></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><span>On Wed, Jun 1, 2016 at 3:49 PM, Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;</span> wrote:<br></span><div><div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Ipshita,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">This is indeed strange.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">As a quick test, can you try (for the gsn_csm_contour_map plot):</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">res@gsnAddCyclic = False</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Also, are you getting any kind of error messages?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can also try using raster contours:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">  res@cnFillMode = &quot;RasterFill&quot;</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you continue to have problems, we may need to take a look at your data, and in particular the lat/lon arrays on the file. It looks like you are doing everything correctly otherwise.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Jun 1, 2016 at 5:33 PM, Ipshita Majhi <span dir="ltr">&lt;<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>&gt;</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><div dir="ltr">Dear NCL,<div><br></div><div>I have been trying to plot this map correctly for a long time. I don&#39;t know where I am going wrong. I will be grateful if you could guide me on this. I am also attaching the pdf version of the map I am getting.</div><div><br></div><div>Best Regards</div><div>Ipshita</div><div><div>;***************************************************</div><div>;*******************************************</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;</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>a=addfile(&quot;/Users/ipshita/Desktop/GlobsnowV2.0/version02/<a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a>&quot;,&quot;r&quot;)</div><div><br></div><div>swe=a-&gt;swe</div><div><br></div><div>time=a-&gt;time</div><div>lat=a-&gt;lat</div><div>lon=a-&gt;lon</div><div><br></div><div>swe_FillValue=-9.99e+08</div><div><br></div><div>printVarSummary(swe)</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><br></div><div>;*******************************************</div><div>; Create plots</div><div>;*******************************************</div><div>    </div><div>    wks = gsn_open_wks(&quot;pdf&quot;,&quot;ease_map&quot;)  </div><div>    gsn_define_colormap(wks,&quot;amwg&quot;)             ; Change color map</div><div><br></div><div>    res                     = True              ; Plot modes desired.</div><div>    res@gsnMaximize         = True              ; Maximize plot</div><div><br></div><div>;---This resource not needed in V6.1.0</div><div>    res@gsnSpreadColors     = True              ; use full range of colormap</div><div><br></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><br></div><div>    res@cnLevelSelectionMode= &quot;ExplicitLevels&quot; ; set explicit contour levels</div><div>    res@cnLevels            = (/0,1,5,10,25,100,200,300,400/)</div><div>    </div><div>;************************************************</div><div>;No georeferencing: simple contour. Draws faster</div><div>;************************************************</div><div>    res@tiMainString = &quot;gsn_csm_contour&quot;</div><div>    plot = gsn_csm_contour(wks,swe,res)    ; contour, no map</div><div><br></div><div>;*******************************************</div><div>; georeferencing: plot on polar projection</div><div>;*******************************************</div><div>                                             ; georeference</div><div>    swe@lat2d = a-&gt;lat</div><div>    swe@lon2d = a-&gt;lon</div><div><br></div><div>    res@trGridType = &quot;TriangularMesh&quot;        ; allow missing coordinates</div><div><br></div><div>    res@gsnPolar   = &quot;NH&quot;                          ; specify the hemisphere</div><div>    res@mpMinLatF  = 35</div><div>    res@tiMainString = &quot;gsn_csm_contour_map_polar&quot;</div><div><br></div><div>    plot = gsn_csm_contour_map_polar(wks,swe(0,:,:),res)   </div><span><font color="#888888"><div><br></div><div><br></div>-- <br><div 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><a href="tel:%28907%29978-4220" value="+19079784220" target="_blank">(907)978-4220</a> <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</font></span></div></div>
<br></div></div>_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div></div></div><div><div><br><br clear="all"><div><br></div>-- <br><div 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><a href="tel:%28907%29978-4220" value="+19079784220" target="_blank">(907)978-4220</a> <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</div></div></div>
</blockquote></div><br></div>
</div></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>