Context Navigation

chemparam_mod.F90 @ 3094

Last change on this file since 3094 was 2836, checked in by abierjon, 2 years ago

VENUS GCM:

INCLUDING THE IONOSPHERE CODE IN VENUS GCM

ATTENTION: INCREASED MEMORY DEMAND

NEEDS AT LEAST 135 GB of allocated memory

==============================================================
===> LIST OF APPENDED FILES AND INTERNAL ADDITIONS <==========
==============================================================

NEW MODEL SPECIES

deftank/traceur-chemistry-IONOSPHERE.def

Neutrals: 4 Species: N, NO, NO2, N(2D)

Ions: 15 species:

CO2+, CO+, O+, O2+, H+,

N2+, H2O+, OH+, C+, HCO+,

H3O+, HCO2+, N+, NO+, elec

FROM 36 to 55 chemical species

NEW KEYWORD OF PHYSIQ.DEF

deftank/physiq-96x96x90-chemistry-IONOSPHERE.def

ok_ionchem: keyword supposed to activate ion chemistry.

(be careful that n, no, no2, n2d and the ion species are in the deftank/traceur-chemistry-IONOSPHERE.def)

ok_jonline: keyword supposed to activate the online photochemistry

(be careful that n, no, no2 and n2d are in the traceur-chemistry-IONOSPHERE.def)

==============================================================
===> LIST OF MODIFIED PROGRAMS <=========================
==============================================================

nonoro_gwd_ran_mod.F90

Change EPFLUXMAX value from 5.E-3 to 1.E-3

photochemistry_venus.F90 (krates; photolysis_online; indices;)

Import of keywords: ok_ionchem, tuneupperatm & ok_jonline
addition of ion species in order
Forcing electroneutrality
Update of the reactions a001 and a002 with taking into account the other species
Change of formula for a002 with the formulation of Baulch et al., 1976 (confirmed by Smith and Robertson, 2008)

photolysis_mod.F90

Modification of the rdsolarflux subroutine to include interpolation with Atlas1 and Atlas3 in connection with E10.7

concentration2.F90

Added chemical species n2d, no, no2 and n in the conductivity calculation.

The ions have been excluded because their sum is 10^{5 times less dense than the neutrals and

their thermal conductivity is unknown}

iono_h.F90

Addition of the phdisrate routine
replace the electronic temperature of Mars by that of

origin = 1: Theis et al. 1980 (Venus) with bilinear interpolation altitude/cos(SZA)
origin = 2: Theis et al. 1984 (Venus) with the formula of the electronic temperature at high solar activity

addition of an ion temperature model based on VIRA

cleshphys.h

added ok_ionchem & ok_jonline in COMMON/clesphys_l/

conf_phys.f90

add ok_ionchem & ok_jonline as parameters to read from physiq.def file set to .false. by default

chemparam_mod.F90

Add species in M_tr and corresponding i_X. Set all i_X to zero before reading traceur-chemistry-IONOSPHERE.def
Added Type_tr table to differentiate species: 1 == neutral, 2 == ION, 3 == liquid, 10 == others

euvheat.F90; hrtherm.F; jthermcalc_e107.F; param_read_e107.F

Normalization with Mars

A.M

Property svn:executable set to *

File size: 33.6 KB

