Context Navigation

source: trunk/LMDZ.TITAN.old/libf/chimtitan/disso.c @ 1990

Last change on this file since 1990 was 1126, checked in by slebonnois, 11 years ago
SL: update of Titan photochemical module to include computation of chemistry up to 1300 km
File size: 21.9 KB

Line
1	/* disso: photodissociation rates */
2	/* GCCM */
3	/* correspond a chimie_simpnit (version 301105) */
4	/* !!! ATTENTION !!! */
5	/* Doit etre mis a jour en fonction de la chimie utilisee ! */
6
7	#include "titan.h"
8
9	void disso_( double KRPD[][NLRT][RDISS+1][15], int *NLAT )
10	{
11	static double sH2[62] = { /* incertain en dessous de 70 et en dessus de 85... */
12	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
13	0.000e+00,0.000e+00,0.000e+00,0.000e+00,1.000e-18,5.000e-18,1.000e-17,
14	9.000e-18,6.500e-18,1.000e-18,1.000e-19 };
15	static double sCH4[62] = {
16	2.852e-19,7.816e-19,1.534e-18,2.069e-18,2.795e-18,4.088e-18,4.543e-18,
17	4.223e-18,3.314e-18,1.565e-18,8.892e-19,8.760e-19,8.792e-19,9.163e-19,
18	2.069e-18,9.378e-18,2.543e-17,3.785e-17,4.066e-17,3.302e-17,2.840e-17,
19	1.800e-17,1.920e-17,1.820e-17,1.840e-17,1.140e-17,2.656e-18,1.256e-19,
20	7.988e-22,1.366e-23,6.740e-24 };
21	static double sCH3CN[62] = {
22	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
23	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
24	0.000e+00,0.000e+00,0.000e+00,0.000e+00,8.000e-17,8.000e-17,7.500e-17,
25	4.800e-17,3.700e-17,2.700e-17,4.100e-17,1.000e-17,8.600e-18,4.000e-18,
26	1.800e-18,1.100e-18,6.400e-19,3.600e-19,2.000e-19,1.200e-19,5.200e-20};
27	static double sC2H2[62] = {
28	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
29	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
30	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,3.000e-17,
31	3.000e-17,3.200e-17,3.640e-17,4.260e-17,1.040e-16,9.900e-18,4.800e-18,
32	1.720e-17,1.782e-17,9.040e-19,9.600e-19,1.294e-18,1.352e-18,1.130e-18,
33	6.680e-19,3.700e-19,3.900e-19,1.660e-19,5.400e-20 };
34	static double sC2H4[62] = {
35	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
36	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
37	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
38	0.000e+00,.1992E-16,.2734E-16,.2594E-16,.1690E-16,.2258E-16,.8507E-17,
39	.1583E-16,.2227E-16,.3056E-16,.3743E-16,.3788E-16,.2737E-16,.3171E-17,
40	.6033E-18,.1223E-18,.7247E-19,1.000e-20 };
41	static double sC2H3CN[62] = {
42	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
43	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
44	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e-00,0.000e-00,0.000e-00,
45	0.000e-00,0.000e-00,1.000e-17,6.600e-18,5.600e-18,3.600e-18,4.300e-18,
46	4.000e-18,3.000e-18,2.900e-18,2.700e-18,2.700e-18,3.300e-18,4.500e-18,
47	6.000e-18,7.100e-18,7.000e-18,5.000e-18,2.500e-18,6.600e-19,1.000e-19,
48	1.500e-20,1.000e-21};
49	static double sC2H6[62] = {
50	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
51	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
52	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,2.000e-17,
53	2.000e-17,2.000e-17,2.060e-17,2.160e-17,1.540e-17,8.060e-18,3.860e-18,
54	1.484e-18,3.060e-19,9.600e-21 };
55	static double sCH3C2H[62] = {
56	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
57	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
58	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,1.000e-17,
59	2.000e-17,4.580e-17,6.240e-17,5.880e-17,5.920e-17,1.940e-17,2.200e-17,
60	2.820e-17,3.380e-17,1.142e-17,6.600e-18,8.100e-18,7.000e-18,2.800e-18,
61	1.600e-18,2.300e-19,4.411e-19,2.119e-19,1.004e-19,2.934e-20,4.157e-21 };
62	static double sCH2CCH2[62] = {
63	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
64	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
65	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
66	0.000e+00,0.000e+00,2.200e-17,2.200e-17,1.700e-17,1.700e-17,1.800e-17,
67	1.200e-17,1.500e-17,1.000e-17,2.100e-17,3.300e-17,4.000e-17,1.300e-17,
68	5.000e-18,2.900e-18,2.601e-18,1.037e-18,9.046e-19,6.565e-19,4.672e-19,
69	3.047e-19,1.579e-19,5.943e-20,2.261e-20 };
70	static double sC3H6[62] = {
71	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
72	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
73	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,3.400e-17,
74	3.300e-17,2.400e-17,2.500e-17,4.000e-17,3.700e-17,2.300e-17,1.900e-17,
75	2.000e-17,1.500e-17,2.200e-17,2.500e-17,4.400e-17,4.200e-17,2.700e-17,
76	1.200e-17,5.800e-18,1.400e-19,9.100e-21 };
77	static double sC3H8[62] = {
78	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
79	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
80	0.000e+00,0.000e+00,0.000e+00,0.000e+00,4.000e-17,4.000e-17,4.000e-17,
81	4.000e-17,4.000e-17,3.280e-17,3.100e-17,2.680e-17,2.200e-17,1.760e-17,
82	6.440e-18,3.000e-18,9.140e-19,7.000e-20 };
83	static double sC4H2[62] = {
84	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
85	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
86	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
87	0.000e+00,0.000e+00,2.500e-17,6.100e-17,4.780e-17,5.780e-17,1.014e-16,
88	5.200e-17,4.200e-17,1.030e-16,1.770e-16,9.100e-17,1.290e-17,2.380e-18,
89	3.400e-19,2.800e-19,2.600e-19,1.513e-19,2.583e-19,3.353e-19,4.115e-19,
90	4.755e-19,4.990e-19,4.399e-19,5.358e-19,2.485e-19,3.931e-19,1.067e-19,
91	3.761e-20,2.370e-20,2.277e-20 };
92	static double sC4H4[62] = {
93	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
94	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
95	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
96	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
97	0.000e+00,0.000e+00,0.000e+00,4.200e-17,5.000e-17,2.000e-17,6.700e-18,
98	6.100e-18,8.000e-18,1.000e-17,1.400e-17,1.900e-17,2.200e-17,2.600e-17,
99	1.600e-17,1.600e-17,4.000e-18,5.700e-19,1.400e-19 };
100	static double sC4H6[62] = {
101	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
102	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
103	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
104	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
105	0.000e+00,0.000e+00,0.000e+00,6.000e-18,1.200e-17,1.500e-17,1.300e-17,
106	1.800e-17,3.100e-17,3.800e-17,6.600e-17,9.600e-17,1.060e-16,8.500e-17,
107	6.700e-17,1.100e-17,1.700e-18,3.300e-19,6.500e-20 };
108	static double sC4H10[62] = {
109	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
110	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
111	0.000e+00,0.000e+00,0.000e+00,0.000e+00,6.000e-17,6.000e-17,6.000e-17,
112	6.000e-17,6.000e-17,5.500e-17,4.400e-17,4.400e-17,3.800e-17,3.100e-17,
113	1.900e-17,4.000e-18,1.300e-18,3.200e-19,2.000e-20 };
114	static double sC6H6[62] = {
115	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
116	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
117	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,7.000e-17,
118	8.000e-17,5.000e-17,3.500e-17,3.500e-17,3.500e-17,2.000e-17,1.500e-17,
119	1.500e-17,1.500e-17,2.000e-17,2.500e-17,4.000e-17,9.500e-17,2.200e-16,
120	1.000e-16,2.000e-17,2.000e-17,2.000e-17,2.000e-17,5.000e-18,1.000e-20,
121	1.000e-20,1.000e-20,1.000e-20,1.000e-20,1.000e-19,2.000e-19,2.500e-19,
122	4.000e-19,2.000e-19,2.000e-19,1.000e-19 };
123	static double sN2[62] = {
124	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
125	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,4.898e-18,1.097e-17,
126	2.192e-17,2.214e-17,2.336e-17,1.679e-17,1.893e-17 };
127	static double sHCN[62] = {
128	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
129	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
130	0.000e+00,0.000e+00,0.000e+00,2.800e-17,3.300e-17,2.800e-17,3.500e-17,
131	4.800e-17,2.800e-17,1.700e-17,3.600e-18,2.500e-18,4.700e-18,2.700e-18,
132	8.600e-19,4.400e-19,2.000e-19,1.429e-19,1.145e-19,7.482e-20,3.852e-20,
133	1.009e-20 };
134	static double sHC3N[62] = {
135	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
136	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
137	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,3.402e-17,
138	3.647e-17,4.317e-17,3.759e-17,7.927e-17,3.796e-17,9.565e-17,1.716e-16,
139	1.247e-16,2.360e-17,8.411e-18,4.400e-18,8.600e-19,7.400e-19,6.200e-19,
140	4.899e-19,3.307e-19,2.128e-19,2.561e-19,2.621e-19,2.737e-19,3.601e-19,
141	1.564e-19,1.816e-19,8.427e-20 };
142	static double sC2N2[62] = {
143	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
144	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
145	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,2.635e-17,
146	3.126e-17,4.392e-17,9.788e-17,1.682e-16,1.064e-16,6.513e-18,2.039e-18,
147	2.828e-18,5.136e-18,8.188e-18,8.857e-18,1.489e-18,9.000e-20,6.200e-20,
148	3.800e-20,4.483e-20,8.618e-20,1.008e-19,7.579e-20,6.666e-20,2.907e-20,
149	2.476e-20,1.142e-20 };
150	static double sH2O[62] = {
151	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
152	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
153	0.000e+00,0.000e+00,0.000e+00,0.000e+00,1.927e-17,9.633e-18,7.440e-18,
154	1.087e-18,8.600e-18,8.850e-18,7.300e-18,4.100e-18,1.240e-18,5.600e-19,
155	7.977e-19,1.888e-18,3.333e-18,4.729e-18,4.963e-18,3.513e-18,1.440e-18,
156	1.563e-19 };
157	static double sCO[62] = {
158	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
159	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
160	0.000e+00,0.000e+00,7.300e-17,2.050e-16,9.000e-17,7.180e-18,3.700e-17};
161	static double sCO2[62] = {
162	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
163	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
164	0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,
165	0.000e+00,1.022e-19,1.070e-19,4.220e-19,8.220e-19,6.640e-19,5.650e-19,
166	5.850e-19,4.400e-19,2.100e-19,8.700e-20,4.000e-20,2.130e-20,3.700e-21,
167	6.930e-22,4.040e-22,2.470e-23,4.800e-24,9.500e-25 };
168	static double sol[62] = {
169	9.45e+08,2.84e+08,4.24e+09,3.11e+09,1.22e+10,6.83e+09,2.54e+09,8.20e+08,
170	1.82e+09,1.30e+09,1.05e+09,1.37e+09,3.87e+08,4.29e+08,2.45e+09,4.09e+09,
171	1.06e+10,1.03e+10,4.08e+09,2.98e+10,7.00e+09,1.94e+09,7.40e+09,4.12e+11,
172	1.36e+10,4.92e+10,3.04e+10,3.81e+10,5.49e+10,1.00e+11,1.26e+11,1.67e+11,
173	2.97e+11,4.92e+11,7.17e+11,9.43e+11,1.34e+12,1.96e+12,2.96e+12,4.34e+12,
174	7.25e+12,1.76e+13,2.21e+13,3.10e+13,2.80e+13,2.97e+13,3.03e+13,3.77e+13,
175	3.44e+13,3.42e+13,7.08e+13,9.58e+13,1.63e+14,1.38e+14,1.26e+14,1.60e+14,
176	2.23e+14,3.70e+14,3.31e+14,3.26e+14,3.59e+14,4.05e+14};
177
178	int i,j,l,s,lat,x;
179	double f,**flux;
180	double flact;
181	char name[60],dir[24];
182	FILE fp,out;
183
184	flux = dm2d(0,NLRT,0,14);
185
186	/* lecture des flux actiniques:
187	- suppose que l'executable est dans $LMDGCM/RUN/xxx/
188	- et les moyennes dans $LMDGCM/INPUT/PHOT(NLAT)/Moy_(lat:1 a NLAT)
189	*/
190	strcpy( dir, "../../INPUT/PHOT" );
191	if( (*NLAT) < 10 )
192	{
193	strcat( dir, "0" );
194	strcat( dir, (const char )ecvt((float)(NLAT),1,&x,&x) );
195	}
196	else
197	strcat( dir, (const char )ecvt((float)(NLAT),2,&x,&x) );
198	strcat( dir, "x/Moy_" );
199	printf( "Directories for actinic fluxes: %s \n", dir );
200
201	for( lat = 0; lat <= (NLAT)-1; lat++ ) / Old array is set equal to 0. */
202	for( j = 0; j <= NLRT-1; j++ )
203	for( i = 0; i <= RDISS; i++ )
204	for( s = 0; s <= 14; s++ )
205	KRPD[lat][j][i][s] = 0.0e0;
206
207	for( i = 0; i <= 13; i++ ) sCH4[i] = 0.0e0;
208	for( i = 0; i <= 13; i++ ) sC2H2[i] = 0.0e0;
209	for( i = 0; i <= 16; i++ ) sC2H4[i] = 0.0e0;
210	for( i = 0; i <= 16; i++ ) sC2H6[i] = 0.0e0;
211
212	for( lat = 0; lat <= (NLAT)-1; lat++ ) / Main loop on latitude */
213	for( i = 10; i <= 310; i += 5 ) /* Main loop on wavelength. */
214	{
215	strcpy( name, dir );
216	if( (lat+1) < 10 )
217	{
218	strcat( name, "0" );
219	strcat( name, (const char *)ecvt((float)(lat+1),1,&x,&x) );
220	}
221	else
222	strcat( name, (const char *)ecvt((float)(lat+1),2,&x,&x) );
223	if( i < 160 ) strcat( name, "/photmoy3a" );
224	else strcat( name, "/photmoy3a" );
225	if( i < 100 )
226	{
227	strcat( name, ".0" );
228	strcat( name, (const char *)ecvt((float)i,2,&x,&x) );
229	}
230	else
231	{
232	strcat( name, "." );
233	strcat( name, (const char *)ecvt((float)i,3,&x,&x) );
234	}
235	if( !( fp = fopen( name, "r" ) ) )
236	{
237	out = fopen( "err.log", "a" );
238	fprintf( out, "I cannot open %s\n", name );
239	fclose( out );
240	exit(0);
241	}
242	for( j = 0; j <= NLRT-1; j++ )
243	{
244	fscanf( fp,"%d ",&l );
245	for( s = 0; s < 15; s++ )
246	{
247	fscanf( fp,"%lg ", &flact );
248	flux[j][s] = flact;
249	}
250	}
251	fclose(fp);
252
253	l = i / 5 - 2; /* Pointer on wavelength. */
254
255	/* taux de photodissociations */
256
257	for( s = 0; s <= 14; s++ )
258	for( j = 0; j <= NLRT-1; j++ )
259	{
260	f = flux[j][s] * sol[l] / ( 9.5e0 * 9.5e0 ); /* !! # de reac de 0 a RDISS-1 !! */
261	if( i == 220 ) KRPD[lat][j][1][s] += 4.4e-17 * f; /* CH3 -> 1CH2 + H */
262	KRPD[lat][j][ 0][s] += sH2[l] * f * 1.00; /* H2 -> H + H */
263	KRPD[lat][j][ 7][s] += sC2H4[l] * f * 0.51; /* C2H4 -> C2H2 + H2 */
264	KRPD[lat][j][ 8][s] += sC2H4[l] * f * 0.49; /* C2H4 -> C2H2 + 2H */
265	KRPD[lat][j][15][s] += sCH2CCH2[l]* f * 0.89; /* CH2CCH2 -> C3H3 + H Jackson 91 */
266	KRPD[lat][j][16][s] += sCH2CCH2[l]* f * 0.11; /* CH2CCH2 -> C3H2 + H2 Jackson 91 */
267	KRPD[lat][j][17][s] += sCH3C2H[l] * f * 0.89; /* CH3C2H -> C3H3 + H Jackson 91 */
268	KRPD[lat][j][18][s] += sCH3C2H[l] * f * 0.11; /* CH3C2H -> C3H2 + H2 Jackson 91 */
269	KRPD[lat][j][19][s] += sC3H6[l] * f * 0.33; /* C3H6 -> CH2CCH2 + H2 */
270	KRPD[lat][j][20][s] += sC3H6[l] * f * 0.17; /* C3H6 -> CH3CCH + H2 */
271	KRPD[lat][j][21][s] += sC3H6[l] * f * 0.03; /* C3H6 -> C2H4 + 3CH2 */
272	KRPD[lat][j][22][s] += sC3H6[l] * f * 0.35; /* C3H6 -> C2H3 + CH3 */
273	KRPD[lat][j][23][s] += sC3H6[l] * f * 0.05; /* C3H6 -> C2H2 + CH4 */
274	KRPD[lat][j][32][s] += sC4H4[l] * f * 0.80; /* C4H4 -> C4H2 + H2 Gladstone 96 */
275	KRPD[lat][j][33][s] += sC4H4[l] * f * 0.20; /* C4H4 -> C2H2 + C2H2 Gladstone 96 */
276	KRPD[lat][j][34][s] += sC4H6[l] * f * 0.04; /* C4H6 -> C4H4 + H2 */
277	KRPD[lat][j][35][s] += sC4H6[l] * f * 0.27; /* C4H6 -> C2H4 + C2H2 */
278	KRPD[lat][j][36][s] += sC4H6[l] * f * 0.69; /* C4H6 -> CH3 + C3H3 */
279	KRPD[lat][j][45][s] += sC6H6[l] * f * 0.04; /* AC6H6 -> C5H3 (prod...) + CH3 */
280	KRPD[lat][j][46][s] += sC6H6[l] * f * 0.96; /* AC6H6 -> AC6H5 + H */
281	KRPD[lat][j][47][s] += sN2[l] * f; /* N2 -> 2N2d */
282	KRPD[lat][j][48][s] += sHCN[l] * f; /* HCN -> H + CN */
283	KRPD[lat][j][51][s] += sC2N2[l] * f * 0.3; /* C2N2 -> 2CN */
284	KRPD[lat][j][52][s] += sCH3CN[l] * f * 1.0; /* CH3CN -> CH3 + CN */
285	KRPD[lat][j][53][s] += sC2N2[l] * f * 0.3; /* C4N2 -> C3N + CN */
286
287	if( i != 125 ) /* Not Lyman alpha */
288	{
289	KRPD[lat][j][ 2][s] += sCH4[l] * f; /* CH4 -> 1CH2 + H2 */
290	KRPD[lat][j][ 9][s] += sC2H6[l] * f * 0.56; /* C2H6 -> C2H4 + H2 */
291	KRPD[lat][j][10][s] += sC2H6[l] * f * 0.14; /* C2H6 -> C2H4 + 2H */
292	KRPD[lat][j][11][s] += sC2H6[l] * f * 0.27; /* C2H6 -> C2H2 + 2H2 */
293	KRPD[lat][j][12][s] += sC2H6[l] * f * 0.02; /* C2H6 -> CH4 + 3CH2 */
294	KRPD[lat][j][13][s] += sC2H6[l] * f * 0.01; /* C2H6 -> 2CH3 */
295	KRPD[lat][j][24][s] += sC3H8[l] * f * 0.94; /* C3H8 -> C3H6 + H2 */
296	KRPD[lat][j][27][s] += sC3H8[l] * f * 0.06; /* C3H8 -> C2H4 + CH4 */
297	}
298	else /* Lyman alpha */
299	{
300	KRPD[lat][j][ 2][s] += sCH4[l] * f * 0.64; /* CH4 -> 1CH2 + H2 */
301	KRPD[lat][j][ 3][s] += sCH4[l] * f * 0.07; /* CH4 -> CH + H2 + H */
302	KRPD[lat][j][ 4][s] += sCH4[l] * f * 0.29; /* CH4 -> CH3 + H */
303	KRPD[lat][j][ 9][s] += sC2H6[l] * f * 0.13; /* C2H6 -> C2H4 + H2 */
304	KRPD[lat][j][10][s] += sC2H6[l] * f * 0.3; /* C2H6 -> C2H4 + 2H */
305	KRPD[lat][j][11][s] += sC2H6[l] * f * 0.25; /* C2H6 -> C2H2 + 2H2 */
306	KRPD[lat][j][12][s] += sC2H6[l] * f * 0.25; /* C2H6 -> CH4 + 3CH2 */
307	KRPD[lat][j][13][s] += sC2H6[l] * f * 0.08; /* C2H6 -> 2CH3 */
308	KRPD[lat][j][24][s] += sC3H8[l] * f * 0.33; /* C3H8 -> C3H6 + H2 */
309	KRPD[lat][j][25][s] += sC3H8[l] * f * 0.09; /* C3H8 -> C2H6 + 3CH2 */
310	KRPD[lat][j][26][s] += sC3H8[l] * f * 0.39; /* C3H8 -> C2H5 + CH3 */
311	KRPD[lat][j][27][s] += sC3H8[l] * f * 0.2; /* C3H8 -> C2H4 + CH4 */
312	}
313	if( i < 145 ) /* C4H10: a revoir avec Jackson & Lias, 1974... */
314	{
315	KRPD[lat][j][37][s] += sC4H10[l]* f * 0.18; /* C4H10 -> C4H8(ieC3H5+CH3)+H2 */
316	KRPD[lat][j][38][s] += sC4H10[l]* f * 0.20; /* C4H10 -> 2 C2H4 + H2 */
317	KRPD[lat][j][39][s] += sC4H10[l]* f * 0.03; /* C4H10 -> C3H6 + CH4 */
318	KRPD[lat][j][40][s] += sC4H10[l]* f * 0.07; /* C4H10 -> C3H6 + CH3 + H */
319	KRPD[lat][j][41][s] += sC4H10[l]* f * 0.00; /* C4H10 -> C2H6 + C2H4 */
320	KRPD[lat][j][42][s] += sC4H10[l]* f * 0.15; /* C4H10 -> C2H6 + C2H2 + H2 */
321	KRPD[lat][j][43][s] += sC4H10[l]* f * 0.27; /* C4H10 -> CH3 + C3H7 */
322	KRPD[lat][j][44][s] += sC4H10[l]* f * 0.10; /* C4H10 -> 2 C2H5 */
323	}
324	else
325	{
326	KRPD[lat][j][37][s] += sC4H10[l]* f * 0.41; /* C4H10 -> C4H8(ieC3H5+CH3)+H2 */
327	KRPD[lat][j][38][s] += sC4H10[l]* f * 0.12; /* C4H10 -> 2 C2H4 + H2 */
328	KRPD[lat][j][39][s] += sC4H10[l]* f * 0.01; /* C4H10 -> C3H6 + CH4 */
329	KRPD[lat][j][40][s] += sC4H10[l]* f * 0.07; /* C4H10 -> C3H6 + CH3 + H */
330	KRPD[lat][j][41][s] += sC4H10[l]* f * 0.02; /* C4H10 -> C2H6 + C2H4 */
331	KRPD[lat][j][42][s] += sC4H10[l]* f * 0.06; /* C4H10 -> C2H6 + C2H2 + H2 */
332	KRPD[lat][j][43][s] += sC4H10[l]* f * 0.24; /* C4H10 -> CH3 + C3H7 */
333	KRPD[lat][j][44][s] += sC4H10[l]* f * 0.07; /* C4H10 -> 2 C2H5 */
334	}
335	if( i < 150 )
336	{
337	KRPD[lat][j][ 5][s] += sC2H2[l] * f * 0.3; /* C2H2 -> C2H + H */
338	KRPD[lat][j][ 6][s] += sC2H2[l] * f * 0.1; /* C2H2 -> C2 + H2 */
339	KRPD[lat][j][49][s] += sHC3N[l] * f * 0.3; /* HC3N -> C2H + CN */
340	KRPD[lat][j][50][s] += sHC3N[l] * f * 0.09; /* HC3N -> H + C3N */
341	}
342	else if( i < 205 )
343	{
344	KRPD[lat][j][ 5][s] += sC2H2[l] * f * 0.08; /* C2H2 -> C2H + H */
345	KRPD[lat][j][ 6][s] += sC2H2[l] * f * 0.1; /* C2H2 -> C2 + H2 */
346	KRPD[lat][j][49][s] += sHC3N[l] * f * 0.05; /* HC3N -> C2H + CN */
347	KRPD[lat][j][50][s] += sHC3N[l] * f * 0.09; /* HC3N -> H + C3N */
348	}
349	else if( i < 245 )
350	{
351	KRPD[lat][j][50][s] += sHC3N[l] * f * 0.09; /* HC3N -> H + C3N */
352	}
353	if( i < 165 )
354	{
355	KRPD[lat][j][28][s] += sC4H2[l] * f * 0.2; /* C4H2 -> C4H + H */
356	KRPD[lat][j][29][s] += sC4H2[l] * f * 0.03; /* C4H2 -> 2C2H */
357	KRPD[lat][j][30][s] += sC4H2[l] * f * 0.1; /* C4H2 -> C2H2 + C2 */
358	KRPD[lat][j][31][s] += sC4H2[l] * f * 0.67; /* C4H2 -> C4H2* */
359	}
360	else if( i < 205 )
361	{
362	KRPD[lat][j][29][s] += sC4H2[l] * f * 0.01; /* C4H2 -> 2C2H */
363	KRPD[lat][j][30][s] += sC4H2[l] * f * 0.06; /* C4H2 -> C2H2 + C2 */
364	KRPD[lat][j][31][s] += sC4H2[l] * f * 0.93; /* C4H2 -> C4H2* */
365	}
366	else
367	{
368	KRPD[lat][j][31][s] += sC4H2[l] * f * 1.00; /* C4H2 -> C4H2* */
369	}
370	if( i < 190 ) KRPD[lat][j][14][s] += 4.e-17 * f; /* C3H3 -> C3H2 + H */
371	}
372	}
373
374	/* taux de dissociation de N2 ( e- et GCR ) */
375
376	for( lat = 0; lat <= (*NLAT)-1; lat++ )
377	for( s = 0; s <= 14; s++ )
378	{
379	for( j = 99; j <= NLRT-1; j++ ) /* level 100 = 200 km */
380	KRPD[lat][j][RDISS][s] = 1.0e-16;
381	for( j = 49; j <= 98; j++ ) /* level 50 = 100 km */
382	KRPD[lat][j][RDISS][s] = 1.0e-17+1.8e-18*(j-49);
383	for( j = 34; j <= 48; j++ ) /* level 35 = 70 km */
384	KRPD[lat][j][RDISS][s] = pow(10.,(-23+0.4*(j-34)));
385	for( j = 0; j <= 33; j++ )
386	KRPD[lat][j][RDISS][s] = 0.0e0;
387	}
388
389	fdm2d( flux, 0, NLEV, 0 );
390	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: