Index: trunk/LMDZ.MARS/libf/aeronomars/calchim.F
===================================================================
--- trunk/LMDZ.MARS/libf/aeronomars/calchim.F	(revision 473)
+++ trunk/LMDZ.MARS/libf/aeronomars/calchim.F	(revision 476)
@@ -493,6 +493,8 @@
             call writediagfi(ngridmx,'jo3','j o3->o1d',
      $                       's-1',3,jo3_3d(1,1))
+           if (callstats) then
             call wstats(ngridmx,'jo3','j o3->o1d',
      $                       's-1',3,jo3_3d(1,1))
+           endif
          end if ! of if (ngridmx.gt.1)
       endif ! of if (output)
Index: trunk/LMDZ.MARS/libf/aeronomars/surfacearea.F
===================================================================
--- trunk/LMDZ.MARS/libf/aeronomars/surfacearea.F	(revision 473)
+++ trunk/LMDZ.MARS/libf/aeronomars/surfacearea.F	(revision 476)
@@ -101,11 +101,15 @@
 
 ! write diagnostics in micron2/cm3
-
-      call wstats(ngrid,"surfdust", "Dust surface area",
+      
+      if (callstats) then
+        call wstats(ngrid,"surfdust", "Dust surface area",
      $            "micron2 cm-3",3,surfdust*1.e6)
+      endif
       call writediagfi(ngrid,"surfdust", "Dust cloud surface area",
      $            "micron2 cm-3",3,surfdust*1.e6)
-      call wstats(ngrid,"surfice", "Ice cloud surface area",
+      if (callstats) then
+        call wstats(ngrid,"surfice", "Ice cloud surface area",
      $            "micron2 cm-3",3,surfice*1.e6)
+      endif
       call writediagfi(ngrid,"surfice", "Ice cloud surface area",
      $            "micron2 cm-3",3,surfice*1.e6)