Changeset 3682 for trunk/LMDZ.TITAN/libf/muphytitan

Timestamp:

Mar 14, 2025, 2:27:06 PM (5 months ago)

Author:

lrosset

Message:

Titan microphysics and physics : Outputs the nucleation rates and growth rates for the condensible species. (LR)

Location:

trunk/LMDZ.TITAN/libf/muphytitan

Files:

• mm_clouds.f90 (modified) (9 diffs)
• mm_globals.f90 (modified) (3 diffs)
• mm_microphysic.f90 (modified) (5 diffs)

trunk/LMDZ.TITAN/libf/muphytitan/mm_clouds.f90

@@ -101 +101 @@
     IF (mm_w_cloud_nucond) THEN
       ! Calls condensation/nucleation (and update saturation ratio diagnostic)
-      call mm_cloud_nucond(dm0a,dm3a,dm0n,dm3n,dm3i,dgazs,mm_gazs_sat)
+      ! ADDED : Extraction of nucleation and growth rates
+      call mm_cloud_nucond(dm0a,dm3a,dm0n,dm3n,dm3i,dgazs,mm_gazs_sat,mm_nrate,mm_grate)
     ENDIF
 
@@ -134 +135 @@
   !-----------------------------------------------------------------------------
 
-  SUBROUTINE mm_cloud_nucond(dm0a,dm3a,dm0n,dm3n,dm3i,dgazs,gazsat)
+  SUBROUTINE mm_cloud_nucond(dm0a,dm3a,dm0n,dm3n,dm3i,dgazs,gazsat,nrate,grate)
     !! Get moments tendencies through nucleation/condensation/evaporation.
     !!
@@ -166 +167 @@
     REAL(kind=mm_wp), DIMENSION(:,:), INTENT(out) :: gazsat
     !! Saturation ratio of each condensible specie.
+
+    ! ADDED : Extraction of the nucleation and growth rate
+    REAL(kind=mm_wp), DIMENSION(:,:),INTENT(out) :: nrate
+    !! Nucleation rate values of each condensible species (\(m^{-2}.s^{-1}\)).
+    REAL(kind=mm_wp), DIMENSION(:,:),INTENT(out) :: grate
+    !! Growth rate values of each condensible species (\(m^{2}.s^{-1}\)).
+
     INTEGER                                       :: i,idx
     REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE :: zdm0a,zdm3a,zdm0n,zdm3n
@@ -175 +183 @@
     DO i = 1, mm_nesp
       call nc_esp(mm_xESPS(i),mm_gazs(:,i),mm_m3ice(:,i),dgazs(:,i),dm3i(:,i), &
-                  zdm0a(:,i),zdm3a(:,i),zdm0n(:,i),zdm3n(:,i),gazsat(:,i))
+                  zdm0a(:,i),zdm3a(:,i),zdm0n(:,i),zdm3n(:,i),gazsat(:,i),nrate(:,i),grate(:,i))
+      ! ADDED : Extraction of nucleation and growth rates
     ENDDO
 
@@ -200 +209 @@
   END SUBROUTINE mm_cloud_nucond
 
-  SUBROUTINE nc_esp(xESP,vapX,m3iX,dvapX,dm3iX,dm0aer,dm3aer,dm0ccn,dm3ccn,Xsat)
+  SUBROUTINE nc_esp(xESP,vapX,m3iX,dvapX,dm3iX,dm0aer,dm3aer,dm0ccn,dm3ccn,Xsat,Xnrate,Xgrate) ! ADDED : arguments nrate, grate
     !! Get moments tendencies through nucleation/condensation/evaporation of a given condensible specie.
     !!
@@ -230 +239 @@
     REAL(kind=mm_wp),INTENT(out), DIMENSION(:) :: Xsat
     !! Saturation ratio values on the vertical grid (--).
+
+    ! ADDED : Extraction of the nucleation and growth rate
+    REAL(kind=mm_wp),INTENT(out), DIMENSION(:) :: Xnrate
+    !! Nucleation rate values on the vertical grid for the species X (\(m^{-2}.s^{-1}\)).
+    REAL(kind=mm_wp),INTENT(out), DIMENSION(:) :: Xgrate
+    !! Growth rate values on the vertical grid for the species X (\(m^{2}.s^{-1}\)).
+
+
     INTEGER                                 :: i
     REAL(kind=mm_wp)                        :: bef,aft
@@ -281 +298 @@
     ! Then, from eq. 2:
     ! Mn(k)[t+dt] = Mn(k)[t] + CST_M(k)/(1+CST_M(k))*Ma(k)[t]            (4)
+    
+    ! Copies the nucleation rate into an output variable
+    Xnrate = nucr
+    
     cm0 = 4._mm_wp*mm_pi*nucr/mm_rm*mm_alpha_f(3._mm_wp)*mm_rcf**3*mm_dt
     cm3 = 4._mm_wp*mm_pi*nucr/mm_rm*mm_alpha_f(6._mm_wp)/mm_alpha_f(3._mm_wp)*mm_rcf**3*mm_dt
@@ -317 +338 @@
     ! gets "true" growth rate
     grate = grate * (newvap/qsat - seq)
+    ! Copies the growth rate into an output variable
+    Xgrate = grate
 
     ! computes tendencies through condensation
@@ -377 +400 @@
     REAL(kind=mm_wp), INTENT(in), DIMENSION(:)  :: sat  !! Saturation ratio of given specie (--).
     TYPE(mm_esp), INTENT(in)                    :: xESP !! X specie properties (--).
