<div dir="ltr"><span class="" style="border-collapse:collapse;font-size:13px">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" style="color:rgb(17,85,204)">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><div class="" style="outline-style:none;outline-width:initial;outline-color:initial;padding-top:10px;padding-right:0px;padding-bottom:10px;padding-left:0px;width:22px;margin-top:2px;margin-right:0px;margin-bottom:0px;margin-left:0px"></div></div></span><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>