2D shadow plot
Plotting a 2D field with shading
draw_2D_shad(ncfiles, values, varn)
  ncfiles= ';' list of files to use as [filen]@[varn]@[dimvals]
    [filen]: name of the file
    [varn]: name of the variable in the file to be shaded
    [dimvals]: ',' list of [dimname]|[value] telling at which dimension of the
      variable a given value is required:
        * [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
  values=[vnamefs]:[dimxvn]:[dimyvn]:[dimxyfmt]:[colorbarvals]:[sminv],[smaxv]:
    [kindfig]:[reverse]:[mapv]:[labels]:[Nrow]:[Ncol]:[globaltitle]:[close]
    [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)
    [kindfig]: kind of figure output file: png, pdf, ...
    [reverse]: Transformation of the values
      * 'transpose': reverse the axes (x-->y, y-->x)
      * 'flip'@[x/y]: flip the axis x or y
    [mapv]: map characteristics: [proj],[res]
      see full documentation: http://matplotlib.org/basemap
      [proj]: projection
        * 'cyl', cilindric
        * 'lcc', lambert conformal
      [res]: resolution:
        * 'c', crude
        * 'l', low
        * 'i', intermediate
        * 'h', high
        * 'f', full
    [labels]: ',' list of consecutive labels for each individual panel ('!' for spaces)     [Nrow]: number of rows in the figure     [Ncol]: number of columns in the figure     [globaltitle]: global title of the figure ('!' for spaces)     [close]: Whether figure should be finished or not
python ${pyHOME}/drawing.py -o draw_multi_2D_shad -f '../PY/wrfout_d01_1995-01-01_00:00:00@T2@west_east|
  -1,south_north|-1,Time|0;../PY/wrfout_d01_1995-01-01_00:00:00@T2@west_east|-1,south_north|-1,Time|1;
  ../PY/wrfout_d01_1995-01-01_00:00:00@T2@west_east|-1,south_north|-1,Time|2;../PY/wrfout_d01_1995-01
  -01_00:00:00@T2@west_east|-1,south_north|-1,Time|3' -S 'tas:XLONG:XLAT:auto:rainbow,auto,auto:
  Srange,Srange:png:None:cyl,l:0!UTC,1!UTC,2!UTC,3!UTC:2:2:tas!at!2001-11-11:True'