2D shadow line time
Plotting a 2D field with shading and a 1D field with a line where one of the axis has temporal values
Plot of a shading 2D field and a 1D variable evolving in time draw_2D_shad_line(ncfile, values, varn)
  ncfile= [ncfiles],[ncfilel] file to use for the shading and for the line
  values= [vnamefs],[vnamefl]:[dimvals]:[dimxvn]:[dimyvn]:[dimxyfmt]:[colorbarvals]:[smin/axv]:[linevalues]:[figt]:[kindfig]:[reverse]:[mapv]:[close]
    [vnamefs]: Name in the figure of the variable to be shaded
    [vnamefl]: Name in the figure of the variable to be lined
    [dimvals]: ',' list of [dimname]|[value] telling at which dimension of the
      variable a given value is required (-1, all the length)
    [dimx/yvn]: name of the variables with the values of the final dimensions (x,y)
    [dimxyfmt]=[dxs],[dxf],[Ndx],[ordx],[dys],[dyf],[Ndy],[ordx]: format of the values at each axis (or 'auto')
      [dxs]: style of x-axis ('auto' for 'pretty')
        'Nfix', values computed at even 'Ndx'
        '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 ('auto' for '%5g')
      [Ndx]: Number of ticks at the x-axis ('auto' for 5)
      [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 ('auto' for '%5g')
      [Ndy]: Number of ticks at the y-axis ('auto' for 5)
      [ordy]: angle of orientation of ticks at the y-axis ('auto' for horizontal)
    [colorbarvals]=[colbarn],[fmtcolorbar],[orientation]
      [colorbarn]: name of the color bar
      [fmtcolorbar]: format of the numbers in the color bar 'C'-like ('auto' for %6g)
      [orientation]: orientation of the colorbar ('vertical' (default, by 'auto'), 'horizontal')
    [smin/axv]: minimum and maximum value for the shading or:
      'Srange': for full range
      'Saroundmean@val': for mean-xtrm,mean+xtrm where xtrm = np.min(mean-min@val,max@val-mean)
      'Saroundminmax@val': for min*val,max*val
      'Saroundpercentile@val': for median-xtrm,median+xtrm where xtrm = np.min(median-percentile_(val),
        percentile_(100-val)-median)
      'Smean@val': for -xtrm,xtrm where xtrm = np.min(mean-min*@val,max*@val-mean)
      'Smedian@val': for -xtrm,xtrm where xtrm = np.min(median-min@val,max@val-median)
      'Spercentile@val': for -xtrm,xtrm where xtrm = np.min(median-percentile_(val),
        percentile_(100-val)-median)
    [linevalues]=[colline],[sminl],[smaxl],[dls],[dlf],[Ndl],[ordl]
      [colline]: name of the color for the line
      [smin/axv]: minimum and maximum value for the line or:
        'Srange': for full range
        'Saroundmean@val': for mean-xtrm,mean+xtrm where xtrm = np.min(mean-min@val,max@val-mean)
        'Saroundminmax@val': for min*val,max*val
        'Saroundpercentile@val': for median-xtrm,median+xtrm where xtrm = np.min(median-percentile_(val),
          percentile_(100-val)-median)
        'Smean@val': for -xtrm,xtrm where xtrm = np.min(mean-min*@val,max*@val-mean)
        'Smedian@val': for -xtrm,xtrm where xtrm = np.min(median-min@val,max@val-median)
        'Spercentile@val': for -xtrm,xtrm where xtrm = np.min(median-percentile_(val),
          percentile_(100-val)-median)
      [dls]: style of line-axis ('auto' for 'pretty')
      [dlf]: format of the labels at the line-axis ('auto' for '%5g')
      [Ndl]: Number of ticks at the line-axis ('auto' for 5)
      [ordl]: angle of orientation of ticks at the line-axis ('auto' for horizontal)
    [smin/axv]: minimum and maximum value for the shading
    [figt]: title of the figure ('|' for spaces)
    [kindfig]: kind of figure
    [reverse]: Transformation of the values
      * 'transpose': reverse the axes (x-->y, y-->x)
      * 'flip'@[x/y]: flip the axis x or y
  [timevals]: [timen]|[units]|[kind]|[tfmt]|[label]|[timeaxis] 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     [label]; label at the graph ('!' for spaces)     [close]: Whether figure should be finished or not
  varn= [varsn],[varnl] name of the variable to plot with shading and with line
$ python ${pyHOME}/drawing.py -o draw_2D_shad_line_time -f wrfout_d01_2001-11-11_00:00:00,wrfout_d01_2001-11-11_00:00:00 -S 'hus,ps;west_east|10,south_north|96,Time|-1,bottom_top|-1;WRFtime;ZNU;auto;rainbow,auto,horizontal;Srange,Srange;k,Srange,Srange,auto,auto,auto,45.;vert.|sec.|hus|&|ps|at|x=10,y=96|temporal|evolution;png;flip@y;WRFtime|hours!since!1949-12-01!00:00:00|Nval,7|%H$^{%M}$|time!($[HH]^{[MM]}$);True' -v QVAPOR,PSFC