Changeset 3806

Timestamp:

Jun 16, 2025, 10:11:09 AM (10 hours ago)

Author:

debatzbr

Message:

Pluto PCM: Clean optical outputs
BBT

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -739 +739 @@
       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 +2232 @@
       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 +2247 @@
       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 +2388 @@
             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