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

Rev	Line
[1305]	1	MODULE chemparam_mod
	2
[1661]	3	!MODULE qui definit les indices des traceurs et leurs masses molaires.
[1305]	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, &
[2836]	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
[1305]	21
[1661]	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
[2836]	32	REAL, DIMENSION(:), SAVE, ALLOCATABLE :: M_tr ! Molecular Mass of tracers
	33	REAL, DIMENSION(:), SAVE, ALLOCATABLE :: type_tr ! type of tracer
[1305]	34
[2795]	35	REAL, DIMENSION(:,:),SAVE, ALLOCATABLE :: no_emission
	36	REAL, DIMENSION(:,:),SAVE, ALLOCATABLE :: o2_emission
[1305]	37	!----------------------------------------------------------------------------
[1661]	38	! DEF FOR CL_SCHEME = 1 (AURELIEN)
	39
[1442]	40	! number of clouds mode modelized
	41	INTEGER, PARAMETER :: nbr_mode = 3
[1305]	42	INTEGER :: i_cloud
	43	INTEGER, SAVE :: cloudmax
	44	INTEGER, SAVE :: cloudmin
[1442]	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
[1305]	52	REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: WH2SO4
	53	REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: rho_droplet
	54	!----------------------------------------------------------------------------
[1661]	55	! DEF FOR CL_SCHEME = 2 (FULL MICROPHYS)
[1305]	56
[1661]	57	!----------------------------------------------------------------------------
	58
[1305]	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
[1442]	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))
[1305]	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
[1442]	82	PRINT*,'nbr_mode',nbr_mode
	83
	84	NBRTOT(:,:,:) = 0.0E+0
	85	WH2SO4(:,:) = 0.0E+0
	86	rho_droplet(:,:) = 0.0E+0
	87
[1305]	88	!=============================================================
	89	! Initialisation cloud layer 1
	90	!=============================================================
	91	! cloudmin et cloudmax niveaux du GCM
[1442]	92	cloudmin= 18
[1305]	93	cloudmax= 50
	94
	95	! radius R_MEDIAN en m (donc *e-6 pour microns)
	96
[1691]	97	R_MEDIAN(:,:,:)=0.0E+0 ! Geometric Average Radius
[1442]	98	STDDEV(:,:,:)=0.0E+0 ! Geometric Std Deviation
	99	K_MASS(:,:,:)=0.0E+0 ! Coeff multimodal
[1305]	100
	101	! ===============================================
	102	! Knollenberg & Hunten, 1980 and James et al 1997
	103	! ===============================================
	104
[1442]	105	! ===============================================
	106	! Initialisation UNIMODALE
	107	! ===============================================
	108
[1305]	109	! Lower Haze: mode 1
[1442]	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
[1305]	120	! DO i_lev=21,23
[1442]	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)
[1305]	127	! END DO
	128
[1442]	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
[1305]	213	! Lower Cloud: mode 3
[1442]	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
[1305]	232
	233	! Middle Cloud: mode 2 prime
	234	! DO i_lev=24,28
[1442]	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)
[1305]	241	! END DO
	242
[1442]	243	! Middle Cloud: mode 3
[1305]	244	! DO i_lev=24,28
[1442]	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)
[1305]	251	! END DO
	252
[1442]	253
	254	! Upper Cloud: mode 1
[1305]	255	! DO i_lev=29,35
[1442]	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)
[1305]	262	! END DO
[1442]	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
[1305]	273
	274	! Upper Cloud: mode 3
	275	! DO i_lev=29,35
[1442]	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)
[1305]	282	! END DO
[1442]	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	! ===============================================
[2048]	316	! Initialisation TRIMODALE Knollenberg
[1442]	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
[2048]	439	! option: upper haze remplacee par extension upper cloud
	440	! => 35 remplace par cloudmax et upper haze commentee
	441	! ===============================================
[1442]	442
	443	! Upper Cloud: mode 1
[2048]	444	DO i_lev=29,35 !cloudmax
[1442]	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
[1305]	453	! Upper Cloud: mode 2
[2048]	454	DO i_lev=29,35 !cloudmax
[1442]	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)
[1305]	461	END DO
[1442]	462
	463	! Upper Cloud: mode 3
