[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