[ncl-talk] White streaks on plot

[Michael Notaro]

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/37cd8a3c/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: less.ps
Type: application/postscript
Size: 310126 bytes
Desc: less.ps
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180920/37cd8a3c/attachment-0001.ps>