<html><head></head><body dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div class="ApplePlainTextBody">Hi, I am plotting snow water equivalent using NCL. The following is my script:<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br><br>begin<br><br>f=addfile("/Users/xueyanzhang/NOHRSC_SNODAS-master/July_mean/201807.nc", "r")<br>swe=f->SWE(0,:,:)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;plot section;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><br> wks = gsn_open_wks("png", "201807_swe_mean")<br> gsn_define_colormap(wks, "BlAqGrYeOrReVi200")<br><br> setvalues NhlGetWorkspaceObjectId()<br> "wsMaximumSize" : 300000000<br> end setvalues<br><br> res = True<br><br> res@gsnAddCyclic = False<br> res@gsnDraw = False <br> res@gsnFrame = False <br> res@gsnMaximize = True<br> res@gsnRightString = "Snow water equivalent (mm)"<br> res@gsnLeftString = "201807 mean"<br> res@gsnStringFontHeightF = 0.01<br><br><br> res@cnFillOn = True<br> res@cnLinesOn = False<br> res@cnLineLabelsOn = False <br> res@cnFillDrawOrder = "PreDraw"<br> res@cnFillMode = "RasterFill" <br> res@cnLineLabelsOn = False<br> res@cnLevelSelectionMode = "ManualLevels"<br> res@cnMinLevelValF = 0<br> res@cnMaxLevelValF = 1000<br> res@cnLevelSpacingF = 50<br><br> res@mpFillOn = False ; turn off map fill<br> res@mpDataBaseVersion = "MediumRes"<br> res@mpOutlineBoundarySets = "USStates"<br> res@mpProjection = "CylindricalEquidistant"<br> res@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last<br> res@mpGeophysicalLineColor = "Black"<br> res@mpNationalLineColor = "Black"<br> res@mpLimbLineColor = "Black"<br> res@mpPerimLineColor = "Black"<br> res@mpUSStateLineColor = "black"<br> res@mpGeophysicalLineThicknessF = 2.<br> res@mpLandFillColor = "white"<br> res@mpInlandWaterFillColor = "white"<br> res@mpOceanFillColor = "white"<br><br> res@mpLimitMode = "LatLon"<br> res@mpMaxLatF = 53<br> res@mpMinLatF = 23<br> res@mpMaxLonF = -65<br> res@mpMinLonF = -125<br> ;res@mpCenterLatF = (53+23)/2<br> ;res@mpCenterLonF = (-65-125)/2<br><br> plot = gsn_csm_contour_map(wks, swe, res)<br><br> draw(plot)<br> frame(wks)<br><br><br>end<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>But the plot is not consistent with the map NCL provides.<br><object type="application/x-apple-msg-attachment" data="cid:402020C7-3B03-4A8B-91E0-99721066CE67@usu.edu" apple-inline="no" id="8DF487CE-5EE3-4064-84B4-BAC3F509A3DF" height="102" width="149" apple-width="yes" apple-height="yes"></object><br>I searched and found it may be because of WGS84, this datum. But I don’t know how to solve it. Can you help me?<br><br>Xueyan</div></div></body></html>