Changeset 3962 for trunk/LMDZ.PLUTO/libf/muphypluto

Timestamp:

Nov 17, 2025, 11:02:59 AM (5 weeks ago)

Author:

debatzbr

Message:

Pluto PCM: Calculates the latent heat released by condensation from the microphysical module.
BBT

File:

• trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_microphysics.F90 (modified) (10 diffs)

trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_microphysics.F90

@@ -22 +22 @@
     USE MP2M_CLOUDS
     USE MP2M_METHODS
+    USE MP2M_CLOUDS_METHODS
     IMPLICIT NONE
     
@@ -35 +36 @@
     CONTAINS
   
-    FUNCTION muphys_all(m3as_prod,dm0as,dm3as,dm0af,dm3af,dm0ccn,dm3ccn,dm3ices,dmugases) RESULT(ret)
+    FUNCTION muphys_all(m3as_prod,dm0as,dm3as,dm0af,dm3af,dm0ccn,dm3ccn,dm3ices,dmugases,dtlc) RESULT(ret)
         !! Compute the evolution of moments tracers through haze microphysics processes.
         !!
@@ -48 +49 @@
         !! before the latter are called to initialize a new step.
         !!
-
         ! Production of the 3rd order moment of the spherical mode distribution (m3.m-2).
         REAL(kind=mm_wp), INTENT(in), DIMENSION(:) :: m3as_prod
@@ -67 +67 @@
         ! Tendencies of each condensible gaz species (mol.mol-1).
         REAL(kind=mm_wp), INTENT(inout), DIMENSION(:,:) :: dmugases
+        ! Latent heat of condensation (J.kg-1).
+        REAL(kind=mm_wp), INTENT(inout), DIMENSION(:) :: dtlc
 
         ! .true. on succes (i.e. model has been initialized at least once previously), .false. otherwise.
         LOGICAL :: ret
         
-        ! Local variables.
+        ! Local variables:
+        !~~~~~~~~~~~~~~~~~
         INTEGER :: i
         ! Production of the spherical aerosols (m3.m-3).
@@ -82 +85 @@
         ALLOCATE(m3a_s_prod(mm_nla))
 
-        ! Sanity check for initialization
+        ! Sanity check for initialization:
+        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
         ret = (mm_ini_col.AND.mm_ini_aer)
         if (.NOT.ret) then
@@ -95 +99 @@
         m3a_s_prod = m3as_prod(mm_nla:1:-1) / mm_dzlev(:)
 
-        ! Calls haze microphysics (/!\ tendencies in X/m-3)
+        ! Initialize latent heat
+        dtlc(:) = 0._mm_wp
+
+        ! Calls haze microphysics (/!\ tendencies in X/m-3):
+        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
         call mm_haze_microphysics(m3a_s_prod,Hdm0as,Hdm3as,Hdm0af,Hdm3af)
 
-        ! Calls cloud microphysics (/!\ tendencies in X/m-3)
+        ! Calls cloud microphysics (/!\ tendencies in X/m-3):
+        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
         if (mm_call_clouds) then
             call mm_cloud_microphysics(Hdm0as,Hdm3as,Hdm0af,Hdm3af,&
@@ -110 +119 @@
                 dm3ices(:,i)  = dm3ices(mm_nla:1:-1,i) * mm_dzlev(mm_nla:1:-1)
                 dmugases(:,i) = dmugases(mm_nla:1:-1,i)
+
+                ! Compute condensation heating rate
+                dtlc(:) = dtlc(:) + (dmugases(:,i) * mm_xESPS(i)%fmol2fmas * mm_LheatX(mm_temp,mm_xESPS(i)))
             enddo
           
         else
-            Cdm0as(:) = 0._mm_wp ; Cdm3as(:) = 0._mm_wp ; Cdm0af(:) = 0._mm_wp ; Cdm3af(:) = 0._mm_wp
+            Cdm0as(:) = 0._mm_wp ; Cdm3as(:) = 0._mm_wp ; Cdm0af(:)    = 0._mm_wp ; Cdm3af(:)     = 0._mm_wp
             dm0ccn(:) = 0._mm_wp ; dm3ccn(:) = 0._mm_wp ; dm3ices(:,:) = 0._mm_wp ; dmugases(:,:) = 0._mm_wp
+            dtlc(:)   = 0._mm_wp
         endif ! end of mm_call_clouds
         
@@ -140 +153 @@
         !! before the latter are called to initialize a new step.
         !!
-
         ! Production of the 3rd order moment of the spherical mode distribution (m3.m-2).
         REAL(kind=mm_wp), INTENT(in), DIMENSION(:) :: m3as_prod
@@ -155 +167 @@
         LOGICAL :: ret
         
-        ! Local variables.
+        ! Local variables:
+        !~~~~~~~~~~~~~~~~~
         ! Production of the spherical aerosols (m3.m-3).
         REAL(kind=mm_wp), DIMENSION(:), ALLOCATABLE :: m3a_s_prod
         ALLOCATE(m3a_s_prod(mm_nla))
         
-        ! Sanity check for initialization
+        ! Sanity check for initialization:
+        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
         ret = (mm_ini_col.AND.mm_ini_aer)
         if (.NOT.ret) then
@@ -173 +187 @@
         m3a_s_prod = m3as_prod(mm_nla:1:-1) / mm_dzlev(:)
 
-        ! Calls haze microphysics
+        ! Calls haze microphysics (/!\ tendencies in X/m-3):
+        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
         call mm_haze_microphysics(m3a_s_prod,dm0as,dm3as,dm0af,dm3af)