Changeset 4054 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Feb 6, 2026, 5:21:02 PM (10 days ago)

Author:

tbertrand

Message:

Mars PCM:
Adding coagulation of dust particles (cf Bertrand et al., 2022, Icarus)
TB

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

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

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

@@ -73 +73 @@
 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)
+integer,save :: coal_kg ! mode for coalescence 0,1, or 2 
+!$OMP THREADPRIVATE(coal_kg)
 logical,save :: latentheat_surfwater ! latent heat release from ground water ice sublimation/condensation
 !$OMP THREADPRIVATE(latentheat_surfwater)

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

@@ -70 +70 @@
      &                        tf_injection, ti_injection, thermochem,
      &                        thermoswater, tituscap, tke_heat_flux,
-     &                        topflows, water
+     &                        topflows, water, coagulation, kernel_b,
+     &                        kernel_g,kernel_de,kernel_ti,coal_kg,
+     &                        fullcoag
 
       IMPLICIT NONE
@@ -379 +381 @@
          call getin_p("topflows",topflows)
          write(*,*)" topflows = ",topflows
+
+! 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
+         endif
 
 ! latent heat release from ground water ice sublimation/condensation

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 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
 
       IMPLICIT NONE
@@ -544 +546 @@
 c entrainment by mountain top dust flows above sub-grid scale topography
       REAL pdqtop(ngrid,nlayer,nq) ! tendency for dust after topmons
+c dust coagulation tendency
+      REAL pdqcoag(ngrid,nlayer,nq) ! tendency for dust coagulation
 
 c when no startfi file is asked for init
@@ -791 +795 @@
 c        ~~~~~~~~~~~~~~~
         if (topflows) call topmons_setup(ngrid)
+
+c        Initialize coagulation parameters
+c        ~~~~~~~~~~~~~~~
+        if (coagulation) call coagul_init()
 
 c        Parameterization of the ATKE routine
@@ -1447 +1455 @@
       ENDIF ! end of if (topflows)
 
-c     3.3 Dust injection from the surface
+c     3.3 Dust coagulation
+c    -------------------------------------------
+      IF (coagulation) THEN
+         pdqcoag(:,:,:)=0.
+         CALL coagul_main(ngrid,nlayer,nq,ptime,ptimestep,
+     &                     pq,pdq,pt,pdt,pplay,pplev,
+     &                     pdqcoag)
+         !Update tendency :
+         DO l=1,nlayer
+          DO ig=1,ngrid
+           pdq(ig,l,igcm_dust_mass)=
+     & pdq(ig,l,igcm_dust_mass)+ pdqcoag(ig,l,igcm_dust_mass)
+           pdq(ig,l,igcm_dust_number)=
+     & pdq(ig,l,igcm_dust_number)+pdqcoag(ig,l,igcm_dust_number)
+          ENDDO
+         ENDDO
+
+         IF (topflows) THEN
+           DO l=1,nlayer
+            DO ig=1,ngrid
+              pdq(ig,l,igcm_topdust_mass)=
+     & pdq(ig,l,igcm_topdust_mass)+
+     &                      pdqcoag(ig,l,igcm_topdust_mass)
+              pdq(ig,l,igcm_topdust_number)=
+     & pdq(ig,l,igcm_topdust_number)+
+     &                      pdqcoag(ig,l,igcm_topdust_number)
+            ENDDO
+           ENDDO
+         ENDIF ! end of if (topflows)
+
+         IF (rdstorm) THEN
+           DO l=1,nlayer
+            DO ig=1,ngrid
+           pdq(ig,l,igcm_stormdust_mass)=
+     & pdq(ig,l,igcm_stormdust_mass)+
+     &                      pdqcoag(ig,l,igcm_stormdust_mass)
+           pdq(ig,l,igcm_stormdust_number)=
+     & pdq(ig,l,igcm_stormdust_number)+
+     &                      pdqcoag(ig,l,igcm_stormdust_number)
+            ENDDO
+           ENDDO
+         ENDIF ! end of if (rdstorm)
+
+      ENDIF ! end of if (coagulation)
+
+c     3.4 Dust injection from the surface
 c    -------------------------------------------
            if (dustinjection.gt.0) then
@@ -3744 +3797 @@
            endif ! (topflows)
 
+           if (coagulation) then
+             call write_output('zdqcoag_dustm',
+     &                     'coagulation tendency',
+     &                     'kg kg-1 s-1',pdqcoag(:,:,igcm_dust_mass))
+             call write_output('zdqcoag_dustn',
+     &                     'coagulation tendency',
+     &                     'nbp kg-1 s-1',pdqcoag(:,:,igcm_dust_number))
+           endif  
+
            if (dustscaling_mode==2) then
              call write_output("dust_rad_adjust",