Changeset 3949 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

Nov 4, 2025, 8:51:03 AM (7 weeks ago)

Author:

debatzbr

Message:

Pluto PCM: Add condensable gas tracers through muphi
BBT

File:

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

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

@@ -407 +407 @@
       ! ---------------------------------
       real gzlat(ngrid,nlayer)           ! Altitude-Latitude-dependent gravity (this should be stored elsewhere...).
+      real pdqmuchem(ngrid,nlayer,nq)    ! Condensable gases through muphi tendency (kg/kg_of_air/s).
       real pdqmufi(ngrid,nlayer,nq)      ! Microphysical tendency (X/kg_of_air/s).
       real pdqmufi_prod(ngrid,nlayer,nq) ! Aerosols production tendency (kg/kg_of_air/s).
@@ -1519 +1520 @@
       if (tracer) then
 
-!      ---------------------------------------
-!      Methane ice condensation in the atmosphere
-!      ----------------------------------------
+         !-------------------------------------------------
+         ! Methane (CH4) ice condensation in the atmosphere
+         !-------------------------------------------------
          rice_ch4(:,:)=0 ! initialization needed for callsedim
          zdqch4cloud(:,:,:)=0.
@@ -1553 +1554 @@
          end if
 
-!      ---------------------------------------
-!      CO ice condensation in the atmosphere
-!      ----------------------------------------
+         !--------------------------------------
+         ! CO ice condensation in the atmosphere
+         !--------------------------------------
          zdqcocloud(:,:,:)=0.
          IF ((carbox).and.(monoxcloud).and.(.not.fast)) THEN
@@ -1587 +1588 @@
          rice_co(:,:)=0 ! initialization needed for callsedim
          END IF  ! of IF (carbox)
+
+         !--------------------------------------------------
+         ! Condensable gases through muphi in the atmosphere
+         !--------------------------------------------------
+         pdqmuchem(:,:,:) = 0.
+         IF (callmufi) THEN
+            call mugas_prof(ngrid,nlayer,nq,zzlay,zzlev,pplay,pt,pdqmuchem)
+         ENDIF
+         pdq(:,:,:) = pdq(:,:,:) + pdqmuchem(:,:,:)
 
          ! ----------------------------------------