[ncl-talk] problems when using the "res at cnFillPalette = "BlGrYeOrReVi200""
Rick Brownrigg
brownrig at ucar.edu
Mon May 10 07:45:21 MDT 2021
Hi,
I looked at the code where the error message arises, and this is indeed a
bit of a mystery. If you might send me one of your WRF files, I'll try to
run the script in the debugger to get a better idea of whats going on.
Rick
On Mon, May 10, 2021 at 12:10 AM ZUO Ziping via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> Dear ncl-team,
>
> The other day when I was trying to use the colors in the "
> *BlGrYeOrReVi200*", I use the *res at cnFillPalette = "BlGrYeOrReVi200"*.
> However, it kept reminding me "*warning:ContourPlotSetValues: color index
> (204) exceeds size of palette, defaulting to foreground color for entry
> (12)*".
> But when I changed the instruction to
> *gsn_define_colormap(wks,"BlGrYeOrReVi200")*, it works fine. I'd like to
> ask the reason for these two instructions. The following is my full script.
>
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
>
> DATADir = "/disk/wrfout/"
> FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d04_2020-06* ")
> numFILES = dimsizes(FILES)
> print("numFILES = " + numFILES)
> print(FILES)
> print (" ")
>
> type = "pdf"
> wks = gsn_open_wks(type,"average_temp")
>
> res = True
> res at MainTitle = "REAL-TIME WRF"
> res at Footer = False
>
> pltres = True
> pltres at PanelPlot = True
> pltres at FramePlot = False
>
> mpres = True
> mpres at mpDataBaseVersion = "HighRes"
> mpres at mpGeophysicalLineColor = "Black"
> mpres at mpGridLineColor = "Black"
> mpres at mpDataResolution = "Finest"
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> a = addfiles(FILES+".nc","r")
>
> times = wrf_user_getvar(a,"times",-1) ; get all times in the file
> ntimes = dimsizes(times) ; number of times in the file
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
>
> do k = 0,23,1
> meantemp = new((/162,213/),float)
> meantemp = 0
> temp = new((/ntimes,162,213/),float)
>
> do it = k,ntimes-1,24 ; TIME LOOP
> print("Working on time: " + times(it) )
> fi = addfile(FILES(it),"r")
> T_2 = fi ->T2 ; T2 in Kelvin
> tc2 = T_2-273.16 ; T2 in C
> temp(it,:,:) = tc2
> meantemp(:,:) = meantemp(:,:)+temp(it,:,:)
> end do
>
> meantemp(:,:) = meantemp(:,:)/30
>
> opts = res
> opts at cnFillOn = True
> opts at ContourParameters = (/23,37,1/)
> opts at cnFillPalette = "BlGrYeOrReVi200"
> contour_FLX = wrf_contour(a[k],wks,meantemp,opts)
> delete(opts)
>
> plot = wrf_map_overlays(a[k],wks,contour_FLX,pltres,mpres)
>
> end do ; END OF TIME LOOP
>
> end
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
>
> Thank you for your time!
>
> Best regards
> Zoe
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210510/b79ea53d/attachment.html>
More information about the ncl-talk
mailing list