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

Timestamp:

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

Author:

lrosset

Message:

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

File:

• trunk/LMDZ.TITAN/libf/muphytitan/mm_clouds.f90 (modified) (9 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