Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/2dplot.py
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/2dplot.py	(revision 296)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/2dplot.py	(revision 296)
@@ -0,0 +1,51 @@
+#! /usr/bin/env python
+
+from netCDF4 import Dataset
+import matplotlib.pyplot as mpl
+import numpy as np
+from myplot import reducefield,getfield,getcoorddef,calculate_bounds,bounds,fmtvar,ptitle,makeplotres
+from matplotlib.pyplot import contourf,colorbar,show,xlabel,ylabel
+from matplotlib.cm import get_cmap
+
+name = "wrfout_d01_9999-09-09_09:00:00_z"
+itime = 12
+#itime = 1
+ndiv = 10
+zey = 0
+var = "W"
+vmin = -1.
+vmax = 1.
+title = "Vertical velocity"
+var = "Um"
+vmin = -2.
+vmax = 18.
+title = "Horizontal velocity"
+#var = "tk"
+#vmin = 150.
+#vmax = 170.
+#title = "Atmospheric temperature"
+
+nc = Dataset(name)
+
+what_I_plot, error = reducefield( getfield(nc,var), d4=itime, d2=zey )
+
+y = nc.variables["vert"][:]
+
+horinp = len(what_I_plot[0,:])
+x = np.linspace(0.,horinp*500.,horinp) / 1000.
+
+zevmin, zevmax = calculate_bounds(what_I_plot,vmin=vmin,vmax=vmax)
+#if colorb in ["def","nobar"]:   palette = get_cmap(name=defcolorb(fvar))
+#else:                           palette = get_cmap(name=colorb)
+palette = get_cmap(name="jet")
+what_I_plot = bounds(what_I_plot,zevmin,zevmax)
+zelevels = np.linspace(zevmin,zevmax)
+contourf( x, y, what_I_plot, zelevels, cmap = palette )
+colorbar(fraction=0.05,pad=0.03,format=fmtvar(var),\
+                       ticks=np.linspace(zevmin,zevmax,ndiv+1),\
+                       extend='neither',spacing='proportional')
+ptitle(title)
+xlabel("Horizontal coordinate (km)")
+ylabel("Altitude (m)")
+makeplotres(var+str(itime),res=200.,disp=False)
+#show()
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/myplot.py
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/myplot.py	(revision 295)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/myplot.py	(revision 296)
@@ -490,4 +490,5 @@
              "WMAX_TH":      "%.1f",\
              "QSURFICE":     "%.0f",\
+             "Um":           "%.0f",\
              "ALBBARE":      "%.2f",\
                     }
@@ -566,4 +567,6 @@
                 "nice":        "http://users.info.unicaen.fr/~karczma/TEACH/InfoGeo/Images/Planets/earthmap1k.jpg",\
                 "blue":        "http://eoimages.gsfc.nasa.gov/ve/2430/land_ocean_ice_2048.jpg",\
+                "blueclouds":  "http://eoimages.gsfc.nasa.gov/ve/2431/land_ocean_ice_cloud_2048.jpg",\
+                "justclouds":  "http://eoimages.gsfc.nasa.gov/ve/2432/cloud_combined_2048.jpg",\
 			}
         ### see http://www.mmedia.is/~bjj/planetary_maps.html
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM	(revision 295)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM	(revision 296)
@@ -107,4 +107,9 @@
 state  real  TAU_ICE   ij   misc  1  -  rd   "TAU_ICE"   "CLOUD OD at 825 cm-1 TES"        ""        #SAVEMARS2 tauTES
 state  real  PDTZ      ikj  misc  1  -  rd   "PDT"       "TEMP TENDENCY"                   "K s-1"   #SAVEMARS3 pdt
+state  real  LMAX_TH   ij   misc  1  -  rd   "LMAX_TH"   "MAXIMUM LEVEL REACHED IN TH"     ""        #SAVEMARS2 lmax_th_out
+state  real  HFMAX_TH  ij   misc  1  -  rd   "HFMAX_TH"  "MAXIMUM TH HEAT FLUX"            "m.K/s"   #SAVEMARS2 hfmax_th
+state  real  WMAX_TH   ij   misc  1  -  rd   "WMAX_TH"   "MAXIMUM TH VERTICAL VELOCITY"    "m/s"     #SAVEMARS2 wmax_th
+
+
 ####
 ####