Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phys_output_mod.F90
===================================================================
--- LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phys_output_mod.F90	(revision 1472)
+++ LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phys_output_mod.F90	(revision 1473)
@@ -427,4 +427,6 @@
   type(ctrl_out),save :: o_pres         = ctrl_out((/ 2, 3, 10, 10, 1 /),'pres')
   type(ctrl_out),save :: o_paprs        = ctrl_out((/ 2, 3, 10, 10, 1 /),'paprs')
+  type(ctrl_out),save :: o_mass        = ctrl_out((/ 2, 3, 10, 10, 1 /),'mass')
+
   type(ctrl_out),save :: o_rneb         = ctrl_out((/ 2, 5, 10, 10, 1 /),'rneb')
   type(ctrl_out),save :: o_rnebcon      = ctrl_out((/ 2, 5, 10, 10, 1 /),'rnebcon')
@@ -1108,4 +1110,5 @@
  CALL histdef3d(iff,o_pres%flag,o_pres%name, "Air pressure", "Pa" )
  CALL histdef3d(iff,o_paprs%flag,o_paprs%name, "Air pressure Inter-Couches", "Pa" )
+ CALL histdef3d(iff,o_mass%flag,o_mass%name, "Masse Couches", "kg/m2" )
  CALL histdef3d(iff,o_rneb%flag,o_rneb%name, "Cloud fraction", "-")
  CALL histdef3d(iff,o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-")
Index: LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phys_output_write.h
===================================================================
--- LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phys_output_write.h	(revision 1472)
+++ LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phys_output_write.h	(revision 1473)
@@ -104,4 +104,9 @@
      s                   o_psol%name,itau_w,zx_tmp_fi2d)
        ENDIF
+
+       IF (o_mass%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),o_mass%name,itau_w,zmasse)
+        ENDIF
+
 
        IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN