Index: trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3585)
+++ trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3590)
@@ -1515,5 +1515,5 @@
    
             pdq(:,:,:) = pdq(:,:,:) + pdqmufi(:,:,:)
-         
+            
          ELSE
             IF (haze) THEN
@@ -2329,5 +2329,22 @@
            endif
 
-         endif
+         endif ! end haze
+
+         if (callmufi) then
+            ! Tracers:
+            call write_output("m0as","Spherical mode 0th order moment","m-3",zq(:,:,micro_indx(1))*int2ext(:,:))
+            call write_output("m3as","Spherical mode 3rd order moment","m3.m-3",zq(:,:,micro_indx(2))*int2ext(:,:))
+            call write_output("m0af","Fractal mode 0th order moment","m-3",zq(:,:,micro_indx(3))*int2ext(:,:))
+            call write_output("m3af","Fractal mode 3rd order moment","m3.m-3",zq(:,:,micro_indx(4))*int2ext(:,:))
+
+            ! Diagnostics:
+            call write_output("rcs","Spherical mode characteristic radius","m",mp2m_rc_sph(:,:))
+            call write_output("rcf","Fractal mode characteristic radius","m",mp2m_rc_fra(:,:))
+            
+           if (optichaze) then
+              call write_output("tau_col",&
+                 "Total aerosol optical depth","opacity",tau_col)
+           endif ! end optichaze
+         endif ! end callmufi
 
       endif ! end of 'tracer'