<div dir="ltr"><div dir="ltr">Hi,<div><br></div><div>I can reproduce the issue. It really appears as though some stretching transformation has been applied during the creation of the NetCDF file -- How/where was the file obtained? I don't know for certain that WGS84 is not the issue, but I would think it should not matter much at the scale of a region as large as CONUS.</div><div><br></div><div>FWIW -- I found this document that describes to original SNODAS:</div><div><br></div><div> <a href="https://nsidc.org/data/g02158">https://nsidc.org/data/g02158</a><br></div><div><br></div><div>Rick</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 5, 2019 at 9:05 AM 张雪岩 <<a href="mailto:xueyanzhang520@gmail.com">xueyanzhang520@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">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/<a href="http://201807.nc" rel="noreferrer" target="_blank">201807.nc</a>", "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>
<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_______________________________________________<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>
</blockquote></div>