<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">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><div><br></div><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></div>