[2048]	464	DO i_lev=29,35 !cloudmax
[1442]	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
[1305]	472
	473	! Upper Haze: mode 1
	474	DO i_lev=36, cloudmax
[1442]	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)
[1305]	481	END DO
[1442]	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
[1305]	503	!=============================================================
	504
[1442]	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
[2048]	689
[1442]	690	!=============================================================
[1305]	691	PRINT*,'==============================='
	692	PRINT*,'FIN Initialisation cloud layer'
	693	PRINT*,'==============================='
	694
[1442]	695	END SUBROUTINE cloud_ini
[1305]	696
[1661]	697	! ===========================================================
	698
[1305]	699	SUBROUTINE chemparam_ini
[1621]	700	USE infotrac_phy, ONLY: nqtot, tname
[1305]	701	IMPLICIT NONE
	702	INTEGER :: i
	703
[2836]	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
[1305]	768	DO i=1, nqtot
	769
	770	PRINT*,'i',i
	771	PRINT*,'tname(i)',tname(i)
	772
	773	SELECT CASE(tname(i))
[2836]	774	! NEUTRAL TRACERS
[1305]	775	CASE('co2')
	776	i_co2=i
	777	PRINT*,'co2',i_co2
[2836]	778	M_tr(i_co2) = 44.0095
	779	type_tr(i_co2) = 1
[1305]	780	CASE('co')
	781	i_co=i
	782	PRINT*,'co',i_co
[2836]	783	M_tr(i_co)=28.0101
	784	type_tr(i_co) = 1
[1305]	785	CASE('h2')
	786	i_h2=i
	787	PRINT*,'h2',i_h2
[1468]	788	M_tr(i_h2)= 2.01588
[2836]	789	type_tr(i_h2) = 1
[1305]	790	CASE('h2o')
	791	i_h2o=i
	792	PRINT*,'h2o',i_h2o
	793	M_tr(i_h2o)=18.0153
[2836]	794	type_tr(i_h2o) = 1
[1305]	795	CASE('o1d')
	796	i_o1d=i
	797	PRINT*,'o1d',i_o1d
[2836]	798	M_tr(i_o1d)=15.994
	799	type_tr(i_o1d) = 1
[1305]	800	CASE('o')
	801	i_o=i
	802	PRINT*,'o',i_o
[2836]	803	M_tr(i_o)=15.994
	804	type_tr(i_o) = 1
[1305]	805	CASE('o2')
	806	i_o2=i
	807	PRINT*,'o2',i_o2
	808	M_tr(i_o2)=31.9988
[2836]	809	type_tr(i_o2) = 1
[1305]	810	CASE('o2dg')
	811	i_o2dg=i
	812	PRINT*,'o2dg',i_o2dg
	813	M_tr(i_o2dg)=31.9988
[2836]	814	type_tr(i_o2dg) = 1
[1305]	815	CASE('o3')
	816	i_o3=i
	817	PRINT*,'o3',i_o3
[2836]	818	M_tr(i_o3)= 47.9982
	819	type_tr(i_o3) = 1
[1305]	820	CASE('h')
	821	i_h=i
	822	PRINT*,'h',i_h
[2836]	823	M_tr(i_h)= 1.00794
	824	type_tr(i_h) = 1
[1305]	825	CASE('oh')
	826	i_oh=i
	827	PRINT*,'oh',i_oh
	828	M_tr(i_oh)=17.0073
[2836]	829	type_tr(i_oh) = 1
[1305]	830	CASE('ho2')
	831	i_ho2=i
	832	PRINT*,'ho2',i_ho2
	833	M_tr(i_ho2)=33.0067
[2836]	834	type_tr(i_ho2) = 1
[1305]	835	CASE('h2o2')
	836	i_h2o2=i
	837	PRINT*,'h2o2',i_h2o2
[2836]	838	M_tr(i_h2o2)=34.0147
	839	type_tr(i_h2o2) = 1
[1305]	840	CASE('cl')
	841	i_cl=i
	842	PRINT*,'cl',i_cl
	843	M_tr(i_cl)=35.453
[2836]	844	type_tr(i_cl) = 1
[1305]	845	CASE('clo')
	846	i_clo=i
	847	PRINT*,'clo',i_clo
	848	M_tr(i_clo)=51.452
