Context Navigation

chemparam_mod.F90 @ 1624

Last change on this file since 1624 was 1621, checked in by emillour, 9 years ago

Further work on full dynamics/physics separation.

LMDZ.COMMON:

added phy_common/vertical_layers_mod.F90 to store information on vertical grid. This is where routines in the physics should get the information.
The contents of vertical_layers_mod intialized via dynphy_lonlat/inigeomphy_mod.F90.

LMDZ.MARS:

physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
created an "ini_tracer_mod" routine in module "tracer_mod" for a cleaner initialization of the later.
removed some purely dynamics-related outputs (etot0, zoom parameters, etc.) from diagfi.nc and stats.nc outputs as these informations are not available in the physics.

LMDZ.GENERIC:

physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
added nqtot to tracer_h.F90.
removed some purely dynamics-related outputs (etot0, zoom parameters, etc.) from diagfi.nc and stats.nc outputs as these informations are not available in the physics.

LMDZ.VENUS:

physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
added infotrac_phy.F90 to store information on tracers in the physics. Initialized via iniphysiq.
added cpdet_phy_mod.F90 to store t2tpot etc. functions to be used in the physics. Initialized via iniphysiq. IMPORTANT: there are some hard-coded constants! These should match what is in cpdet_mod.F90 in the dynamics.
got rid of references to moyzon_mod module within the physics. The required variables (tmoy, plevmoy) are passed to the physics as arguments to physiq.

LMDZ.TITAN:

added infotrac_phy.F90 to store information on tracers in the physics. Initialized via iniphysiq.
added cpdet_phy_mod.F90 to store t2tpot etc. functions to be used in the physics.
Extra work required to completely decouple physics and dynamics: moyzon_mod should be cleaned up and information passed from dynamics to physics as as arguments. Likewise moyzon_ch and moyzon_mu should not be queried from logic_mod (which is in the dynamics).

Property svn:executable set to *

File size: 25.9 KB

Rev	Line
[1305]	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, &
[1442]	15	i_clso2, i_oscl, i_h2oliq, i_h2so4liq, &
	16	i_n2
[1305]	17
	18	REAL, DIMENSION(:), SAVE, ALLOCATABLE :: M_tr
	19
	20
	21	!----------------------------------------------------------------------------
[1442]	22	! number of clouds mode modelized
	23	INTEGER, PARAMETER :: nbr_mode = 3
[1305]	24	INTEGER :: i_cloud
	25	INTEGER, SAVE :: cloudmax
	26	INTEGER, SAVE :: cloudmin
[1442]	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
[1305]	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
[1442]	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))
[1305]	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
[1442]	61	PRINT*,'nbr_mode',nbr_mode
	62
	63	NBRTOT(:,:,:) = 0.0E+0
	64	WH2SO4(:,:) = 0.0E+0
	65	rho_droplet(:,:) = 0.0E+0
	66
[1305]	67	!=============================================================
	68	! Initialisation cloud layer 1
	69	!=============================================================
	70	! cloudmin et cloudmax niveaux du GCM
[1442]	71	cloudmin= 18
[1305]	72	cloudmax= 50
	73
	74	! radius R_MEDIAN en m (donc *e-6 pour microns)
	75
[1442]	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
[1305]	79
	80	! ===============================================
	81	! Knollenberg & Hunten, 1980 and James et al 1997
	82	! ===============================================
	83
[1442]	84	! ===============================================
	85	! Initialisation UNIMODALE
	86	! ===============================================
	87
[1305]	88	! Lower Haze: mode 1
[1442]	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
[1305]	99	! DO i_lev=21,23
[1442]	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)
[1305]	106	! END DO
	107
[1442]	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
[1305]	192	! Lower Cloud: mode 3
[1442]	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
[1305]	211
	212	! Middle Cloud: mode 2 prime
	213	! DO i_lev=24,28
[1442]	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)
[1305]	220	! END DO
	221
[1442]	222	! Middle Cloud: mode 3
[1305]	223	! DO i_lev=24,28
[1442]	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)
[1305]	230	! END DO
	231
[1442]	232
	233	! Upper Cloud: mode 1
[1305]	234	! DO i_lev=29,35
[1442]	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)
[1305]	241	! END DO
[1442]	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
[1305]	252
	253	! Upper Cloud: mode 3
	254	! DO i_lev=29,35
[1442]	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)
[1305]	261	! END DO
[1442]	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
[1305]	429	! Upper Cloud: mode 2
	430	DO i_lev=29,35
[1442]	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)
[1305]	437	END DO
[1442]	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
[1305]	448
	449	! Upper Haze: mode 1
	450	DO i_lev=36, cloudmax
[1442]	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)
[1305]	457	END DO
[1442]	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
[1305]	479	!=============================================================
	480
[1442]	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	!=============================================================
[1305]	666	PRINT*,'==============================='
	667	PRINT*,'FIN Initialisation cloud layer'
	668	PRINT*,'==============================='
	669
[1442]	670	END SUBROUTINE cloud_ini
[1305]	671
	672	SUBROUTINE chemparam_ini
[1621]	673	USE infotrac_phy, ONLY: nqtot, tname
[1305]	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
[1468]	696	M_tr(i_h2)= 2.01588
[1305]	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
[1442]	821	CASE('n2')
	822	i_n2=i
	823	M_tr(i_n2)=28.013
[1305]	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