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

Last change on this file since 937 was 705, checked in by emillour, 13 years ago

Mars GCM:

  • Added possibility to run with a varying EUV cycle following real one. The flag solvarmod=1 triggers this behaviour, with companion flag solvaryear=## , where ## is the Mars Year (from 23 to 30). Setting solvarmod=0 reverts to 'old' behaviour, where there is a constant EUV forcing throughout the run, set by the "solarcondate" flag.
  • Needs corresponding input data files ("param_v6" subdirectory of "EUV" subdirectory in "datadir").
  • Added files jthermcalc_e107.F and param_read_e107.F in "aeronomars", modified files euvheat.F90, hrtherm.F, chemthermos.F90, param_v4.h and param_read.F in "aeronomars" and files inifis.F, physiq.F and callkeys.h in "phymars".

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