Changeset 3882 for trunk/LMDZ.PLUTO/libf

Timestamp:

Aug 7, 2025, 10:12:58 AM (4 months ago)

Author:

debatzbr

Message:

Pluto PCM: Clean outputs for photolysis and mufi
BBT

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (6 diffs)

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

@@ -465 +465 @@
             endwhere
             nueffrad(:,:,1) = exp(sig**2) - 1
-            !if (exp(sig**2) - 1 > 0.1) then
-            !   nueffrad(:,:,1) = exp(sig**2) - 1
-            !else
-            !   nueffrad(:,:,1) = 0.1
-            !endif
             ! Fractal aerosols
             sig = 0.35
@@ -478 +473 @@
             endwhere
             nueffrad(:,:,2) = exp(sig**2) - 1
-            !if (exp(sig**2) - 1 > 0.1) then
-            !   nueffrad(:,:,2) = exp(sig**2) - 1
-            !else
-            !   nueffrad(:,:,2) = 0.1
-            !endif
 
          else
@@ -669 +659 @@
                 do nw=1,L_NSPECTI
                     if(qsiaer(k,nw,iaer).gt.1.05*qxiaer(k,nw,iaer))then
-                    message='Serious problems with qsvaer values'
+                    message='Serious problems with qsiaer values'
                     call abort_physic(subname,message,1)
                     endif

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

@@ -1968 +1968 @@
       zrho_haze(:,:)=0.
       zdqrho_photprec(:,:)=0.
-      IF (haze.and.optichaze) then
+      IF (haze.or.callmufi) then
          DO ig=1,ngrid
             DO l=1,nlayer
@@ -2251 +2251 @@
       endif
 
+      ! Optical diagnostics.
+      call write_output("tau_col",&
+            "Aerosol surface opacity at reference visible wavelength","",tau_col)
       ! 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))
+      call write_output('dtauv_01','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,1))  ! 5.398 um (17x27)
+      call write_output('dtauv_23','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,23)) ! 0.941 um (17x27)
+      call write_output('dtauv_24','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,24)) ! 0.700 um (17x27)
+      call write_output('dtauv_27','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,27)) ! 0.119 um (17x27)
+      call write_output('dtaui_01','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,1))  ! 550.0 um (17x27)
+      call write_output('dtaui_17','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,17)) ! 3.531 um (17x27)
 
       if (calllott) then
@@ -2345 +2350 @@
          endif
 
-         ! Haze
-         if (haze) then
+         ! Haze or callmufi
+         if (haze.or.callmufi) then
 
             if (haze_radproffix)then
@@ -2356 +2361 @@
 
             ! 3D Tendencies
-            call write_output("zdqrho_photprec","zdqrho_photprec",&
+            call write_output("zdqrho_photprec","Photolysis rate",&
                         "kg.m-3.s-1",zdqrho_photprec(:,:))
-            call write_output("zdqphot_prec","zdqphot_prec","",&
+            
+            !call write_output("zdqhaze_ch4","Photolysis of CH4 tendancy","kg.kg-1.s-1",&
+            !         zdqhaze(:,:,igcm_ch4_gas))
+            !call write_output("zdqhaze_prec","Precursor haze tendancy","kg.kg-1.s-1",&
+            !         zdqhaze(:,:,igcm_prec_haze))
+
+            call write_output("zdqphot_ch4","Photolysis of CH4 tendancy","kg.kg-1.s-1",&
+                                                zdqphot_ch4(:,:))
+            call write_output("zdqphot_prec","Precursor haze tendancy","kg.kg-1.s-1",&
                                                 zdqphot_prec(:,:))
-            call write_output("zdqhaze_ch4","zdqhaze_ch4","",&
-                     zdqhaze(:,:,igcm_ch4_gas))
-            call write_output("zdqhaze_prec","zdqhaze_prec","",&
-                     zdqhaze(:,:,igcm_prec_haze))
-            call write_output("zdqphot_ch4","zdqphot_ch4","",&
-                                                zdqphot_ch4(:,:))
-            call write_output("zdqconv_prec","zdqconv_prec","",&
+            call write_output("zdqconv_prec","Conversion of prec_haze into haze tendancy","kg.kg-1.s-1",&
                                                 zdqconv_prec(:,:))
 
@@ -2377 +2384 @@
                endif
             endif
-
-            if (optichaze) then
-               call write_output("tau_col",&
-               "Total aerosol optical depth","opacity",tau_col)
-            endif
-
-         endif ! end haze
+         endif ! end Haze or callmufi
 
          if (callmufi) then
@@ -2393 +2394 @@
 
             ! Diagnostics:
+            call write_output("aers_prec","Spherical aerosols precipitations","kg.m-2.s-1",mp2m_aer_s_prec(:))
+            call write_output("aerf_prec","Fractal aerosols precipitations","kg.m-2.s-1",mp2m_aer_f_prec(:))
+            call write_output("aers_w","Spherical aerosol settling velocity","m.s-1",mp2m_aer_s_w(:,:))
+            call write_output("aerf_w","Fractal aerosol settling velocity","m.s-1",mp2m_aer_f_w(:,:))
             call write_output("rcs","Characteristic radius of spherical aerosols","m",mp2m_rc_sph(:,:))
             call write_output("rcf","Characteristic radius of fractal aerosols","m",mp2m_rc_fra(:,:))
-
-            if (optichaze) then
-               call write_output("tau_col",&
-               "Aerosol surface opacity at reference visible wavelength","",tau_col)
-            endif ! end optichaze
          endif ! end callmufi