Changeset 408 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Nov 22, 2011, 2:26:02 PM (13 years ago)

Author:

emillour

Message:

Mars GCM:
Added FL corrections to aeronomars/photochemistry.F :

• Bug correction on the chemistry time step (now set to 1/3rd of the physics time step)
• Updated kinetics coefficients (according to JPL 2011)

EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F

@@ -86 +86 @@
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 c     ctimestep : chemistry timestep (s)                         c
-c                 10 minutes                                     c
-cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-c
-      ctimestep = 10.*60.
-c
-      phychemrat = int(ptimestep/ctimestep)
+c                 1/3 of physical timestep                       c
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+c
+      phychemrat = 3
+c
+      ctimestep = ptimestep/real(phychemrat)
 c
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -1521 +1521 @@
 ccc      b003: o(1d) + h2  -> oh + h
 c
-c        jpl 2003
-c
-         b003(l) = 1.1e-10
-c
-c        nair et al., 1994
-c
-c        b003(l) = 1.0e-10
+c        jpl 2011
+c
+         b003(l) = 1.2e-10
 c    
 ccc      b004: o(1d) + o2  -> o + o2
@@ -1581 +1577 @@
 ccc      c002: o + oh -> o2 + h
 c
-c        jpl 2003
+c        jpl 2011
 c 
-         c002(l) = 2.2e-11*exp(120./t(l))
+         c002(l) = 1.8e-11*exp(180./t(l))
 c
 c        robertson and smith, j. chem. phys. a 110, 6673, 2006
@@ -1667 +1663 @@
 ccc      c011: h + o2 + co2 -> ho2 + co2
 c
-c        jpl 2003
-c 
-c        ak0 = 2.5*5.7e-32*(t(l)/300.)**(-1.6)
-c        ak1 = 7.5e-11*(t(l)/300.)**(0.0)
-c
-c        jpl 2006
+c        jpl 2011
 c
          ak0 = 2.5*4.4e-32*(t(l)/300.)**(-1.3)
-         ak1 = 4.7e-11*(t(l)/300.)**(-0.2)
+         ak1 = 7.5e-11*(t(l)/300.)**(0.2)
 c
          rate = (ak0*dens(l))/(1. + ak0*dens(l)/ak1)
@@ -1727 +1718 @@
 ccc      c016: ho2 + ho2 + co2 -> h2o2 + o2 + co2
 c
-c        jpl 2003
-c
-         c016(l) = 2.5*1.7e-33
-     $              *exp(1000./t(l))*dens(l)
+c        jpl 2011
+c
+         c016(l) = 2.5*2.1e-33
+     $              *exp(920./t(l))*dens(l)
 c
 ccc      c017: oh + oh + co2 -> h2o2 + co2
@@ -1785 +1776 @@
 ccc      d003: no + ho2 -> no2 + oh
 c
-c        jpl 2003
-c
-         d003(l) = 3.5e-12*exp(250./t(l))
+c        jpl 2011
+c
+         d003(l) = 3.3e-12*exp(270./t(l))
 c
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -2006 +1997 @@
       return
       end
+