Changeset 4047 for trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

Timestamp:

Feb 5, 2026, 5:50:37 PM (3 weeks ago)

Author:

aslmd

Message:

Titan PCM:

• add possibility to release latent heat to the atmosphere due to gases condensation.
• correct sign error on temperature tendency due to latent heat release in phytitan/cond_muphy.F90
• add latent heat release flag in callphys.def ("latent_heat" : .true. or .false.).

EMo

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -196 +196 @@
       real,intent(in) :: pv(ngrid,nlayer)      ! Meridional wind component (m/s).
       real,intent(in) :: pt(ngrid,nlayer)      ! Temperature (K).
-      real,intent(in) :: pq(ngrid,nlayer,nq)   ! Tracers (kg/kg_of_air).
+      real,intent(in) :: pq(ngrid,nlayer,nq)   ! Tracers (X/kg_of_air). (X is kg or nothing (number of particles))
       real,intent(in) :: flxw(ngrid,nlayer)    ! Vertical mass flux (ks/s) at lower boundary of layer
 
@@ -431 +431 @@
 
         ! Allocate saved arrays (except for 1D model, where this has already been done)
+        ! For mesoscale it is done in the interface
 #ifndef MESOSCALE
         if (ngrid>1) call phys_state_var_init(nq)
@@ -1177 +1178 @@
                enddo
                call cond_muphy(ngrid,nlayer,pt,dmuficond,zdtlc)
-               !pdt(:,:) = pdt(:,:) + zdtlc(:,:)
+               if (latent_heat) then
+                  pdt(:,:) = pdt(:,:) + zdtlc(:,:)
+               endif
             endif
          endif ! callmufi