Changeset 4047 for trunk

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

Location:

trunk/LMDZ.TITAN

Files:

• deftank/callphys.def.allmuparams (modified) (5 diffs)
• libf/phytitan/callkeys_mod.F90 (modified) (1 diff)
• libf/phytitan/cond_muphy.F90 (modified) (2 diffs)
• libf/phytitan/inifis_mod.F90 (modified) (1 diff)
• libf/phytitan/physiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.TITAN/deftank/callphys.def.allmuparams

@@ -104 +104 @@
 # If yes, number of ices ? (must be compatible with traceur.def AND microphysical model)
 nices      = 1
+# Apply condensation heating rate?
+latent_heat = .true.
 # Use new optics for clouds ?
 opt4clouds = .true.
@@ -137 +139 @@
 rho_aer                = 600.
 # Enable/disable Haze production process 
-haze_production        = T
+haze_production        = .true.
 # Enable/disable Haze coagulation process 
-haze_coagulation       = T
+haze_coagulation       = .true.
 # Coagulation interactions, a combination of:
 #    0 - no interactions (same as haze_coagulation == F)
@@ -148 +150 @@
 haze_coag_interactions = 7
 # Enable/disable Haze sedimentation process 
-haze_sedimentation     = T
+haze_sedimentation     = .true.
 # Disable Fiadero correction for sedimentation process
-no_fiadero             = T
+no_fiadero             = .true.
 # Fiadero correction minimum ratio threshold
 fiadero_min_ratio      = 0.1
@@ -156 +158 @@
 fiadero_max_ratio      = 10.
 # Force settling velocity to M0
-wsed_m0                = T
+wsed_m0                = .true.
 # Force settling velocity to M3
-wsed_m3                = F
+wsed_m3                = .false.
 # Enable/disable clouds sedimentation process
 # (automatically set to F if clouds microphysics is not enabled)
-clouds_sedimentation   = T
+clouds_sedimentation   = .true.
 # Enable/disable clouds nucleation and condensation processes
 # (automatically set to F if clouds microphysics is not enabled)
-clouds_nuc_cond        = T
+clouds_nuc_cond        = .true.
 # Condensible species configuration file
 # (not needed if clouds microphysics is not enabled)
 specie_cfg             = mp2m_species.cfg
 # Enable/disable spherical mode transfert probability
-transfert_probability = T
+transfert_probability = .true.
 # Spherical mode transfert probability look-up tables file
 # (optional if 'transfert_probability' is False) 
@@ -175 +177 @@
 # Electric charging coagulation correction
 # If set to .false. then no correction is assumed
-electric_charging     = T 
+electric_charging     = .true.
 # File for the electric charging correction factor.
 # (optional if 'electric_charging' is False) 

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

@@ -21 +21 @@
       logical,save :: callchim, callmufi, callclouds
 !$OMP THREADPRIVATE(callchim,callmufi,callclouds)
+      logical,save :: latent_heat
+!$OMP THREADPRIVATE(latent_heat)
       logical,save :: global1d
 !$OMP THREADPRIVATE(global1d)

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

@@ -69 +69 @@
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 do iq = 1, size(ices_indx)
-    dtlc(:,:) = dtlc(:,:) + (dqmuficond(:,:,iq) * Lc(:,:,iq))
+    dtlc(:,:) = dtlc(:,:) - (dqmuficond(:,:,iq) * Lc(:,:,iq))
 enddo
 
@@ -75 +75 @@
 ! Condensation heating rate :
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! If ice formation  : dqmuficond > 0 --> dtlc > 0
-! Else vaporisation : dqmuficond < 0 --> dtlc < 0
+! If ice formation  : dqmuficond < 0 --> dtlc > 0
+! Else vaporisation : dqmuficond > 0 --> dtlc < 0
 dtlc(:,:) = dtlc(:,:) / cpp ! [K.s-1]
 

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

@@ -490 +490 @@
      endif
 
+     write(*,*) "Apply condensation heating rate?"
+     latent_heat=.false. ! default value
+     call getin_p("latent_heat",latent_heat)
+     write(*,*)" latent_heat = ",latent_heat
+
      write(*,*) "Disable the coupling of microphysics within rad. transf. ?"
      write(*,*) "If disabled we will assume a planetwide vert. profile of extinction ..."

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