Index: /trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 =================================================================== --- /trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 (revision 2842) +++ /trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 (revision 2843) @@ -2700,8 +2700,6 @@ !--- b003: o(1d) + h2 -> oh + h -! jpl 2003 -! b003(:) = 1.1E-10 - -! jpl 2016 +! jpl 2011 + b003(:) = 1.2E-10 @@ -2751,53 +2749,51 @@ !--- c002: o + oh -> o2 + h +! jpl 2011 + + c002(:) = 1.8E-11*exp(180./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = c002(:) + +!--- c003: h + o3 -> oh + o2 + ! jpl 2003 -! c002(:) = 2.2E-11*exp(120./t(:)) - -! jpl 2016 - c002(:) = 1.8E-11*exp(180./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = c002(:) - -!--- c003: h + o3 -> oh + o2 + + c003(:) = 1.4E-10*exp(-470./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = c003(:) + +!--- c004: h + ho2 -> oh + oh + +! jpl 2006 + + c004(:) = 7.2E-11 + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = c004(:) + +!--- c005: h + ho2 -> h2 + o2 + +! jpl 2006 + + c005(:) = 6.9E-12 + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = c005(:) + +!--- c006: h + ho2 -> h2o + o + +! jpl 2006 + + c006(:) = 1.6E-12 + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = c006(:) + +!--- c007: oh + ho2 -> h2o + o2 ! jpl 2003 - c003(:) = 1.4E-10*exp(-470./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = c003(:) - -!--- c004: h + ho2 -> oh + oh - -! jpl 2006 - - c004(:) = 7.2E-11 - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = c004(:) - -!--- c005: h + ho2 -> h2 + o2 - -! jpl 2006 - - c005(:) = 6.9E-12 - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = c005(:) - -!--- c006: h + ho2 -> h2o + o - -! jpl 2006 - - c006(:) = 1.6E-12 - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = c006(:) - -!--- c007: oh + ho2 -> h2o + o2 - -! jpl 2003 - c007(:) = 4.8E-11*exp(250./t(:)) @@ -2807,13 +2803,11 @@ !--- c008: ho2 + ho2 -> h2o2 + o2 -! jpl 2006 - -! c008(:) = 3.5E-13*exp(430./t(:)) +! jpl 2015 + + c008(:) = 3.0E-13*exp(460./t(:)) ! christensen et al., grl, 13, 2002 -! c008(:) = 1.5E-12*exp(19./t(:)) - -! jpl 2016 - c008(:) = 3.0E-13*exp(460./t(:)) + +! c008(:) = 1.5E-12*exp(19./t(:)) nb_reaction_3 = nb_reaction_3 + 1 @@ -2892,8 +2886,6 @@ !--- c016: ho2 + ho2 + co2 -> h2o2 + o2 + co2 -! jpl 2003 -! c016(:) = 2.5*1.7E-33*exp(1000./t(:))*conc(:) - -! jpl 2016 +! jpl 2011 + c016(:) = 2.5*2.1E-33*exp(920./t(:))*conc(:) @@ -3135,53 +3127,51 @@ !--- f004: hcl + oh -> h2o + cl +! jpl 2009 + + f004(:) = 1.8E-12*exp(-250./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = f004(:) + +!--- f005: clo + o -> cl + o2 + ! jpl 2006 -! f004(:) = 2.6E-12*exp(-350./t(:)) - -! jpl 2016 - f004(:) = 1.8E-12*exp(-250./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = f004(:) - -!--- f005: clo + o -> cl + o2 + + f005(:) = 2.8E-11*exp(85./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = f005(:) + +!--- f006: clo + oh -> cl + ho2 ! jpl 2006 - f005(:) = 2.8E-11*exp(85./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = f005(:) - -!--- f006: clo + oh -> cl + ho2 + f006(:) = 7.4E-12*exp(270./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = f006(:) + +!--- f007: clo + oh -> hcl + o2 ! jpl 2006 - f006(:) = 7.4E-12*exp(270./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = f006(:) - -!--- f007: clo + oh -> hcl + o2 + f007(:) = 6.0E-13*exp(230./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = f007(:) + +!--- f008: cl + h2 -> hcl + h ! jpl 2006 - f007(:) = 6.0E-13*exp(230./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = f007(:) - -!--- f008: cl + h2 -> hcl + h + f008(:) = 3.05E-11*exp(-2270./t(:)) + + nb_reaction_4 = nb_reaction_4 + 1 + v_4(:,nb_reaction_4) = f008(:) + +!--- f009: cl + o3 -> clo + o2 ! jpl 2006 - f008(:) = 3.05E-11*exp(-2270./t(:)) - - nb_reaction_4 = nb_reaction_4 + 1 - v_4(:,nb_reaction_4) = f008(:) - -!--- f009: cl + o3 -> clo + o2 - -! jpl 2006 - f009(:) = 2.3E-11*exp(-200./t(:)) @@ -3191,8 +3181,6 @@ !--- f010: cl + ho2 -> clo + oh -! jpl 2006 -! f010(:) = 4.1E-11*exp(-450./t(:)) - -! jpl 2016 +! jpl 2009 + f010(:) = 3.6E-11*exp(-375./t(:)) @@ -3202,8 +3190,6 @@ !--- f011: cl + ho2 -> hcl + o2 -! jpl 2006 -! f011(:) = 1.8E-11*exp(170./t(:)) - -! jpl 2016 +! jpl 2009 + f011(:) = 1.4E-11*exp(270./t(:)) @@ -3493,7 +3479,6 @@ !--- g001: s + o2 -> so + o -! g001(:) = 2.3E-12 - -! jpl 2016 +! jpl 2015 + g001(:) = 1.6E-12*exp(100./t(:)) @@ -3503,4 +3488,6 @@ !--- g002: s + o3 -> so + o2 +! jpl 2015 + g002(:) = 1.2E-11 @@ -3510,8 +3497,7 @@ !--- g003: so + o2 -> so2 + o -! g003(:) = 1.25E-13*exp(-2190./t(:)) - -! jpl 2016 - g003(:) = 1.6E-13*exp(-2280./t(:)) +! jpl 2015 + + g003(:) = 1.6E-13*exp(-2280./t(:)) nb_reaction_4 = nb_reaction_4 + 1 @@ -3520,4 +3506,6 @@ !--- g004: so + o3 -> so2 + o2 +! jpl 2015 + g004(:) = 3.4E-12*exp(-1100./t(:)) @@ -3527,4 +3515,6 @@ !--- g005: so + oh -> so2 + h +! jpl 2015 + g005(:) = 2.7E-11*exp(335./t(:)) @@ -3533,4 +3523,6 @@ !--- g006: s + oh -> so + h + +! jpl 2015 g006(:) = 6.6E-11