Changeset 3597 for trunk/LMDZ.VENUS/libf/phyvenus/photolysis_mod.F90

Timestamp:

Jan 22, 2025, 3:47:23 PM (38 hours ago)

Author:

flefevre

Message:

Updated S2 ultraviolet cross-section from Frandsen (personal communication, 2024).
The cross-section file can be downloaded from

​https://owncloud.latmos.ipsl.fr/index.php/s/QfazKRX8qfEz97F

and must be put in the INPUT/cross_sections directory.

File:

• trunk/LMDZ.VENUS/libf/phyvenus/photolysis_mod.F90 (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/photolysis_mod.F90

@@ -1571 +1571 @@
       REAL y1(kdata)
       INTEGER ierr, dummy
-      INTEGER i, n
-
+      INTEGER i, n, iopt
       CHARACTER*100 fil
       integer :: kin, kout ! input/output logical units
@@ -1578 +1577 @@
       kin = 10
 
-      fil = 'cross_sections/s2_millsBestEst_1560_5059.txt'
-      print*, 'section efficace S2: ', fil
+!---------------------------------------------------------------------------------
+!     iopt = 1       Frank Mills' best estimate (1998)
+!     iopt = 2       Benjamin Frandsen calculation (personal communication, 2024)
+!---------------------------------------------------------------------------------
+
+      iopt = 2
+
+      if (iopt == 1) then
+         fil = 'cross_sections/s2_millsBestEst_1560_5059.txt'
+      else if (iopt == 2) then
+         fil = 'cross_sections/s2_cross_sections_frandsen.txt'
+      end if
 
       if(is_master) then
 
+      print*, 'S2 cross-section : ', fil
       OPEN(UNIT=kin,FILE=fil,STATUS='old') 
-
-      n = 1203 
-      DO i = 1, n
-        READ(kin,*) x1(i), y1(i)
-        x1(i) = x1(i)/10
-      END DO
+      
+      if (iopt == 1) then
+         n = 1203 
+         do i = 1, n
+            read(kin,*) x1(i), y1(i)
+            x1(i) = x1(i)/10.
+         end do
+      else if (iopt == 2) then
+         n = 406 
+         do i = 1, 4
+            read(kin,*)
+         end do
+         do i = 1, n
+            read(kin,*) x1(i), y1(i)
+         end do
+      end if
+
       CLOSE (kin)