Line
1	MODULE chemparam_mod
2
3	!MODULE qui definit 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_n2, i_he, i_n, i_no, &
16	i_no2, i_n2d, &
17	i_co2plus, i_coplus, i_oplus, i_o2plus, &
18	i_n2plus, i_hplus, i_h2oplus, i_nplus, &
19	i_ohplus, i_cplus, i_noplus, i_h3oplus, &
20	i_hcoplus, i_hco2plus, i_elec
21
22	INTEGER, SAVE :: i_h2oliq, i_h2so4liq
23
24	INTEGER, SAVE :: i_m0_aer, i_m3_aer, &
25	i_m0_mode1drop, i_m0_mode1ccn, &
26	i_m3_mode1sa, i_m3_mode1w, i_m3_mode1ccn, &
27	i_m0_mode2drop, i_m0_mode2ccn, &
28	i_m3_mode2sa, i_m3_mode2w, i_m3_mode2ccn
29
30	INTEGER, SAVE :: nmicro ! number of microphysical tracers
31
32	REAL, DIMENSION(:), SAVE, ALLOCATABLE :: M_tr ! Molecular Mass of tracers
33	REAL, DIMENSION(:), SAVE, ALLOCATABLE :: type_tr ! type of tracer
34
35	REAL, DIMENSION(:,:),SAVE, ALLOCATABLE :: no_emission
36	REAL, DIMENSION(:,:),SAVE, ALLOCATABLE :: o2_emission
37	!----------------------------------------------------------------------------
38	! DEF FOR CL_SCHEME = 1 (AURELIEN)
39
40	! number of clouds mode modelized
41	INTEGER, PARAMETER :: nbr_mode = 3
42	INTEGER :: i_cloud
43	INTEGER, SAVE :: cloudmax
44	INTEGER, SAVE :: cloudmin
45	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: R_MEDIAN
46	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: STDDEV
47
48	! K_MASS coefficient correspondant à la partie condensee de chaque mode
49	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: K_MASS
50
51	REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: NBRTOT
52	REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: WH2SO4
53	REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: rho_droplet
54	!----------------------------------------------------------------------------
55	! DEF FOR CL_SCHEME = 2 (FULL MICROPHYS)
56
57	!----------------------------------------------------------------------------
58
59	!----------------------------------------------------------------------------
60	CONTAINS
61	!----------------------------------------------------------------------------
62
63	SUBROUTINE cloud_ini(nbr_lon,nbr_lev)
64
65	!=============================================================
66	! cloud_ini definit le champ 3D des caracteristiques du nuage
67
68	INTEGER :: nbr_lon,nbr_lev,i_lev
69
70	ALLOCATE(NBRTOT(nbr_lon,nbr_lev,nbr_mode))
71	ALLOCATE(R_MEDIAN(nbr_lon,nbr_lev,nbr_mode))
72	ALLOCATE(K_MASS(nbr_lon,nbr_lev,nbr_mode))
73	ALLOCATE(STDDEV(nbr_lon,nbr_lev,nbr_mode))
74	ALLOCATE(WH2SO4(nbr_lon,nbr_lev))
75	ALLOCATE(rho_droplet(nbr_lon,nbr_lev))
76
77	PRINT*,'=========================='
78	PRINT*,'Initialisation cloud layer'
79	PRINT*,'=========================='
80	PRINT*,'nbr_lon',nbr_lon
81	PRINT*,'nbr_lev',nbr_lev
82	PRINT*,'nbr_mode',nbr_mode
83
84	NBRTOT(:,:,:) = 0.0E+0
85	WH2SO4(:,:) = 0.0E+0
86	rho_droplet(:,:) = 0.0E+0
87
88	!=============================================================
89	! Initialisation cloud layer 1
90	!=============================================================
91	! cloudmin et cloudmax niveaux du GCM
92	cloudmin= 18
93	cloudmax= 50
94
95	! radius R_MEDIAN en m (donc *e-6 pour microns)
96
97	R_MEDIAN(:,:,:)=0.0E+0 ! Geometric Average Radius
98	STDDEV(:,:,:)=0.0E+0 ! Geometric Std Deviation
99	K_MASS(:,:,:)=0.0E+0 ! Coeff multimodal
100
101	! ===============================================
102	! Knollenberg & Hunten, 1980 and James et al 1997
103	! ===============================================
104
105	! ===============================================
106	! Initialisation UNIMODALE
107	! ===============================================
108
109	! Lower Haze: mode 1
110	! DO i_lev=cloudmin,20
111	! R_MEDIAN(:,i_lev,1)=0.2e-6
112	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
113	! STDDEV(:,i_lev,1)=1.56
114	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
115	! K_MASS(:,i_lev,1)=1.0
116	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
117	! END DO
118
119	! Lower Cloud: mode 3
120	! DO i_lev=21,23
121	! R_MEDIAN(:,i_lev,1)=3.65e-6
122	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
123	! STDDEV(:,i_lev,1)=1.28
124	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
125	! K_MASS(:,i_lev,1)=1.0
126	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
127	! END DO
128
129	! Middle Cloud: mode 2 prime
130	! DO i_lev=24,28
131	! R_MEDIAN(:,i_lev,1)=1.4e-6
132	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
133	! STDDEV(:,i_lev,1)=1.23
134	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
135	! K_MASS(:,i_lev,1)=1.0
136	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
137	! END DO
138
139	! Upper Cloud: mode 2
140	! DO i_lev=29,35
141	! R_MEDIAN(:,i_lev,1)=1.0e-6
142	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
143	! STDDEV(:,i_lev,1)=1.29
144	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
145	! K_MASS(:,i_lev,1)=1.0
146	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
147	! END DO
148
149	! Upper Haze: mode 1
150	! DO i_lev=36, cloudmax
151	! R_MEDIAN(:,i_lev,1)=0.2e-6
152	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
153	! STDDEV(:,i_lev,1)=2.16
154	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
155	! K_MASS(:,i_lev,1)=1.0
156	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
157	! END DO
158
159	! ===============================================
160	! Initialisation TRIMODALE
161	! ===============================================
162
163	! Lower Haze: mode 1
164	! DO i_lev=cloudmin,20
165	! R_MEDIAN(:,i_lev,1)=0.3e-6
166	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
167	! STDDEV(:,i_lev,1)=1.56
168	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
169	! K_MASS(:,i_lev,1)=1.0
170	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
171	! END DO
172
173	! Lower Haze: mode 2
174	! DO i_lev=cloudmin,20
175	! R_MEDIAN(:,i_lev,2)=1.4e-6
176	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
177	! STDDEV(:,i_lev,2)=1.23
178	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
179	! K_MASS(:,i_lev,2)=0.0
180	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
181	! END DO
182
183	! Lower Haze: mode 3
184	! DO i_lev=cloudmin,20
185	! R_MEDIAN(:,i_lev,3)=3.65e-6
186	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
187	! STDDEV(:,i_lev,3)=1.28
188	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
189	! K_MASS(:,i_lev,3)=0.
190	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
191	! END DO
192
193	! Lower Cloud: mode 1
194	! DO i_lev=21,23
195	! R_MEDIAN(:,i_lev,1)=0.3e-6
196	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
197	! STDDEV(:,i_lev,1)=1.56
198	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
199	! K_MASS(:,i_lev,1)=0.1
200	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
201	! END DO
202
203	! Lower Cloud: mode 2 prime
204	! DO i_lev=21,23
205	! R_MEDIAN(:,i_lev,2)=1.4e-6
206	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
207	! STDDEV(:,i_lev,2)=1.23
208	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
209	! K_MASS(:,i_lev,2)=0.4
210	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
211	! END DO
212
213	! Lower Cloud: mode 3
214	! DO i_lev=21,23
215	! R_MEDIAN(:,i_lev,3)=3.65e-6
216	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
217	! STDDEV(:,i_lev,3)=1.28
218	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
219	! K_MASS(:,i_lev,3)=0.5
220	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
221	! END DO
222
223	! Middle Cloud: mode 1
224	! DO i_lev=24,28
225	! R_MEDIAN(:,i_lev,1)=0.3e-6
226	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
227	! STDDEV(:,i_lev,1)=1.56
228	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
229	! K_MASS(:,i_lev,1)=0.0
230	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
231	! END DO
232
233	! Middle Cloud: mode 2 prime
234	! DO i_lev=24,28
235	! R_MEDIAN(:,i_lev,2)=1.4e-6
236	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
237	! STDDEV(:,i_lev,2)=1.23
238	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
239	! K_MASS(:,i_lev,2)=0.8
240	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
241	! END DO
242
243	! Middle Cloud: mode 3
244	! DO i_lev=24,28
245	! R_MEDIAN(:,i_lev,3)=3.65e-6
246	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
247	! STDDEV(:,i_lev,3)=1.28
248	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
249	! K_MASS(:,i_lev,3)=0.2
250	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
251	! END DO
252
253
254	! Upper Cloud: mode 1
255	! DO i_lev=29,35
256	! R_MEDIAN(:,i_lev,1)=0.3e-6
257	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
258	! STDDEV(:,i_lev,1)=1.56
259	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
260	! K_MASS(:,i_lev,1)=0.15
261	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
262	! END DO
263
264	! Upper Cloud: mode 2
265	! DO i_lev=29,35
266	! R_MEDIAN(:,i_lev,2)=1.0e-6
267	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
268	! STDDEV(:,i_lev,2)=1.29
269	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
270	! K_MASS(:,i_lev,2)=0.85
271	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
272	! END DO
273
274	! Upper Cloud: mode 3
275	! DO i_lev=29,35
276	! R_MEDIAN(:,i_lev,3)=3.65e-6
277	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
278	! STDDEV(:,i_lev,3)=1.28
279	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
280	! K_MASS(:,i_lev,3)=0.0
281	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
282	! END DO
283
284	! Upper Haze: mode 1
285	! DO i_lev=36, cloudmax
286	! R_MEDIAN(:,i_lev,1)=0.3e-6
287	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
288	! STDDEV(:,i_lev,1)=1.56
289	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
290	! K_MASS(:,i_lev,1)=1.0
291	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
292	! END DO
293
294	! Upper Haze: mode 2
295	! DO i_lev=36, cloudmax
296	! R_MEDIAN(:,i_lev,2)=1.e-6
297	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
298	! STDDEV(:,i_lev,2)=1.29
299	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
300	! K_MASS(:,i_lev,2)=0.0
301	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
302	! END DO
303
304	! Upper Haze: mode 3
305	! DO i_lev=36, cloudmax
306	! R_MEDIAN(:,i_lev,3)=3.65e-6
307	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
308	! STDDEV(:,i_lev,3)=2.16
309	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
310	! K_MASS(:,i_lev,3)=0.0
311	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
312	! END DO
313	!=============================================================
314
315	! ===============================================
316	! Initialisation TRIMODALE Knollenberg
317	! ===============================================
318
319	! Lower Haze: mode 1
320	DO i_lev=cloudmin,22
321	R_MEDIAN(:,i_lev,1)=0.1e-6
322	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
323	STDDEV(:,i_lev,1)=1.57
324	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
325	K_MASS(:,i_lev,1)=1.0
326	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
327	END DO
328
329	! Lower Haze: mode 2
330	DO i_lev=cloudmin,22
331	R_MEDIAN(:,i_lev,2)=1.4e-6
332	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
333	STDDEV(:,i_lev,2)=1.23
334	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
335	K_MASS(:,i_lev,2)=0.0
336	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
337	END DO
338
339	! Lower Haze: mode 3
340	DO i_lev=cloudmin,22
341	R_MEDIAN(:,i_lev,3)=3.65e-6
342	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
343	STDDEV(:,i_lev,3)=1.28
344	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
345	K_MASS(:,i_lev,3)=0.0
346	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
347	END DO
348
349	! Pre Cloud: mode 1
350	DO i_lev=23,23
351	R_MEDIAN(:,i_lev,1)=0.15e-6
352	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
353	STDDEV(:,i_lev,1)=1.8
354	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
355	K_MASS(:,i_lev,1)=0.04
356	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
357	END DO
358
359	! Pre Cloud: mode 2
360	DO i_lev=23,23
361	R_MEDIAN(:,i_lev,2)=1.0e-6
362	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
363	STDDEV(:,i_lev,2)=1.29
364	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
365	K_MASS(:,i_lev,2)=0.96
366	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
367	END DO
368
369	! Pre Cloud: mode 3
370	DO i_lev=23,23
371	R_MEDIAN(:,i_lev,3)=3.65e-6
372	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
373	STDDEV(:,i_lev,3)=1.28
374	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
375	K_MASS(:,i_lev,3)=0.0
376	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
377	END DO
378
379	! Lower Cloud: mode 1
380	DO i_lev=24,24
381	R_MEDIAN(:,i_lev,1)=0.2e-6
382	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
383	STDDEV(:,i_lev,1)=1.8
384	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
385	K_MASS(:,i_lev,1)=0.014
386	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
387	END DO
388
389	! Lower Cloud: mode 2
390	DO i_lev=24,24
391	R_MEDIAN(:,i_lev,2)=1.0e-6
392	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
393	STDDEV(:,i_lev,2)=1.29
394	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
395	K_MASS(:,i_lev,2)=0.02
396	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
397	END DO
398
399	! Lower Cloud: mode 3
400	DO i_lev=24,24
401	R_MEDIAN(:,i_lev,3)=3.65e-6
402	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
403	STDDEV(:,i_lev,3)=1.28
404	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
405	K_MASS(:,i_lev,3)=0.966
406	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
407	END DO
408
409	! Middle Cloud: mode 1
410	DO i_lev=25,28
411	R_MEDIAN(:,i_lev,1)=0.15e-6
412	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
413	STDDEV(:,i_lev,1)=1.9
414	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
415	K_MASS(:,i_lev,1)=0.0084
416	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
417	END DO
418
419	! Middle Cloud: mode 2 prime
420	DO i_lev=25,28
421	R_MEDIAN(:,i_lev,2)=1.4e-6
422	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
423	STDDEV(:,i_lev,2)=1.23
424	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
425	K_MASS(:,i_lev,2)=0.21
426	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
427	END DO
428
429	! Middle Cloud: mode 3
430	DO i_lev=25,28
431	R_MEDIAN(:,i_lev,3)=3.65e-6
432	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
433	STDDEV(:,i_lev,3)=1.28
434	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
435	K_MASS(:,i_lev,3)=0.7816
436	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
437	END DO
438
439	! option: upper haze remplacee par extension upper cloud
440	! => 35 remplace par cloudmax et upper haze commentee
441	! ===============================================
442
443	! Upper Cloud: mode 1
444	DO i_lev=29,35 !cloudmax
445	R_MEDIAN(:,i_lev,1)=0.2e-6
446	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
447	STDDEV(:,i_lev,1)=2.16
448	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
449	K_MASS(:,i_lev,1)=0.72
450	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
451	END DO
452
453	! Upper Cloud: mode 2
454	DO i_lev=29,35 !cloudmax
455	R_MEDIAN(:,i_lev,2)=1.0e-6
456	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
457	STDDEV(:,i_lev,2)=1.29
458	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
459	K_MASS(:,i_lev,2)=0.28
460	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
461	END DO
462
463	! Upper Cloud: mode 3
464	DO i_lev=29,35 !cloudmax
465	R_MEDIAN(:,i_lev,3)=3.65e-6
466	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
467	STDDEV(:,i_lev,3)=1.28
468	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
469	K_MASS(:,i_lev,3)=0.0
470	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
471	END DO
472
473	! Upper Haze: mode 1
474	DO i_lev=36, cloudmax
475	R_MEDIAN(:,i_lev,1)=0.2e-6
476	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
477	STDDEV(:,i_lev,1)=2.16
478	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
479	K_MASS(:,i_lev,1)=1.0
480	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
481	END DO
482
483	! Upper Haze: mode 2
484	DO i_lev=36, cloudmax
485	R_MEDIAN(:,i_lev,2)=1.e-6
486	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
487	STDDEV(:,i_lev,2)=1.29
488	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
489	K_MASS(:,i_lev,2)=0.0
490	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
491	END DO
492
493	! Upper Haze: mode 3
494	DO i_lev=36, cloudmax
495	R_MEDIAN(:,i_lev,3)=3.65e-6
496	PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
497	STDDEV(:,i_lev,3)=2.16
498	PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
499	K_MASS(:,i_lev,3)=0.0
500	PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
501	END DO
502
503	!=============================================================
504
505	! ===============================================================
506	! Initialisation TRIMODALE "Knollenberg" sans Mode3, Mode2 etendu
507	! ===============================================================
508
509	! Lower Haze: mode 1
510	! DO i_lev=cloudmin,22
511	! R_MEDIAN(:,i_lev,1)=0.1e-6
512	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
513	! STDDEV(:,i_lev,1)=1.57
514	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
515	! K_MASS(:,i_lev,1)=1.0
516	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
517	! END DO
518
519	! Lower Haze: mode 2
520	! DO i_lev=cloudmin,22
521	! R_MEDIAN(:,i_lev,2)=1.4e-6
522	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
523	! STDDEV(:,i_lev,2)=1.23
524	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
525	! K_MASS(:,i_lev,2)=0.0
526	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
527	! END DO
528
529	! Lower Haze: mode 3
530	! DO i_lev=cloudmin,22
531	! R_MEDIAN(:,i_lev,3)=3.65e-6
532	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
533	! STDDEV(:,i_lev,3)=1.28
534	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
535	! K_MASS(:,i_lev,3)=0.0
536	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
537	! END DO
538
539	! Pre Cloud: mode 1
540	! DO i_lev=23,23
541	! R_MEDIAN(:,i_lev,1)=0.15e-6
542	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
543	! STDDEV(:,i_lev,1)=1.8
544	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
545	! K_MASS(:,i_lev,1)=0.04
546	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
547	! END DO
548
549	! Pre Cloud: mode 2
550	! DO i_lev=23,23
551	! R_MEDIAN(:,i_lev,2)=1.0e-6
552	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
553	! STDDEV(:,i_lev,2)=1.29
554	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
555	! K_MASS(:,i_lev,2)=0.96
556	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
557	! END DO
558
559	! Pre Cloud: mode 3
560	! DO i_lev=23,23
561	! R_MEDIAN(:,i_lev,3)=3.65e-6
562	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
563	! STDDEV(:,i_lev,3)=1.28
564	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
565	! K_MASS(:,i_lev,3)=0.0
566	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
567	! END DO
568
569	! Lower Cloud: mode 1
570	! DO i_lev=24,24
571	! R_MEDIAN(:,i_lev,1)=0.2e-6
572	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
573	! STDDEV(:,i_lev,1)=1.8
574	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
575	! K_MASS(:,i_lev,1)=0.014
576	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
577	! END DO
578
579	! Lower Cloud: mode 2
580	! DO i_lev=24,24
581	! R_MEDIAN(:,i_lev,2)=1.0e-6
582	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
583	! STDDEV(:,i_lev,2)=1.6
584	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
585	! K_MASS(:,i_lev,2)=0.986
586	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
587	! END DO
588
589	! Lower Cloud: mode 3
590	! DO i_lev=24,24
591	! R_MEDIAN(:,i_lev,3)=3.65e-6
592	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
593	! STDDEV(:,i_lev,3)=1.28
594	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
595	! K_MASS(:,i_lev,3)=0.
596	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
597	! END DO
598
599	! Middle Cloud: mode 1
600	! DO i_lev=25,28
601	! R_MEDIAN(:,i_lev,1)=0.15e-6
602	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
603	! STDDEV(:,i_lev,1)=1.9
604	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
605	! K_MASS(:,i_lev,1)=0.0084
606	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
607	! END DO
608
609	! Middle Cloud: mode 2 prime
610	! DO i_lev=25,28
611	! R_MEDIAN(:,i_lev,2)=1.4e-6
612	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
613	! STDDEV(:,i_lev,2)=1.6
614	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
615	! K_MASS(:,i_lev,2)=0.9916
616	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
617	! END DO
618
619	! Middle Cloud: mode 3
620	! DO i_lev=25,28
621	! R_MEDIAN(:,i_lev,3)=3.65e-6
622	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
623	! STDDEV(:,i_lev,3)=1.28
624	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
625	! K_MASS(:,i_lev,3)=0.0
626	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
627	! END DO
628
629
630	! Upper Cloud: mode 1
631	! DO i_lev=29,35
632	! R_MEDIAN(:,i_lev,1)=0.2e-6
633	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
634	! STDDEV(:,i_lev,1)=2.16
635	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
636	! K_MASS(:,i_lev,1)=0.72
637	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
638	! END DO
639
640	! Upper Cloud: mode 2
641	! DO i_lev=29,35
642	! R_MEDIAN(:,i_lev,2)=1.0e-6
643	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
644	! STDDEV(:,i_lev,2)=1.29
645	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
646	! K_MASS(:,i_lev,2)=0.28
647	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
648	! END DO
649
650	! Upper Cloud: mode 3
651	! DO i_lev=29,35
652	! R_MEDIAN(:,i_lev,3)=3.65e-6
653	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
654	! STDDEV(:,i_lev,3)=1.28
655	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
656	! K_MASS(:,i_lev,3)=0.0
657	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
658	! END DO
659
660	! Upper Haze: mode 1
661	! DO i_lev=36, cloudmax
662	! R_MEDIAN(:,i_lev,1)=0.2e-6
663	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
664	! STDDEV(:,i_lev,1)=2.16
665	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
666	! K_MASS(:,i_lev,1)=1.0
667	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
668	! END DO
669
670	! Upper Haze: mode 2
671	! DO i_lev=36, cloudmax
672	! R_MEDIAN(:,i_lev,2)=1.e-6
673	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
674	! STDDEV(:,i_lev,2)=1.29
675	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
676	! K_MASS(:,i_lev,2)=0.0
677	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
678	! END DO
679
680	! Upper Haze: mode 3
681	! DO i_lev=36, cloudmax
682	! R_MEDIAN(:,i_lev,3)=3.65e-6
683	! PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
684	! STDDEV(:,i_lev,3)=2.16
685	! PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
686	! K_MASS(:,i_lev,3)=0.0
687	! PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
688	! END DO
689
690	!=============================================================
691	PRINT*,'==============================='
692	PRINT*,'FIN Initialisation cloud layer'
693	PRINT*,'==============================='
694
695	END SUBROUTINE cloud_ini
696
697	! ===========================================================
698
699	SUBROUTINE chemparam_ini
700	USE infotrac_phy, ONLY: nqtot, tname
701	IMPLICIT NONE
702	INTEGER :: i
703
704	ALLOCATE(M_tr(nqtot)) ! Molecular Mass of Chemistry Tracers
705	ALLOCATE(type_tr(nqtot)) ! Type of chemistry tracers 1: neutral, 2: ion, 3: liquide
706
707	! Initialisation index chemistry tracers
708
709	! Neutral gas Tracer
710	i_co2 = 0
711	i_co = 0
712	i_h2 = 0
713	i_h2o = 0
714	i_o1d = 0
715	i_o = 0
716	i_o2 = 0
717	i_o2dg = 0
718	i_o3 = 0
719	i_h = 0
720	i_oh = 0
721	i_ho2 = 0
722	i_h2o2 = 0
723	i_cl = 0
724	i_clo = 0
725	i_cl2 = 0
726	i_hcl = 0
727	i_hocl = 0
728	i_clco = 0
729	i_clco3 = 0
730	i_cocl2 = 0
731	i_s = 0
732	i_so = 0
733	i_so2 = 0
734	i_so3 = 0
735	i_s2o2 = 0
736	i_ocs = 0
737	i_hso3 = 0
738	i_h2so4 = 0
739	i_s2 = 0
740	i_clso2 = 0
741	i_oscl = 0
742	i_n2 = 0
743	i_he = 0
744	i_n2d = 0
745	i_n = 0
746	i_no = 0
747	i_no2 = 0
748	! ION TRACERS FOR IONCHEM = TRUE
749	i_co2plus = 0
750	i_coplus = 0
751	i_oplus = 0
752	i_o2plus = 0
753	i_n2plus = 0
754	i_hplus = 0
755	i_h2oplus = 0
756	i_nplus = 0
757	i_ohplus = 0
758	i_cplus = 0
759	i_noplus = 0
760	i_h3oplus = 0
761	i_hcoplus = 0
762	i_hco2plus = 0
763	i_elec = 0
764	! MICROPHYSICAL TRACERS FOR CL_SCHEME=1
765	i_h2oliq = 0
766	i_h2so4liq = 0
767
768	DO i=1, nqtot
769
770	PRINT*,'i',i
771	PRINT*,'tname(i)',tname(i)
772
773	SELECT CASE(tname(i))
774	! NEUTRAL TRACERS
775	CASE('co2')
776	i_co2=i
777	PRINT*,'co2',i_co2
778	M_tr(i_co2) = 44.0095
779	type_tr(i_co2) = 1
780	CASE('co')
781	i_co=i
782	PRINT*,'co',i_co
783	M_tr(i_co)=28.0101
784	type_tr(i_co) = 1
785	CASE('h2')
786	i_h2=i
787	PRINT*,'h2',i_h2
788	M_tr(i_h2)= 2.01588
789	type_tr(i_h2) = 1
790	CASE('h2o')
791	i_h2o=i
792	PRINT*,'h2o',i_h2o
793	M_tr(i_h2o)=18.0153
794	type_tr(i_h2o) = 1
795	CASE('o1d')
796	i_o1d=i
797	PRINT*,'o1d',i_o1d
798	M_tr(i_o1d)=15.994
799	type_tr(i_o1d) = 1
800	CASE('o')
801	i_o=i
802	PRINT*,'o',i_o
803	M_tr(i_o)=15.994
804	type_tr(i_o) = 1
805	CASE('o2')
806	i_o2=i
807	PRINT*,'o2',i_o2
808	M_tr(i_o2)=31.9988
809	type_tr(i_o2) = 1
810	CASE('o2dg')
811	i_o2dg=i
812	PRINT*,'o2dg',i_o2dg
813	M_tr(i_o2dg)=31.9988
814	type_tr(i_o2dg) = 1
815	CASE('o3')
816	i_o3=i
817	PRINT*,'o3',i_o3
818	M_tr(i_o3)= 47.9982
819	type_tr(i_o3) = 1
820	CASE('h')
821	i_h=i
822	PRINT*,'h',i_h
823	M_tr(i_h)= 1.00794
824	type_tr(i_h) = 1
825	CASE('oh')
826	i_oh=i
827	PRINT*,'oh',i_oh
828	M_tr(i_oh)=17.0073
829	type_tr(i_oh) = 1
830	CASE('ho2')
831	i_ho2=i
832	PRINT*,'ho2',i_ho2
833	M_tr(i_ho2)=33.0067
834	type_tr(i_ho2) = 1
835	CASE('h2o2')
836	i_h2o2=i
837	PRINT*,'h2o2',i_h2o2
838	M_tr(i_h2o2)=34.0147
839	type_tr(i_h2o2) = 1
840	CASE('cl')
841	i_cl=i
842	PRINT*,'cl',i_cl
843	M_tr(i_cl)=35.453
844	type_tr(i_cl) = 1
845	CASE('clo')
846	i_clo=i
847	PRINT*,'clo',i_clo
848	M_tr(i_clo)=51.452
849	type_tr(i_clo) = 1
850	CASE('cl2')
851	i_cl2=i
852	PRINT*,'cl2',i_cl2
853	M_tr(i_cl2)=70.906
854	type_tr(i_cl2) = 1
855	CASE('hcl')
856	i_hcl=i
857	PRINT*,'hcl',i_hcl
858	M_tr(i_hcl)=36.461
859	type_tr(i_hcl) = 1
860	CASE('hocl')
861	i_hocl=i
862	PRINT*,'hocl',i_hocl
863	M_tr(i_hocl)=52.46
864	type_tr(i_hocl) = 1
865	CASE('clco')
866	i_clco=i
867	PRINT*,'clco',i_clco
868	M_tr(i_clco)=63.463
869	type_tr(i_clco) = 1
870	CASE('clco3')
871	i_clco3=i
872	PRINT*,'clco3',i_clco3
873	M_tr(i_clco3)=95.462
874	type_tr(i_clco3) = 1
875	CASE('cocl2')
876	i_cocl2=i
877	PRINT*,'cocl2',i_cocl2
878	M_tr(i_cocl2)=98.916
879	type_tr(i_cocl2) = 1
880	CASE('s')
881	i_s=i
882	PRINT*,'s',i_s
883	M_tr(i_s)=32.065
884	type_tr(i_s) = 1
885	CASE('so')
886	i_so=i
887	PRINT*,'so',i_so
888	M_tr(i_so)=48.0644
889	type_tr(i_so) = 1
890	CASE('so2')
891	i_so2=i
892	PRINT*,'so2',i_so2
893	M_tr(i_so2)=64.064
894	type_tr(i_so2) = 1
895	CASE('so3')
896	i_so3=i
897	PRINT*,'so3',i_so3
898	M_tr(i_so3)=80.063
899	type_tr(i_so3) = 1
900	CASE('s2o2')
901	i_s2o2=i
902	PRINT*,'s2o2',i_s2o2
903	M_tr(i_s2o2)= 96.1288
904	type_tr(i_s2o2) = 1
905	CASE('ocs')
906	i_ocs=i
907	PRINT*,'ocs',i_ocs
908	M_tr(i_ocs)=60.0751
909	type_tr(i_ocs) = 1
910	CASE('hso3')
911	i_hso3=i
912	PRINT*,'hso3',i_hso3
913	M_tr(i_hso3)=81.071
914	type_tr(i_hso3) = 1
915	CASE('h2so4')
916	i_h2so4=i
917	PRINT*,'h2so4',i_h2so4
918	M_tr(i_h2so4)=98.078
919	type_tr(i_h2so4) = 1
920	CASE('s2')
921	i_s2=i
922	PRINT*,'s2',i_s2
923	M_tr(i_s2)=64.13
924	type_tr(i_s2) = 1
925	CASE('clso2')
926	i_clso2=i
927	PRINT*,'clso2',i_clso2
928	M_tr(i_clso2)=99.517
929	type_tr(i_clso2) = 1
930	CASE('oscl')
931	i_oscl=i
932	PRINT*,'oscl',i_oscl
933	M_tr(i_oscl)=83.517
934	type_tr(i_oscl) = 1
935	CASE('n2')
936	i_n2=i
937	PRINT*,'n2',i_n2
938	M_tr(i_n2)=28.013
939	type_tr(i_n2) = 1
940	CASE('he')
941	i_he=i
942	PRINT*,'he',i_he
943	M_tr(i_he)=4.0026
944	type_tr(i_he) = 1
945	CASE('n2d')
946	i_n2d=i
947	PRINT*,'n2d',i_n2d
948	M_tr(i_n2d)=14.0067
949	type_tr(i_n2d) = 1
950	CASE('n')
951	i_n=i
952	PRINT*,'n',i_n
953	M_tr(i_n)=14.0067
954	type_tr(i_n) = 1
955	CASE('no')
956	i_no=i
957	PRINT*,'no',i_no
958	M_tr(i_no)=30.0061
959	type_tr(i_no) = 1
960	CASE('no2')
961	i_no2=i
962	PRINT*,'no2',i_no2
963	M_tr(i_no2)=46.0055
964	type_tr(i_no2) = 1
965	! ION TRACERS FOR OK_IONCHEM = TRUE
966	CASE('co2plus')
967	i_co2plus=i
968	PRINT*,'co2plus',i_co2plus
969	M_tr(i_co2plus)=44.0095
970	type_tr(i_co2plus) = 2
971	CASE('coplus')
972	i_coplus=i
973	PRINT*,'coplus',i_coplus
974	M_tr(i_coplus)=28.0101
975	type_tr(i_coplus) = 2
976	CASE('oplus')
977	i_oplus=i
978	PRINT*,'oplus',i_oplus
979	M_tr(i_oplus)=15.994
980	type_tr(i_oplus) = 2
981	CASE('o2plus')
982	i_o2plus=i
983	PRINT*,'o2plus',i_o2plus
984	M_tr(i_o2plus)=31.9988
985	type_tr(i_o2plus) = 2
986	CASE('n2plus')
987	i_n2plus=i
988	PRINT*,'n2plus',i_n2plus
989	M_tr(i_n2plus)=28.013
990	type_tr(i_n2plus) = 2
991	CASE('hplus')
992	i_hplus=i
993	PRINT*,'hplus',i_hplus
994	M_tr(i_hplus)=1.00794
995	type_tr(i_hplus) = 2
996	CASE('h2oplus')
997	i_h2oplus=i
998	PRINT*,'h2oplus',i_h2oplus
999	M_tr(i_h2oplus)=18.0153
1000	type_tr(i_h2oplus) = 2
1001	CASE('nplus')
1002	i_nplus=i
1003	PRINT*,'nplus',i_nplus
1004	M_tr(i_nplus)=14.0067
1005	type_tr(i_nplus) = 2
1006	CASE('ohplus')
1007	i_ohplus=i
1008	PRINT*,'ohplus',i_ohplus
1009	M_tr(i_ohplus)=17.0073
1010	type_tr(i_ohplus) = 2
1011	CASE('cplus')
1012	i_cplus=i
1013	PRINT*,'cplus',i_cplus
1014	M_tr(i_cplus)=12.011
1015	type_tr(i_cplus) = 2
1016	CASE('noplus')
1017	i_noplus=i
1018	PRINT*,'noplus',i_noplus
1019	M_tr(i_noplus)=30.0061
1020	type_tr(i_noplus) = 2
1021	CASE('h3oplus')
1022	i_h3oplus=i
1023	PRINT*,'h3oplus',i_h3oplus
1024	M_tr(i_h3oplus)=19.0232
1025	type_tr(i_h3oplus) = 2
1026	CASE('hcoplus')
1027	i_hcoplus=i
1028	PRINT*,'hcoplus',i_hcoplus
1029	M_tr(i_hcoplus)=29.0180
1030	type_tr(i_hcoplus) = 2
1031	CASE('hco2plus')
1032	i_hco2plus=i
1033	PRINT*,'hco2plus',i_hco2plus
1034	M_tr(i_hco2plus)=45.
1035	type_tr(i_hco2plus) = 2
1036	CASE('elec')
1037	i_elec=i
1038	PRINT*,'elec',i_elec
1039	M_tr(i_elec)=1./1822.89
1040	type_tr(i_elec) = 2
1041	! MICROPHYSICAL TRACERS FOR CL_SCHEME=1
1042	CASE('h2oliq')
1043	i_h2oliq=i
1044	PRINT*,'h2oliq',i_h2oliq
1045	M_tr(i_h2oliq)=18.0153
1046	type_tr(i_h2oliq) = 3
1047	CASE('h2so4liq')
1048	i_h2so4liq=i
1049	PRINT*,'h2so4liq',i_h2so4liq
1050	M_tr(i_h2so4liq)=98.078
1051	type_tr(i_h2so4liq) = 3
1052	! MICROPHYSICAL TRACERS FOR CL_SCHEME=2
1053	CASE('M0_aer')
1054	i_m0_aer=i
1055	type_tr(i_m0_aer) = 10
1056	PRINT*,'M0_aer',i_m0_aer
1057	CASE('M3_aer')
1058	i_m3_aer=i
1059	type_tr(i_m3_aer) = 10
1060	PRINT*,'M3_aer',i_m3_aer
1061	CASE('M0_m1drop')
1062	i_m0_mode1drop=i
1063	type_tr(i_m0_mode1drop) = 10
1064	PRINT*,'M0_m1drop',i_m0_mode1drop
1065	CASE('M0_m1ccn')
1066	i_m0_mode1ccn=i
1067	type_tr(i_m0_mode1ccn) = 10
1068	PRINT*,'M0_m1ccn',i_m0_mode1ccn
1069	CASE('M3_m1sa')
1070	i_m3_mode1sa=i
1071	type_tr(i_m3_mode1sa) = 10
1072	PRINT*,'M3_m1sa',i_m3_mode1sa
1073	CASE('M3_m1w')
1074	i_m3_mode1w=i
1075	type_tr(i_m3_mode1w) = 10
1076	PRINT*,'M3_m1w',i_m3_mode1w
1077	CASE('M3_m1ccn')
1078	i_m3_mode1ccn=i
1079	type_tr(i_m3_mode1ccn) = 10
1080	PRINT*,'M3_m1ccn',i_m3_mode1ccn
1081	CASE('M0_m2drop')
1082	i_m0_mode2drop=i
1083	type_tr(i_m0_mode2drop) = 10
1084	PRINT*,'M0_m2drop',i_m0_mode2drop
1085	CASE('M0_m2ccn')
1086	i_m0_mode2ccn=i
1087	type_tr(i_m0_mode2ccn) = 10
1088	PRINT*,'M0_m2ccn',i_m0_mode2ccn
1089	CASE('M3_m2sa')
1090	i_m3_mode2sa=i
1091	type_tr(i_m3_mode2sa) = 10
1092	PRINT*,'M3_m2sa',i_m3_mode2sa
1093	CASE('M3_m2w')
1094	i_m3_mode2w=i
1095	type_tr(i_m3_mode2w) = 10
1096	PRINT*,'M3_m2w',i_m3_mode2w
1097	CASE('M3_m2ccn')
1098	i_m3_mode2ccn=i
1099	type_tr(i_m3_mode2ccn) = 10
1100	PRINT*,'M3_m2ccn',i_m3_mode2ccn
1101	END SELECT
1102
1103	! PRINT*,'M_tr(i)',M_tr(i)
1104	END DO
1105
1106	END SUBROUTINE chemparam_ini
1107
1108	! ===========================================================
1109
1110	SUBROUTINE vapors4muphy_ini(nlon,nlev,trac)
1111	USE infotrac_phy, ONLY: nqtot, tname
1112	IMPLICIT NONE
1113
1114	integer :: nlon, nlev
1115	real :: trac(nlon,nlev,nqtot) ! traceur ( en vmr)
1116
1117	! integer :: i
1118	! real :: trac1d(nlev,2) ! traceur lu ( en vmr)
1119
1120	! lecture d'un fichier texte contenant les profils de trac1d(:1) = H2O et trac1d(:,2) = H2SO4
1121	! DO i=1,nlon
1122	! trac(i,:,i_h2o) = trac1d(:,1)
1123	! trac(i,:,i_h2so4) = trac1d(:,2)
1124	! ENDDO
1125
1126	! intitialisation profils altitude H2O et H2SO4
1127	! profil H2O initial vap+liq == que vap
1128	trac(:,1:24,i_h2o) = 30.E-6 !
1129	trac(:,25:50,i_h2o) = 1.E-6 !
1130
1131	trac(:,:,i_h2so4) = 3.E-9 ! Limite sup Sandor 2012
1132	trac(:,23:50,i_h2so4) = 2.E-6 ! Profil H2SO4 initial => vap+liq
1133
1134	END SUBROUTINE vapors4muphy_ini
1135
1136	END MODULE chemparam_mod
1137

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

Download in other formats:

Original Format