[2836]	849	type_tr(i_clo) = 1
[1305]	850	CASE('cl2')
	851	i_cl2=i
	852	PRINT*,'cl2',i_cl2
	853	M_tr(i_cl2)=70.906
[2836]	854	type_tr(i_cl2) = 1
[1305]	855	CASE('hcl')
	856	i_hcl=i
	857	PRINT*,'hcl',i_hcl
	858	M_tr(i_hcl)=36.461
[2836]	859	type_tr(i_hcl) = 1
[1305]	860	CASE('hocl')
	861	i_hocl=i
	862	PRINT*,'hocl',i_hocl
	863	M_tr(i_hocl)=52.46
[2836]	864	type_tr(i_hocl) = 1
	865	CASE('clco')
[1305]	866	i_clco=i
	867	PRINT*,'clco',i_clco
	868	M_tr(i_clco)=63.463
[2836]	869	type_tr(i_clco) = 1
	870	CASE('clco3')
[1305]	871	i_clco3=i
	872	PRINT*,'clco3',i_clco3
	873	M_tr(i_clco3)=95.462
[2836]	874	type_tr(i_clco3) = 1
	875	CASE('cocl2')
[1305]	876	i_cocl2=i
	877	PRINT*,'cocl2',i_cocl2
	878	M_tr(i_cocl2)=98.916
[2836]	879	type_tr(i_cocl2) = 1
	880	CASE('s')
[1305]	881	i_s=i
	882	PRINT*,'s',i_s
	883	M_tr(i_s)=32.065
[2836]	884	type_tr(i_s) = 1
	885	CASE('so')
[1305]	886	i_so=i
	887	PRINT*,'so',i_so
	888	M_tr(i_so)=48.0644
[2836]	889	type_tr(i_so) = 1
	890	CASE('so2')
[1305]	891	i_so2=i
	892	PRINT*,'so2',i_so2
	893	M_tr(i_so2)=64.064
[2836]	894	type_tr(i_so2) = 1
	895	CASE('so3')
[1305]	896	i_so3=i
	897	PRINT*,'so3',i_so3
	898	M_tr(i_so3)=80.063
[2836]	899	type_tr(i_so3) = 1
	900	CASE('s2o2')
[1305]	901	i_s2o2=i
	902	PRINT*,'s2o2',i_s2o2
	903	M_tr(i_s2o2)= 96.1288
[2836]	904	type_tr(i_s2o2) = 1
	905	CASE('ocs')
[1305]	906	i_ocs=i
	907	PRINT*,'ocs',i_ocs
	908	M_tr(i_ocs)=60.0751
[2836]	909	type_tr(i_ocs) = 1
	910	CASE('hso3')
[1305]	911	i_hso3=i
	912	PRINT*,'hso3',i_hso3
	913	M_tr(i_hso3)=81.071
[2836]	914	type_tr(i_hso3) = 1
	915	CASE('h2so4')
[1305]	916	i_h2so4=i
	917	PRINT*,'h2so4',i_h2so4
	918	M_tr(i_h2so4)=98.078
[2836]	919	type_tr(i_h2so4) = 1
	920	CASE('s2')
[1305]	921	i_s2=i
	922	PRINT*,'s2',i_s2
	923	M_tr(i_s2)=64.13
[2836]	924	type_tr(i_s2) = 1
	925	CASE('clso2')
[1305]	926	i_clso2=i
	927	PRINT*,'clso2',i_clso2
	928	M_tr(i_clso2)=99.517
[2836]	929	type_tr(i_clso2) = 1
	930	CASE('oscl')
[1305]	931	i_oscl=i
	932	PRINT*,'oscl',i_oscl
	933	M_tr(i_oscl)=83.517
[2836]	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
[1661]	1041	! MICROPHYSICAL TRACERS FOR CL_SCHEME=1
[1305]	1042	CASE('h2oliq')
	1043	i_h2oliq=i
	1044	PRINT*,'h2oliq',i_h2oliq
	1045	M_tr(i_h2oliq)=18.0153
[2836]	1046	type_tr(i_h2oliq) = 3
[1305]	1047	CASE('h2so4liq')
	1048	i_h2so4liq=i
	1049	PRINT*,'h2so4liq',i_h2so4liq
	1050	M_tr(i_h2so4liq)=98.078
