Changeset 3497 for trunk/LMDZ.TITAN/libf/phytitan/evapCH4.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/evapCH4.F90 (modified) (4 diffs)

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

@@ -79 +79 @@
 real, parameter :: humCH4 = 0.4 ! Imposed surface humidity for CH4 [-]
 
-real, parameter :: Flnp = 0.08  ! Fraction occupied by lakes (North Pole)
-real, parameter :: Flsp = 0.02  ! Fraction occupied by lakes (South Pole)
-real, parameter :: Flml = 0.02  ! Fraction not infiltrated into the ground (Mid latitudes)
+real, parameter :: Flnp = 0.10  ! Fraction occupied by lakes (North Pole)
+real, parameter :: Flsp = 0.01  ! Fraction occupied by lakes (South Pole)
+real, parameter :: Flml = 1.    ! Fraction not infiltrated into the ground (Mid latitudes)
 
 real, parameter :: mmolair = 28.e-3 ! Molar mass of air [kg.mol-1]
 real, parameter :: mmolCH4 = 16.e-3 ! Molar mass of CH4 [kg.mol-1]
-real, parameter :: rhoiCH4 = 425.   ! Density of ice of CH4 [kg.m-3]
+real, parameter :: rholCH4 = 425.   ! Density of liquid CH4 [kg.m-3]
 
 real, parameter :: TcCH4  = 190.56  ! Critical point of CH4 [K]
@@ -131 +131 @@
   Tlake = tsurf(ig) - 7 ! Lakes are 2-7K less than surface.
   qsatCH4 = (1.0e5 / pplev(ig,1)) * exp(1.051e1 - 1.110e3/Tlake - 4.341e3/Tlake**2 + 1.035e5/Tlake**3 - 7.910e5/Tlake**4)
-  ! CH4 : 0.80 * qsat because of dissolution in N2
-  qsatCH4 = 0.80 * qsatCH4
+  ! CH4 : 0.70 * qsat because of dissolution in N2
+  qsatCH4 = 0.70 * qsatCH4
 
   ! Flux at the surface [kg.m-2.s-1] :
@@ -171 +171 @@
 
   ! Trend of CH4's tank [m]
-  dtankCH4 = - (newpqCH4 - pqCH4(ig,1)) * rhoair * (zzlev(ig,2) - zzlev(ig,1)) * mmolCH4 / mmolair / rhoiCH4
+  dtankCH4 = - (newpqCH4 - pqCH4(ig,1)) * rhoair * (zzlev(ig,2) - zzlev(ig,1)) * mmolCH4 / mmolair / rholCH4
 
   ! New tank depth of CH4 [m]
@@ -177 +177 @@
     tankCH4(ig) = tankCH4(ig) + dtankCH4
   else
-    !write(*,*) 'Evaporation CH4 : Empty lakes...'
-    newpqCH4 = tankCH4(ig) / (rhoair * (zzlev(ig,2) - zzlev(ig,1)) * mmolCH4 / mmolair / rhoiCH4) + pqCH4(ig,1)
+    newpqCH4 = tankCH4(ig) / (rhoair * (zzlev(ig,2) - zzlev(ig,1)) * mmolCH4 / mmolair / rholCH4) + pqCH4(ig,1)
     tankCH4(ig) = 1.e-30
   endif