Context Navigation

chemparam_mod.F90 @ 1518

Last change on this file since 1518 was 1468, checked in by slebonnois, 10 years ago
SL: bug corrections
Property svn:executable set to ``*
File size: 25.9 KB

Line
1	MODULE chemparam_mod
2
3	!MODULE qui définit les indices des traceurs et leurs masses molaires.
4	! utilise aussi pour variables communes nuages/photochimie
5
6	IMPLICIT NONE
7
8	!----------------------------------------------------------------------------
9	INTEGER, SAVE :: i_co2, i_co, i_h2, i_h2o, i_o1d, &
10	i_o, i_o2, i_o2dg, i_o3, i_h, &
11	i_oh, i_ho2, i_h2o2, i_cl, i_clo, &
12	i_cl2, i_hcl, i_hocl, i_clco, i_clco3, &
13	i_cocl2, i_s, i_so, i_so2, i_so3, &
14	i_s2o2, i_ocs, i_hso3, i_h2so4, i_s2, &
15	i_clso2, i_oscl, i_h2oliq, i_h2so4liq, &
16	i_n2
17
18	REAL, DIMENSION(:), SAVE, ALLOCATABLE :: M_tr
19
20
21	!----------------------------------------------------------------------------
22	! number of clouds mode modelized
23	INTEGER, PARAMETER :: nbr_mode = 3
24	INTEGER :: i_cloud
25	INTEGER, SAVE :: cloudmax
26	INTEGER, SAVE :: cloudmin
27	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: R_MEDIAN
28	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: STDDEV
29
30	! K_MASS coefficient correspondant à la partie condensee de chaque mode
31	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: K_MASS
32
33	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: NBRTOT
34	REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: WH2SO4
35	REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: rho_droplet
36	!----------------------------------------------------------------------------
37
38	!----------------------------------------------------------------------------
39	CONTAINS
40	!----------------------------------------------------------------------------
41
42	SUBROUTINE cloud_ini(nbr_lon,nbr_lev)
43
44	!=============================================================
45	! cloud_ini definit le champ 3D des caracteristiques du nuage
46
47	INTEGER :: nbr_lon,nbr_lev,i_lev
48
49	ALLOCATE(NBRTOT(nbr_lon,nbr_lev,nbr_mode))
50	ALLOCATE(R_MEDIAN(nbr_lon,nbr_lev,nbr_mode))
51	ALLOCATE(K_MASS(nbr_lon,nbr_lev,nbr_mode))
52	ALLOCATE(STDDEV(nbr_lon,nbr_lev,nbr_mode))
53	ALLOCATE(WH2SO4(nbr_lon,nbr_lev))
54	ALLOCATE(rho_droplet(nbr_lon,nbr_lev))
55
56	PRINT*,'=========================='
57	PRINT*,'Initialisation cloud layer'
58	PRINT*,'=========================='
59	PRINT*,'nbr_lon',nbr_lon
60	PRINT*,'nbr_lev',nbr_lev
61	PRINT*,'nbr_mode',nbr_mode
62
63	NBRTOT(:,:,:) = 0.0E+0
64	WH2SO4(:,:) = 0.0E+0
65	rho_droplet(:,:) = 0.0E+0
66
67	!=============================================================
68	! Initialisation cloud layer 1
69	!=============================================================
70	! cloudmin et cloudmax niveaux du GCM
71	cloudmin= 18
72	cloudmax= 50
73
74	! radius R_MEDIAN en m (donc *e-6 pour microns)
75
76	R_MEDIAN(:,:,:)=0.0E+0 ! Geometric Average Radius
77	STDDEV(:,:,:)=0.0E+0 ! Geometric Std Deviation
78	K_MASS(:,:,:)=0.0E+0 ! Coeff multimodal
79
80	! ===============================================
81	! Knollenberg & Hunten, 1980 and James et al 1997
82	! ===============================================
83
84	! ===============================================
85	! Initialisation UNIMODALE
86	! ===============================================
87
88	! Lower Haze: mode 1
89	! DO i_lev=cloudmin,20
90	! R_MEDIAN(:,i_lev,1)=0.2e-6
91	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
92	! STDDEV(:,i_lev,1)=1.56
93	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
94	! K_MASS(:,i_lev,1)=1.0
95	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
96	! END DO
97
98	! Lower Cloud: mode 3
99	! DO i_lev=21,23
100	! R_MEDIAN(:,i_lev,1)=3.65e-6
101	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
102	! STDDEV(:,i_lev,1)=1.28
103	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
104	! K_MASS(:,i_lev,1)=1.0
105	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
106	! END DO
107
108	! Middle Cloud: mode 2 prime
109	! DO i_lev=24,28
110	! R_MEDIAN(:,i_lev,1)=1.4e-6
111	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
112	! STDDEV(:,i_lev,1)=1.23
113	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
114	! K_MASS(:,i_lev,1)=1.0
115	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
116	! END DO
117
118	! Upper Cloud: mode 2
119	! DO i_lev=29,35
120	! R_MEDIAN(:,i_lev,1)=1.0e-6
121	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
122	! STDDEV(:,i_lev,1)=1.29
123	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
124	! K_MASS(:,i_lev,1)=1.0
125	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
126	! END DO
127
128	! Upper Haze: mode 1
129	! DO i_lev=36, cloudmax
130	! R_MEDIAN(:,i_lev,1)=0.2e-6
131	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
132	! STDDEV(:,i_lev,1)=2.16
133	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
134	! K_MASS(:,i_lev,1)=1.0
135	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
136	! END DO
137
138	! ===============================================
139	! Initialisation TRIMODALE
140	! ===============================================
141
142	! Lower Haze: mode 1
143	! DO i_lev=cloudmin,20
144	! R_MEDIAN(:,i_lev,1)=0.3e-6
145	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
146	! STDDEV(:,i_lev,1)=1.56
147	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
148	! K_MASS(:,i_lev,1)=1.0
149	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
150	! END DO
151
152	! Lower Haze: mode 2
153	! DO i_lev=cloudmin,20
154	! R_MEDIAN(:,i_lev,2)=1.4e-6
155	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
156	! STDDEV(:,i_lev,2)=1.23
157	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
158	! K_MASS(:,i_lev,2)=0.0
159	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
160	! END DO
161
162	! Lower Haze: mode 3
163	! DO i_lev=cloudmin,20
164	! R_MEDIAN(:,i_lev,3)=3.65e-6
165	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
166	! STDDEV(:,i_lev,3)=1.28
167	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
168	! K_MASS(:,i_lev,3)=0.
169	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
170	! END DO
171
172	! Lower Cloud: mode 1
173	! DO i_lev=21,23
174	! R_MEDIAN(:,i_lev,1)=0.3e-6
175	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
176	! STDDEV(:,i_lev,1)=1.56
177	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
178	! K_MASS(:,i_lev,1)=0.1
179	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
180	! END DO
181
182	! Lower Cloud: mode 2 prime
183	! DO i_lev=21,23
184	! R_MEDIAN(:,i_lev,2)=1.4e-6
185	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
186	! STDDEV(:,i_lev,2)=1.23
187	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
188	! K_MASS(:,i_lev,2)=0.4
189	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
190	! END DO
191
192	! Lower Cloud: mode 3
193	! DO i_lev=21,23
194	! R_MEDIAN(:,i_lev,3)=3.65e-6
195	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
196	! STDDEV(:,i_lev,3)=1.28
197	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
198	! K_MASS(:,i_lev,3)=0.5
199	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
200	! END DO
201
202	! Middle Cloud: mode 1
203	! DO i_lev=24,28
204	! R_MEDIAN(:,i_lev,1)=0.3e-6
205	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
206	! STDDEV(:,i_lev,1)=1.56
207	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
208	! K_MASS(:,i_lev,1)=0.0
209	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
210	! END DO
211
212	! Middle Cloud: mode 2 prime
213	! DO i_lev=24,28
214	! R_MEDIAN(:,i_lev,2)=1.4e-6
215	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
216	! STDDEV(:,i_lev,2)=1.23
217	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
218	! K_MASS(:,i_lev,2)=0.8
219	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
220	! END DO
221
222	! Middle Cloud: mode 3
223	! DO i_lev=24,28
224	! R_MEDIAN(:,i_lev,3)=3.65e-6
225	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
226	! STDDEV(:,i_lev,3)=1.28
227	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
228	! K_MASS(:,i_lev,3)=0.2
229	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
230	! END DO
231
232
233	! Upper Cloud: mode 1
234	! DO i_lev=29,35
235	! R_MEDIAN(:,i_lev,1)=0.3e-6
236	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
237	! STDDEV(:,i_lev,1)=1.56
238	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
239	! K_MASS(:,i_lev,1)=0.15
240	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
241	! END DO
242
243	! Upper Cloud: mode 2
244	! DO i_lev=29,35
245	! R_MEDIAN(:,i_lev,2)=1.0e-6
246	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
247	! STDDEV(:,i_lev,2)=1.29
248	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
249	! K_MASS(:,i_lev,2)=0.85
250	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
251	! END DO
252
253	! Upper Cloud: mode 3
254	! DO i_lev=29,35
255	! R_MEDIAN(:,i_lev,3)=3.65e-6
256	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
257	! STDDEV(:,i_lev,3)=1.28
258	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
259	! K_MASS(:,i_lev,3)=0.0
260	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
261	! END DO
262
263	! Upper Haze: mode 1
264	! DO i_lev=36, cloudmax
265	! R_MEDIAN(:,i_lev,1)=0.3e-6
266	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
267	! STDDEV(:,i_lev,1)=1.56
268	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
269	! K_MASS(:,i_lev,1)=1.0
270	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
271	! END DO
272
273	! Upper Haze: mode 2
274	! DO i_lev=36, cloudmax
275	! R_MEDIAN(:,i_lev,2)=1.e-6
276	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
277	! STDDEV(:,i_lev,2)=1.29
278	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
279	! K_MASS(:,i_lev,2)=0.0
280	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
281	! END DO
282
283	! Upper Haze: mode 3
284	! DO i_lev=36, cloudmax
285	! R_MEDIAN(:,i_lev,3)=3.65e-6
286	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
287	! STDDEV(:,i_lev,3)=2.16
288	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
289	! K_MASS(:,i_lev,3)=0.0
290	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
291	! END DO
292	!=============================================================
293
294	! ===============================================
295	! Initialisation TRIMODALE Knollenberg
296	! ===============================================
297
298	! Lower Haze: mode 1
299	DO i_lev=cloudmin,22
300	R_MEDIAN(:,i_lev,1)=0.1e-6
301	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
302	STDDEV(:,i_lev,1)=1.57
303	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
304	K_MASS(:,i_lev,1)=1.0
305	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
306	END DO
307
308	! Lower Haze: mode 2
309	DO i_lev=cloudmin,22
310	R_MEDIAN(:,i_lev,2)=1.4e-6
311	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
312	STDDEV(:,i_lev,2)=1.23
313	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
314	K_MASS(:,i_lev,2)=0.0
315	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
316	END DO
317
318	! Lower Haze: mode 3
319	DO i_lev=cloudmin,22
320	R_MEDIAN(:,i_lev,3)=3.65e-6
321	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
322	STDDEV(:,i_lev,3)=1.28
323	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
324	K_MASS(:,i_lev,3)=0.0
325	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
326	END DO
327
328	! Pre Cloud: mode 1
329	DO i_lev=23,23
330	R_MEDIAN(:,i_lev,1)=0.15e-6
331	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
332	STDDEV(:,i_lev,1)=1.8
333	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
334	K_MASS(:,i_lev,1)=0.04
335	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
336	END DO
337
338	! Pre Cloud: mode 2
339	DO i_lev=23,23
340	R_MEDIAN(:,i_lev,2)=1.0e-6
341	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
342	STDDEV(:,i_lev,2)=1.29
343	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
344	K_MASS(:,i_lev,2)=0.96
345	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
346	END DO
347
348	! Pre Cloud: mode 3
349	DO i_lev=23,23
350	R_MEDIAN(:,i_lev,3)=3.65e-6
351	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
352	STDDEV(:,i_lev,3)=1.28
353	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
354	K_MASS(:,i_lev,3)=0.0
355	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
356	END DO
357
358	! Lower Cloud: mode 1
359	DO i_lev=24,24
360	R_MEDIAN(:,i_lev,1)=0.2e-6
361	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
362	STDDEV(:,i_lev,1)=1.8
363	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
364	K_MASS(:,i_lev,1)=0.014
365	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
366	END DO
367
368	! Lower Cloud: mode 2
369	DO i_lev=24,24
370	R_MEDIAN(:,i_lev,2)=1.0e-6
371	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
372	STDDEV(:,i_lev,2)=1.29
373	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
374	K_MASS(:,i_lev,2)=0.02
375	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
376	END DO
377
378	! Lower Cloud: mode 3
379	DO i_lev=24,24
380	R_MEDIAN(:,i_lev,3)=3.65e-6
381	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
382	STDDEV(:,i_lev,3)=1.28
383	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
384	K_MASS(:,i_lev,3)=0.966
385	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
386	END DO
387
388	! Middle Cloud: mode 1
389	DO i_lev=25,28
390	R_MEDIAN(:,i_lev,1)=0.15e-6
391	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
392	STDDEV(:,i_lev,1)=1.9
393	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
394	K_MASS(:,i_lev,1)=0.0084
395	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
396	END DO
397
398	! Middle Cloud: mode 2 prime
399	DO i_lev=25,28
400	R_MEDIAN(:,i_lev,2)=1.4e-6
401	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
402	STDDEV(:,i_lev,2)=1.23
403	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
404	K_MASS(:,i_lev,2)=0.21
405	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
406	END DO
407
408	! Middle Cloud: mode 3
409	DO i_lev=25,28
410	R_MEDIAN(:,i_lev,3)=3.65e-6
411	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
412	STDDEV(:,i_lev,3)=1.28
413	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
414	K_MASS(:,i_lev,3)=0.7816
415	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
416	END DO
417
418
419	! Upper Cloud: mode 1
420	DO i_lev=29,35
421	R_MEDIAN(:,i_lev,1)=0.2e-6
422	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
423	STDDEV(:,i_lev,1)=2.16
424	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
425	K_MASS(:,i_lev,1)=0.72
426	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
427	END DO
428
429	! Upper Cloud: mode 2
430	DO i_lev=29,35
431	R_MEDIAN(:,i_lev,2)=1.0e-6
432	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
433	STDDEV(:,i_lev,2)=1.29
434	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
435	K_MASS(:,i_lev,2)=0.28
436	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
437	END DO
438
439	! Upper Cloud: mode 3
440	DO i_lev=29,35
441	R_MEDIAN(:,i_lev,3)=3.65e-6
442	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
443	STDDEV(:,i_lev,3)=1.28
444	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
445	K_MASS(:,i_lev,3)=0.0
446	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
447	END DO
448
449	! Upper Haze: mode 1
450	DO i_lev=36, cloudmax
451	R_MEDIAN(:,i_lev,1)=0.2e-6
452	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
453	STDDEV(:,i_lev,1)=2.16
454	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
455	K_MASS(:,i_lev,1)=1.0
456	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
457	END DO
458
459	! Upper Haze: mode 2
460	DO i_lev=36, cloudmax
461	R_MEDIAN(:,i_lev,2)=1.e-6
462	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
463	STDDEV(:,i_lev,2)=1.29
464	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
465	K_MASS(:,i_lev,2)=0.0
466	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
467	END DO
468
469	! Upper Haze: mode 3
470	DO i_lev=36, cloudmax
471	R_MEDIAN(:,i_lev,3)=3.65e-6
472	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
473	STDDEV(:,i_lev,3)=2.16
474	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
475	K_MASS(:,i_lev,3)=0.0
476	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
477	END DO
478
479	!=============================================================
480
481	! ===============================================================
482	! Initialisation TRIMODALE "Knollenberg" sans Mode3, Mode2 etendu
483	! ===============================================================
484
485	! Lower Haze: mode 1
486	! DO i_lev=cloudmin,22
487	! R_MEDIAN(:,i_lev,1)=0.1e-6
488	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
489	! STDDEV(:,i_lev,1)=1.57
490	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
491	! K_MASS(:,i_lev,1)=1.0
492	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
493	! END DO
494
495	! Lower Haze: mode 2
496	! DO i_lev=cloudmin,22
497	! R_MEDIAN(:,i_lev,2)=1.4e-6
498	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
499	! STDDEV(:,i_lev,2)=1.23
500	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
501	! K_MASS(:,i_lev,2)=0.0
502	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
503	! END DO
504
505	! Lower Haze: mode 3
506	! DO i_lev=cloudmin,22
507	! R_MEDIAN(:,i_lev,3)=3.65e-6
508	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
509	! STDDEV(:,i_lev,3)=1.28
510	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
511	! K_MASS(:,i_lev,3)=0.0
512	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
513	! END DO
514
515	! Pre Cloud: mode 1
516	! DO i_lev=23,23
517	! R_MEDIAN(:,i_lev,1)=0.15e-6
518	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
519	! STDDEV(:,i_lev,1)=1.8
520	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
521	! K_MASS(:,i_lev,1)=0.04
522	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
523	! END DO
524
525	! Pre Cloud: mode 2
526	! DO i_lev=23,23
527	! R_MEDIAN(:,i_lev,2)=1.0e-6
528	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
529	! STDDEV(:,i_lev,2)=1.29
530	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
531	! K_MASS(:,i_lev,2)=0.96
532	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
533	! END DO
534
535	! Pre Cloud: mode 3
536	! DO i_lev=23,23
537	! R_MEDIAN(:,i_lev,3)=3.65e-6
538	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
539	! STDDEV(:,i_lev,3)=1.28
540	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
541	! K_MASS(:,i_lev,3)=0.0
542	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
543	! END DO
544
545	! Lower Cloud: mode 1
546	! DO i_lev=24,24
547	! R_MEDIAN(:,i_lev,1)=0.2e-6
548	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
549	! STDDEV(:,i_lev,1)=1.8
550	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
551	! K_MASS(:,i_lev,1)=0.014
552	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
553	! END DO
554
555	! Lower Cloud: mode 2
556	! DO i_lev=24,24
557	! R_MEDIAN(:,i_lev,2)=1.0e-6
558	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
559	! STDDEV(:,i_lev,2)=1.6
560	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
561	! K_MASS(:,i_lev,2)=0.986
562	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
563	! END DO
564
565	! Lower Cloud: mode 3
566	! DO i_lev=24,24
567	! R_MEDIAN(:,i_lev,3)=3.65e-6
568	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
569	! STDDEV(:,i_lev,3)=1.28
570	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
571	! K_MASS(:,i_lev,3)=0.
572	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
573	! END DO
574
575	! Middle Cloud: mode 1
576	! DO i_lev=25,28
577	! R_MEDIAN(:,i_lev,1)=0.15e-6
578	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
579	! STDDEV(:,i_lev,1)=1.9
580	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
581	! K_MASS(:,i_lev,1)=0.0084
582	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
583	! END DO
584
585	! Middle Cloud: mode 2 prime
586	! DO i_lev=25,28
587	! R_MEDIAN(:,i_lev,2)=1.4e-6
588	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
589	! STDDEV(:,i_lev,2)=1.6
590	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
591	! K_MASS(:,i_lev,2)=0.9916
592	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
593	! END DO
594
595	! Middle Cloud: mode 3
596	! DO i_lev=25,28
597	! R_MEDIAN(:,i_lev,3)=3.65e-6
598	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
599	! STDDEV(:,i_lev,3)=1.28
600	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
601	! K_MASS(:,i_lev,3)=0.0
602	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
603	! END DO
604
605
606	! Upper Cloud: mode 1
607	! DO i_lev=29,35
608	! R_MEDIAN(:,i_lev,1)=0.2e-6
609	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
610	! STDDEV(:,i_lev,1)=2.16
611	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
612	! K_MASS(:,i_lev,1)=0.72
613	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
614	! END DO
615
616	! Upper Cloud: mode 2
617	! DO i_lev=29,35
618	! R_MEDIAN(:,i_lev,2)=1.0e-6
619	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
620	! STDDEV(:,i_lev,2)=1.29
621	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
622	! K_MASS(:,i_lev,2)=0.28
623	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
624	! END DO
625
626	! Upper Cloud: mode 3
627	! DO i_lev=29,35
628	! R_MEDIAN(:,i_lev,3)=3.65e-6
629	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
630	! STDDEV(:,i_lev,3)=1.28
631	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
632	! K_MASS(:,i_lev,3)=0.0
633	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
634	! END DO
635
636	! Upper Haze: mode 1
637	! DO i_lev=36, cloudmax
638	! R_MEDIAN(:,i_lev,1)=0.2e-6
639	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
640	! STDDEV(:,i_lev,1)=2.16
641	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
642	! K_MASS(:,i_lev,1)=1.0
643	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
644	! END DO
645
646	! Upper Haze: mode 2
647	! DO i_lev=36, cloudmax
648	! R_MEDIAN(:,i_lev,2)=1.e-6
649	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
650	! STDDEV(:,i_lev,2)=1.29
651	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
652	! K_MASS(:,i_lev,2)=0.0
653	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
654	! END DO
655
656	! Upper Haze: mode 3
657	! DO i_lev=36, cloudmax
658	! R_MEDIAN(:,i_lev,3)=3.65e-6
659	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
660	! STDDEV(:,i_lev,3)=2.16
661	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
662	! K_MASS(:,i_lev,3)=0.0
663	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
664	! END DO
665	!=============================================================
666	PRINT*,'==============================='
667	PRINT*,'FIN Initialisation cloud layer'
668	PRINT*,'==============================='
669
670	END SUBROUTINE cloud_ini
671
672	SUBROUTINE chemparam_ini
673	USE infotrac
674	IMPLICIT NONE
675	INTEGER :: i
676
677	ALLOCATE(M_tr(nqtot))
678
679	DO i=1, nqtot
680
681	PRINT*,'i',i
682	PRINT*,'tname(i)',tname(i)
683
684	SELECT CASE(tname(i))
685	CASE('co2')
686	i_co2=i
687	PRINT*,'co2',i_co2
688	M_tr(i_co2) = 44.0095
689	CASE('co')
690	i_co=i
691	PRINT*,'co',i_co
692	M_tr(i_co)=28.0101
693	CASE('h2')
694	i_h2=i
695	PRINT*,'h2',i_h2
696	M_tr(i_h2)= 2.01588
697	CASE('h2o')
698	i_h2o=i
699	PRINT*,'h2o',i_h2o
700	M_tr(i_h2o)=18.0153
701	CASE('o1d')
702	i_o1d=i
703	PRINT*,'o1d',i_o1d
704	M_tr(i_o1d)=15.994
705	CASE('o')
706	i_o=i
707	PRINT*,'o',i_o
708	M_tr(i_o)=15.994
709	CASE('o2')
710	i_o2=i
711	PRINT*,'o2',i_o2
712	M_tr(i_o2)=31.9988
713	CASE('o2dg')
714	i_o2dg=i
715	PRINT*,'o2dg',i_o2dg
716	M_tr(i_o2dg)=31.9988
717	CASE('o3')
718	i_o3=i
719	PRINT*,'o3',i_o3
720	M_tr(i_o3)= 47.9982
721	CASE('h')
722	i_h=i
723	PRINT*,'h',i_h
724	M_tr(i_h)= 1.00794
725	CASE('oh')
726	i_oh=i
727	PRINT*,'oh',i_oh
728	M_tr(i_oh)=17.0073
729	CASE('ho2')
730	i_ho2=i
731	PRINT*,'ho2',i_ho2
732	M_tr(i_ho2)=33.0067
733	CASE('h2o2')
734	i_h2o2=i
735	PRINT*,'h2o2',i_h2o2
736	M_tr(i_h2o2)=34.0147
737	CASE('cl')
738	i_cl=i
739	PRINT*,'cl',i_cl
740	M_tr(i_cl)=35.453
741	CASE('clo')
742	i_clo=i
743	PRINT*,'clo',i_clo
744	M_tr(i_clo)=51.452
745	CASE('cl2')
746	i_cl2=i
747	PRINT*,'cl2',i_cl2
748	M_tr(i_cl2)=70.906
749	CASE('hcl')
750	i_hcl=i
751	PRINT*,'hcl',i_hcl
752	M_tr(i_hcl)=36.461
753	CASE('hocl')
754	i_hocl=i
755	PRINT*,'hocl',i_hocl
756	M_tr(i_hocl)=52.46
757	CASE('clco')
758	i_clco=i
759	PRINT*,'clco',i_clco
760	M_tr(i_clco)=63.463
761	CASE('clco3')
762	i_clco3=i
763	PRINT*,'clco3',i_clco3
764	M_tr(i_clco3)=95.462
765	CASE('cocl2')
766	i_cocl2=i
767	PRINT*,'cocl2',i_cocl2
768	M_tr(i_cocl2)=98.916
769	CASE('s')
770	i_s=i
771	PRINT*,'s',i_s
772	M_tr(i_s)=32.065
773	CASE('so')
774	i_so=i
775	PRINT*,'so',i_so
776	M_tr(i_so)=48.0644
777	CASE('so2')
778	i_so2=i
779	PRINT*,'so2',i_so2
780	M_tr(i_so2)=64.064
781	CASE('so3')
782	i_so3=i
783	PRINT*,'so3',i_so3
784	M_tr(i_so3)=80.063
785	CASE('s2o2')
786	i_s2o2=i
787	PRINT*,'s2o2',i_s2o2
788	M_tr(i_s2o2)= 96.1288
789	CASE('ocs')
790	i_ocs=i
791	PRINT*,'ocs',i_ocs
792	M_tr(i_ocs)=60.0751
793	CASE('hso3')
794	i_hso3=i
795	PRINT*,'hso3',i_hso3
796	M_tr(i_hso3)=81.071
797	CASE('h2so4')
798	i_h2so4=i
799	PRINT*,'h2so4',i_h2so4
800	M_tr(i_h2so4)=98.078
801	CASE('s2')
802	i_s2=i
803	PRINT*,'s2',i_s2
804	M_tr(i_s2)=64.13
805	CASE('clso2')
806	i_clso2=i
807	PRINT*,'clso2',i_clso2
808	M_tr(i_clso2)=99.517
809	CASE('oscl')
810	i_oscl=i
811	PRINT*,'oscl',i_oscl
812	M_tr(i_oscl)=83.517
813	CASE('h2oliq')
814	i_h2oliq=i
815	PRINT*,'h2oliq',i_h2oliq
816	M_tr(i_h2oliq)=18.0153
817	CASE('h2so4liq')
818	i_h2so4liq=i
819	PRINT*,'h2so4liq',i_h2so4liq
820	M_tr(i_h2so4liq)=98.078
821	CASE('n2')
822	i_n2=i
823	M_tr(i_n2)=28.013
824	END SELECT
825
826	! PRINT*,'M_tr(i)',M_tr(i)
827	END DO
828
829	END SUBROUTINE chemparam_ini
830
831	END MODULE chemparam_mod
832

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

Download in other formats:

Original Format