[2836]	1051	type_tr(i_h2so4liq) = 3
[1661]	1052	! MICROPHYSICAL TRACERS FOR CL_SCHEME=2
[2836]	1053	CASE('M0_aer')
[1661]	1054	i_m0_aer=i
[2836]	1055	type_tr(i_m0_aer) = 10
[1661]	1056	PRINT*,'M0_aer',i_m0_aer
	1057	CASE('M3_aer')
	1058	i_m3_aer=i
[2836]	1059	type_tr(i_m3_aer) = 10
[1661]	1060	PRINT*,'M3_aer',i_m3_aer
[2836]	1061	CASE('M0_m1drop')
[1661]	1062	i_m0_mode1drop=i
[2836]	1063	type_tr(i_m0_mode1drop) = 10
[1661]	1064	PRINT*,'M0_m1drop',i_m0_mode1drop
[2836]	1065	CASE('M0_m1ccn')
[1661]	1066	i_m0_mode1ccn=i
[2836]	1067	type_tr(i_m0_mode1ccn) = 10
[1661]	1068	PRINT*,'M0_m1ccn',i_m0_mode1ccn
	1069	CASE('M3_m1sa')
	1070	i_m3_mode1sa=i
[2836]	1071	type_tr(i_m3_mode1sa) = 10
[1661]	1072	PRINT*,'M3_m1sa',i_m3_mode1sa
	1073	CASE('M3_m1w')
	1074	i_m3_mode1w=i
[2836]	1075	type_tr(i_m3_mode1w) = 10
[1661]	1076	PRINT*,'M3_m1w',i_m3_mode1w
	1077	CASE('M3_m1ccn')
	1078	i_m3_mode1ccn=i
[2836]	1079	type_tr(i_m3_mode1ccn) = 10
[1661]	1080	PRINT*,'M3_m1ccn',i_m3_mode1ccn
[2836]	1081	CASE('M0_m2drop')
[1661]	1082	i_m0_mode2drop=i
[2836]	1083	type_tr(i_m0_mode2drop) = 10
[1661]	1084	PRINT*,'M0_m2drop',i_m0_mode2drop
[2836]	1085	CASE('M0_m2ccn')
[1661]	1086	i_m0_mode2ccn=i
[2836]	1087	type_tr(i_m0_mode2ccn) = 10
[1661]	1088	PRINT*,'M0_m2ccn',i_m0_mode2ccn
	1089	CASE('M3_m2sa')
	1090	i_m3_mode2sa=i
[2836]	1091	type_tr(i_m3_mode2sa) = 10
[1661]	1092	PRINT*,'M3_m2sa',i_m3_mode2sa
	1093	CASE('M3_m2w')
	1094	i_m3_mode2w=i
[2836]	1095	type_tr(i_m3_mode2w) = 10
[1661]	1096	PRINT*,'M3_m2w',i_m3_mode2w
	1097	CASE('M3_m2ccn')
	1098	i_m3_mode2ccn=i
[2836]	1099	type_tr(i_m3_mode2ccn) = 10
[1661]	1100	PRINT*,'M3_m2ccn',i_m3_mode2ccn
[1305]	1101	END SELECT
	1102
	1103	! PRINT*,'M_tr(i)',M_tr(i)
	1104	END DO
	1105
	1106	END SUBROUTINE chemparam_ini
[1661]	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
[1691]	1117	! integer :: i
	1118	! real :: trac1d(nlev,2) ! traceur lu ( en vmr)
[1661]	1119
	1120	! lecture d'un fichier texte contenant les profils de trac1d(:1) = H2O et trac1d(:,2) = H2SO4
[1691]	1121	! DO i=1,nlon
	1122	! trac(i,:,i_h2o) = trac1d(:,1)
	1123	! trac(i,:,i_h2so4) = trac1d(:,2)
	1124	! ENDDO
[1661]	1125
[1691]	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 !
[1661]	1130
[1691]	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
[1661]	1133
	1134	END SUBROUTINE vapors4muphy_ini
	1135
[1305]	1136	END MODULE chemparam_mod
	1137

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

Download in other formats:

Original Format