Changeset 3962

Timestamp:

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

Author:

debatzbr

Message:

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

Location:

trunk/LMDZ.PLUTO/libf

Files:

• muphypluto/mp2m_microphysics.F90 (modified) (10 diffs)
• phypluto/mp2m_calmufi.F90 (modified) (9 diffs)
• phypluto/physiq_mod.F90 (modified) (3 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)
         

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

@@ -1 +1 @@
 MODULE mp2m_calmufi
     use tracer_h
-    use comcstfi_mod, only : mugaz
+    use comcstfi_mod, only : mugaz, cpp
     use callkeys_mod, only : call_haze_prod_pCH4, haze_rho,&
                              callmuclouds
@@ -44 +44 @@
     CONTAINS
     
-    SUBROUTINE calmufi(dt, plev, zlev, play, zlay, g3d, temp, pq, zdqfi, zdqmufi_prod, zdqmufi)
+    SUBROUTINE calmufi(dt, plev, zlev, play, zlay, g3d, temp, pq, zdqfi, zdqmufi_prod, zdqmufi, zdtcond)
         !! Interface subroutine to YAMMS model for LMD PCM.
         !!
@@ -62 +62 @@
         REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: zdqmufi_prod ! Aerosols production tendency (kg/kg_of_air/s).
         REAL(kind=8), DIMENSION(:,:,:), INTENT(OUT) :: zdqmufi      ! Microphysical tendency for tracers (X.m-2 --> X.kg-1.s-1).
+        REAL(kind=8), DIMENSION(:,:),   INTENT(OUT) :: zdtcond      ! Condensation heating rate (K.s-1).
         
         ! Local tracers:
@@ -90 +91 @@
         REAL(kind=8), DIMENSION(:,:), ALLOCATABLE :: dm3ices  ! Tendencies of the 3rd order moments of each ice components (m3.m-2). 
         REAL(kind=8), DIMENSION(:,:), ALLOCATABLE :: dmugases ! Tendencies of each condensible gas species (mol.mol-1).
+        REAL(kind=8), DIMENSION(:), ALLOCATABLE   :: dtlc     ! Latent heat of condensation (J.kg-1).
         
         ! Local variables:
@@ -129 +131 @@
         ALLOCATE(dm3ices(nlay,nmicro_ices))
         ALLOCATE(dmugases(nlay,nmicro_ices))
+        ALLOCATE(dtlc(nlay))
 
         ALLOCATE(int2ext(nlon,nlay))  
@@ -138 +141 @@
         ! Initialization of zdqmufi here since intent=out and no action performed on every tracers
         zdqmufi(:,:,:) = 0.D0
+
+        ! Initialization of zdtcond here since intent=out
+        zdtcond(:,:) = 0.D0
         
         ! Initialize tracers updated with former processes from physics
@@ -191 +197 @@
             
             ! Initializes tendencies
-            dm0as(:) = 0._mm_wp ; dm3as(:) = 0._mm_wp ; dm0af(:) = 0._mm_wp ; dm3af(:) = 0._mm_wp
+            dm0as(:)  = 0._mm_wp ; dm3as(:)  = 0._mm_wp ; dm0af(:)     = 0._mm_wp ; dm3af(:)      = 0._mm_wp
             dm0ccn(:) = 0._mm_wp ; dm3ccn(:) = 0._mm_wp ; dm3ices(:,:) = 0._mm_wp ; dmugases(:,:) = 0._mm_wp
+            dtlc(:)   = 0._mm_wp
             
         !----------------------------
@@ -200 +207 @@
             ! Call microphysics
             if (callmuclouds) then
-                if(.NOT.mm_muphys(m3as_prod,dm0as,dm3as,dm0af,dm3af,dm0ccn,dm3ccn,dm3ices,dmugases)) then
+                if(.NOT.mm_muphys(m3as_prod,dm0as,dm3as,dm0af,dm3af,dm0ccn,dm3ccn,dm3ices,dmugases,dtlc)) then
                     call abort_program(error("mm_muphys (clouds) aborted -> initialization not done !",-1))
                 endif
@@ -232 +239 @@
                     zdqmufi(ilon,:,micro_gas_indx(i)) = dmugases(:,i) * (mmol(micro_gas_indx(i))/mugaz)
                 enddo
+
+                ! Compute condensation heating rate in K.s-1
+                zdtcond(ilon,:) = dtlc(:) / cpp / dt
             endif ! End of callmuclouds
             

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

@@ -411 +411 @@
       real pdqmufi_prod(ngrid,nlayer,nq) ! Aerosols production tendency (kg/kg_of_air/s).
       real int2ext(ngrid,nlayer)         ! Intensive to extensive factor (kg_air/m3: X/kg_air --> X/m3).
+      real zdtcond(ngrid,nlayer)         ! Condensation heating rate (K.s-1).
 
 ! Local variables for LOCAL CALCULATIONS:
@@ -1617 +1618 @@
             pdqmufi(:,:,:) = 0.
 
-            call calmufi(ptimestep,pplev,zzlev,pplay,zzlay,gzlat,pt,pq,pdq,pdqmufi_prod,pdqmufi)
+            call calmufi(ptimestep,pplev,zzlev,pplay,zzlay,gzlat,pt,pq,pdq,pdqmufi_prod,pdqmufi,zdtcond)
 
             pdq(:,:,:) = pdq(:,:,:) + pdqmufi(:,:,:)
@@ -2536 +2537 @@
                   call write_output("m3"//TRIM(str(6:)),"Volume of "//TRIM(str(6:))//" ice","m3.m-3",zq(:,:,micro_ice_indx(iq))*int2ext(:,:))
                enddo
+               call write_output("dtcond","Condensation heating rate","K.s-1",zdtcond(:,:))
 
                ! Diagnostics: