Changeset 2326 for trunk/LMDZ.TITAN/libf/phytitan/calc_ysat.F90

Timestamp:

May 18, 2020, 4:33:35 PM (5 years ago)

Author:

jvatant

Message:

Update Titan reference photochemistry (reaction constants,branching ratios, condensation rates) according to Vuitton et al 2019.
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/calc_ysat.F90 (modified) (4 diffs)

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

@@ -10 +10 @@
   !     S. Lebonnois 
   !     -> inicondens.F in old Titan, with T,P in planetary average
-  !     J. Vatant d'Ollone (2017)
-  !     -> Adapt to new physics, move to F90 and compute every grid point 
+  !     J. Vatant d'Ollone
+  !       -> 2017 : Adapt to new physics, move to F90 and compute every grid point 
+  !       -> 2018 : Update saturation profiles from recent litterature ( Vuitton 2019 ) 
   !     ==============================================================================
 
@@ -62 +63 @@
      else if(trim(cnames(ic)).eq."C2H2") then
 
-        ysat(:,:,ic) = 10.0**(6.09748e0-1644.1e0/temp(:,:)+7.42346e0                        &
-             * alog10(1.0e3/temp(:,:)) ) / press(:,:)*1013.25e0/760.0
+!        ysat(:,:,ic) = 10.0**(6.09748e0-1644.1e0/temp(:,:)+7.42346e0                        &
+!             * alog10(1.0e3/temp(:,:)) ) / press(:,:)*1013.25e0/760.0
+
+         ysat(:,:,ic) = 1.E5 * exp(13.4-2536./temp(:,:)) / press(:,:) ! Fray and Schmidt (2009)
 
      else if(trim(cnames(ic)).eq."C2H4") then
 
-        where (temp(:,:).lt.89.0)
+        where (temp(:,:).lt.89.0) ! not far from Fray and Schmidt, 2009
            ysat(:,:,ic) = 10.0**(1.5477e0 + (1.0e0/temp(:,:) - 0.011e0)                     &
                 * (16537.0e0*(1.0e0/temp(:,:) - 0.011e0) - 1038.1e0))                       &
                 / press(:,:) * 1.01325e0 / 760.0
-        elsewhere (temp(:,:).lt.104.0)
+        elsewhere (temp(:,:).lt.104.0) ! not far from Fray and Schmidt, 2009
            ysat(:,:,ic) = 10.0**(8.724e0 - 901.6e0/(temp(:,:) - 2.555e0) )                  &
                 / press(:,:) * 1013.25e0 / 760.0
@@ -85 +88 @@
 
         where (temp(:,:).lt.90.)
-           ysat(:,:,ic) = 10.0**(10.01e0-1085.0e0/(temp(:,:)-0.561e0) )                     &
-                / press(:,:) * 1013.25e0 / 760.0e0
+!           ysat(:,:,ic) = 10.0**(10.01e0-1085.0e0/(temp(:,:)-0.561e0) )                     &
+!                / press(:,:) * 1013.25e0 / 760.0e0
+           ysat(:,:,ic) = 1.E5 * exp ( 15.11 - 2207./temp(:,:) - 24110./(temp(:,:)**2)       &
+                        + 7.744E5/(temp(:,:)**3) - 1.161E7/(temp(:,:)**4)                    &
+                        + 6.763E7/(temp(:,:)**5) ) / press(:,:) ! Fray and Schmidt (2009)
         elsewhere
            ysat(:,:,ic) = 10.0**(5.9366e0 - 1086.17e0/temp(:,:) + 3.83464e0                 &
@@ -135 +141 @@
      else if(trim(cnames(ic)).eq."AC6H6")  then
 
-        x = 1.0e0 - temp(:,:) / 562.2e0
-        ysat(:,:,ic)= 48.9e3 * exp( ( 1.33213 * x**1.5 - 6.98273 * x                        &
-             - x**3 * (2.62863 + 3.33399 * x**3) ) * 562.2e0/temp(:,:) ) / press(:,:)
+        !x = 1.0e0 - temp(:,:) / 562.2e0
+        !ysat(:,:,ic)= 48.9e3 * exp( ( 1.33213 * x**1.5 - 6.98273 * x                        &
+        !     - x**3 * (2.62863 + 3.33399 * x**3) ) * 562.2e0/temp(:,:) ) / press(:,:)
 
+        ysat(:,:,ic) = 1.E5 * exp (17.35-5663./temp(:,:)) / press(:,:) ! Fray and Schmidt (2009)
+        
      else if(trim(cnames(ic)).eq."HCN")  then
 
-        ysat(:,:,ic)= 10.0**(8.6165e0 - 1516.5e0/(temp(:,:) - 26.2e0) ) / press(:,:) * 1013.25e0 / 760.0e0
+        !ysat(:,:,ic)= 10.0**(8.6165e0 - 1516.5e0/(temp(:,:) - 26.2e0) ) / press(:,:) * 1013.25e0 / 760.0e0
+        
+        ysat(:,:,ic) = 1.E5 * exp ( 13.93 - 3624./temp(:,:) - 1.325E5/(temp(:,:)**2)       &
+                     + 6.314E6/(temp(:,:)**3) - 1.128E8/(temp(:,:)**4) ) / press(:,:)  ! Fray and Schmidt (2009)
 
      else if(trim(cnames(ic)).eq."CH3CN")  then