-    REAL(kind=mm_wp), INTENT(out), DIMENSION(:) :: rate !! The nucleation rate (\(m^-{2}.s^{-1}\)).
+    REAL(kind=mm_wp), INTENT(out), DIMENSION(:) :: rate !! The nucleation rate (\(m^{-2}.s^{-1}\)).
     INTEGER          :: i
     REAL(kind=mm_wp) :: r,t,s,sig,nX,rstar,gstar,x,zeldov,deltaf,fsh,fstar

trunk/LMDZ.TITAN/libf/muphytitan/mm_globals.f90

@@ -522 +522 @@
   REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE, SAVE :: mm_ice_fluxes
 
+  !> Condensible components nucleation rates.
+  !!
+  !! It is a 2D-array with the vertical layers in first dimension and the number of condensible
+  !! species in the second.
+  !! It is updated in [[mm_clouds(module):mm_cloud_microphysics(subroutine)]].
+  REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE, SAVE :: mm_nrate
+
+  !> Condensible components growth rates.
+  !!
+  !! It is a 2D-array with the vertical layers in first dimension and the number of condensible
+  !! species in the second.
+  !! It is updated in [[mm_clouds(module):mm_cloud_microphysics(subroutine)]].
+  REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE, SAVE :: mm_grate
+
   !> Condensible species saturation ratio (--).
   !!
@@ -583 +597 @@
   !$OMP THREADPRIVATE(mm_rcs,mm_rcf,mm_drad,mm_drho)
   !$OMP THREADPRIVATE(mm_m0as_vsed,mm_m3as_vsed,mm_m0af_vsed,mm_m3af_vsed)
-  !$OMP THREADPRIVATE(mm_aer_s_flux,mm_aer_f_flux,mm_ccn_vsed,mm_ccn_flux,mm_ice_prec,mm_ice_fluxes,mm_gazs_sat)
+  !$OMP THREADPRIVATE(mm_aer_s_flux,mm_aer_f_flux,mm_ccn_vsed,mm_ccn_flux,mm_ice_prec,mm_ice_fluxes,mm_gazs_sat, mm_nrate, mm_grate)
   !$OMP THREADPRIVATE(mm_m0as_min,mm_m3as_min,mm_rcs_min,mm_m0af_min,mm_m3af_min,mm_rcf_min,mm_m0n_min,mm_m3cld_min)
   !$OMP THREADPRIVATE(mm_nla,mm_nle)
@@ -1251 +1265 @@
       ALLOCATE(mm_gazs_sat(mm_nla,mm_nesp)) ; mm_gazs_sat(:,:) = 0._mm_wp
     ENDIF
+    ! ADDED : Allocate memory for nucleation and growth rate
+    IF (.NOT.ALLOCATED(mm_nrate)) THEN
+      ALLOCATE(mm_nrate(mm_nla,mm_nesp)) ; mm_nrate(:,:) = 0._mm_wp
+    ENDIF
+      IF (.NOT.ALLOCATED(mm_grate)) THEN
+      ALLOCATE(mm_grate(mm_nla,mm_nesp)) ; mm_grate(:,:) = 0._mm_wp
+    ENDIF
 
     ! note mm_dzlev already from top to ground

trunk/LMDZ.TITAN/libf/muphytitan/mm_microphysic.f90

@@ -163 +163 @@
   END FUNCTION muphys_nocld
 
-  SUBROUTINE mm_diagnostics(dt,aer_prec,aer_s_w,aer_f_w,aer_s_flux,aer_f_flux,ccn_prec,ccn_w,ccn_flux,ice_prec,ice_fluxes,gazs_sat)
+  SUBROUTINE mm_diagnostics(dt,aer_prec,aer_s_w,aer_f_w,aer_s_flux,aer_f_flux,ccn_prec,ccn_w,ccn_flux,ice_prec,ice_fluxes,gazs_sat,nrate,grate)
     !! Get various diagnostic fields of the microphysics.
     !!
@@ -171 +171 @@
     !! - Settling velocity (aerosols -total-, CCN and ices)
     !! - Precipitations (aerosols -total-, CCN and ices)
-    !! - condensible gazes saturation ratio
+    !! - Condensible gazes saturation ratio
+    !! - Nucleation and growth rates
     !!
     !! @note
@@ -197 +198 @@
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:) :: gazs_sat   !! Condensible gaz saturation ratios (--).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: ice_prec   !! Ice precipitations (\(kg.m^{-2}.s^{-1}\)).
+    ! Nucleation and growth rates
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:)   :: nrate     !! Nucleation rate (\(m^{-2}.s^{-1}\)).
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:)   :: grate     !! Growth rate (\(m^{2}.s^{-1}\)).
 
     IF (PRESENT(aer_prec))   aer_prec   = ABS(mm_aer_prec) / dt
@@ -211 +215 @@
       IF (PRESENT(ice_fluxes)) ice_fluxes = mm_ice_fluxes(mm_nla:1:-1,:)
       IF (PRESENT(gazs_sat))   gazs_sat   = mm_gazs_sat(mm_nla:1:-1,:)
+      ! Nucleation and growth rates
+      IF (PRESENT(nrate))      nrate      = mm_nrate(mm_nla:1:-1,:)
+      IF (PRESENT(grate))      grate      = mm_grate(mm_nla:1:-1,:)
     ELSE
       IF (PRESENT(ccn_prec))   ccn_prec   = 0._mm_wp
@@ -218 +225 @@
       IF (PRESENT(ice_fluxes)) ice_fluxes = 0._mm_wp
       IF (PRESENT(gazs_sat))   gazs_sat   = 0._mm_wp
+      ! Nucleation and growth rates
+      IF (PRESENT(nrate))      nrate      = 0._mm_wp
+      IF (PRESENT(grate))      grate      = 0._mm_wp
     ENDIF
   END SUBROUTINE mm_diagnostics