Changeset 4118


Ignore:
Timestamp:
Mar 10, 2026, 5:51:25 PM (5 weeks ago)
Author:
mlefevre
Message:

Correction to the polysulfur chemical network in the Venus PCM. ML

Location:
trunk/LMDZ.VENUS/libf/phyvenus
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90

    r4111 r4118  
    20132013
    20142014!===========================================================
    2015 !      g063: SO3 + CO -> CO2 + SO2
    2016 !===========================================================
    2017 
    2018 nb_reaction_4 = nb_reaction_4 + 1
    2019 
    2020 indice_4(nb_reaction_4) = z4spec(1.0, i_so3, 1.0, i_co, 1.0, i_co2, 1.0, i_so2)
    2021 
    2022 !===========================================================
    2023 !      g064: SO2 + CO -> CO2 + SO
    2024 !===========================================================
    2025 
    2026 nb_reaction_4 = nb_reaction_4 + 1
    2027 
    2028 indice_4(nb_reaction_4) = z4spec(1.0, i_so2, 1.0, i_co, 1.0, i_co2, 1.0, i_so)
    2029 
    2030 !===========================================================
    2031 !      g065: S2 + CO -> OCS + S
     2015!      g063: SO3 + CO -> SO2 + CO2
     2016!===========================================================
     2017
     2018nb_reaction_4 = nb_reaction_4 + 1
     2019
     2020indice_4(nb_reaction_4) = z4spec(1.0, i_so3, 1.0, i_co, 1.0, i_so2, 1.0, i_co2)
     2021
     2022!===========================================================
     2023!      g064: S2 + CO -> OCS + S
    20322024!===========================================================
    20332025
    … …  
    20372029
    20382030!===========================================================
    2039 !      g066: H + OCS -> CO + HS
     2031!      g065: H + OCS -> CO + HS
    20402032!===========================================================
    20412033
    … …  
    20452037
    20462038!===========================================================
    2047 !      g067: CO + HS -> H + OCS
    2048 !===========================================================
    2049 
    2050 nb_reaction_4 = nb_reaction_4 + 1
    2051 
    2052 indice_4(nb_reaction_4) = z4spec(1.0, i_co, 1.0, i_hs, 1.0, i_h, 1.0, i_ocs)
    2053 
    2054 !===========================================================
    2055 !      g068: S + SO2 -> SO + SO
     2039!      g066: S + SO2 -> SO + SO
    20562040!===========================================================
    20572041
    … …  
    20612045
    20622046!===========================================================
    2063 !      g069: HS + CO -> OCS + H
     2047!      g067: HS + CO -> OCS + H
    20642048!===========================================================
    20652049
    … …  
    20692053
    20702054!===========================================================
    2071 !      g070: SO2 + H -> OH + SO
     2055!      g068: SO2 + H -> OH + SO
    20722056!===========================================================
    20732057
    … …  
    20772061
    20782062!===========================================================
    2079 !      g071: OCS + OH -> CO2 + HS
     2063!      g069: OCS + OH -> CO2 + HS
    20802064!===========================================================
    20812065
    … …  
    20852069
    20862070!===========================================================
    2087 !      g072: S + S4 -> S3 + S2
     2071!      g070: S + S4 -> S3 + S2
    20882072!===========================================================
    20892073
    … …  
    20932077
    20942078!===========================================================
    2095 !      g073: H2 + HS -> H + H2S
     2079!      g071: H2 + HS -> H + H2S
    20962080!===========================================================
    20972081
    … …  
    21012085
    21022086!===========================================================
    2103 !      g074: H + S3 -> HS + S2
    2104 !===========================================================
    2105 
    2106 nb_reaction_4 = nb_reaction_4 + 1
    2107 
    2108 indice_4(nb_reaction_4) = z4spec(1.0, i_h, 1.0, i_s3, 1.0, i_hs, 1.0, i_s2)
    2109 
    2110 !===========================================================
    2111 !      g075: SO3 + CO -> CO2 + SO2
    2112 !===========================================================
    2113 
    2114 nb_reaction_4 = nb_reaction_4 + 1
    2115 
    2116 indice_4(nb_reaction_4) = z4spec(1.0, i_so3, 1.0, i_co, 1.0, i_co2, 1.0,i_so2)
    2117 
    2118 !===========================================================
    2119 !      g076: HS + CO2 -> OCS + OH
     2087!      g072: HS + CO2 -> OCS + OH
    21202088!===========================================================
    21212089
    … …  
    21252093
    21262094!===========================================================
    2127 !      g077: S4 + H -> S3 + HS
     2095!      g073: S4 + H -> S3 + HS
    21282096!===========================================================
    21292097
    … …  
    21332101
    21342102!===========================================================
    2135 !      g078: HS + HS -> H2 + S2
     2103!      g074: HS + HS -> H2 + S2
    21362104!===========================================================
    21372105
    … …  
    30082976                       g057, g058, g059, g060, g061, g062, g063,   &
    30092977                       g064, g065, g066, g067, g068, g069, g070,   &
    3010                        g071, g072, g073, g074, g075, g076, g077,   &
    3011                        g078,                                       &
     2978                       g071, g072, g073, g074,                     &
    30122979                       h001, h002, h003,                           &
    30132980                       i001, i002, i003, i004, i005, i006,         &
    … …  
    44294396!  Atkinson et al 2004 Atmospheric Chemistry & Physics (KIDA) 200-430K
    44304397
    4431       g047(:) = 3.7E-11*exp(280.0/t(:))
     4398      g047(:) = 3.7E-11*exp(208.0/t(:))
    44324399
    44334400      nb_reaction_4 = nb_reaction_4 + 1
    … …  
    45814548      v_4(:,nb_reaction_4) = g062(:)
    45824549
    4583 !---  g063: SO3 + CO -> SO + CO2
     4550!---  g063: SO3 + CO -> SO2 + CO2
    45844551!     Krasnopolsky 2007 Icarus
    45854552
    … …  
    45894556      v_4(:,nb_reaction_4) = g063(:)     
    45904557
    4591 !---  g064: SO2 + CO -> CO2 + SO
     4558!---  g064: S2 + CO -> OCS + S
    45924559!     Krasnopolsky 2007 Icarus
    45934560
    4594       g064(:) = 4.5E-12*exp(-24300./t(:))
     4561      g064(:) = 1.E-12*exp(-17460/t(:))
    45954562
    45964563      nb_reaction_4 = nb_reaction_4 + 1
    45974564      v_4(:,nb_reaction_4) = g064(:)     
    45984565
    4599 !---  g065: S2 + CO -> OCS + S
     4566!---  g065: H + OCS -> CO + HS
    46004567!     Krasnopolsky 2007 Icarus
    46014568
    4602       g065(:) = 1.E-12*exp(-17460/t(:))
    4603 
    4604       nb_reaction_4 = nb_reaction_4 + 1
    4605       v_4(:,nb_reaction_4) = g065(:)     
    4606 
    4607 !---  g066: H + OCS -> CO + HS
     4569      g065(:) = 1.2E-11*exp(-1950/t(:))
     4570
     4571      nb_reaction_4 = nb_reaction_4 + 1
     4572      v_4(:,nb_reaction_4) = g065(:)
     4573
     4574!---  g066: S + SO2 -> SO + SO
    46084575!     Krasnopolsky 2007 Icarus
    46094576
    4610       g066(:) = 1.2E-11*exp(-1950/t(:))
    4611 
    4612       nb_reaction_4 = nb_reaction_4 + 1
    4613       v_4(:,nb_reaction_4) = g066(:)
    4614 
    4615 !---  g067: CO + HS -> H + OCS
    4616 !     Harada, N. et al 2010 ApJ (KIDA)
    4617 
    4618       g067(:) = 5.95E-14*(t(:)/300)**1.12 * exp(-8330/t(:))
     4577      g066(:) = 2.3E-11*exp(-5200/t(:))
     4578
     4579      nb_reaction_4 = nb_reaction_4 + 1
     4580      v_4(:,nb_reaction_4) = g066(:)     
     4581
     4582!---  g067: HS + CO -> OCS + H
     4583!     Kurbanov & Mamedov (1995) / Zhanle et al 2016 ApJ
     4584
     4585      g067(:) = 4.2e-14*exp(-7660/t(:))
    46194586
    46204587      nb_reaction_4 = nb_reaction_4 + 1
    46214588      v_4(:,nb_reaction_4) = g067(:)     
    46224589
    4623 !---  g068: S + SO2 -> SO + SO
     4590!---  g068: SO2 + H -> SO + OH
    46244591!     Krasnopolsky 2007 Icarus
    46254592
    4626       g068(:) = 2.3E-11*exp(-5200/t(:))
    4627 
    4628       nb_reaction_4 = nb_reaction_4 + 1
    4629       v_4(:,nb_reaction_4) = g068(:)     
    4630 
    4631 !---  g069: HS + CO -> OCS + H
    4632 !     Kurbanov & Mamedov (1995) / Zhanle et al 2016 ApJ
    4633 
    4634       g069(:) = 4.2e-14*exp(-7660/t(:))
    4635 
    4636       nb_reaction_4 = nb_reaction_4 + 1
    4637       v_4(:,nb_reaction_4) = g069(:)     
    4638 
    4639 !---  g070: SO2 + H -> SO + OH
    4640 !     Krasnopolsky 2007 Icarus
    4641 
    4642       g070(:) = 3.7e-9*exp(-14350/t(:))
    4643 
    4644       nb_reaction_4 = nb_reaction_4 + 1
    4645       v_4(:,nb_reaction_4) = g070(:)
    4646 
    4647 !---  g071: OCS + OH -> CO2 + HS
     4593      g068(:) = 3.7e-9*exp(-14350/t(:))
     4594
     4595      nb_reaction_4 = nb_reaction_4 + 1
     4596      v_4(:,nb_reaction_4) = g068(:)
     4597
     4598!---  g069: OCS + OH -> CO2 + HS
    46484599!     Atkinson 2004 / Zhanle et al 2016 ApJ
    46494600
    4650       g071(:) = 1.1e-13*exp(-1200/t(:))
    4651 
    4652       nb_reaction_4 = nb_reaction_4 + 1
    4653       v_4(:,nb_reaction_4) = g071(:)
    4654 
    4655 !---  g072: S + S4 -> S2 + S3
     4601      g069(:) = 1.1e-13*exp(-1200/t(:))
     4602
     4603      nb_reaction_4 = nb_reaction_4 + 1
     4604      v_4(:,nb_reaction_4) = g069(:)
     4605
     4606!---  g070: S + S4 -> S2 + S3
    46564607!     xMoses 2002 / Hobbs et al. 2020 MNRAS
    46574608!     Zhanle et al 2016 ApJ
    46584609
    4659       g072(:) = 4e-11*exp(-500/t(:))
     4610      g070(:) = 4e-11*exp(-500/t(:))
     4611
     4612      nb_reaction_4 = nb_reaction_4 + 1
     4613      v_4(:,nb_reaction_4) = g070(:)     
     4614
     4615!---  g071: H2 + HS -> H + H2S
     4616!     Krasnopolsky 2007
     4617
     4618!     Leen, T.M. et al 1988 ApJ / KIDA
     4619      g071(:) = 6.52E-12 * (t(:)/300)**(9.00E-2) * exp(-8.05E+3/t(:))
     4620
     4621      nb_reaction_4 = nb_reaction_4 + 1
     4622      v_4(:,nb_reaction_4) = g071(:)     
     4623
     4624!---  g072: HS + CO2 -> OCS + OH
     4625!     Krasnopolsky 2007
     4626
     4627      g072(:) = 2.6e-13*exp(-19360/t(:))
    46604628
    46614629      nb_reaction_4 = nb_reaction_4 + 1
    46624630      v_4(:,nb_reaction_4) = g072(:)     
    46634631
    4664 !---  g073: H2 + HS -> H + H2S
    4665 !     Krasnopolsky 2007
    4666 !      g061(:) = 3e-11*exp(-7930/t(:))
    4667 
    4668 !     Leen, T.M. et al 1988 ApJ / KIDA
    4669       g073(:) = 6.52E-12 * (t(:)/300)**(9.00E-2) * exp(-8.05E+3/t(:))
    4670 
    4671       nb_reaction_4 = nb_reaction_4 + 1
    4672       v_4(:,nb_reaction_4) = g073(:)     
    4673 
    4674 !---  g074: H + S3 -> HS + S2
    4675 !     Krasnopolsky 2007
    4676 
    4677       g074(:) = 1.2e-10*exp(-1950/t(:))
    4678 
    4679       nb_reaction_4 = nb_reaction_4 + 1
    4680       v_4(:,nb_reaction_4) = g074(:)     
    4681 
    4682 !---  g075: SO3 + CO -> CO2 + SO2
    4683 !     Krasnopolsky 2007
    4684 
    4685       g075(:) = 1.e-11*exp(-13000/t(:))
    4686 
    4687       nb_reaction_4 = nb_reaction_4 + 1
    4688       v_4(:,nb_reaction_4) = g075(:)     
    4689 
    4690 !---  g076: HS + CO2 -> OCS + OH
    4691 !     Krasnopolsky 2007
    4692 
    4693       g076(:) = 2.6e-13*exp(-19360/t(:))
    4694 
    4695       nb_reaction_4 = nb_reaction_4 + 1
    4696       v_4(:,nb_reaction_4) = g076(:)     
    4697 
    4698 !---  g077: S4 + H -> S3 + HS
     4632!---  g073: S4 + H -> S3 + HS
    46994633!     Zhanle et al 2016 ApJ
    47004634
    4701       g077(:) = 5.e-11*exp(-500/t(:))
    4702 
    4703       nb_reaction_4 = nb_reaction_4 + 1
    4704       v_4(:,nb_reaction_4) = g077(:)
    4705 
    4706 !      g078: HS + HS -> H2 + S2
     4635      g073(:) = 5.e-11*exp(-500/t(:))
     4636
     4637      nb_reaction_4 = nb_reaction_4 + 1
     4638      v_4(:,nb_reaction_4) = g073(:)
     4639
     4640!      g074: HS + HS -> H2 + S2
    47074641!      Zanle et al 2016 ApJ
    47084642
    4709       g078(:) = 3.e-11*exp(-20600/t(:))
     4643      g074(:) = 3.e-11*exp(-20600/t(:))
    47104644
    47114645      nb_reaction_3 = nb_reaction_3 + 1
    4712       v_3(:,nb_reaction_3) = g078(:)         
     4646      v_3(:,nb_reaction_3) = g074(:)         
    47134647
    47144648!----------------------------------------------------------------------

  • trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F

    r4108 r4118  
    157157           ! default numbers if no ion and no deuterium chemistry
    158158   
    159            nb_reaction_4_max = 128     ! set number of bimolecular reactions
     159           nb_reaction_4_max = 124     ! set number of bimolecular reactions
    160160           nb_reaction_3_max = 19     ! set number of quadratic reactions
    161161           nquench           = 16     ! set number of first-order reactions:
Note: See TracChangeset for help on using the changeset viewer.