Changeset 4063 for trunk

Timestamp:

Feb 10, 2026, 9:47:13 PM (3 days ago)

Author:

tbertrand

Message:

MARS PCM:
dust coagulation : just renaming some flags so that it is more human readable, following suggestion by Grand Manitou Ehouarn
TB

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callkeys_mod.F90 (modified) (1 diff)
• coagulation_mod.F90 (deleted)
• conf_phys.F (modified) (3 diffs)
• dust_coagulation_mod.F90 (added)
• physiq_mod.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/phymars/callkeys_mod.F90

@@ -64 +64 @@
 logical,save :: topflows ! entrainment by mountain top dust flows parametrization
 !$OMP THREADPRIVATE(topflows)
-logical,save :: coagulation ! coagulation of dust particles
-logical,save :: fullcoag ! coagulation of dust particles using full equations instead of lookup tables
-logical,save :: kernel_b, kernel_g, kernel_de, kernel_ti ! kernels for different types of dust coagulation
-!$OMP THREADPRIVATE(coagulation,fullcoag,kernel_b,kernel_g,kernel_de,kernel_ti)
+logical,save :: dust_coagulation ! coagulation of dust particles
+!$OMP THREADPRIVATE(coagulation)
 integer,save :: coal_kg ! mode for coalescence 0,1, or 2 
 !$OMP THREADPRIVATE(coal_kg)

trunk/LMDZ.MARS/libf/phymars/conf_phys.F

@@ -46 +46 @@
      &                     ads_massive_ice
       use nonoro_gwd_mix_mod, only: calljliu_gwimix
+      use dust_coagulation_mod, only: dust_coag_kernel_b,
+     &      dust_coag_kernel_g,dust_coag_kernel_de,dust_coag_kernel_ti,
+     &      coal_coeff_mode,full_coag_equations
       use lwxn_mod, only: lwxn_linear, lwxn_alphan, lwxn_ncouche
       use callkeys_mod, only: startphy_file, activice, activeco2ice,
@@ -72 +75 @@
      &                        tf_injection, ti_injection, thermochem,
      &                        thermoswater, tituscap, tke_heat_flux,
-     &                        topflows, water, coagulation, kernel_b,
-     &                        kernel_g,kernel_de,kernel_ti,coal_kg,
-     &                        fullcoag
+     &                        topflows, water, dust_coagulation
       use write_output_mod, only: output_diagfi
 
@@ -389 +390 @@
 ! dust particle coagulation
          write(*,*)"call coagulation of dust"
-         coagulation=.false. ! default value
-         call getin_p("coagulation",coagulation)
-         write(*,*)" coagulation = ",coagulation
-         fullcoag=.false. ! default value 
-         kernel_b=.true. ! default value (if coagulation=true)
-         kernel_g=.false. ! default value
-         kernel_de=.false. ! default value
-         kernel_ti=.false. ! default value
-         coal_kg=0 ! default value
-         if (coagulation) then
-           write(*,*)" coagulation fullcoag and mode= ",fullcoag,coal_kg
-           write(*,*)" coagulation kernels= ",kernel_b,kernel_g,
-     &                                        kernel_de,kernel_ti
+         dust_coagulation=.false. ! default value
+         call getin_p("dust_coagulation",dust_coagulation)
+         write(*,*)" dust_coagulation = ",dust_coagulation
+         full_coag_equations=.false. ! default value 
+         dust_coag_kernel_b=.true. ! default value (if coagulation=true)
+         dust_coag_kernel_g=.false. ! default value
+         dust_coag_kernel_de=.false. ! default value
+         dust_coag_kernel_ti=.false. ! default value
+         coal_coeff_mode=0 ! default value
+         if (dust_coagulation) then
+           write(*,*)" coagulation full_coag_equations: ",
+     &                                            full_coag_equations
+           write(*,*)" coagulation coal_coeff_mode: ",coal_coeff_mode
+           write(*,*)" coagulation kernels= ",dust_coag_kernel_b,
+     &     dust_coag_kernel_g,dust_coag_kernel_de,dust_coag_kernel_ti
          endif
 

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -25 +25 @@
       use calcstormfract_mod, only: calcstormfract
       use topmons_mod, only: topmons,topmons_setup
-      use coagulation_mod, only: coagul_main,coagul_init
+      use dust_coagulation_mod, only: dust_coagulation_main,
+     &                                dust_coagulation_init
       use nltecool_mod, only: nltecool
       use nlte_tcool_mod, only: nlte_tcool
@@ -150 +151 @@
       use callkeys_mod, only: photochem, callthermos
       use callkeys_mod, only: startphy_file
-      use callkeys_mod, only: coagulation
+      use callkeys_mod, only: dust_coagulation
 
       IMPLICIT NONE
@@ -794 +795 @@
         if (topflows) call topmons_setup(ngrid)
 
-c        Initialize coagulation parameters
+c        Initialize dust coagulation parameters
 c        ~~~~~~~~~~~~~~~
-        if (coagulation) call coagul_init()
+        if (dust_coagulation) call dust_coagulation_init()
 
 c        Parameterization of the ATKE routine
@@ -1455 +1456 @@
 c     3.3 Dust coagulation
 c    -------------------------------------------
-      IF (coagulation) THEN
+      IF (dust_coagulation) THEN
          pdqcoag(:,:,:)=0.
-         CALL coagul_main(ngrid,nlayer,nq,ptime,ptimestep,
+         CALL dust_coagulation_main(ngrid,nlayer,nq,ptime,ptimestep,
      &                     pq,pdq,pt,pdt,pplay,pplev,
      &                     pdqcoag)
@@ -1496 +1497 @@
          ENDIF ! end of if (rdstorm)
 
-      ENDIF ! end of if (coagulation)
+      ENDIF ! end of if (dust coagulation)
 
 c     3.4 Dust injection from the surface
@@ -3791 +3792 @@
            endif ! (topflows)
 
-           if (coagulation) then
+           if (dust_coagulation) then
              call write_output('zdqcoag_dustm',
-     &                     'coagulation tendency',
+     &                     'dust coagulation tendency',
      &                     'kg kg-1 s-1',pdqcoag(:,:,igcm_dust_mass))
              call write_output('zdqcoag_dustn',
-     &                     'coagulation tendency',
+     &                     'dust coagulation tendency',
      &                     'nbp kg-1 s-1',pdqcoag(:,:,igcm_dust_number))
            endif