source: trunk/LMDZ.MARS/libf/aeronomars/param_v4.h @ 657

Last change on this file since 657 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: 9.2 KB
Line 
1!**********************************************************************
2
3!       param.cmn
4       
5!       fgg     mar 02          first version
6!**********************************************************************
7
8
9      common/uv/ crscabsi2,c1_16,c17_24,c25_29,c30_31,c32,                         &
10     &   c33,c34,c35,c36,jfotsout,co2crsc195,co2crsc295,                           &
11     &   t0,fluxtop,freccen,jabsifotsintpar
12     
13     
14
15      real crscabsi2(nabs,16)    !cross section
16      real c1_16(nz2,16)   !Col. int. 1 (o2+o+h2+n) (cm^-2)
17      real c17_24(nz2)     !Col. int. 17-24 (co2+o2+n2+no+co+no2) (cm^-2)
18      real c25_29(nz2)     !Col. int. 25-29 (co2+o2+h2o+h2o2+no+co+no2) (cm^-2)
19      real c30_31(nz2)     !Col. int. 30-31 (co2+o2+h2o+h2o2+no+no2)
20      real c32(nz2)        !col. int. 32 (co2+h2o2+o2+no+no2) (cm^-2)
21      real c33(nz2)        !col. int. 33 (h2o2+o2+no2) (cm^-2)
22      real c34(nz2)        !col. int. 34 (h2o2+o2+o3+no2) (cm^-2)
23      real c35(nz2)        !col. int. 35 (h2o2+o3+no2) (cm^-2)
24      real c36(nz2)        !col. int. 36 (o3+no2) (cm^-2)
25      real jfotsout(ninter,nabs,nlayermx)
26      real co2crsc195(9)
27      real co2crsc295(9)
28      real t0(nz2)
29      real fluxtop(ninter)
30      real freccen(ninter)          !representative wavelenght
31      real jabsifotsintpar(ninter,nabs,nz2)
32
33
34      common/fotoquimica/  ch2, ch3, ch4, ch5, ch7,ch9,ch10,ch11,                &
35     &                     ch13,ch14,ch15,ch18,ch19,ch20,ch21,ch22,              &
36     &                     ch23,ch24,ch30,ch31,ch32,ch33,ch34,ch35,              &
37     &                     ch36,ch37,ch38,ch39,ch40,ch41,ch42,ch43,              &
38     &                     ch45,ch46,ch47,ch48,ch49,ch50,ch55,ch56,              &
39     &                     ch57,ch58,ch59,ch62,ch63,ch64,ch65,                   &
40     &                     ch66,ch67,ch68,ch69,ch70,ch71,ch72,                   &
41     &                     ch73,ch74,ch75,ch76,ch85,ch86,ch87,                   &
42     &                     jdistot,jdistot_b,jion
43
44      real*8 ch2, ch3, ch4, ch5, ch7,ch9,ch10,ch11,ch13,ch14,ch15,ch18
45      real*8 ch19,ch20,ch21,ch22,ch23,ch24,ch30,ch31,ch32,ch33,ch34
46      real*8 ch35,ch36,ch37,ch38,ch39,ch40,ch41,ch42,ch43,ch45
47      real*8 ch46,ch47,ch48,ch49,ch50,ch55,ch56,ch57,ch58,ch59,ch62
48      real*8 ch63,ch64,ch65,ch66,ch67,ch68,ch69,ch70,ch71
49      real*8 ch72,ch73,ch74,ch75,ch76,ch85,ch86,ch87
50
51                                                  !reaction rates
52      real jdistot(nabs,nlayermx)    !photodissociation rates
53      real jdistot_b(nabs,nlayermx)
54     
55      real jion(nabs,nlayermx,4)
56
57
58
59      common/prod/ Pco2,Po2,Po3p,Pco,Ph,Poh,Pho2,Ph2,Ph2o,Po1d,Ph2o2,          &
60     &             Po3,Pn,Pno,Pno2,Pn2,Pn2d,                                   &
61     &             Pco2plus,Poplus,Po2plus,Pelect,Pcoplus,Pcplus,              &
62     &             Pnplus, Pnoplus, Pn2plus, Phplus, Phco2plus,                &
63     &             Pco2tot,Po2tot,Po3ptot,Pcotot,Phtot,Pohtot,                 &
64     &             Pho2tot,Ph2tot,Ph2otot,Po1dtot,Ph2o2tot,                    &
65     &             Po3tot,Pntot,Pnotot,Pno2tot,Pn2tot,Pn2dtot,                 &
66     &             Pcoplustot,Pcplustot,                                       &
67     &             Pco2plustot,Poplustot,Po2plustot,Pelecttot,                 &
68     &             Pnplustot,Pnoplustot,Pn2plustot,Phplustot,                  &
69     &             Phco2plustot 
70
71       
72      real*8 Pco2(nlayermx,nreact),Po2(nlayermx,nreact)
73      real*8 Po3p(nlayermx,nreact)
74      real*8 Pco(nlayermx,nreact),Ph(nlayermx,nreact)
75      real*8 Poh(nlayermx,nreact)
76      real*8 Pho2(nlayermx,nreact),Ph2(nlayermx,nreact)
77      real*8 Ph2o(nlayermx,nreact)
78      real*8 Po1d(nlayermx,nreact),Ph2o2(nlayermx,nreact)
79      real*8 Po3(nlayermx,nreact),Pn(nlayermx,nreact)
80      real*8 Pno(nlayermx,nreact)
81      real*8 Pno2(nlayermx,nreact),Pn2(nlayermx,nreact)
82      real*8 Pn2d(nlayermx,nreact)
83      real*8 Pco2plus(nlayermx,nreact),Poplus(nlayermx,nreact)
84      real*8 Po2plus(nlayermx,nreact), Pelect(nlayermx,nreact)
85      real*8 Pcoplus(nlayermx,nreact),Pcplus(nlayermx,nreact)
86      real*8 Pnplus(nlayermx,nreact),Pnoplus(nlayermx,nreact)
87      real*8 Pn2plus(nlayermx,nreact),Phplus(nlayermx,nreact)
88      real*8 Phco2plus(nlayermx,nreact)
89      real*8 Pco2tot(nlayermx),Po2tot(nlayermx)
90      real*8 Po3ptot(nlayermx),Pcotot(nlayermx)
91      real*8 Phtot(nlayermx)
92      real*8 Pohtot(nlayermx),Pho2tot(nlayermx)
93      real*8 Ph2tot(nlayermx),Ph2otot(nlayermx)
94      real*8 Po1dtot(nlayermx),Ph2o2tot(nlayermx)
95      real*8 Po3tot(nlayermx),Pntot(nlayermx),Pnotot(nlayermx)
96      real*8 Pno2tot(nlayermx),Pn2tot(nlayermx),Pn2dtot(nlayermx)
97      real*8 Pco2plustot(nlayermx), Poplustot(nlayermx)
98      real*8 Po2plustot(nlayermx),Pelecttot(nlayermx)
99      real*8 Pcoplustot(nlayermx), Pcplustot(nlayermx)
100      real*8 Pnplustot(nlayermx), Pnoplustot(nlayermx)
101      real*8 Pn2plustot(nlayermx),Phplustot(nlayermx)
102      real*8 Phco2plustot(nlayermx)
103
104
105      common/loss/ Lco2,Lo2,Lo3p,Lco,Lh,Loh,Lho2,Lh2,Lh2o,Lo1d,Lh2o2,         &
106     &             Lo3,Ln,Lno,Lno2,Ln2,Ln2d,Lelect,                           &
107     &             Lco2plus,Lo2plus,Loplus, Lcoplus, Lcplus,                  &
108     &             Lnplus,Lnoplus,Ln2plus, Lhplus, Lhco2plus,                 &
109     &             Lco2tot,Lo2tot,Lo3ptot,Lcotot,Lhtot,Lohtot,                &
110     &             Lho2tot,Lh2tot,Lh2otot,Lo1dtot,Lh2o2tot,                   &
111     &             Lo3tot,Lntot,Lnotot,Lno2tot,Ln2tot,Ln2dtot,                &
112     &             Lcoplustot,Lcplustot,Lco2plustot,Loplustot,                &
113     &             Lo2plustot,Lelecttot,Lnplustot,Lnoplustot,                 &
114     &             Ln2plustot,Lhplustot,Lhco2plustot
115
116       
117      real*8 Lco2(nlayermx,nreact),Lo2(nlayermx,nreact)
118      real*8 Lo3p(nlayermx,nreact)
119      real*8 Lco(nlayermx,nreact),Lh(nlayermx,nreact)
120      real*8 Loh(nlayermx,nreact)
121      real*8 Lho2(nlayermx,nreact),Lh2(nlayermx,nreact)
122      real*8 Lh2o(nlayermx,nreact)
123      real*8 Lo1d(nlayermx,nreact),Lh2o2(nlayermx,nreact)
124      real*8 Lo3(nlayermx,nreact),Ln(nlayermx,nreact)
125      real*8 Lno(nlayermx,nreact)
126      real*8 Lno2(nlayermx,nreact),Ln2(nlayermx,nreact)
127      real*8 Ln2d(nlayermx,nreact)
128      real*8 Lco2plus(nlayermx,nreact), Loplus(nlayermx,nreact)
129      real*8 Lo2plus(nlayermx,nreact), Lelect(nlayermx,nreact)
130      real*8 Lcoplus(nlayermx,nreact), Lcplus(nlayermx,nreact)
131      real*8 Lnplus(nlayermx,nreact), Lnoplus(nlayermx,nreact)
132      real*8 Ln2plus(nlayermx,nreact), Lhplus(nlayermx,nreact)
133      real*8 Lhco2plus(nlayermx,nreact)
134      real*8 Lco2tot(nlayermx),Lo2tot(nlayermx)
135      real*8 Lo3ptot(nlayermx),Lcotot(nlayermx)
136      real*8 Lhtot(nlayermx),Lohtot(nlayermx)
137      real*8 Lho2tot(nlayermx),Lh2tot(nlayermx)
138      real*8 Lh2otot(nlayermx),Lo1dtot(nlayermx)
139      real*8 Lh2o2tot(nlayermx)
140      real*8 Lo3tot(nlayermx),Lntot(nlayermx),Lnotot(nlayermx)
141      real*8 Lno2tot(nlayermx),Ln2tot(nlayermx)
142      real*8 Ln2dtot(nlayermx)
143      real*8 Lco2plustot(nlayermx), Loplustot(nlayermx)
144      real*8 Lo2plustot(nlayermx), Lelecttot(nlayermx)
145      real*8 Lcoplustot(nlayermx), Lcplustot(nlayermx)
146      real*8 Lnplustot(nlayermx), Lnoplustot(nlayermx)
147      real*8 Ln2plustot(nlayermx), Lhplustot(nlayermx)
148      real*8 Lhco2plustot(nlayermx)
149
150
151
152      common/vida/ tminco2,tmino2,tmino3p,tminco,tminh,tminoh,tminho2,        &
153     &               tminh2,tminh2o,tmino1d,tminh2o2,tmino3,tminn,            &
154     &               tminno,tminno2,tminn2,tminn2d,                           &
155     &               tminco2plus,tminoplus,tmino2plus,tmincoplus,             &
156     &               tmincplus,tminnplus,tminnoplus,                          &
157     &               tminn2plus,tminhplus,tminhco2plus
158
159
160
161        real*8 tminco2(nlayermx),tmino2(nlayermx)
162        real*8 tmino3p(nlayermx),tminco(nlayermx)
163        real*8 tminh(nlayermx),tminoh(nlayermx)
164        real*8 tminho2(nlayermx),tminh2(nlayermx)
165        real*8 tminh2o(nlayermx),tmino1d(nlayermx),tminh2o2(nlayermx)
166        real*8 tmino3(nlayermx),tminn(nlayermx),tminno(nlayermx)
167        real*8 tminno2(nlayermx),tminn2(nlayermx),tminn2d(nlayermx)
168        real*8 tminco2plus(nlayermx),tminoplus(nlayermx)
169        real*8 tmino2plus(nlayermx)
170        real*8 tmincoplus(nlayermx)
171        real*8 tmincplus(nlayermx)
172        real*8 tminnplus(nlayermx),tminnoplus(nlayermx)
173        real*8 tminn2plus(nlayermx)
174        real*8 tminhplus(nlayermx),tminhcoplus(nlayermx)
175        real*8 tminhco2plus(nlayermx)
176
177
178      common/flujoctes/ fluxtophr,ct1,p1,ct2,p2   
179 
180        real fluxtophr(ninter)
181        real ct1(ninter),p1(ninter),ct2(ninter),p2(ninter)       
182
183      common/phdisef/ efdisco2, efdiso2, efdish2, efdish2o, efdish2o2,        &
184     &                efdiso3,efdiso,efdisn,efdish,efdisno,efdisn2,           &
185     &                efdisno2,efdisco,efionco2,efionn2,efionco,              &
186     &                efiono3p,efionn,efionno,efionh
187 
188
189        real efdisco2(ninter), efdiso2(ninter), efdish2o(ninter)
190        real efdish2o2(ninter), efdish2(ninter), efdiso3(ninter)
191        real efdiso(ninter), efdisn(ninter), efdish(ninter)
192        real efdisno(ninter), efdisn2(ninter), efdisno2(ninter)
193        real efdisco(ninter)
194        real efionco2(ninter,4)
195        real efionn2(ninter,2)
196        real efionco(ninter,3)
197        real efiono3p(ninter),efionn(ninter)
198        real efionno(ninter),efionh(ninter)
199
200
201
202
203
204
Note: See TracBrowser for help on using the repository browser.