Function to plot an ensembles of data following an axis-time
draw_ensembles_time(ncfile, values)
ncfile = '#' list of files [filen];[slice];[tvarn];[varn];[kinddata]
[filen]: name of the file
[dimslice]: '|' list of [dimn],[dimrange]
[dimn]: name of the dimension
[dimrange]: range for the slice
* [integer]: which value of the dimension
* -1: all along the dimension
* -9: last value of the dimension
* [beg]@[end]@[inc] slice from [beg] to [end] every [inc]
* NOTE, no dim name all the dimension size
[tvarn]: name of the variable time ('WRFtime', for times from WRF)
[varn]: name of the variable to use for the values
[kinddata]: which kind of use should be done:
'data': to be plot as data
'ensemble': to be used as a member of the ensemble
values = [dimxyfmt];[valaxis];[axislab];[axisunits];[rangevals];[datavals];[timevals];[valskind];
[labels];[figenskind];[graphtit];[legvalues];[kindfig];[fign];[outfile];[close]
[dimxyfmt]=[dxs],[dxf],[Ndx],[ordx],[dys],[dyf],[Ndy],[ordy]: format of the values at each
axis (or single 'auto' for 'auto','auto')
[dxs]: style of x-axis ('auto' for 'pretty')
'Nfix', values computed at even 'Ndx'
'Stringlist', a given list of strings
'Vfix', values computed at even 'Ndx' increments
'pretty', values computed following aprox. 'Ndx' at 'pretty' intervals (2.,2.5,4,5,10)
[dxf]: format of the labels at the x-axis
[Ndx]: Number of ticks at the x-axis (for 'Stringlist', '@' list of labels)
[ordx]: angle of orientation of ticks at the x-axis ('auto' for horizontal)
[dys]: style of y-axis ('auto' for 'pretty')
[dyf]: format of the labels at the y-axis
[Ndy]: Number of ticks at the y-axis
[ordy]: angle of orientation of ticks at the y-axis ('auto' for horizontal)
[valaxis]: which axes should be used for the data
[axislab]: Label for the values-axis ('!' for spaces)
[axisunits]: units of the values-label
[rangevals]: Range of the axis with the values (single 'None' for 'auto','auto')
[vmin],[vmax]: minimum and maximum values
'auto': the computed minimumm or maximum of the values
[datavals]: [color]:[linetype]:[linewidht]:[markertype]:[markersize] values for the data-lines
[color]: '@' list of colors of the lines (single value for the same, or 'auto')
[linetype]: '@' list of type of the lines (single value for the same, or 'auto')
[linewidht]: '@' list of width of the lines (single value for the same, or 'auto')
[markertype]: '@' list of type of the markers (single value for the same, or 'auto')
[markersize]: '@' list of size of the markers (single value for the same, or 'auto')
[timevals]: [timen]|[units]|[kind]|[tfmt] time labels characteristics
[timen]; name of the time variable
[units]; units string according to CF conventions ([tunits] since
[YYYY]-[MM]-[DD] [[HH]:[MI]:[SS]], '!' for spaces)
[kind]; kind of output
'Nval': according to a given number of values as 'Nval',[Nval]
'exct': according to an exact time unit as 'exct',[tunit];
tunit= [Nunits],[tu]; [tu]= 'c': centuries, 'y': year, 'm': month,
'w': week, 'd': day, 'h': hour, 'i': minute, 's': second,
'l': milisecond
[tfmt]; desired format
[valskind]: way to compute the values to be plotted
'FromFile': from the files with the 'data' flag
'FromEns',[usedims],[stats]: from all the files with the 'ensemble' flag with statistics [stats],
throughout dimensions [usedims]
[usedims]= '@' list of dimensions to take to compute the statistics
[stats]= statistics to use
'direct': without statistics
'mean': mean value of all files with ensemble flag
'median': median value of all files with ensemble flag
'weightedmean,[stats],[weights]': weighted mean value all files with ensemble flag
[stats] = statistics to use from each member: 'direct', 'mean', 'median'
[weights] = '@' list of weights [0,1]
'FromIndivEns',[usedims],[stats]: from each files with the 'ensemble' flag with statistics [stats],
throughout dimensions [usedims]
[usedims]= '@' list of dimensions to take to compute the statistics
[stats]= statistics to use:
'mean': mean value of each file with ensemble flag
'median': median value of each file with ensemble flag
[labels]: ',' list of labels for the legend ('None' for no-labels, '!' for spaces)
[figenskind]=[kindname],[valueskind] kind of ensemble plot
'cloud',[color],[limits]: ensembles values as a cloud
[color]: name of the color to use ('auto' for '#DDDDDD')
[limits]: which stadistics to use for the plot ('auto' for extremes)
'5quantile': 5quantile,95quantile of the values
'direct': below and up values are tacking directly from the provided data
'extremes': min,max of the values
'stddev': standard deviation fo the values (only available if there is one line)
'spaghetti',[color],[width],[type]: a line for each member
[color]: name of the color to use for all lines ('auto' for '#DDDDDD')
[width]: width of the line to use for all lines ('auto' for 0.5)
[type]: type of the line to use for all lines ('auto' for '-')
[graphtit]: title of the figure '!' for spaces
[legvalues]=[locleg]|[fontsize]:
[locleg]: location of the legend (0, autmoatic)
1: 'upper right', 2: 'upper left', 3: 'lower left', 4: 'lower right',
5: 'right', 6: 'center left', 7: 'center right', 8: 'lower center',
9: 'upper center', 10: 'center'
[fontsize]: font size for the legend (auto for 12)
[kindfig]: kind of output figure (png, pdf, ps)
[fign]: name of the figure (without extension)
[outfile]: wehter or not a file with the result of the computation of the ensemble should be
written or not (labelled 'ensemble_ofile.nc')
[close]: whether figure should be close or not
$ python ${pyHOME}/drawing.py -o draw_ensembles_time -f 'Forcing_rain_SESA.nc;time_counter,-1|stats,0;time_centered;rain_stats;ensemble#Forcing_rain_SESA.nc;time_counter,-1|stats,1;time_centered;rain_stats;ensemble#Forcing_rain_SESA.nc;time_counter,-1|stats,2;time_centered;rain_stats;data' -S 'auto;y;tas;mm/d;auto,auto;#0000FF:-:2.:,:2.;year|seconds!since!1958-01-01!00:00:00|exct,5,y|%Y;FromFile;$pr_{space-mean}^{SESA}$;cloud,#AAAADD,auto;$pr_{space-mean}$!evolution!at!SESA;0|10;png;Forcing_pr_spacemean;no;yes'
$ python drawing.py -o ${pyHOME}/draw_ensembles_time -f 'wrfout_d01_1995-01-01_00:00:00;Time,-1|south_north,0@44@1;WRFtime;T2;ensemble#wrfout_d01_1995-01-01_00:00:00;Time,-1|south_north,45@99@1;WRFtime;T2;ensemble' -S 'auto;y;tas;K;auto,auto;#FF0000@#0000FF:-:2.:,:2.;day|seconds!since!1949-12-01!00:00:00|exct,2,d|%D;FromIndivEns,south_north@west_east,mean;$T2_{space-mean}^{Shalf}$,$T2_{space-mean}^{Nhalf}$;cloud,#DD6464,auto;$T2_{space-mean}$ evolution;0|10;png;T2_spaceensembleIndiv;yes;yes'