Changeset 4104 for trunk/LMDZ.VENUS/libf

Timestamp:

Mar 6, 2026, 2:34:42 PM (4 weeks ago)

Author:

mlefevre

Message:

Addition of the polysulfur species as tracers in the Venus PCM. ML

File:

• trunk/LMDZ.VENUS/libf/phyvenus/chemparam_mod.F90 (modified) (5 diffs)

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

@@ -22 +22 @@
                  i_ocs, i_hso3, i_h2so4, i_s2,           &
                  i_clso2, i_cl2so2, i_oscl, i_n2, i_he,  &
-                 i_n, i_no, i_no2, i_n2d,                &
+                 i_n, i_no, i_no2, i_n2d, i_h2s, i_hs,   & 
+                 i_s3, i_s4, i_s8,                       &
                  i_co2plus, i_coplus, i_oplus, i_o2plus, &
                  i_n2plus, i_hplus, i_h2oplus, i_nplus,  &
@@ -28 +29 @@
                  i_hcoplus, i_hco2plus, i_elec
 
-integer, save :: i_h2oliq, i_h2so4liq
+integer, save :: i_h2oliq, i_h2so4liq, i_s2con,i_s3con, i_s4con,i_s8con
 
 integer, save :: i_m0_aer, i_m3_aer,                       &
@@ -1554 +1555 @@
       i_no       = 0
       i_no2      = 0
+      i_h2s      = 0
+      i_hs       = 0
+      i_s3       = 0
+      i_s4       = 0
+      i_s8       = 0      
 
       ! ions
@@ -1577 +1583 @@
       i_h2oliq   = 0
       i_h2so4liq = 0
+      i_s2con    = 0
+      i_s3con    = 0
+      i_s4con    = 0
+      i_s8con    = 0      
 
       ! MICROPHYSICAL TRACERS FOR CL_SCHEME=2
@@ -1803 +1813 @@
                m_tr(i_no2) = 46.0055
                type_tr(i_no2) = 1
+            case('h2s')
+                i_h2s = i
+                print*,'h2s',i_h2s
+                m_tr(i_h2s) = 34.082
+                type_tr(i_h2s) = 1
+            case('hs')
+                i_hs = i
+                print*,'hs',i_hs
+                m_tr(i_hs) = 33.073
+                type_tr(i_hs) = 1
+            case('s3')
+                i_s3=i
+                PRINT*,'s3',i_s3
+                m_tr(i_s3) = 96.198
+                type_tr(i_s3) = 1
+            case('s4')
+                i_s4 = i
+                PRINT*,'s4',i_s4
+                m_tr(i_s4) = 128.26
+                type_tr(i_s4) = 1                  
+            case('s8')
+                i_s8 = i
+                PRINT*,'s8',i_s8
+                m_tr(i_s8) = 256.52
+                type_tr(i_s8) = 1                
+            case('s2con')
+                i_s2con = i
+                PRINT*,'s2con',i_s2con
+                m_tr(i_s2con) = 64.13
+                type_tr(i_s2con) = 1
+            case('s3con')
+                i_s3con=i
+                PRINT*,'s3con',i_s3con
+                m_tr(i_s3con) = 96.198
+                type_tr(i_s3con) = 1
+            case('s4con')
+                i_s4con=i
+                PRINT*,'s4con',i_s4con
+                m_tr(i_s4con) = 128.26
+                type_tr(i_s4con) = 1
+            case('s8con')
+                i_s8con = i
+                PRINT*,'s8con',i_s8con
+                m_tr(i_s8con) = 256.52
+                type_tr(i_s8con) = 1                 
 
             ! ions