Index: trunk/tools/documentation/plotting/2Dshad_contdisc.html
===================================================================
--- trunk/tools/documentation/plotting/2Dshad_contdisc.html	(revision 2128)
+++ trunk/tools/documentation/plotting/2Dshad_contdisc.html	(revision 2128)
@@ -0,0 +1,107 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      2D shadow plot with discontinous one
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting one continuous fields with shading and another discrete one with points
+      <IMG WIDTH="100%" SRC="figs/obs-sim_sndmap_20121027000000_PRW_50lev_assigned.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+    draw_2D_shad_contdisc(ncfiles, values)<BR>
+    &nbsp;&nbsp;  ncfiles= [contfilen];[contvarn];[dimxvarn];[dimvvarn];[dimvals]@[discfilen];<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[discvarn];[dimxvarn];[dimyvarn];[dimvals] files and variables to use<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[contfilen]: name of the file with the continuous varible<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[contvarn]: name of the continuos variable<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimxvarn]: name of the variable with values for the x-dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvvarn]: name of the variable with values for the y-dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvals]: ',' list of [dimname]|[value] values to slice the variable:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [integer]: which value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -1: all along the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -9: last value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [beg]@[end]@[inc] slice from [beg] to [end] every [inc]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE, no dim name all the dimension size<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discfilen]: name of the file with the discrete varible<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discvarn]: name of the discrete variable<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;* NOTE: limits of the graph will be computed from the continuous variable<BR>
+    &nbsp;&nbsp;values=[vnamefs]:[dvarxn],[dvaryn]:[dimxyfmt]:[colorbarvals]:[sminv],[smaxv]:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discvals]:[figt]:[kindfig]:[reverse]:[mapv]:[plotrange]:[close]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[vnamefs]: Name in the figure of the variable to be shaded<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dvarxn],[dvaryn]: name of the dimensions for the final x-axis and y-axis at <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;the figure (from contfilen)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimxyfmt]=[dxs],[dxf],[Ndx],[ordx],[dys],[dyf],[Ndy],[ordx]: format of the <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values at each axis (or 'auto')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dxs]: style of x-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nfix', values computed at even 'Ndx'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Vfix', values computed at even 'Ndx' increments<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'pretty', values computed following aprox. 'Ndx' at 'pretty' intervals (2.,2.5,4,5,10)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dxf]: format of the labels at the x-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndx]: Number of ticks at the x-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordx]: angle of orientation of ticks at the x-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dys]: style of y-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dyf]: format of the labels at the y-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndy]: Number of ticks at the y-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordy]: angle of orientation of ticks at the y-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[colorbarvals]=[colbarn],[fmtcolorbar],[orientation]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colorbarn]: name of the color bar<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[fmtcolorbar]: format of the numbers in the color bar 'C'-like ('auto' for %6g)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[orientation]: orientation of the colorbar ('vertical' ['auto'], 'horizontal')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE: single 'auto' for 'rainbow,%6g,vertical'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[smin/axv]: minimum and maximum values for the shading or string for each for:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Srange': for full range<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundmean@val': for mean-xtrm,mean+xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;xtrm = np.min(mean-min@val,max@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundminmax@val': for min*val,max*val<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundpercentile@val': for median-xtrm,median+xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;xtrm = np.min(median-percentile_(val),percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smean@val': for -xtrm,xtrm where xtrm = np.min(mean-min*@val,max*@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smedian@val': for -xtrm,xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xtrm = np.min(median-min@val,max@val-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Spercentile@val': for -xtrm,xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xtrm = np.min(median-percentile_(val),percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discvals]= [type],[size],[lwidth],[lcol] characteristics of the points for the discrete field<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[type]: type of point. Any marker from matoplib must be filled !<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[size]: size of point <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lwidth]: width of the line around the point<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lcol]: color of the line around the point<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'auto': for [type]='o', [size]=5, [lwdith]=0.25, [lcol]='#000000'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[figt]: title of the figure ('!' for spaces)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[kindfig]: kind of figure output (ps, png, pdf)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[reverse]: Transformation of the values to plot<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'transpose': reverse the axes (x-->y, y-->x)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'flip'@[x/y]: flip the axis x or y<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[mapv]: map characteristics: [proj],[res]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;see full documentation: http://matplotlib.org/basemap/<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[proj]: projection<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'cyl', cilindric<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'lcc', lambert conformal<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[res]: resolution of the coastaline data:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'c', crude<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'l', low<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'i', intermediate<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'h', high<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'f', full<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[plotrange]: range of the plot<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'strict': map covers only the minimum and maximum lon,lats from the <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;locations of the discrtete points<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'sponge,'[dlon],[dlat]: map covers an extended [dlon],[dlat] from the <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;minimum and maximum lon,lats from the locations of the discrtete points<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'fullcontinuous': map covers all the shadding area<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'lonlatbox,[lonSW],[latSW],[lonNE],[latNE]': plotted map only covers a lon,lat box<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[close]: Whether figure should be finished or not<BR>
+    </DIV>
+    <DIV CLASS="valins">
+      $ python ${pyHOME}/drawing.py -o draw_2D_shad_contdisc -f 'simcdx_vars_cape_50lev_assigned.nc;PRW;LON;LAT;<BR>
+&nbsp;&nbsp;Time|23,time|23,west_east|-1,south_north|-1@all_sounding_1D.nc;prw;stslon;stslat;time|14,lon|-1,lat|-1'<BR>
+&nbsp;&nbsp;-S 'PRW:west_east,south_north:auto:Blues,auto,auto:0.,35.:auto:obs!,!50leva!PRW!on!2012/10/27!00!UTC:png:<BR>
+&nbsp;&nbsp;None:cyl,i:lonlatbox,6.,44.,12.,48.5:yes'
+    </DIV>
+  </BODY>
+</HTML>
+
Index: trunk/tools/documentation/plotting/2Dshad_contdisc_time.html
===================================================================
--- trunk/tools/documentation/plotting/2Dshad_contdisc_time.html	(revision 2128)
+++ trunk/tools/documentation/plotting/2Dshad_contdisc_time.html	(revision 2128)
@@ -0,0 +1,108 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      2D shadow plot
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting one continuous fields with shading and another discrete one with points with a time-axis
+      <IMG WIDTH="100%" SRC="figs/2Dshad_obs-sim_comparison_time.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+    draw_2D_shad_contdisc_time(ncfiles, values)<BR>
+    &nbsp;&nbsp;  ncfiles= [contfilen];[contvarn];[dimtvarn];[dimvvarn];[dimvals]@[discfilen];<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[discvarn];[dimtvarn];[dimvvarn];[dimvals] files and variables to use<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[contfilen]: name of the file with the continuous varible<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[contvarn]: name of the continuos variable<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimtvarn]: name of the variable with values for the time-dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvvarn]: name of the variable with values for the variable-dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvals]: ',' list of [dimname]|[value] values to slice the variable:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [integer]: which value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -1: all along the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -9: last value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [beg]@[end]@[inc] slice from [beg] to [end] every [inc]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE, no dim name all the dimension size<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discfilen]: name of the file with the discrete varible<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discvarn]: name of the discrete variable<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;* NOTE: limits of the graph will be computed from the continuous variable<BR>
+    &nbsp;&nbsp;  values=[vnamefs];[varaxis];[timevals];[dimvfmt];[colorbarvals];[sminv],[smaxv];<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;   [discvals];[figt];[kindfig];[reverse];[mapv];[close]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[vnamefs]: Name in the figure of the variable to be shaded<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[varaxis]: axis in the figure for the variable-axis ('x' or 'y')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[timevals]: [units]|[kind]|[tfmt]|[label] time labels characteristics<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[units]: units string according to CF conventions ([tunits] since <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[YYYY]-[MM]-[DD] [[HH]:[MI]:[SS]], '!' for spaces), <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'auto' for minutes!since!1949-12-01!00:00:00<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[kind]: kind of output<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nval': according to a given number of values as 'Nval',[Nval]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'exct': according to an exact time unit as 'exct',[tunit]; <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;  tunit= [Nunits],[tu]; [tu]= 'c': centuries, 'y': year, 'm': month,<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;   'w': week, 'd': day, 'h': hour, 'i': minute, 's': second, <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;   'l': milisecond<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[tfmt]: desired format (combination of 'C'-style values and LaTeX)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[label]: label at the graph ('!' for spaces, combination of 'C'-style <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values and LaTeX)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvfmt]=[[dvs],[dvf],[Ndv],[ordv]: format of the values at variable-axis <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(or 'auto')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dvs]: style of variable-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nfix', values computed at even 'Ndx'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Vfix', values computed at even 'Ndx' increments<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'pretty', values computed following aprox. 'Ndx' at 'pretty' intervals (2.,2.5,4,5,10)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dvf]: format of the labels at the var-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndv]: Number of ticks at the var-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordv]: angle of orientation of ticks at the var-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[colorbarvals]=[colbarn],[fmtcolorbar],[orientation]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colorbarn]: name of the color bar<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[fmtcolorbar]: format of the numbers in the color bar 'C'-like ('auto' for %6g)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[orientation]: orientation of the colorbar ('vertical' ['auto'], 'horizontal')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE: single 'auto' for 'rainbow,%6g,vertical'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[smin/axv]: minimum and maximum values for the shading or string for each for:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Srange': for full range<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundmean@val': for mean-xtrm,mean+xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;xtrm = np.min(mean-min@val,max@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundminmax@val': for min*val,max*val<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundpercentile@val': for median-xtrm,median+xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;xtrm = np.min(median-percentile_(val),percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smean@val': for -xtrm,xtrm where xtrm = np.min(mean-min*@val,max*@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smedian@val': for -xtrm,xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xtrm = np.min(median-min@val,max@val-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Spercentile@val': for -xtrm,xtrm where <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xtrm = np.min(median-percentile_(val),percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[discvals]= [type],[size],[lwidth],[lcol] characteristics of the points for the discrete field<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[type]: type of point. Any marker from matoplib must be filled !<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[size]: size of point <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lwidth]: width of the line around the point<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lcol]: color of the line around the point<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'auto': for [type]='o', [size]=5, [lwdith]=0.25, [lcol]='#000000'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[figt]: title of the figure ('!' for spaces)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[kindfig]: kind of figure output (ps, png, pdf)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[reverse]: Transformation of the values to plot<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'transpose': reverse the axes (x-->y, y-->x)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'flip'@[x/y]: flip the axis x or y<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[mapv]: map characteristics: [proj],[res]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;see full documentation: http://matplotlib.org/basemap/<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[proj]: projection<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'cyl', cilindric<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'lcc', lambert conformal<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[res]: resolution of the coastaline data:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'c', crude<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'l', low<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'i', intermediate<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'h', high<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'f', full<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[close]: Whether figure should be finished or not<BR>
+    </DIV>
+    <DIV CLASS="valins">
+      $ python ${pyHOME}/drawing.py -o draw_2D_shad_contdisc_time -f 'simout_snddiags.nc;ta;time;pres;<BR>
+&nbsp;&nbsp;time|-1,pres|-1@UWyoming_snd_87576.nc;ta;time;pres;time|-1,pres|-1' -S 'ta;y;auto|exct,12,h|%d$^{%H}$|<BR>
+&nbsp;&nbsp;time!($[DD]^{[HH]}$);Vfix,auto,50.,auto;auto;Srange,Srange;auto;obs!&!sim!Ezeiza!airport!sounding;png;<BR>
+&nbsp;&nbsp;flip@y;None;yes'
+    </DIV>
+  </BODY>
+</HTML>
+
Index: trunk/tools/documentation/plotting/2Dshad_map.html
===================================================================
--- trunk/tools/documentation/plotting/2Dshad_map.html	(revision 2128)
+++ trunk/tools/documentation/plotting/2Dshad_map.html	(revision 2128)
@@ -0,0 +1,106 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      2D shadow map
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting a 2D maksed field with shading on top of a map
+      <IMG WIDTH="100%" SRC="figs/shad-map_convini_Case6_era5_1a.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+    draw_2Dshad_map(ncfile, values, varn)<BR>
+    plotting a shadow field with a background map<BR>
+    &nbsp;&nbsp;ncfile= file to use<BR>
+    &nbsp;&nbsp;values=[vnamemap];[vnamefs];[dimvals];[dimxvn];[dimyvn];[dimxyfmt];<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colorbarvals];[sminv],[smaxv];[figt];[kindfig];[mapkind];[mapv];[close]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[vnamemap]: Name of the variable to be used for map<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[vnamefs]: Name in the figure of the variable to be shaded<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvals]: ',' list of [dimname]|[value] telling at which dimension of the<BR> 
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;variable a given value is required:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [integer]: which value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -1: all along the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -9: last value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [beg]@[end]@[inc] slice from [beg] to [end] every [inc]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE, no dim name all the dimension size<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimx/yvn]: name of the variables with the values of the final dimensions (x,y)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimxyfmt]=[dxs],[dxf],[Ndx],[ordx],[dys],[dyf],[Ndy],[ordx]: format of the values at each axis (or 'auto')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dxs]: style of x-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nfix', values computed at even 'Ndx'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Vfix', values computed at even 'Ndx' increments<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'pretty', values computed following aprox. 'Ndx' at 'pretty' intervals (2.,2.5,4,5,10)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dxf]: format of the labels at the x-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndx]: Number of ticks at the x-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordx]: angle of orientation of ticks at the x-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dys]: style of y-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dyf]: format of the labels at the y-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndy]: Number of ticks at the y-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordy]: angle of orientation of ticks at the y-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[colorbarvals]=[colbarn]#[fmtcolorbar]#[orientation]#[colorbarticks]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colorbarn]: name of the color bar for the values to plot<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[fmtcolorbar]: format of the numbers in the color bar 'C'-like ('auto' for %6g)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[orientation]: orientation of the colorbar ('vertical' (default, by 'auto'), 'horizontal')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colorbarticks]: kind of colorbarticks<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'direct': direct values<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'time'@[units]|[kind]|[tfmt]|[label] time labels characteristics<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[units]: units string according to CF conventions ([tunits] since <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[YYYY]-[MM]-[DD] [[HH]:[MI]:[SS]], '!' for spaces), <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'auto' for minutes!since!1949-12-01!00:00:00<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[kind]: kind of output<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nval': according to a given number of values as 'Nval',[Nval]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'exct': according to an exact time unit as 'exct',[tunit]; <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;tunit= [Nunits],[tu]; [tu]= 'c': centuries, 'y': year, 'm': month,<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'w': week, 'd': day, 'h': hour, 'i': minute, 's': second, <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'l': milisecond<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[tfmt]: desired format (combination of 'C'-style values and LaTeX)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[label]: label at the graph ('!' for spaces, combination of 'C'-style <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values and LaTeX)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[smin/axv]: minimum and maximum value for the shading or:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Srange': for full range<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundmean@val': for mean-xtrm,mean+xtrm where xtrm = np.min(mean-min@val,max@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundminmax@val': for min*val,max*val<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundpercentile@val': for median-xtrm,median+xtrm where xtrm = np.min(median-percentile_(val),<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smean@val': for -xtrm,xtrm where xtrm = np.min(mean-min*@val,max*@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smedian@val': for -xtrm,xtrm where xtrm = np.min(median-min@val,max@val-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Spercentile@val': for -xtrm,xtrm where xtrm = np.min(median-percentile_(val),<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[figt]: title of the figure ('|' for spaces)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[kindfig]: kind of file for the figure output (ps, eps, png, pdf, ...)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[mapfmt] = [cbarmap]|[mapkind]|[lonlatbox] values for the map <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;cbarmap: name of the colorbar for the map<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;mapkind: kind of map to use in the plot<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'direct': values are they come<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'shadow',[pos],[enhance]: pseudo-shadding from a given location of the sun<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[pos]: 'N', 'NW'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[enhance]: enhance factor for the shading ('auto' for 1./5.)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;lonlatbox: [lonSW],[latSW],[lonNE],[latNE] to plot only a lon,lat box<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'full': for the whole domain<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[mapv]: map characteristics: [proj],[res]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;see full documentation: <A HREF="http://matplotlib.org/basemap/" CLASS="la" TARGET="_blank">http://matplotlib.org/basemap/</A><BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[proj]: projection<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'cyl', cilindric<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'lcc', lambert conformal<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[res]: resolution:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'c', crude<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'l', low<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'i', intermediate<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'h', high<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'f', full<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[close]: Whether figure should be finished or not<BR>
+    &nbsp;&nbsp;varn= [varsn] name of the variable to plot with shading<BR>
+    </DIV>
+    <DIV CLASS="valins">
+      $ python ${pyHOME}/drawing.py -o draw_2Dshad_map -S 'orog;$conv^{ini}$;lat|-1,lon|-1;lon;lat;<BR>
+    &nbsp;&nbsp;auto;rainbow#auto#auto#time@minutes!since!1949-12-01!00:00:00|exct,6,h|%d$^{%H}$|<BR>
+    &nbsp;&nbsp;date!([DD]$^{[HH]}$);Srange,Srange;convini|Case6|era5_1a;png;binary|shadow,NW,auto|<BR>
+    &nbsp;&nbsp;-70.,-36,-62.,-30.;cyl,c;yes' -v convini -f UBA_ERA5_1a_2D.nc
+    </DIV>
+  </BODY>
+</HTML>
+
Index: trunk/tools/documentation/plotting/Nlines_time.html
===================================================================
--- trunk/tools/documentation/plotting/Nlines_time.html	(revision 1948)
+++ trunk/tools/documentation/plotting/Nlines_time.html	(revision 2128)
@@ -57,4 +57,5 @@
     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;9: 'upper center', 10: 'center'<BR>
     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[fontsize]: font size for the legend (auto for 12)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ncols]: number of columns in the legend<BR>
     &nbsp;&nbsp;&nbsp;&nbsp;[graphk]: kind of the graphic<BR>
     &nbsp;&nbsp;&nbsp;&nbsp;[lines]: ',' list of type of lines, None for automatic, single value all the same<BR>
