Index: trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3805)
+++ trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3806)
@@ -739,5 +739,7 @@
       if (callmufi.or.haze) then
          call lymalpha(zls,zfluxuv)
-         print*, 'Haze lyman-alpha zls,zfluxuv=',zls,zfluxuv
+         if (is_master) then
+            print*, 'Haze lyman-alpha zls,zfluxuv=',zls,zfluxuv
+         endif
       end if
 
@@ -2230,19 +2232,4 @@
       endif ! end of 'enertest'
 
-      ! Diagnostics of optical thickness
-      ! Warning this is exp(-tau), I let you postproc with -log to have tau itself
-      !do nw=1,L_NSPECTV
-      !   write(str2,'(i2.2)') nw
-      !   call write_output('dtauv'//str2,'Layer optical thickness attenuation in VI band '//str2,'',int_dtauv(:,nlayer:1:-1,nw))
-      !enddo
-      !do nw=1,L_NSPECTI
-      !   write(str2,'(i2.2)') nw
-      !   call write_output('dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',int_dtaui(:,nlayer:1:-1,nw))
-      !enddo
-      call write_output('dtauv_01','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,1))
-      call write_output('dtauv_23','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,23))
-      call write_output('dtaui_01','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,1))
-      call write_output('dtaui_17','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,17))
-
       ! Temporary inclusions for heating diagnostics.
       if (.not.fast) then
@@ -2260,4 +2247,11 @@
       endif
 
+      ! Diagnostics of optical thickness (dtau = dtau_gas + dtau_rayaer + dtau_cont).
+      ! Warning this is exp(-tau), I let you postproc with -log to have tau itself
+      call write_output('dtauv_01','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,1))
+      call write_output('dtauv_23','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,23))
+      call write_output('dtaui_01','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,1))
+      call write_output('dtaui_17','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,17))
+
       ! Output tracers.
       if (tracer) then
@@ -2394,5 +2388,5 @@
             if (optichaze) then
                call write_output("tau_col",&
-               "Total aerosol optical depth","opacity",tau_col)
+               "Aerosol surface opacity at reference visible wavelength","",tau_col)
             endif ! end optichaze
          endif ! end callmufi