[ncl-talk] empty plots in panel-plot

Noelia otero noeli1680 at gmail.com
Tue Apr 11 01:43:29 MDT 2017


Hi Mary,

Thanks a lot. I was trying to set min and max values, but I wasn't taking a
good value for levelSpacing ..,I follow your suggestion, although I had to
increase considerably the number of levels. Now it seems to work.

Many thanks again,

Best

Noelia

2017-04-11 7:47 GMT+02:00 Mary Haley <haley at ucar.edu>:

> 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/20170411/5238169b/attachment.html 


More information about the ncl-talk mailing list