Changeset 4047 for trunk/LMDZ.TITAN/libf/phytitan/cond_muphy.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/cond_muphy.F90 (modified) (2 diffs)

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]