Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 4140)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 4141)
@@ -3138,10 +3138,39 @@
                     call wstats(ngrid,"vmr_"//trim(noms(iq)),
      $                        "Volume mixing ratio","mol/mol",3,vmr)
-                    if ((noms(iq).eq."o")
-     $             .or. (noms(iq).eq."co2")
+                    if ((noms(iq).eq."co2")
+     $             .or. (noms(iq).eq."co")
+     $             .or. (noms(iq).eq."o")
+     $             .or. (noms(iq).eq."o1d")
+     $             .or. (noms(iq).eq."o2")
      $             .or. (noms(iq).eq."o3")
+     $             .or. (noms(iq).eq."h")
+     $             .or. (noms(iq).eq."h2")
+     $             .or. (noms(iq).eq."oh")
+     $             .or. (noms(iq).eq."ho2")
+     $             .or. (noms(iq).eq."h2o2")
+     $             .or. (noms(iq).eq."n2")
      $             .or. (noms(iq).eq."ar")
-     $             .or. (noms(iq).eq."o2")
-     $             .or. (noms(iq).eq."h2o_vap") ) then
+     $             .or. (noms(iq).eq."h2o_ice")
+     $             .or. (noms(iq).eq."h2o_vap")
+     $             .or. (noms(iq).eq."n")
+     $             .or. (noms(iq).eq."no")
+     $             .or. (noms(iq).eq."no2")
+     $             .or. (noms(iq).eq."n2d")
+     $             .or. (noms(iq).eq."he")
+     $             .or. (noms(iq).eq."co2plus")
+     $             .or. (noms(iq).eq."oplus")
+     $             .or. (noms(iq).eq."o2plus")
+     $             .or. (noms(iq).eq."coplus")
+     $             .or. (noms(iq).eq."cplus")
+     $             .or. (noms(iq).eq."nplus")
+     $             .or. (noms(iq).eq."noplus")
+     $             .or. (noms(iq).eq."n2plus")
+     $             .or. (noms(iq).eq."hplus")
+     $             .or. (noms(iq).eq."hco2plus")
+     $             .or. (noms(iq).eq."hcoplus")
+     $             .or. (noms(iq).eq."h2oplus")
+     $             .or. (noms(iq).eq."h3oplus")
+     $             .or. (noms(iq).eq."ohplus")
+     $             .or. (noms(iq).eq."elec") ) then
                       call write_output("vmr_"//trim(noms(iq)),
      $                         "Volume mixing ratio","mol/mol",vmr(:,:))