source: trunk/LMDZ.MARS/libf/aeronomars/chimiedata.h @ 1448

Last change on this file since 1448 was 1430, checked in by flefevre, 10 years ago

aeronomars : preparation de l'implementation du solveur chimique ASIS

1) extension de la table de photodissociation "jmars" a la thermosphere

et restriction aux seules especes necessaires.

2) modernisation de la subroutine d'interpolation dans cette

table de photodissociation.

3) chimiedata.h compatible avec solveur ASIS

ATTENTION!! la nouvelle table de photodissociation

jmars.20140930

doit etre presente dans /datadir sinon crash modele...

File size: 1.5 KB
Line 
1!--------------------------------------------
2!     data for photochemistry
3!--------------------------------------------
4
5!--------------------------------------------
6!     dimensions of photolysis lookup table
7!--------------------------------------------
8
9      integer, parameter :: nd    = 13  ! species
10      integer, parameter :: nz    = 143 ! altitude
11      integer, parameter :: nozo  = 7   ! ozone
12      integer, parameter :: nsza  = 27  ! solar zenith angle
13      integer, parameter :: ntemp = 4   ! temperature
14      integer, parameter :: ntau  = 8   ! dust
15
16!--------------------------------------------
17
18      real, parameter :: kb = 1.3806e-23
19
20      common/chimiedata/jphot,colairtab,table_ozo
21
22      real jphot(ntemp,nsza,nz,nozo,ntau,nd)
23      real colairtab(nz)
24      real szatab(nsza)
25      real table_ozo(nozo)
26      real tautab(ntau)
27
28      data szatab/0.,  5., 10., 15., 20., 25.,                          &
29     &            30., 35., 40., 45., 50., 55.,                         &
30     &            60., 65., 70., 75., 80., 82.,                         &
31     &            84., 86., 88., 90., 91., 92.,                         &
32     &            93., 94., 95./
33
34      data tautab/0., 0.2, 0.4, 0.6, 0.8, 1., 2., 4./
35
36!--------------------------------------------
37!     number of reactions in ASIS solver
38!--------------------------------------------
39
40      integer, parameter :: nb_phot_max       = 18
41      integer, parameter :: nb_reaction_3_max = 6
42      integer, parameter :: nb_reaction_4_max = 30
Note: See TracBrowser for help on using the repository browser.