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

Timestamp:

Apr 26, 2024, 4:27:26 PM (7 months ago)

Author:

slebonnois

Message:

Titan PCM update : optics + microphysics

File:

• trunk/LMDZ.TITAN/libf/phytitan/calc_ysat.F90 (modified) (1 diff)

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

@@ -82 +82 @@
       else if(trim(cnames(ic)).eq."C2H4") then
          ! not far from Fray and Schmidt (2009)
-         where (temp(:,:).lt.89.0)
-            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)
-            ysat(:,:,ic) = 10.0**(8.724e0 - 901.6e0/(temp(:,:) - 2.555e0) )                  &
-                  / press(:,:) * 1013.25e0 / 760.0
-         elsewhere (temp(:,:).lt.120.0)
-            ysat(:,:,ic) = 10.0**(50.79e0 - 1703.0e0/temp(:,:) - 17.141e0                    &
-                  * alog10(temp(:,:)) ) / press(:,:) * 1013.25e0 / 760.0
-         elsewhere (temp(:,:).lt.155.0)
-            ysat(:,:,ic) = 10.0**(6.74756e0 - 585.0e0/(temp(:,:) - 18.16e0) )                &
-                  / press(:,:) * 1013.25e0 / 760.0
-         endwhere
+         !where (temp(:,:).lt.89.0)
+         !   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)
+         !   ysat(:,:,ic) = 10.0**(8.724e0 - 901.6e0/(temp(:,:) - 2.555e0) )                  &
+         !         / press(:,:) * 1013.25e0 / 760.0
+         !elsewhere (temp(:,:).lt.120.0)
+         !   ysat(:,:,ic) = 10.0**(50.79e0 - 1703.0e0/temp(:,:) - 17.141e0                    &
+         !         * alog10(temp(:,:)) ) / press(:,:) * 1013.25e0 / 760.0
+         !elsewhere (temp(:,:).lt.155.0)
+         !   ysat(:,:,ic) = 10.0**(6.74756e0 - 585.0e0/(temp(:,:) - 18.16e0) )                &
+         !         / press(:,:) * 1013.25e0 / 760.0
+         !endwhere
+
+         ! Fray and Schmidt (2009)
+         ysat(:,:,ic) = (1.e3 / press(:,:)) * exp(1.540e1 - 2.206e3/temp(:,:)                &
+                         - 1.216e4/temp(:,:)**2 + 2.843e5/temp(:,:)**3 - 2.203e6/temp(:,:)**4)
       
       ! C2H6 :