@@ -67,5 +68,7 @@
     &nbsp;&nbsp;&nbsp;&nbsp;[pointsizes]: ',' list of widths for the lines, None for automatic, single<BR>
     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;value all the same<BR>
-    &nbsp;&nbsp;&nbsp;&nbsp;[pointfreq]: frequency of point plotting, 'all' for all time steps<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[pointfreq]: frequency of point plotting (every given number of values),<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'all' to plot all values using all time steps<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'runmean',[Nsteps]: computing a running mean of [Nsteps] along all values
     &nbsp;&nbsp;&nbsp;&nbsp;[period]: which period to plot<BR>
     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'-1': all period<BR>
Index: trunk/tools/documentation/plotting/WRFeta_levels.html
===================================================================
--- trunk/tools/documentation/plotting/WRFeta_levels.html	(revision 2128)
+++ trunk/tools/documentation/plotting/WRFeta_levels.html	(revision 2128)
@@ -0,0 +1,48 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      WRF eta levels
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting WRF eta levels distribution
+      <IMG WIDTH="100%" SRC="figs/WRFeta_levels.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+    draw_WRFeta_levels(filenames, values)<BR>
+    Function to plot vertical levels from WRF wrfinput_d[nn] file<BR>
+    &nbsp;&nbsp;filenames= ',' list of name of files to use<BR>
+    &nbsp;&nbsp;values = [labs]:[colors]:[markers]:[legvals]:[imgtit]:[imgkind]:[newfile]:[close]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[labs]= ',' list of labels for the plot (LaTeX like)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[colors]= ',' list of colors for the lines ('auto' for automatic)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[markers]= '@' list of markers for the lines ('auto' for automatic)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[legvals]=[loclegend]|[fonstsize]|[ncol] values for the legend<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[locleg]: location of the legend (0, automatic)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;1: 'upper right', 2: 'upper left', 3: 'lower left', 4: 'lower right',<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;5: 'right', 6: 'center left', 7: 'center right', 8: 'lower center',<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;9: 'upper center', 10: 'center'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[fontsize]: font size for the legend ('auto' for 12)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ncol]: number of columns ('auto' for 1)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[imgtit]= title of the image ('!' for spaces)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[imgkind]= kind of file output of the image (ps, pns, pdf, ...)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[newfile]= whether the possible pre-existing file with the required data be <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;should removed or not ('yes/no')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[close]= whether figure should be closed or not<BR>
+    &nbsp;&nbsp;* NOTE: A multi-panel plot will be drawn called `WRFeta_levels.png', and also<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;an individual plot for each panel will be produced
+    </DIV>
+    <DIV CLASS="valins">
+      $ python $pyHOME/drawing.py -o draw_WRFeta_levels -f '120lev_cdxwrf2/wrfinput_d02,<BR>
+      &nbsp;&nbsp;120lev_assigned/wrfinput_d02,80lev_cdxwrf2/wrfinput_d02,80lev_assigned/wrfinput_d02,<BR>
+      &nbsp;&nbsp;50lev_cdxwrf2/wrfinput_d02,50lev_assigned/wrfinput_d02,38lev_cdxwrf2/wrfinput_d02,<BR>
+      &nbsp;&nbsp;38lev_assigned/wrfinput_d02' -S '120lev,120leva,80lev,8leva,50lev,50leva,38lev,<BR>
+      &nbsp;&nbsp;38leva:auto:auto:0|4|2:FPS!Alps!vertical!levels:png:no:true'
+    </DIV>
+  </BODY>
+</HTML>
+
Index: trunk/tools/documentation/plotting/gallery.html
===================================================================
--- trunk/tools/documentation/plotting/gallery.html	(revision 1948)
+++ trunk/tools/documentation/plotting/gallery.html	(revision 2128)
@@ -15,9 +15,11 @@
       <A HREF="2Dshad_time.html"><IMG SRC="figs/ta_valley_fixpixel-grid.png" ALT="2D shadow time" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="2Dshad_cont.html"><IMG SRC="figs/huss-tas_Blues_shad-cont.png" ALT="2D shadow contour" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
-      <A HREF="2Dshad_contdisc_time.html"><IMG SRC="figs/2Dshad_obs-sim_comparison_time.png" ALT="2D shadow continuous and discrete" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
+      <A HREF="2Dshad_contdisc.html"><IMG SRC="figs/obs-sim_sndmap_20121027000000_PRW_50lev_assigned.png" ALT="2D shadow continuous and discrete" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
+      <A HREF="2Dshad_contdisc_time.html"><IMG SRC="figs/2Dshad_obs-sim_comparison_time.png" ALT="2D shadow continuous and discrete with time" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="2Dshad_2cont.html"><IMG SRC="figs/2Dfields_shadow-2contour.png" ALT="2D shadow 2 contour" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="2Dshad_cont_time.html"><IMG SRC="figs/huss-tas_shad_cont_time_BuPu.png" ALT="2D shadow contour time" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="2Dshad_line.html"><IMG SRC="figs/hus-hgt_vert-sec_shad-line.png" ALT="2D shadow line" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="2Dshad_line_time.html"><IMG SRC="figs/hus-ps_vert-sec_shad-line_time.png" ALT="2D shadow line time" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
+      <A HREF="2Dshad_map.html"><IMG SRC="figs/shad-map_convini_Case6_era5_1a.png" ALT="2D maksed shadow on top of a map" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="draw_bar.html"><IMG SRC="figs/bar_rain_CRU-NCEP_norm.png" ALT="bar char plot" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="bar_line.html"><IMG SRC="figs/bar-line_rain_CRU-NCEP_norm.png" ALT="bar char and lines plot" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
@@ -33,4 +35,5 @@
       <A HREF="movie_2Dshad.html"><IMG SRC="figs/movie_2D_shadow.gif" ALT="movie 2D shadow" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="multi_2Dshad.html"><IMG SRC="figs/multi_2Dfields_shadow.png" ALT="multi 2D shadow" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
+      <A HREF="multiWindRose.html"><IMG SRC="figs/WindRose_obs-sim_step_Case6_87418_20151118000000.png" ALT="multi Wind Rose" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="Neighbourgh_evol.html"><IMG SRC="figs/Neighbourghood_evol.png" ALT="Neighbourghood evolution" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="Nlines.html"><IMG SRC="figs/Nlines_hurs_wrfout_tturb-xmean-last.png" ALT="Nlines" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
@@ -42,4 +45,5 @@
       <A HREF="pointsZval.html"><IMG SRC="figs/pointsZval.png" ALT="plot points by Z value" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="river_desc.html"><IMG SRC="figs/or_rivers.png" ALT="plot rivers from ORCHIDEE's description" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>     
+      <A HREF="stations_map.html"><IMG SRC="figs/stations_map.png" ALT="plot stations map" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>     
       <A HREF="subbasin.html"><IMG SRC="figs/Caceres_subbasin.png" ALT="plot sub-basin from ORCHIDEE's routing" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>     
       <A HREF="Taylor.html"><IMG SRC="figs/Taylor.png" ALT="Taylor diagram" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
@@ -58,4 +62,5 @@
       <A HREF="WindRose.html"><IMG SRC="figs/WindRose_anglespeedfreq.png" ALT="WindRose with frequencies of angle and wind speed" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="WindRose.html"><IMG SRC="figs/WindRose_linepoint-multicoltime.png" ALT="WindRose of line-mark as consecutive values changing color" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
+      <A HREF="WRFeta_levels.html"><IMG SRC="figs/WRFeta_levels.png" ALT="WRF eta-levels" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
     </DIV>
   </BODY>
Index: trunk/tools/documentation/plotting/multiWindRose.html
===================================================================
--- trunk/tools/documentation/plotting/multiWindRose.html	(revision 2128)
+++ trunk/tools/documentation/plotting/multiWindRose.html	(revision 2128)
@@ -0,0 +1,82 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      2D shadow plot
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting multiple WindRose at the same time
+      <IMG WIDTH="100%" SRC="figs/WindRose_obs-sim_step_Case6_87418_20151118000000.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+      draw_multiWindRose(ncfiles, values)<BR>
+      Function to plot multiple wind rose (from where the dinw blows)<BR>
+      &nbsp;&nbsp;ncfiles='#' separated list of [filen]@[dimvals]@[uvarn],[vvarn]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;[filen]: name of the file<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;[dimvals]: ';' list of [dimn]|[dvalue] dimension and slice along dimension to retrieve the winds<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dimn]: name of the dimension<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dvalue]: value for the slice in the given dimension<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [integer]: which value of the dimension<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -1: all along the dimension<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -9: last value of the dimension<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [beg],[end],[freq] slice from [beg] to [end] every [freq]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE, no dim name all the dimension size<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;No value takes all the range of the dimension<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;[uvarn],[vvarn]: name of the x and y wind components<BR>
+      &nbsp;&nbsp;values=[kindRose]:[labels]:[imgtit]:[imgkind]:[kindlabelsangle]:[close]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;[kindRose]: [kind];[value1];[...[valueN]] Kind of rose to plot and values of the kind<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'linepoint': consecutive (time, height, level, ...) line-point angle and speed values. Three different species<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'multicol': line-marker color changing according to a third variable [extravarn]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values: [extravarn];[lines];[markers];[freqmarkers];[colbar];[Nang]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'multicoltime': line-marker color changing according to a temporal variable [extravarn]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values: [extravarn];[lines];[markers];[colbar];[Nang];[timekind];[timefmt];[timelabel]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'singlecol': same color for the line-marker<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values: [lines];[markers];[cols];[Nang]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'scatter': a marker for each wind at different values (time, height, level, ...). Three different species<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'multicol':marker color changing according to a third variable [extravarn]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values: [extravarn];[markers];[colbar];[Nang]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'multicoltime': marker color changing according to a temporal variable [extravarn]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values: [extravarn];[lines];[markers];[colbar];[Nang];[timekind];[timefmt];[timelabel]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'singlecol': same color for all the markers<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;values: [markers];[cols];[Nang]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;meaning (where apply):<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[extravarn]: name of the extra variable <BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lines]: '@' separated list of type of lines (matplotlib ref)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[markers]: '@' separated list of type of markers to use (matplotlib ref)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[freqmarkers]: '@' separated list of frequency of drawing markers ('auto' for every 10)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[cols] = '@' separated list of colors ('#[RR][GG][BB])<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colbar]: name of the colorbar ('auto' for 'spectral_r')<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Nang]: number of angles to divide the rose ('auto' for 8)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Nspeed]: number of speeds to divide the wind speed distribution ('auto' for 8)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[maxspeed]: maximum wind speed used to compute the frequency of distributions ('auto' for 40.)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[timekind]; time computation of ticks<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nval': according to a given number of values as 'Nval',[Nval]<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'exct': according to an exact time unit as 'exct',[tunit]; <BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;tunit= [Nunits],[tu]; [tu]= 'c': centuries, 'y': year, 'm': month,<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'w': week, 'd': day, 'h': hour, 'i': minute, 's': second, <BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'l': milisecond<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[timefmt]; desired format of time labels (C-like)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[timelabel]; label of time colorbar at the graph ('!' for spaces)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;labels: ',' separated list of labels for the legend ('None' for no label)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;imgtit: title of the image ('!' for spaces)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;imgkind: kind of file output of the image (ps, pns, pdf, ...)<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;kindlabelsangle: kind of labels for the angles of the wind Rose<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'cardianals': Following combinations of 'N', 'E', 'S', 'W' according to Nang<BR>
+      &nbsp;&nbsp;&nbsp;&nbsp;close: whether figure should be closed or not<BR>
+    </DIV>
+    <DIV CLASS="valins">
+    $ python ${pyHOME}/drawing.py -o draw_multiWindRose -f 'sounding_uava_87418.nc@pres|1,1401,1;time|<BR>
+    &nbsp;&nbsp;33@ua,va#UBA_ERA-I_1a_snd_mendoza.nc@pres|-1;time|0@ua,va#UBA_ERA5_1a_snd_mendoza.nc@pres|-1;<BR>
+    &nbsp;&nbsp;time|0@ua,va#UBA_morrison_MYNN_snd_mendoza.nc@pres|-1;time|0@ua,va#SMN_snd_mendoza.nc@pres|-1;<BR>
+    &nbsp;&nbsp;time|0@ua,va#UBA_ERA-I_smnphys_snd_mendoza.nc@pres|-1;time|0@ua,va' -S 'linepoint;multicol;<BR>
+    &nbsp;&nbsp;pres;-;x@+@>@s@*@<;1@10@10@10@10@10;spectral;auto:obs,1a,era5_1a,sens,SMN,smnphys:WindRose!<BR>
+    &nbsp;&nbsp;obs!,!sim!comparison!at!87418!on!2015!Nov!18!00!UTC:png:auto:True'
+    </DIV>
+  </BODY>
+</HTML>
+
Index: trunk/tools/documentation/plotting/multi_2Dshad.html
===================================================================
--- trunk/tools/documentation/plotting/multi_2Dshad.html	(revision 2128)
+++ trunk/tools/documentation/plotting/multi_2Dshad.html	(revision 2128)
@@ -0,0 +1,87 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      2D shadow plot
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting a 2D field with shading
+      <IMG WIDTH="100%" SRC="figs/multi_2Dfields_shadow.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+    draw_2D_shad(ncfiles, values, varn)<BR>
+    &nbsp;&nbsp;ncfiles= ';' list of files to use as [filen]@[varn]@[dimvals]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[filen]: name of the file<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[varn]: name of the variable in the file to be shaded<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimvals]: ',' list of [dimname]|[value] telling at which dimension of the <BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;variable a given value is required:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [integer]: which value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -1: all along the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* -9: last value of the dimension<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* [beg]@[end]@[inc] slice from [beg] to [end] every [inc]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* NOTE, no dim name all the dimension size<BR>
+    &nbsp;&nbsp;values=[vnamefs]:[dimxvn]:[dimyvn]:[dimxyfmt]:[colorbarvals]:[sminv],[smaxv]:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[kindfig]:[reverse]:[mapv]:[labels]:[Nrow]:[Ncol]:[globaltitle]:[close]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimx/yvn]: name of the variables with the values of the final dimensions (x,y)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[dimxyfmt]=[dxs],[dxf],[Ndx],[ordx],[dys],[dyf],[Ndy],[ordx]: format of the values at each axis (or 'auto')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dxs]: style of x-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Nfix', values computed at even 'Ndx'<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Vfix', values computed at even 'Ndx' increments<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'pretty', values computed following aprox. 'Ndx' at 'pretty' intervals (2.,2.5,4,5,10)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dxf]: format of the labels at the x-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndx]: Number of ticks at the x-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordx]: angle of orientation of ticks at the x-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dys]: style of y-axis ('auto' for 'pretty')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[dyf]: format of the labels at the y-axis ('auto' for '%5g')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[Ndy]: Number of ticks at the y-axis ('auto' for 5)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[ordy]: angle of orientation of ticks at the y-axis ('auto' for horizontal)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[colorbarvals]=[colbarn],[fmtcolorbar],[orientation]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[colorbarn]: name of the color bar<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[fmtcolorbar]: format of the numbers in the color bar 'C'-like ('auto' for %6g)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[orientation]: orientation of the colorbar ('vertical' (default, by 'auto'), 'horizontal')<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[smin/axv]: minimum and maximum value for the shading or:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Srange': for full range<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundmean@val': for mean-xtrm,mean+xtrm where xtrm = np.min(mean-min@val,max@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundminmax@val': for min*val,max*val<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Saroundpercentile@val': for median-xtrm,median+xtrm where xtrm = np.min(median-percentile_(val),<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smean@val': for -xtrm,xtrm where xtrm = np.min(mean-min*@val,max*@val-mean)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Smedian@val': for -xtrm,xtrm where xtrm = np.min(median-min@val,max@val-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'Spercentile@val': for -xtrm,xtrm where xtrm = np.min(median-percentile_(val),<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;percentile_(100-val)-median)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[kindfig]: kind of figure output file: png, pdf, ...<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[reverse]: Transformation of the values<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'transpose': reverse the axes (x-->y, y-->x)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'flip'@[x/y]: flip the axis x or y<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[mapv]: map characteristics: [proj],[res]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;see full documentation: <HREF CLASS="lc" HREF="http://matplotlib.org/basemap" TARGET="_blank">http://matplotlib.org/basemap</A><BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[proj]: projection<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'cyl', cilindric<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'lcc', lambert conformal<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[res]: resolution:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'c', crude<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'l', low<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'i', intermediate<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'h', high<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'f', full<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[labels]: ',' list of consecutive labels for each individual panel ('!' for spaces)
+    &nbsp;&nbsp;&nbsp;&nbsp;[Nrow]: number of rows in the figure
+    &nbsp;&nbsp;&nbsp;&nbsp;[Ncol]: number of columns in the figure
+    &nbsp;&nbsp;&nbsp;&nbsp;[globaltitle]: global title of the figure ('!' for spaces)
+    &nbsp;&nbsp;&nbsp;&nbsp;[close]: Whether figure should be finished or not<BR>
+    </DIV>
+    <DIV CLASS="valins">
+    python ${pyHOME}/drawing.py -o draw_multi_2D_shad -f '../PY/wrfout_d01_1995-01-01_00:00:00@T2@west_east|<BR>
+    &nbsp;&nbsp;-1,south_north|-1,Time|0;../PY/wrfout_d01_1995-01-01_00:00:00@T2@west_east|-1,south_north|-1,Time|1;<BR>
+    &nbsp;&nbsp;../PY/wrfout_d01_1995-01-01_00:00:00@T2@west_east|-1,south_north|-1,Time|2;../PY/wrfout_d01_1995-01<BR>
+    &nbsp;&nbsp;-01_00:00:00@T2@west_east|-1,south_north|-1,Time|3' -S 'tas:XLONG:XLAT:auto:rainbow,auto,auto:<BR>
+    &nbsp;&nbsp;Srange,Srange:png:None:cyl,l:0!UTC,1!UTC,2!UTC,3!UTC:2:2:tas!at!2001-11-11:True'
+    </DIV>
+  </BODY>
+</HTML>
+
Index: trunk/tools/documentation/plotting/plotting.html
===================================================================
--- trunk/tools/documentation/plotting/plotting.html	(revision 1948)
+++ trunk/tools/documentation/plotting/plotting.html	(revision 2128)
@@ -10,4 +10,5 @@
       <A CLASS="lc" HREF="2Dshad_time.html" TARGET="value">2D shadow time</A><BR>
       <A CLASS="lc" HREF="2Dshad_cont.html" TARGET="value">2D shadow contour</A><BR>
+      <A CLASS="lc" HREF="2Dshad_contdisc.html" TARGET="value">2D shadow cont. & disc.</A><BR>
       <A CLASS="lc" HREF="2Dshad_contdisc_time.html" TARGET="value">2D shadow cont. & disc. time</A><BR>
       <A CLASS="lc" HREF="2Dshad_2cont.html" TARGET="value">2D shadow 2 contour</A><BR>
@@ -15,4 +16,5 @@
       <A CLASS="lc" HREF="2Dshad_line.html" TARGET="value">2D shadow line</A><BR>
       <A CLASS="lc" HREF="2Dshad_line_time.html" TARGET="value">2D shadow line time</A><BR>
+      <A CLASS="lc" HREF="2Dshad_map.html" TARGET="value">2D masked shadow on map</A><BR>
       <A CLASS="lc" HREF="draw_bar.html" TARGET="value">bar</A><BR>
       <A CLASS="lc" HREF="bar_line.html" TARGET="value">bar & line</A><BR>
@@ -25,4 +27,5 @@
       <A CLASS="lc" HREF="movie_2Dshad.html" TARGET="value">movie 2D shadow</A><BR>
       <A CLASS="lc" HREF="multi_2Dshad.html" TARGET="value">multi 2D shadow</A><BR>
+      <A CLASS="lc" HREF="multiWindRose.html" TARGET="value">multi Wind Rose</A><BR>
       <A CLASS="lc" HREF="Neighbourgh_evol.html" TARGET="value">Neighbourghood evolution</A><BR>
       <A CLASS="lc" HREF="Nlines.html" TARGET="value">Nlines</A><BR>
@@ -33,4 +36,5 @@
       <A CLASS="lc" HREF="pointsZval.html" TARGET="value">points Zval</A><BR>
       <A CLASS="lc" HREF="river_desc.html" TARGET="value">ORCHIDEE's river description</A><BR>
+      <A CLASS="lc" HREF="stations_map.html" TARGET="value">Map of stations</A><BR>
       <A CLASS="lc" HREF="subbasin.html" TARGET="value">ORCHIDEE's sub-basin</A><BR>
       <A CLASS="lc" HREF="Taylor.html" TARGET="value">Taylor diagram</A><BR>
@@ -44,4 +48,5 @@
       <A CLASS="lc" HREF="vertical_levels.html" TARGET="value">vertical levels</A><BR>
       <A CLASS="lc" HREF="WindRose.html" TARGET="value">Wind Rose</A><BR>
+      <A CLASS="lc" HREF="WRFeta_levels.html" TARGET="value">WRF eta-levels</A><BR>
      <BR>
      ------- ------ ----- ---- --- -- -<BR>
Index: trunk/tools/documentation/plotting/stations_map.html
===================================================================
--- trunk/tools/documentation/plotting/stations_map.html	(revision 2128)
+++ trunk/tools/documentation/plotting/stations_map.html	(revision 2128)
@@ -0,0 +1,55 @@
+<HTML>
+  <HEAD>
+    <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=utf-8"/>
+    <LINK REL="stylesheet" TYPE="text/css" HREF="../PyNCplot.css"/>
+  </HEAD>
+
+  </HEAD>
+  <BODY>
+    <DIV CLASS="valheader">
+      stations_map
+    </DIV>
+    <DIV CLASS="valimg">
+      Plotting a map with locations of stations
+      <IMG WIDTH="100%" SRC="figs/stations_map.png"></IMG>
+    </DIV>
+    <DIV CLASS="valhelp">
+    draw_stations_map(filename, values)<BR>
+    Function to plot a map with the stations<BR>
+    &nbsp;&nbsp;filename= name of ASCII file with the values of the stations as: ('#' comments)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;[sttype],[label],[lon],[lat],[height]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[sttype]: type of the station<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[label]: label of the station (LaTeX format, 'None' for no label)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lon]: longitude of the station<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[lat]: latitude of the station<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[height]: height of the station<BR>
+    &nbsp;&nbsp;values=[stypevals]:[mapvals]:[SWNEbox]:[imgtit]:[imgkind]:[close]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;stypevals='@' separated list of [sttype]|[marker]|[color]|[size]|[fontsize]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[sttype]: label of the type of the marker<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[marker]: type of the marker<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[color]: color of the marker<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[size]: size of the marker<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;mapvals= map characteristics: [proj],[res]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;see full documentation: http://matplotlib.org/basemap/<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[proj]: projection<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'cyl', cilindric<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'lcc', lambert-conformal<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[res]: resolution:<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'c', crude<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'l', low<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'i', intermediate<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'h', high<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;* 'f', full<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;SWNEbox: ',' list of the vertexs of the map [SWlon, NElon, SWlat, NElat]<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'lonlatxtrms': to use the extremes of the lon,lat stations<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;imgtit: title of the image ('!' for spaces)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;imgkind: kind of file output of the image (ps, pns, pdf, ...)<BR>
+    &nbsp;&nbsp;&nbsp;&nbsp;close: whether figure should be closed or not<BR>
+    </DIV>
+    <DIV CLASS="valins">
+      $ python ${pyHOME}/drawing.py -o draw_stations_map -f stations.inf -S 'snd|x|#FFAAAA|2|8@sfc|<BR>
+      &nbsp;&nbsp;x|#AAFFAA|2|8:cyl,f:lonlatxtrms:FPS!Alps!stations:png:true'
+    </DIV>
+  </BODY>
+</HTML>
+
