source: trunk/LMDZ.MARS/libf/aeronomars/iono.h @ 643

Last change on this file since 643 was 635, checked in by emillour, 13 years ago

Mars GCM: Update of the chemistry package, including:

  • 93 reactions are accounted for (instead of 22); tracking 28 species (instead of 11)
  • computation of photoabsorption using raytracing
  • improved time stepping in the photochemistry
  • updated parameters (cross-sections); with this new version input files

are in 'EUV/param_v5' of "datafile" directory.

  • transition between lower and upper atmosphere chemistry set to 0.1 Pa (calchim.F90)
  • Lots of code clean-up: removed obsolete files column.F, param_v3.h, flujo.F, phdisrate.F, ch.F, interpfast.F, paramfoto.F, getch.F Converted chemtermos.F -> chemthermos.F90 and euvheat.F -> euvheat.F90. Added paramfoto_compact.F , param_v4.h and iono.h
  • Upadted surfacearea.F
  • Cleaned initracer.F and callkeys.h (removed use of obsolete "nqchem" and "oldnames" case when initializing tracers).
  • Minor correction in "callsedim": compute "rdust" and/or "rice" only when it makes sense.

FGG+FL+EM
File size: 2.0 KB
Line 
1      common/eqF/ o1d_eq,ho2_eq,oh_eq,h_eq,n2d_eq,no2_eq,
2     $ o3_eq,no_eq,
3     $ cplus_eq,coplus_eq,oplus_eq,n2plus_eq,hplus_eq,co2plus_eq,
4     $ o2plus_eq,noplus_eq,nplus_eq,hco2plus_eq
5
6        character*1 o1d_eq(nlayermx),ho2_eq(nlayermx),oh_eq(nlayermx)
7      character*1 h_eq(nlayermx),n2d_eq(nlayermx),no2_eq(nlayermx)
8      character*1 o3_eq(nlayermx),no_eq(nlayermx),cplus_eq(nlayermx)
9      character*1 coplus_eq(nlayermx),oplus_eq(nlayermx)
10      character*1 n2plus_eq(nlayermx)
11      character*1 hplus_eq(nlayermx),co2plus_eq(nlayermx)
12      character*1 o2plus_eq(nlayermx)
13      character*1 noplus_eq(nlayermx), nplus_eq(nlayermx) 
14      character*1 hco2plus_eq(nlayermx)     
15
16      common/taus/ tauco2,tauo2,tauo3p,tauco,tauh,tauoh,tauho2,tauh2,
17     $ tauh2o,tauo1d,tauh2o2,tauo3,taun,tauno,taun2,taun2d,tauno2,
18     $ tauco2plus,tauoplus,tauo2plus,taucoplus,taucplus,taunplus,
19     $ taunoplus,taun2plus,tauhplus,tauhco2plus
20
21      real*8    tauco2(nlayermx,nreact)
22      real*8    tauo2(nlayermx,nreact)
23      real*8    tauo3p(nlayermx,nreact)
24      real*8    tauco(nlayermx,nreact)
25      real*8    tauh(nlayermx,nreact)
26      real*8    tauoh(nlayermx,nreact)
27      real*8    tauho2(nlayermx,nreact)
28      real*8    tauh2(nlayermx,nreact)
29      real*8    tauh2o(nlayermx,nreact)
30      real*8    tauo1d(nlayermx,nreact)
31      real*8    tauh2o2(nlayermx,nreact)
32      real*8    tauo3(nlayermx,nreact)
33      real*8    taun(nlayermx,nreact)
34      real*8    tauno(nlayermx,nreact)
35      real*8    taun2(nlayermx,nreact)
36      real*8    taun2d(nlayermx,nreact)
37      real*8    tauno2(nlayermx,nreact)
38      real*8    tauco2plus(nlayermx,nreact)
39      real*8    tauoplus(nlayermx,nreact)
40      real*8    tauo2plus(nlayermx, nreact)
41      real*8    taucoplus(nlayermx, nreact)
42      real*8    taucplus(nlayermx, nreact)
43      real*8    taunplus(nlayermx, nreact) 
44      real*8    taunoplus(nlayermx, nreact)
45      real*8    taun2plus(nlayermx, nreact)
46      real*8    tauhplus(nlayermx, nreact)
47      real*8    tauhco2plus(nlayermx, nreact)
Note: See TracBrowser for help on using the repository browser.