[ncl-talk] Question about SNODAS plot

Rick Brownrigg brownrig at ucar.edu
Wed Jun 5 10:26:11 MDT 2019


Hi,

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.

FWIW -- I found this document that describes to original SNODAS:

    https://nsidc.org/data/g02158

Rick

On Wed, Jun 5, 2019 at 9:05 AM 张雪岩 <xueyanzhang520 at gmail.com> wrote:

> Hi, I am plotting snow water equivalent using NCL. The following is my
> script:
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
>
> f=addfile("/Users/xueyanzhang/NOHRSC_SNODAS-master/July_mean/201807.nc",
> "r")
> swe=f->SWE(0,:,:)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;plot
> section;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> wks = gsn_open_wks("png", "201807_swe_mean")
> gsn_define_colormap(wks, "BlAqGrYeOrReVi200")
>
> setvalues NhlGetWorkspaceObjectId()
>   "wsMaximumSize" : 300000000
> end setvalues
>
> res = True
>
> res at gsnAddCyclic           = False
> res at gsnDraw                = False
> res at gsnFrame               = False
> res at gsnMaximize            = True
> res at gsnRightString         = "Snow water equivalent (mm)"
> res at gsnLeftString          = "201807 mean"
> res at gsnStringFontHeightF   = 0.01
>
>
> res at cnFillOn               = True
> res at cnLinesOn              = False
> res at cnLineLabelsOn         = False
> res at cnFillDrawOrder        = "PreDraw"
> res at cnFillMode             = "RasterFill"
> res at cnLineLabelsOn         = False
> res at cnLevelSelectionMode   = "ManualLevels"
> res at cnMinLevelValF         = 0
> res at cnMaxLevelValF         = 1000
> res at cnLevelSpacingF        = 50
>
> res at mpFillOn                = False           ; turn off map fill
> res at mpDataBaseVersion       = "MediumRes"
> res at mpOutlineBoundarySets   = "USStates"
> res at mpProjection     = "CylindricalEquidistant"
> res at mpOutlineDrawOrder      = "PostDraw"  ; draw continental outline last
> res at mpGeophysicalLineColor  = "Black"
> res at mpNationalLineColor     = "Black"
> res at mpLimbLineColor         = "Black"
> res at mpPerimLineColor        = "Black"
> res at mpUSStateLineColor      = "black"
> res at mpGeophysicalLineThicknessF = 2.
> res at mpLandFillColor         = "white"
> res at mpInlandWaterFillColor  = "white"
> res at mpOceanFillColor        = "white"
>
> res at mpLimitMode              = "LatLon"
> res at mpMaxLatF                = 53
> res at mpMinLatF                = 23
> res at mpMaxLonF                = -65
> res at mpMinLonF                = -125
> ;res at mpCenterLatF             = (53+23)/2
> ;res at mpCenterLonF             = (-65-125)/2
>
> plot = gsn_csm_contour_map(wks, swe, res)
>
> draw(plot)
> frame(wks)
>
>
> end
>
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> But the plot is not consistent with the map NCL provides.
>
> 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?
>
> Xueyan_______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190605/6adf67f0/attachment.html>


More information about the ncl-talk mailing list