[ncl-talk] White streaks on plot
Michael Notaro
mnotaro at wisc.edu
Thu Sep 20 09:17:54 MDT 2018
Never mind. Apparently the plot is fine. It is an issue
of the way that ghostscript is displaying it on NASA Discover.
Michael
Michael Notaro
Associate Director
Nelson Institute Center for Climatic Research
University of Wisconsin-Madison
Phone: (608) 261-1503
Email: mnotaro at wisc.edu
________________________________
From: Michael Notaro
Sent: Thursday, September 20, 2018 9:55 AM
To: ncl-talk at ucar.edu
Subject: White streaks on plot
Please see the attached plot. Does anyone know why I am getting white
streaks across my plot and how to get rid of them? I am using NCL 6.4.0.
Thanks, Michael
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/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
periods=(/"nov2014","dec2014","jan2015","feb2015","mar2015"/)
simulated=new((/420,520,5/),float)
do ip=0,4
a=addfile("../nudging/bias_corr_save_albedo_nudging_nldas2_"+periods(ip)+".nc","r")
XLAT=a->XLAT
XLONG=a->XLONG
simulated(:,:,ip)=a->climo_model
end do
model=dim_avg(simulated)
model at _FillValue=1e+35
nlat=420
nlon=520
iy=ispan(0,419,1)
jx=ispan(0,519,1)
iy!0="iy"
iy&iy=iy
jx!0="jx"
jx&jx=jx
lat2d=XLAT
lon2d=XLONG
lat2d!0="iy"
lat2d!1="jx"
lat2d at iy=iy
lat2d at jx=jx
lat2d at units="degrees_north"
lat2d at rcm_mapping="rcm_map"
lat2d at coordinates="iy jx"
lon2d!0="iy"
lon2d!1="jx"
lon2d at iy=iy
lon2d at jx=jx
lon2d at units="degrees_east"
lon2d at rcm_mapping="rcm_map"
lon2d at coordinates="iy jx"
copy_VarCoords(lat2d,model)
model at lat2d=lat2d
model at lon2d=lon2d
wks = gsn_open_wks("ps","less")
gsn_define_colormap(wks,"BkBlAqGrYeOrReViWh200")
res=True
res at gsnAddCyclic = False
res at cnFillOn = True
res at cnRasterModeOn = True
res at cnLinesOn = False
res at mpProjection = "LambertConformal"
res at mpLambertParallel1F = 37.
res at mpLambertParallel2F = 49.
res at mpLambertMeridianF = -98.
res at cnLevelSelectionMode = "ExplicitLevels"
res at cnLevels = (/.1,.15,.2,.25,.3,.35,.4,.45,.5,.55,.6/)
res at cnFillColors=(/130,120,110,100,90,80,70,60,50,40,30,20/)
res at tiMainFontHeightF = 0.02
res at tmYROn = False
res at tmXTOn = False
res at mpOceanFillColor = -1
res at mpLandFillColor = "white"
res at pmLabelBarDisplayMode = True
res at lbOrientation = "vertical"
res at mpOutlineBoundarySets = "AllBoundaries"
res at mpGeophysicalLineThicknessF = 2.5
res at tiMainString = "Albedo"
res at mpLimitMode = "Corners"
res at mpLeftCornerLatF = lat2d(0,0)-.2
res at mpLeftCornerLonF = lon2d(0,0)-.2
res at mpRightCornerLatF = lat2d(nlat-2,nlon-2)+.2
res at mpRightCornerLonF = lon2d(nlat-2,nlon-2)+.2
plot = gsn_csm_contour_map(wks,model,res)
end
Michael Notaro
Associate Director
Nelson Institute Center for Climatic Research
University of Wisconsin-Madison
Phone: (608) 261-1503
Email: mnotaro at wisc.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180920/d70bceab/attachment.html>
More information about the ncl-talk
mailing list