[ncl-talk] empty plots in panel-plot
Mary Haley
haley at ucar.edu
Mon Apr 10 23:47:43 MDT 2017
Dear Noelia,
When you are paneling a series of color contour plots and then paneling
them later, you first want to be sure you use a fixed set of contour levels
across all the plots, so you can then use a single labelbar in the final
panel plot.
The function nice_mnmxintvl
is what you can use to fix your contour levels
mnmxlvl = nice_mnmxintvl(min(data),max(data),16,False)
res at cnLevelSelectionMode = "ManualLevels"
res at cnMinLevelValF = mnmxlvl(0)
res at cnMaxLevelValF = mnmxlvl(1)
res at cnLevelSpacingF = mnmxlvl(2)
You will need to get the min and max of your data across all plots,
and pass that to
The "min(data)" and "max(data)" are from another example. You will need to
plug in the min/max of your own data across all plots.
You now want to turn off the individual labelbar you're getting for each
plot:
res at lbLabelBarOn = False
and then turn on the panel labelbar when you call gsn_panel:
pres = True
pres at gsnPanelLabelBar = True
...
gsn_panel(...)
You may want to see our "panel" examples page:
http://www.ncl.ucar.edu/Applications/panel.shtml
I think panel_35.ncl uses the nice_mnmxintvl function.
--Mary
On Mon, Apr 10, 2017 at 11:37 AM, Noelia otero <noeli1680 at gmail.com> wrote:
> Hi!
>
> I am plotting in a panel plots seasonal means from different models, and I
> got some "empty" plots for some specific models and seasons. My guess is
> that it is due to the large spread of values (or some outliers). Maybe,
> it's a simple question..but I don't know how to set the resources to get a
> good label bar. I attach the final plot that I got with my code (see a part
> of the code below):
>
> ; i= model ,s =season
>
> model = addfile(fils(i), "r")
> c=str_split(fils(i), "/")
> nmod=str_split(c(6), "#")
> mmod=nmod(0)
> print(nmod)
> sfile=systemfunc("ls " + dir2 + mmod +"*"+ season(s)+".nc")
> smodel=addfile(sfile, "r")
>
> svar = smodel ->$var$
> svar!1="lat"
> svar!2="lon"
>
> ;------Getting statistics for each grid point-------
>
> stats_seas_var = new((/nlat,mlon/), typeof(svar),-9e+33)
>
> ;-------get 95th possition 25 of statistics---------
>
>
> opt=False
> do nl=0,nlat-1
> do ml=0,mlon-1
>
> seas_st_var = stat_dispersion( svar(:,nl,ml), opt)
> stats_seas_var(nl,ml) = seas_st_var(0)
>
> ;stat_dispersion: (0) mean, (1) std, (2), min, (14),max
> end do
> end do
>
>
> stats_seas_var at long_name = "mean"
> stats_seas_var!0="lat"
> stats_seas_var!1="lon"
> stats_seas_var&lat = lat
> stats_seas_var&lon = lon
>
> res = True
> res at gsnMaximize = True
> res at cnFillOn = True
> res at cnLinesOn = False
> res at gsnAddCyclic = False
> res at gsnDraw = False
> res at gsnFrame = False
> res at mpMinLonF = min(lon)
> res at mpMinLatF = min(lat)
> res at mpMaxLonF = max(lon)
> res at mpMaxLatF = max(lat)
> res at lbLabelAutoStride=True
> res at cnLevelSelectionMode = "ManualLevels"
> res at cnMinLevelValF = 0.
> res at cnMaxLevelValF = 10.
> res at cnLevelSpacingF = 5.
> res at lbLabelBarOn = False
> res at gsnCenterString =nmod(0)
> res at gsnRightString =season(s)
> plotseas = gsn_csm_contour_map(wks,stats_seas_var,res)
>
> I really appreciate any suggestion to improve it,
>
> Many thanks in advance,
>
> Noelia
>
>
> _______________________________________________
> 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/20170410/0224440a/attachment.html
More information about the ncl-talk
mailing list