Changeset 3497 for trunk/LMDZ.TITAN/libf/phytitan/calc_condlh.F90

Timestamp:

Nov 7, 2024, 11:13:46 AM (2 weeks ago)

Author:

debatzbr

Message:

Add AC6H6 condensation in the microphysics

File:

• trunk/LMDZ.TITAN/libf/phytitan/calc_condlh.F90 (modified) (5 diffs)

trunk/LMDZ.TITAN/libf/phytitan/calc_condlh.F90

@@ -33 +33 @@
 !   C2H6 --> 3
 !   HCN  --> 4
+!   AC6H6--> 5
 !====================================================================    
 
@@ -56 +57 @@
 real, parameter :: mmolC2H6 = 30.e-3 ! Molar mass of C2H6 [kg.mol-1]
 real, parameter :: mmolHCN  = 27.e-3 ! Molar mass of HCN [kg.mol-1]
+real, parameter :: mmolAC6H6 = 78.e-3 ! Molar mass of AC6H6 [kg.mol-1]
 
 real, parameter :: TcCH4  = 190.56 ! Critical point of CH4  [K]
@@ -61 +63 @@
 real, parameter :: TcC2H6 = 305.35 ! Critical point of C2H6 [K]
 real, parameter :: TcHCN  = 456.70 ! Critical point of HCN [K]
+real, parameter :: TcAC6H6 = 562.10 ! Critical point of AC6H6 [K]
 
 ! Local variables :
@@ -111 +114 @@
             enddo ! End of loop on vertical layers
         enddo ! End of main loop on the horizontal grid
-
+    
+    ! AC6H6 :
+    !--------
+    else if(trim(nameOfTracer(gazs_indx(iq))) .eq. "AC6H6") then
+        
+        do ig = 1, ngrid ! Main loop on the horizontal grid
+            do l = 1, nlayer ! Loop on vertical layers
+                ftx = (1. - pt(ig,l)) / TcAC6H6
+                if (ftx .le. 1.e-3) ftx = 1.e-3
+                Lcx(ig,l,iq) = 8.314 * TcAC6H6 * (7.08 * ftx**0.354 + 10.95 * 9.9e-2 * ftx**0.456) / mmolAC6H6
+            enddo ! End of loop on vertical layers
+        enddo ! End of main loop on the horizontal grid
 
     ! HCN :
@@ -127 +141 @@
 
     else
-        write(*,*) 'Species ',trim(nameOfTracer(gazs_indx(iq))),' not considered... Check traceur.def'
+        write(*,*) 'Species ',trim(nameOfTracer(gazs_indx(iq))),' not considered in calc_condlh.F90... Check traceur.def'
         Lcx(ig,l,iq) = 0.0D0
     endif