<div dir="ltr"><div class="gmail_default" style="font-size:small">Ipshita,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks for providing the files, this was helpful!</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The key to plotting this seems to be using raster mode:</div><div class="gmail_default" style="font-size:small">
<p class=""><font face="monospace, monospace"><span class=""> res@</span><span class="">cnFillMode</span><span class=""> = </span>"RasterFill"</font></p><p class="">You also want to make sure you have both of these set:</p>
<p class=""><font face="monospace, monospace"><span class=""> res@</span><span class="">trGridType</span><span class=""> = </span><span class="">"TriangularMesh"</span><span class=""> </span>; allow missing coordinates <br> res@<span class="">gsnAddCyclic</span> = False ; don't wrap longitude point</font></p>
<p class="">A minor thing: you have this line, which doesn't have any purpose:<br></p><p class="">
</p><p class=""><font face="monospace, monospace">swe_FillValue=-9.99e+08</font></p><p class="">I think you meant to type:</p><p class="">
</p><p class=""><font face="monospace, monospace">swe@_FillValue=-9.99e+08</font></p><p class="">but this is not needed because "swe" already has the correct _FillValue attribute attached.</p><p class="">Finally, this is not a problem, but there's a better way to set your colors for the contours. Instead of calling "gsn_define_colormap" and setting "res@gsnSpreadColors" to True, you can simply set "res@cnFillPalette" to "amwg". The color map will be spanned automatically.</p><p class="">Please see the attached script.</p><p class="">--Mary</p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 2, 2016 at 4:40 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">Hi Mary,<div><br></div><div>I transfered the files . They are named <a href="http://GlobSnow_SWE_L3A_20110101_v2.0.nc" target="_blank">GlobSnow_SWE_L3A_20110101_v2.0.nc</a></div><div>and <a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a>.</div><div><br></div><div>Thank you for your help.</div><div>Best Regards</div><span class="HOEnZb"><font color="#888888"><div>Ipshita</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 1:49 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@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 class="gmail_default" style="font-size:small">The only thing I can think of at this point is perhaps there'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><font color="#888888"><br></font></span></div><span><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><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 2, 2016 at 3:30 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">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 = "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 "$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>a=addfile("/Users/ipshita/Desktop/GlobsnowV2.0/version02/<a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a>","r")</div></span><div>b=addfile("/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>","r")</div><div>swe=a->swe</div><div><br></div><div>time=a->time</div><div><br></div><div>lat=b->lat</div><div>lon=b->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("png","ease_map") </div><span><div> gsn_define_colormap(wks,"amwg") ; 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= "ExplicitLevels" ; 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 = "Polar Plot"</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 = "TriangularMesh" ; allow missing coordinates</div><div> res@gsnAddCyclic = False</div><span><div> res@gsnPolar = "NH" ; specify the hemisphere</div><div> res@mpMinLatF = 35</div><div> </div><div> res@tiMainString = "gsn_csm_contour_map_polar"</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"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></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 = "RasterFill"</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"><<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><div dir="ltr">Dear NCL,<div><br></div><div>I have been trying to plot this map correctly for a long time. I don'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 "$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>a=addfile("/Users/ipshita/Desktop/GlobsnowV2.0/version02/<a href="http://monthly_82_13.nc" target="_blank">monthly_82_13.nc</a>","r")</div><div><br></div><div>swe=a->swe</div><div><br></div><div>time=a->time</div><div>lat=a->lat</div><div>lon=a->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("pdf","ease_map") </div><div> gsn_define_colormap(wks,"amwg") ; 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= "ExplicitLevels" ; 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 = "gsn_csm_contour"</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->lat</div><div> swe@lon2d = a->lon</div><div><br></div><div> res@trGridType = "TriangularMesh" ; allow missing coordinates</div><div><br></div><div> res@gsnPolar = "NH" ; specify the hemisphere</div><div> res@mpMinLatF = 35</div><div> res@tiMainString = "gsn_csm_contour_map_polar"</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 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>