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

Timestamp:

Nov 17, 2025, 11:02:59 AM (3 months ago)

Author:

debatzbr

Message:

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

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• mp2m_calmufi.F90 (modified) (9 diffs)
• physiq_mod.F90 (modified) (3 diffs)

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: