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

Last change on this file since 461 was 334, checked in by aslmd, 13 years ago

LMDZ.MARS:

28/10/11 == FL + AS

ADDED THE NEW FRAMEWORK FOR PHOTOCHEMISTRY
This is not the last version. A new rewritten version of calchim.F [using LAPACK] is yet to be committed.
--> A new version for photochemistry routines : *_B.F no longer exists, new routines in aeronomars
D 333 libf/aeronomars/photochemist_B.F
D 333 libf/aeronomars/init_chimie_B.F
A 0 libf/aeronomars/read_phototable.F
M 333 libf/aeronomars/calchim.F
A 0 libf/aeronomars/photochemistry.F
M 333 libf/aeronomars/chimiedata.h
--> Changed calls to calchim and watercloud [surfdust and surfice needed] in physiq.F
--> Left commented code for outputs in physiq.F [search for FL]
--> Added settings which works for 35 levels in inidissip.F according to FL. Commented for the moment.
--> Checked compilation and run, looks fine. Note that 'jmars.20101220' is needed.

File size: 1.2 KB
RevLine 
[38]1c ------------------------------------
2c ---  Data for chemical routines  ---
3c ------------------------------------
4c
[334]5      integer nd, nozo, nr, nsza, ntemp, ntau
[38]6c
[334]7      parameter (nd    = 29)
[38]8      parameter (nozo  = 7)
9      parameter (nr    = nd + 28)
10      parameter (nsza  = 27)
11      parameter (ntemp = 4)
[334]12      parameter (ntau  = 8)
[38]13c
[334]14      real kb
15      parameter (kb = 1.3806e-23)
[38]16c
[334]17      common/chimiedata/jphot,colairtab,table_ozo
[38]18c
[334]19      real jphot(ntemp,nsza,0:200,nozo,ntau,nd)
20      real colairtab(0:200)
21      real szatab(nsza)
22      real table_ozo(nozo)
23      real tautab(ntau)
[38]24c
[334]25      data szatab/0.,  5., 10., 15., 20., 25.,
26     $            30., 35., 40., 45., 50., 55.,
27     $            60., 65., 70., 75., 80., 82.,
28     $            84., 86., 88., 90., 91., 92.,
29     $            93., 94., 95./
[38]30c
[334]31      data tautab/0., 0.2, 0.4, 0.6, 0.8, 1., 2., 4./
[38]32c
33      integer       ncomp
[334]34      parameter    (ncomp = 15)
[38]35      character*10  nomchem(ncomp)
36      data          nomchem/"co2","co","o","o1d","o2","o3","h","h2",
[334]37     $                      "oh","ho2","h2o2", "ch4", "n2", "ar", "h2o"/
[38]38c
[334]39      real mmolchem(ncomp)
40      data mmolchem/44.,28.,16.,16.,32.,48., 1.,
41     $               2.,17.,33.,34.,16.,28.,40.,18./
Note: See TracBrowser for help on using the repository browser.