Changeset 1013


Ignore:
Timestamp:
Jul 26, 2013, 12:54:55 PM (11 years ago)
Author:
emillour
Message:

Mars GCM:

  • Upgrade on the thermospheric photochemical reaction rates. These are now read in from a file "chemthermos_reactionrates.def".

FGG + JYC

Location:
trunk/LMDZ.MARS
Files:
2 added
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/README

    r1009 r1013  
    18851885- Converted run0 and run_mcd scripts to bash.
    18861886
     1887== 26/07/2013 == FGG +JYC
     1888- Upgrade on the thermospheric photochemical reaction rates. These are now
     1889  read in from a file "chemthermos_reactionrates.def".

  • trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

    r710 r1013  
    876876        IMPLICIT NONE
    877877
    878         INTEGER :: nl
    879         INTEGER :: l,ig
    880         REAL,DIMENSION(nl) :: T1,DT1,DT2,DT3
    881         REAL*8,DIMENSION(nl) :: T2
    882         REAL :: dtime
     878        INTEGER,INTENT(IN) :: nl,ig
     879        REAL,INTENT(IN),DIMENSION(nl) :: T1,DT1,DT2,DT3
     880        REAL*8,INTENT(OUT),DIMENSION(nl) :: T2
     881        REAL,INTENT(IN) :: dtime
     882        INTEGER :: l
    883883        DO l=1,nl
    884884        T2(l)=T1(l)*1D0+(DT1(l)*dble(dtime)+                            &
    … …  
    897897#include "dimensions.h"
    898898
    899         INTEGER :: nl,nq
    900         INTEGER :: l,iq,ig
    901         INTEGER,dimension(nq) :: gc
    902         REAL,DIMENSION(nl,nqmx) :: Q1,DQ
    903         REAL*8,DIMENSION(nl,nq) :: Q2
    904         REAL :: dtime
     899        INTEGER,INTENT(IN) :: nl,nq
     900        INTEGER,INTENT(IN) :: ig
     901        INTEGER,INTENT(IN),dimension(nq) :: gc
     902        REAL,INTENT(IN),DIMENSION(nl,nqmx) :: Q1,DQ
     903        REAL*8,INTENT(OUT),DIMENSION(nl,nq) :: Q2
     904        REAL,INTENT(IN) :: dtime
     905        INTEGER :: l,iq
    905906        DO l=1,nl
    906907        DO iq=1,nq

  • trunk/LMDZ.MARS/libf/aeronomars/param_v4.h

    r705 r1013  
    4444     &                     ch57,ch58,ch59,ch62,ch63,ch64,ch65,                   &
    4545     &                     ch66,ch67,ch68,ch69,ch70,ch71,ch72,                   &
    46      &                     ch73,ch74,ch75,ch76,ch85,ch86,ch87,                   &
    47      &                     jdistot,jdistot_b,jion
     46     &                     ch73,ch74,ch75,ch76,ch85,ch86,ch87,rcoef,             &
     47     &                     jdistot,jdistot_b,jion 
    4848
    4949      real*8 ch2, ch3, ch4, ch5, ch7,ch9,ch10,ch11,ch13,ch14,ch15,ch18
    … …  
    5353      real*8 ch63,ch64,ch65,ch66,ch67,ch68,ch69,ch70,ch71
    5454      real*8 ch72,ch73,ch74,ch75,ch76,ch85,ch86,ch87
    55 
     55      real*8 rcoef(61,3)
    5656                                                  !reaction rates
    5757      real jdistot(nabs,nlayermx)    !photodissociation rates
    … …  
    5959     
    6060      real jion(nabs,nlayermx,4)
     61
    6162
    6263

  • trunk/LMDZ.MARS/libf/aeronomars/paramfoto_compact.F

    r690 r1013  
    2727c     local variables
    2828
    29       integer ij0,l0
    3029      real*8  deltat,timefrac_sec
    3130      real*8  co2xnew,o2xnew,o3pxnew,coxnew,hxnew,ohxnew
    … …  
    121120c**********************************************************************
    122121
    123 
    124         ij0=302
    125         l0=25
    126122
    127123C     Start: altitude loop
    … …  
    180176         endif
    181177
    182 
    183178         !Calculation of photodissociation and photoionization rates
    184179         !from photoabsorption rates and ionization-to-dissociation
    … …  
    193188            enddo
    194189         end do
     190         
    195191
    196192         !Reaction rates
    197          call getch( ig, chemthermod,tx8, co2xini, zx(i))
     193         call getch( ig, chemthermod,tx8, zx(i))
    198194                 
    199195         !Lifetimes and temporal integration
    … …  
    213209! JYC criteria added to avoid instabilities in (H) + (O+) <-> (H+) + (O) reactions when H+ is important
    214210
    215         fmargin1=1
     211        fmargin1=5
    216212        !Only if ion chemistry requested
    217213        if(chemthermod.eq.3) then
    … …  
    231227         numpasos = int( timefrac_sec / deltat )
    232228         alfa_laststep = 1.d0 + timefrac_sec/deltat - dble(numpasos)
    233 
    234229         do paso=1,numpasos 
    235230
    … …  
    603598
    604599         end do   
     600         
    605601         !New concentrations to be returned
    606602         rm(i,i_co2)     = real(co2xnew)
    … …  
    642638      end do     
    643639cccccc End altitude loop
    644      
    645640
    646641      return
    … …  
    12721267     $        jfotsout(inter,1,i)*fluxtop(inter)*efionco2(inter,4)
    12731268
     1269
    12741270c     O2
    12751271         jdis(2,inter,i) = jfotsout(inter,2,i) * fluxtop(inter)
    … …  
    12851281         jion(2,i,1)=jion(2,i,1) +
    12861282     $        jfotsout(inter,2,i) * fluxtop(inter) * (1.-efdiso2(inter))
     1283!         if(ig.eq.1.and.i.eq.39)write(*,*)'paramfoto_compact/1355',
     1284!     $        jion(2,i,1),jfotsout(inter,2,i),fluxtop(inter)
    12871285         
    12881286c     O3P
    … …  
    13811379c**********************************************************************
    13821380c***************************************************************************
    1383         subroutine getch (ig,chemthermod,tt,concco2,zkm)
     1381        subroutine getch (ig,chemthermod,tt,zkm)
    13841382
    13851383 
    … …  
    13981396
    13991397c Arguments     
    1400         integer           ig,chemthermod
    1401         real*8          tt      ! Temperature
    1402         real*8            concco2
    1403         real              zkm   !  Altitude in km
     1398        integer       ig,chemthermod
     1399        real*8        tt         ! Temperature
     1400        real          zkm       !  Altitude in km
    14041401
    14051402c local variables:
    1406         double precision        tcte   
    1407         double precision        t_elect ! electronic temperatures
    1408         double precision        val ! valores de alturas corresp a t_elect
    1409         double precision        zhanson(9),tehanson(9)
    1410         double precision        incremento
    1411         integer                 ii, i1, i2
    1412         real                    phelec_fact_ch86
     1403        real*8        tcte   
     1404        real*8        t_elect ! electronic temperatures
     1405        real*8        val ! valores de alturas corresp a t_elect
     1406        real*8        zhanson(9),tehanson(9)
     1407        real*8        incremento
     1408        integer       ii, i1, i2
    14131409       
    14141410c**************************************************************************
    … …  
    15131509        ch85 = 0.d0
    15141510        ch86 = 0.d0
    1515        
     1511
    15161512
    15171513        !Reaction rates
    1518 ch2: h + o2 + co2 --> ho2 + co2
    1519         ch2 = 1.425d-31 * (tcte / 300.)**(-1.6d0)
    1520 
    1521 ch3: o + ho2 --> oh + o2
    1522         ch3 = 3.0d-11 * exp(200.d0 / tcte)
     1514!ch2: h + o2 + co2 --> ho2 + co2
     1515        ! JPL 2003 (low pressure limit)*2.5
     1516!       ch2 = 1.425d-31 * (tcte / 300.)**(-1.6d0)
     1517        ! JPL 2011 (low pressure limit)*2.5
     1518!        ch2 = 1.1e-31 * (tcte / 300.)**(-1.3)
     1519        ch2=rcoef(1,1)*((tcte/300.)**rcoef(1,2))*exp(rcoef(1,3)/tcte)
     1520
     1521!ch3: o + ho2 --> oh + o2
     1522        ! JPL 2011:
     1523!       ch3 = 3.0d-11 * exp(200.d0 / tcte)
     1524        ch3=rcoef(2,1)*((tcte/300.)**rcoef(2,2))*exp(rcoef(2,3)/tcte)
    15231525
    15241526ch4: co + oh --> co2 + h
    1525 !       ch4 = 3.2d-13 * exp(-300.d0 / tcte)
    1526         ch4 = 1.57d-13 + 3.54d-33*concco2
     1527        !Nair et al, 1994:
     1528        !ch4 = 3.2d-13 * exp(-300.d0 / tcte)
     1529        !mccabe et al., grl, 28, 3135, 2001
     1530        !ch4 = 1.57d-13 + 3.54d-33*concco2
     1531        !JPL 2011 (low pressure limit):
     1532!        ch4 = 1.5d-13 * (tcte/300.)**0.6
     1533        ch4=rcoef(3,1)*((tcte/300.)**rcoef(3,2))*exp(rcoef(3,3)/tcte)
    15271534
    15281535ch5: ho2 + ho2 --> h2o2 + o2
    1529         ch5 = 2.3d-13 * exp(600.d0 / tcte)
     1536        !JPL 2003:
     1537        !ch5 = 2.3d-13 * exp(600.d0 / tcte)
     1538        !JPL 2011:
     1539!        ch5 = 3.0d-13 * exp(460.d0 / tcte)
     1540        ch5=rcoef(4,1)*((tcte/300.)**rcoef(4,2))*exp(rcoef(4,3)/tcte)
    15301541
    15311542ch7: oh + ho2 --> h2o + o2
    1532         ch7 = 4.8d-11 * exp(250.d0 / tcte)
     1543        !JPL 2011:
     1544!       ch7 = 4.8d-11 * exp(250.d0 / tcte)
     1545        ch7=rcoef(5,1)*((tcte/300.)**rcoef(5,2))*exp(rcoef(5,3)/tcte)
    15331546
    15341547ch9: o(1d) + h2o --> 2oh
    1535         ch9 = 2.2d-10
     1548        !JPL 2003:
     1549        !ch9 = 2.2d-10
     1550        !JPL 2011:
     1551!        ch9 = 1.63d-10 * exp(60.d0 / tcte)
     1552        ch9=rcoef(6,1)*((tcte/300.)**rcoef(6,2))*exp(rcoef(6,3)/tcte)
    15361553
    15371554ch10: o + o + co2 --> o2 + co2
     1555        !JPL 1990:
    15381556!       ch10 = 1.1d-27 * (tcte **(-2.0d0)) !Estandard en el 1-D
    1539         ch10 = 1.3d-34 * exp(900.d0/tcte)
     1557        !Tsang and Hampson, 1986:
     1558!       ch10 = 1.3d-34 * exp(900.d0/tcte)
     1559        ch10=rcoef(7,1)*((tcte/300.)**rcoef(7,2))*exp(rcoef(7,3)/tcte)
    15401560
    15411561ch11: o + oh --> o2 + h
    1542         ch11 = 2.2d-11 * exp(120.d0 / tcte)
     1562        !JPL 2003:
     1563        !ch11 = 2.2d-11 * exp(120.d0 / tcte)
     1564        !JPL 2011:
     1565!        ch11 = 1.8d-11 * exp(180.d0 /tcte)
     1566        ch11=rcoef(8,1)*((tcte/300.)**rcoef(8,2))*exp(rcoef(8,3)/tcte)
    15431567
    15441568ch13: h + ho2 --> h2 + o2
    1545         ch13 = 6.5d-12
     1569        !JPL 2003:
     1570        !ch13 = 6.5d-12
     1571        !JPL 2011:
     1572!        ch13 = 6.9d-12
     1573        ch13=rcoef(9,1)*((tcte/300.)**rcoef(9,2))*exp(rcoef(9,3)/tcte)
    15461574
    15471575ch14: o(1d) + h2 --> h + oh
    1548         ch14 = 1.1d-10
     1576        !JPL 2003:
     1577        !ch14 = 1.1d-10
     1578        !JPL 2011:
     1579!        ch14 = 1.2d-10
     1580        ch14=rcoef(10,1)*((tcte/300.)**rcoef(10,2))*
     1581     $       exp(rcoef(10,3)/tcte)
    15491582
    15501583ch15: oh + h2 --> h + h2o
    1551         ch15 = 5.5d-12 * exp (-2000.d0 / tcte)
     1584        !JPL 2003:
     1585        !ch15 = 5.5d-12 * exp (-2000.d0 / tcte)
     1586        !JPL 2011:
     1587!        ch15 = 2.8d-12 * exp (-1800.d0 / tcte)
     1588        ch15=rcoef(11,1)*((tcte/300.)**rcoef(11,2))*
     1589     $       exp(rcoef(11,3)/tcte)
    15521590
    15531591ch18: oh + h2o2 --> h2o + ho2
    1554         ch18 = 2.9d-12 * exp (-160.d0 / tcte)
     1592        !JPL 2003:
     1593        !ch18 = 2.9d-12 * exp (-160.d0 / tcte)
     1594        !JPL 2011:
     1595!        ch18 = 1.8d-12
     1596        ch18=rcoef(12,1)*((tcte/300.)**rcoef(12,2))*
     1597     $       exp(rcoef(12,3)/tcte)
    15551598
    15561599ch19: o(1d) + co2 --> o + co2
    1557         ch19 = 7.4d-11 * exp(120.d0 / tcte)
     1600        !JPL 2003:
     1601        !ch19 = 7.4d-11 * exp(120.d0 / tcte)
     1602        !JPL 2011:
     1603!        ch19 = 7.5d-11 * exp(115.d0 / tcte)       
     1604        ch19=rcoef(13,1)*((tcte/300.)**rcoef(13,2))*
     1605     $       exp(rcoef(13,3)/tcte)
    15581606
    15591607ch20: o(1d) + o2 --> o + o2
    1560         ch20 = 3.2d-11 * exp (70.d0 / tcte)
     1608        !JPL 2003:
     1609        !ch20 = 3.2d-11 * exp (70.d0 / tcte)
     1610        !JPL 2011:
     1611!        ch20 = 3.3d-11 * exp(55.d0 / tcte)
     1612        ch20=rcoef(14,1)*((tcte/300.)**rcoef(14,2))*
     1613     $       exp(rcoef(14,3)/tcte)
    15611614
    15621615ch21: o + o2 + co2 --> o3 + co2
    1563         ch21 = 1.5d-33 * ((tcte / 300.d0) ** (-2.4d0))
     1616        !JPL 2011 * 2.5:
     1617!       ch21 = 1.5d-33 * ((tcte / 300.d0) ** (-2.4d0))
     1618        ch21=rcoef(15,1)*((tcte/300.)**rcoef(15,2))*
     1619     $       exp(rcoef(15,3)/tcte)
    15641620       
    15651621
    … …  
    15681624
    15691625ch22: o3 + h --> o2 + oh
    1570            ch22 = 1.4d-10 * exp (-470.d0 / tcte)
     1626           !JPL 2011:
     1627!          ch22 = 1.4d-10 * exp (-470.d0 / tcte)
     1628           ch22=rcoef(16,1)*((tcte/300.)**rcoef(16,2))*
     1629     $       exp(rcoef(16,3)/tcte)
    15711630
    15721631ch23: o3 + oh --> ho2 + o2
    1573            ch23 = 1.7d-12 * exp (-940.d0 / tcte)
     1632           !JPL 2011:
     1633!          ch23 = 1.7d-12 * exp (-940.d0 / tcte)
     1634           ch23=rcoef(17,1)*((tcte/300.)**rcoef(17,2))*
     1635     $       exp(rcoef(17,3)/tcte)
    15741636
    15751637ch24: o3 + ho2 --> oh + 2o2
    1576            ch24 = 1.0d-14 * exp (-490.d0 / tcte)
     1638           !JPL 2011:
     1639!          ch24 = 1.0d-14 * exp (-490.d0 / tcte)
     1640           ch24=rcoef(18,1)*((tcte/300.)**rcoef(18,2))*
     1641     $       exp(rcoef(18,3)/tcte)
    15771642
    15781643        endif
    … …  
    15831648c N chemistry
    15841649ch30: n + no --> n2 + o
    1585            ch30 = 2.1d-11 * exp (100.d0 / tcte)
     1650           !JPL 2011:
     1651!          ch30 = 2.1d-11 * exp (100.d0 / tcte)
     1652           ch30=rcoef(19,1)*((tcte/300.)**rcoef(19,2))*
     1653     $       exp(rcoef(19,3)/tcte)
    15861654
    15871655ch31: n2 + o(1d) --> n2 + o
    1588            ch31 = 2.5d-11 * exp (110.d0 / tcte)
     1656           !JPL 2011:
     1657!          ch31 = 2.15d-11 * exp (110.d0 / tcte)
     1658           ch31=rcoef(20,1)*((tcte/300.)**rcoef(20,2))*
     1659     $       exp(rcoef(20,3)/tcte)
    15891660
    15901661ch32: n + o2 --> no + o
    1591            ch32 = 1.5d-11 * exp (-3600.d0 / tcte)
     1662           !JPL 2011:
     1663!          ch32 = 1.5d-11 * exp (-3600.d0 / tcte)
     1664           ch32=rcoef(21,1)*((tcte/300.)**rcoef(21,2))*
     1665     $       exp(rcoef(21,3)/tcte)
    15921666
    15931667ch33: n + oh --> no + h
    1594            ch33 = 3.8d-11 * exp (85.d0 / tcte)
     1668           !Atkinson et al., 1989 (usado en Nair et al., 1994)
     1669!          ch33 = 3.8d-11 * exp (85.d0 / tcte)
     1670           ch33=rcoef(22,1)*((tcte/300.)**rcoef(22,2))*
     1671     $       exp(rcoef(22,3)/tcte)
    15951672
    15961673ch34: n + o3 --> no + o2
    1597            ch34 = 1.0d-16
     1674           !JPL 2011 (it is an upper limit):
     1675!          ch34 = 1.0d-16
     1676           ch34=rcoef(23,1)*((tcte/300.)**rcoef(23,2))*
     1677     $       exp(rcoef(23,3)/tcte)
    15981678
    15991679ch35: n + ho2 --> no + oh
    1600            ch35 = 2.2d-11
    1601 
    1602 ch36: n(2d) + n2 --> n + n2
    1603            ch36 = 6.9d-13
     1680           !Brune et al., 1983 (from Nair et al., 1994)
     1681!          ch35 = 2.2d-11
     1682           ch35=rcoef(24,1)*((tcte/300.)**rcoef(24,2))*
     1683     $       exp(rcoef(24,3)/tcte)
     1684
     1685ch36: n(2d) + o --> n + o
     1686           !Fell et al., 1990 (from Nair et al., 1994)
     1687           !ch36 = 6.9d-13
     1688           !Herron, 1999:
     1689!           ch36 = 3.3d-12 * exp(-260.d0 / tcte)
     1690           ch36=rcoef(25,1)*((tcte/300.)**rcoef(25,2))*
     1691     $       exp(rcoef(25,3)/tcte)
    16041692
    16051693ch37: n(2d) + n2 --> n + n2
    1606            ch37 = 1.7d-14
     1694           !Herron, 1999:
     1695           !Coincides with Nair et al., 1994:
     1696!          ch37 = 1.7d-14
     1697           ch37=rcoef(26,1)*((tcte/300.)**rcoef(26,2))*
     1698     $       exp(rcoef(26,3)/tcte)
    16071699
    16081700ch38: n(2d) + co2 --> no + co
    1609            ch38 = 3.5d-13
     1701           !Pipper et al., 1989 (from Nair et al., 1994):
     1702           !ch38 = 3.5d-13
     1703           !Herron, 1999:
     1704!           ch38 = 3.6d-13
     1705           ch38=rcoef(27,1)*((tcte/300.)**rcoef(27,2))*
     1706     $       exp(rcoef(27,3)/tcte)
    16101707
    16111708ch39: no + ho2 --> no2+oh
    1612            ch39 = 3.5d-12 * exp (250.d0 / tcte)
     1709           !JPL 2006:
     1710           !ch39 = 3.5d-12 * exp (250.d0 / tcte)
     1711           !JPL 2011:
     1712!           ch39 = 3.3d-12 * exp(270.d0 / tcte)
     1713           ch39=rcoef(28,1)*((tcte/300.)**rcoef(28,2))*
     1714     $       exp(rcoef(28,3)/tcte)
    16131715
    16141716ch40: o + no + co2 --> no2 + co2
    1615            ch40 = 2.5d0 * 9.0d-32 * ((tcte / 300.d0) ** (-1.5d0))
     1717           !JPL 2011 * 2.5 (low pressure limit)
     1718!          ch40 = 2.5d0 * 9.0d-32 * ((tcte / 300.d0) ** (-1.5d0))
     1719           ch40=rcoef(29,1)*((tcte/300.)**rcoef(29,2))*
     1720     $       exp(rcoef(29,3)/tcte)
    16161721
    16171722ch41: o + no2 --> no + o2
    1618            ch41 = 5.1d-12 * exp (210.d0 / tcte)
     1723           !JPL 2011:
     1724!          ch41 = 5.1d-12 * exp (210.d0 / tcte)
     1725           ch41=rcoef(30,1)*((tcte/300.)**rcoef(30,2))*
     1726     $       exp(rcoef(30,3)/tcte)
    16191727
    16201728ch42: no + o3 --> no2 + o2
    1621            ch42 = 3.0d-12 * exp (-1500.d0 / tcte)
     1729           !JPL 2011:
     1730!          ch42 = 3.0d-12 * exp (-1500.d0 / tcte)
     1731           ch42=rcoef(31,1)*((tcte/300.)**rcoef(31,2))*
     1732     $       exp(rcoef(31,3)/tcte)
    16221733
    16231734ch43: h + no2 --> no + oh
    1624            ch43 = 4.0d-10 * exp (-340.d0 / tcte)
     1735           !JPL 2011:
     1736!          ch43 = 4.0d-10 * exp (-340.d0 / tcte)
     1737           ch43=rcoef(32,1)*((tcte/300.)**rcoef(32,2))*
     1738     $       exp(rcoef(32,3)/tcte)
    16251739
    16261740ch45: n + o --> no
    1627            ch45 = 2.8d-17 * ((300.d0 / tcte) ** 0.5)
     1741           !Du and Dalgarno, 1990
     1742!          ch45 = 2.8d-17 * ((300.d0 / tcte) ** 0.5)
     1743           ch45=rcoef(33,1)*((tcte/300.)**rcoef(33,2))*
     1744     $       exp(rcoef(33,3)/tcte)
    16281745           
    16291746        endif    !of if(chemthermod.ge.2)
    … …  
    16351752c Ionosphere
    16361753c
    1637 ch46:  co2+ +  O --> O2+ + CO
    1638            ch46 = 5.0d-11
     1754ch46:  co2+ +  O2 --> O2+ + CO2
     1755           !Moffat et al., 2005 (given by GG):
     1756           !ch46 = 5.0d-11
     1757           !Copp et al., 1982:
     1758           !ch46 = 5.3d-11
     1759           !Aninich 1993 (from Fox and Sung 2001):
     1760!           ch46 = 5.5d-11 * (300.d0/t_elect)**0.82
     1761           ch46=rcoef(34,1)*((t_elect/300.)**rcoef(34,2))*
     1762     $       exp(rcoef(34,3)/t_elect)
     1763           
    16391764
    16401765ch47: CO2+ + O -->  O+  +  CO2
    1641            ch47 = 1.64d-10
     1766           !Original (incorrect) value (corresponds to ch48):
     1767           !ch47 = 1.64d-10
     1768           !Fehsenfeld et al., 1970 (from UMIST,
     1769           !Fox and Sung 2001, Krasnopolsky 2002):
     1770!           ch47 = 9.6d-11
     1771           ch47=rcoef(35,1)*((t_elect/300.)**rcoef(35,2))*
     1772     $       exp(rcoef(35,3)/t_elect)
    16421773
    16431774ch48: CO2+ + O --> O2+  + CO
    1644            ch48 = 9.6d-11
     1775           !Original (incorrect) value (corresponds to ch47):
     1776           !ch48 = 9.6d-11
     1777           !Fehsenfeld et al., 1970 (from UMIST,
     1778           !Fox and Sung 2001, Krasnopolsky 2002):
     1779!           ch48 = 1.64d-10
     1780           ch48=rcoef(36,1)*((t_elect/300.)**rcoef(36,2))*
     1781     $       exp(rcoef(36,3)/t_elect)
    16451782
    16461783ch49: O2+ + elect --> O + O
    1647            ch49 = 2.0d-7*(300.d0/t_elect)**(0.7d0)
     1784           !Alge et al., 1983:
     1785           !Here we do not divide into reaction producing different
     1786           !O atomic states. O + O(1d) seems to be the dominant products
     1787           !(see Fox and Sung 2002). We should consider dividing
     1788           !into two different reactions
     1789!          ch49 = 2.0d-7*(300.d0/t_elect)**(0.7d0)
     1790           ch49=rcoef(37,1)*((t_elect/300.)**rcoef(37,2))*
     1791     $       exp(rcoef(37,3)/t_elect)
    16481792
    16491793ch50: O+  + CO2 --> O2+  + CO
    1650            ch50 = 9.4d-10
     1794           !Adams et al., 1980 (from UMIST):
     1795!          ch50 = 9.4d-10
     1796           !Anicich 1993 (from Fox and Sung 2001):
     1797           !ch50 = 1.1d-9
     1798           ch50=rcoef(38,1)*((t_elect/300.)**rcoef(38,2))*
     1799     $       exp(rcoef(38,3)/t_elect)
    16511800
    16521801ch55: CO2+ +  e ---->   CO  +   O
    1653            ch55 = 3.8d-7*(300.d0/t_elect)**(0.5d0)
     1802           !Mitchell, 1990 (from UMIST):
     1803!          ch55 = 3.8d-7*(300.d0/t_elect)**(0.5d0)
     1804           !Gougousi et al., 1997 (from Fox and Sung 2001):
     1805           !ch55 = 3.5d-7*(300.d0/t_elect)**0.5d0
     1806           ch55=rcoef(39,1)*((t_elect/300.)**rcoef(39,2))*
     1807     $       exp(rcoef(39,3)/t_elect)
    16541808       
    16551809ch56: O+ + CO2  --->  O2   +   CO+
    1656            ch56 = 9.4d-10
     1810           !Original, Kim et al., 1989:
     1811           !ch56 = 9.4d-10
     1812           !It does not appear in any other paper. Its presence in
     1813           !Kim et al., 1989 is probably a confusion with ch50.
     1814!           ch56 = 0.d0
     1815           ch56=rcoef(40,1)*((t_elect/300.)**rcoef(40,2))*
     1816     $       exp(rcoef(40,3)/t_elect)
    16571817
    16581818ch57: CO+ + CO2  ---> CO2+ + CO
    1659            ch57 = 1.0d-9
     1819           !Adams et al., 1978 (from UMIST):
     1820!          ch57 = 1.0d-9
     1821           !Anicich 1993 (from Fox and Sung 2001):
     1822           !ch57 = 1.1d-9
     1823           ch57=rcoef(41,1)*((t_elect/300.)**rcoef(41,2))*
     1824     $       exp(rcoef(41,3)/t_elect)
     1825
    16601826
    16611827ch58: CO+ + O   --->   O+  + CO
    1662            ch58 = 1.4d-10
     1828           !Fenhsenfeld et al. 1970 (from UMIST, F&S2001, K2002):
     1829!          ch58 = 1.4d-10
     1830           ch58=rcoef(42,1)*((t_elect/300.)**rcoef(42,2))*
     1831     $       exp(rcoef(42,3)/t_elect)
    16631832 
    16641833ch59: C+ +  CO2 --->  CO+  + CO    !!!!  NEW  !!!
    1665            ch59 = 1.1d-9
     1834           !Fahey et al., 1981 (from UMIST, F&S2001, K2002):
     1835!          ch59 = 1.1d-9
     1836           ch59=rcoef(43,1)*((t_elect/300.)**rcoef(43,2))*
     1837     $       exp(rcoef(43,3)/t_elect)
    16661838
    16671839ch62:  CO2+  +  NO   -->  NO+  +  CO2
    1668            ch62 = 1.2d-10
     1840           !Copp et al., 1982 (from UMIST):
     1841!          ch62 = 1.2d-10
     1842           !Anicich 1993 (from Fox and Sung 2001):
     1843           !ch62 = 1.23d-10
     1844           ch62=rcoef(44,1)*((t_elect/300.)**rcoef(44,2))*
     1845     $       exp(rcoef(44,3)/t_elect)
    16691846
    16701847ch63:  CO2+  +  N    -->  NO  + CO+
    1671            ch63 = 1.0d-11
     1848           !Kim et al., 1989:
     1849           !ch63 = 1.0d-11
     1850           !Scott et al., 1998 (from Fox and Sung 2001):
     1851!           ch63 = 3.4d-10
     1852           ch63=rcoef(45,1)*((t_elect/300.)**rcoef(45,2))*
     1853     $       exp(rcoef(45,3)/t_elect)
    16721854
    16731855ch64:   O2+  +  NO    -->  NO+  + O2
    1674            ch64 = 4.5d-10
     1856           !Middey and Vigiano 1999 (from Fox and Sung 2001):
     1857!          ch64 = 4.5d-10
     1858           !Aninich 1993 (from UMIST):
     1859           !ch64 = 4.6d-10
     1860           ch64=rcoef(46,1)*((t_elect/300.)**rcoef(46,2))*
     1861     $       exp(rcoef(46,3)/t_elect)
    16751862
    16761863ch65:  O2+  +  N2    -->  NO+  + NO
    1677            ch65 = 1.0d-16
     1864           !Original from GG, Moffat 2005:
     1865           !ch65 = 1.0d-16
     1866           !Ferguson 1973 (from Fox and Sung 2001):
     1867!           ch65 = 1.0d-15
     1868           ch65=rcoef(47,1)*((t_elect/300.)**rcoef(47,2))*
     1869     $       exp(rcoef(47,3)/t_elect)
    16781870
    16791871ch66:  O2+  +  N    -->  NO+  + O
    1680            ch66 = 1.2d-10
     1872           !Kim et al., 1989:
     1873           !ch66 = 1.2d-10
     1874           !Scott et al., 1998 (from Fox and Sung 2001):
     1875!           ch66 = 1.0d-10
     1876           !Goldan et al., 1966 (from UMIST):
     1877           !ch66 = 1.8d-10
     1878           ch66=rcoef(48,1)*((t_elect/300.)**rcoef(48,2))*
     1879     $       exp(rcoef(48,3)/t_elect)
    16811880
    16821881
    16831882ch67:   O+  +  N2    -->  NO+  + N
    1684            ch67 = 1.2d-12 * (300.d0/t_elect)**(0.41d0)
     1883           !Moffat 2005:
     1884           !ch67 = 1.2d-12 * (300.d0/t_elect)**(0.41d0)
     1885           !Hierl et al. 1997 (from Fox and Sung 2001):
     1886!           ch67 = 1.2d-12 * (300.d0/t_elect)**0.45d0
     1887           !Adams et al., 1980 (from UMIST):
     1888           !ch67=2.42d-12 * (300.d0/t_elec)**(-0.21)*exp(44./t_elec)
     1889           ch67=rcoef(49,1)*((t_elect/300.)**rcoef(49,2))*
     1890     $       exp(rcoef(49,3)/t_elect)
    16851891
    16861892ch68:  N2+  +  CO2    -->  CO2+  + N2
    1687            ch68 = 7.7d-10
     1893           !Adams et al. 1980 (from UMIST):
     1894           !ch68 = 7.7d-10
     1895           !Dotan et al. 2000 (from F&S2001):
     1896!           ch68 = 9.0d-10 * (300./t_elect)**0.23
     1897           ch68=rcoef(50,1)*((t_elect/300.)**rcoef(50,2))*
     1898     $       exp(rcoef(50,3)/t_elect)
    16881899
    16891900ch69:   N2+  +  O3p    -->  NO+  + N
    1690            ch69 = 1.3d-10
     1901           !McFarland et al., 1974 (from UMIST):
     1902           !ch69 = 1.3d-10
     1903           !Scott et al. 1999 (from F&S2001):
     1904!           ch69 = 1.33d-10 * (300./t_elect)**0.44
     1905           ch69=rcoef(51,1)*((t_elect/300.)**rcoef(51,2))*
     1906     $       exp(rcoef(51,3)/t_elect)
    16911907
    16921908ch70:  N2+  +  CO    -->  N2  + CO+
    1693            ch70 = 7.4d-11
     1909           !Adams et al., 1980 (from UMIST):
     1910!          ch70 = 7.4d-11
     1911           !Frost et al., 1998 (from F&S2001):
     1912           !ch70 = 7.6d-11
     1913           ch70=rcoef(52,1)*((t_elect/300.)**rcoef(52,2))*
     1914     $       exp(rcoef(52,3)/t_elect)
    16941915
    16951916ch71:  N2+  +  e-    -->  N  + N
    1696            ch71 = 3.5d-7 * (300.d0/t_elect)**(0.5d0)
     1917           !Moffat 2005
     1918           !ch71 = 3.5d-7 * (300.d0/t_elect)**(0.5d0)
     1919           !Peterson et al. 1998 (from UMIST):
     1920!           ch71 = 1.7d-7 * (300.d0/t_elect)**0.3
     1921           !Zipf 1980+Kella et al 1996 (from F&S2001):
     1922           !ch71 = 2.2d-7 * (300.d0/t_elect)**0.39
     1923           ch71=rcoef(53,1)*((t_elect/300.)**rcoef(53,2))*
     1924     $       exp(rcoef(53,3)/t_elect)
    16971925
    16981926
    16991927ch72:  N2+  +  O3p    -->  O+  + N2
    1700            ch72 = 4.1d-10
     1928           !Moffat 2005:
     1929           !ch72 = 4.1d-10
     1930           !McFarland et al. 1974 (From UMIST):
     1931           !ch72 = 1.1d-11
     1932           !Scott et al., 1999 (from F&S2001):
     1933!           ch72 = 7.0d-12 * (300.d0/t_elect)**0.23
     1934           ch72=rcoef(54,1)*((t_elect/300.)**rcoef(54,2))*
     1935     $       exp(rcoef(54,3)/t_elect)
     1936           
    17011937
    17021938ch73  CO+  +  H    -->  H+  + CO
    1703            ch73 = 4.0d-10
     1939           !Scott et al., 1997 (from F&S2001):
     1940!          ch73 = 4.0d-10
     1941           !Federer et al. 1984 (from UMIST):
     1942           !ch73 = 7.5d-10
     1943           ch73=rcoef(55,1)*((t_elect/300.)**rcoef(55,2))*
     1944     $       exp(rcoef(55,3)/t_elect)
    17041945
    17051946
    17061947ch74:   O+  +  H    -->  H+  + O
    1707            ch74 = 5.7d-10 * (tt/300.d0)**(0.36d0)
     1948           !Krasnopolsky 2002:
     1949           !ch74 = 5.7d-10 * (tt/300.d0)**(0.36d0)
     1950           !Stancil et al. 1999 (from UMIST):
     1951!           ch74 = 5.66d-10*(tt/300.)**0.36*exp(8.6/tt)
     1952           !Aninich 1993 (from F&S2001):
     1953           !ch74 = 6.4e-10
     1954           ch74=rcoef(56,1)*((tcte/300.)**rcoef(56,2))*
     1955     $       exp(rcoef(56,3)/tcte)
    17081956
    17091957ch75:  NO+  +  e-  -->  N  +  O
    1710            ch75 = 4.3d-7 * (300.d0/t_elect)**(0.37d0)
    1711 
     1958           !Mitchel 1990 (from UMIST):
     1959!          ch75 = 4.3d-7 * (300.d0/t_elect)**(0.37d0)
     1960           !Vejby-Christensen et al. 1996 (from F&S2001):
     1961           !ch75=4.0d-7 * (300.d0/t_elect)**0.5d0
     1962           ch75=rcoef(57,1)*((t_elect/300.)**rcoef(57,2))*
     1963     $       exp(rcoef(57,3)/t_elect)
    17121964
    17131965ch76:   H+  +  O3p  -->  O+  +  H
    1714            ch76 = 7.3d-10 * (tt/300.d0)**(0.23d0) * exp(-226./tt)
     1966           !Krasnopolsky et al. 2002:
     1967           !ch76 = 7.3d-10 * (tt/300.d0)**(0.23d0) * exp(-226./tt)
     1968           !Stancil et al. 1999 (from UMIST):
     1969!           ch76 = 6.86e-10* (tt/300.)**0.26*exp(-224.3/tt)
     1970           ch76=rcoef(58,1)*((tcte/300.)**rcoef(58,2))*
     1971     $       exp(rcoef(58,3)/tcte)
    17151972
    17161973ch85:  N+  +  CO2    -->  CO2+  + N
    1717            ch85 = 1.d-9
     1974           !Krasnopolsky 2002:
     1975           !ch85 = 1.d-9
     1976           !Adams et al. 1980 (from UMIST):
     1977!           ch85 = 7.5d-10
     1978           !Aninich et al. 1993 (from F&S2001):
     1979           !ch85 = 9.2d-10
     1980           ch85=rcoef(59,1)*((t_elect/300.)**rcoef(59,2))*
     1981     $       exp(rcoef(59,3)/t_elect)
    17181982
    17191983ch86:  H2 + CO2+  --> H + HCO2+
    1720            ch86 = 8.7d-10
    1721 
    1722 ch87:  HCO2+ + e -> CO2 + H
    1723            ch87 = 3.4d-7 * (300.d0/t_elect)**(0.5d0)
     1984           !Scott et al. 1998 (from F&S2001 and K2002):
     1985           !ch86 = 8.7d-10
     1986           !Copp et al. 1982 (from UMIST):
     1987!           ch86 = 9.5d-10
     1988           ch86=rcoef(60,1)*((t_elect/300.)**rcoef(60,2))*
     1989     $       exp(rcoef(60,3)/t_elect)
     1990
     1991
     1992c     h87:  HCO2+ + e -> CO2 + H
     1993           !Krasnopolsky 2002:
     1994!           ch87 = 3.4d-7 * (300.d0/t_elect)**(0.5d0)
     1995           !UMIST 2012: the reactions has 3 different sets of products: CO2+H,
     1996           !CO+O+H (dominante) y CO+OH. Habria que tener esto en cuenta
     1997           ch87=rcoef(61,1)*((t_elect/300.)**rcoef(61,2))*
     1998     $       exp(rcoef(61,3)/t_elect)
    17241999
    17252000           
    … …  
    23952670
    23962671         ! O+
    2397             oplus_eq(i)='N'
    2398 !         oplus_eq(i)='Y'
    2399 !         xn_comp_en_EQ = xn_comp_en_EQ + 1
     2672!            oplus_eq(i)='N'
     2673         oplus_eq(i)='Y'
     2674         xn_comp_en_EQ = xn_comp_en_EQ + 1
    24002675
    24012676         ! N2+
    … …  
    24042679
    24052680         ! H+
    2406             hplus_eq(i)='N'
    2407 
    2408 !         hplus_eq(i)='Y'
    2409 !         xn_comp_en_EQ = xn_comp_en_EQ + 1
     2681!            hplus_eq(i)='N'
     2682
     2683         hplus_eq(i)='Y'
     2684         xn_comp_en_EQ = xn_comp_en_EQ + 1
    24102685
    24112686         ! CO2+

  • trunk/LMDZ.MARS/libf/phymars/physiq.F

    r1005 r1013  
    470470         if(callnlte.and.nltemodel.eq.2) call nlte_setup
    471471         if(callnirco2.and.nircorr.eq.1) call NIR_leedat
     472         if(thermochem) call chemthermos_readini
    472473
    473474        IF (tracer.AND.water.AND.(ngridmx.NE.1)) THEN
    … …  
    17271728               call wstats(ngrid,'albedo',
    17281729     &                         'albedo',
    1729      &                         '',2,albedo(1:ngridmx,1))
     1730     &                         '',2,albedo(1,1))
    17301731               call wstats(ngrid,"mtot",
    17311732     &                    "total mass of water vapor","kg/m2",
    … …  
    20442045            CALL WRITEDIAGFI(ngridmx,'albedo',
    20452046     &                         'albedo',
    2046      &                         '',2,albedo(1:ngridmx,1))
     2047     &                         '',2,albedo(1,1))
    20472048              if (tifeedback) then
    20482049                 call WRITEDIAGSOIL(ngridmx,"soiltemp",
Note: See TracChangeset for help on using the changeset viewer.