Context Navigation

initracer.F @ 2732

Last change on this file since 2732 was 2659, checked in by cmathe, 3 years ago
ticket 100 resolved
File size: 35.3 KB

Rev	Line
[1224]	1	SUBROUTINE initracer(ngrid,nq,qsurf)
[38]	2
[1036]	3	use tracer_mod
[2321]	4	use comcstfi_h, only: pi
[2409]	5	use dust_param_mod, only: doubleq, submicron, dustbin
[38]	6	IMPLICIT NONE
	7	c=======================================================================
	8	c subject:
	9	c --------
	10	c Initialization related to tracer
	11	c (transported dust, water, chemical species, ice...)
	12	c
	13	c Name of the tracer
	14	c
	15	c Test of dimension :
[1036]	16	c Initialize tracer related data in tracer_mod, using tracer names provided
	17	c by the dynamics in "infotrac"
[38]	18	c
	19	c
	20	c author: F.Forget
	21	c ------
	22	c Modifs: Franck Montmessin, Sebastien Lebonnois (june 2003)
	23	c Ehouarn Millour (oct. 2008) identify tracers by their names
	24	c=======================================================================
	25
	26
[1974]	27	include "callkeys.h"
[38]	28
[1036]	29	integer,intent(in) :: ngrid ! number of atmospheric columns
	30	integer,intent(in) :: nq ! number of tracers
	31	real,intent(out) :: qsurf(ngrid,nq) ! tracer on surface (e.g. kg.m-2)
	32
[283]	33	integer iq,ig,count
[38]	34	real r0_lift , reff_lift, nueff_lift
[1974]	35	real r0_storm,reff_storm
[38]	36	c Ratio of small over large dust particles (used when both
	37	c doubleq and the submicron mode are active); In Montmessin
	38	c et al. (2002), a value of 25 has been deduced;
	39	real, parameter :: popratio = 25.
[1974]	40	character(len=30) :: txt ! to store some text
[38]	41
	42	c-----------------------------------------------------------------------
[1036]	43	c radius(nq) ! aerosol particle radius (m)
	44	c rho_q(nq) ! tracer densities (kg.m-3)
	45	c alpha_lift(nq) ! saltation vertical flux/horiz flux ratio (m-1)
	46	c alpha_devil(nq) ! lifting coeeficient by dust devil
[38]	47	c rho_dust ! Mars dust density
	48	c rho_ice ! Water ice density
	49	c doubleq ! if method with mass (iq=1) and number(iq=2) mixing ratio
	50	c varian ! Characteristic variance of log-normal distribution
	51	c-----------------------------------------------------------------------
	52
[1036]	53
[38]	54	c------------------------------------------------------------
	55	c NAME and molar mass of the tracer
	56	c------------------------------------------------------------
	57
	58	! Identify tracers by their names: (and set corresponding values of mmol)
	59	! 0. initialize tracer indexes to zero:
[1036]	60	igcm_dustbin(1:nq)=0
[1617]	61	igcm_co2_ice=0
	62	igcm_ccnco2_mass=0
	63	igcm_ccnco2_number=0
[2589]	64	igcm_ccnco2_meteor_mass=0
	65	igcm_ccnco2_meteor_number=0
[2562]	66	igcm_ccnco2_h2o_mass_ice=0
	67	igcm_ccnco2_h2o_mass_ccn=0
	68	igcm_ccnco2_h2o_number=0
[38]	69	igcm_dust_mass=0
	70	igcm_dust_number=0
[358]	71	igcm_ccn_mass=0
	72	igcm_ccn_number=0
[38]	73	igcm_dust_submicron=0
	74	igcm_h2o_vap=0
	75	igcm_h2o_ice=0
[2312]	76	igcm_hdo_vap=0
	77	igcm_hdo_ice=0
[1974]	78	igcm_stormdust_mass=0
	79	igcm_stormdust_number=0
[2199]	80	igcm_topdust_mass=0
	81	igcm_topdust_number=0
[38]	82	igcm_co2=0
	83	igcm_co=0
	84	igcm_o=0
	85	igcm_o1d=0
	86	igcm_o2=0
	87	igcm_o3=0
	88	igcm_h=0
[2461]	89	igcm_d=0
	90	igcm_hd=0
[38]	91	igcm_h2=0
[2461]	92	igcm_od=0
	93	igcm_do2=0
	94	igcm_hdo2=0
[38]	95	igcm_oh=0
	96	igcm_ho2=0
	97	igcm_h2o2=0
[459]	98	igcm_ch4=0
[38]	99	igcm_n2=0
	100	igcm_ar=0
	101	igcm_ar_n2=0
[171]	102	igcm_n=0
	103	igcm_no=0
	104	igcm_no2=0
	105	igcm_n2d=0
[1660]	106	igcm_he=0
[171]	107	igcm_co2plus=0
	108	igcm_oplus=0
	109	igcm_o2plus=0
	110	igcm_coplus=0
	111	igcm_cplus=0
	112	igcm_nplus=0
	113	igcm_noplus=0
	114	igcm_n2plus=0
	115	igcm_hplus=0
[635]	116	igcm_hco2plus=0
[2284]	117	igcm_hcoplus=0
[2302]	118	igcm_h2oplus=0
[2321]	119	igcm_h3oplus=0
	120	igcm_ohplus=0
[171]	121	igcm_elec=0
[2616]	122
[38]	123	! 1. find dust tracers
	124	count=0
	125	if (dustbin.gt.0) then
[1036]	126	do iq=1,nq
[38]	127	txt=" "
	128	write(txt,'(a4,i2.2)')'dust',count+1
	129	if (noms(iq).eq.txt) then
	130	count=count+1
	131	igcm_dustbin(count)=iq
	132	mmol(iq)=100.
	133	endif
[1036]	134	enddo !do iq=1,nq
[38]	135	endif ! of if (dustbin.gt.0)
	136	if (doubleq) then
[1036]	137	do iq=1,nq
[38]	138	if (noms(iq).eq."dust_mass") then
	139	igcm_dust_mass=iq
	140	count=count+1
	141	endif
	142	if (noms(iq).eq."dust_number") then
	143	igcm_dust_number=iq
	144	count=count+1
	145	endif
	146	enddo
	147	endif ! of if (doubleq)
[740]	148	if (microphys) then
[1036]	149	do iq=1,nq
[358]	150	if (noms(iq).eq."ccn_mass") then
	151	igcm_ccn_mass=iq
	152	count=count+1
	153	endif
	154	if (noms(iq).eq."ccn_number") then
	155	igcm_ccn_number=iq
	156	count=count+1
	157	endif
	158	enddo
[740]	159	endif ! of if (microphys)
[38]	160	if (submicron) then
[1036]	161	do iq=1,nq
[38]	162	if (noms(iq).eq."dust_submicron") then
	163	igcm_dust_submicron=iq
	164	mmol(iq)=100.
	165	count=count+1
	166	endif
	167	enddo
	168	endif ! of if (submicron)
[1974]	169	if (rdstorm) then
	170	do iq=1,nq
	171	if (noms(iq).eq."stormdust_mass") then
	172	igcm_stormdust_mass=iq
	173	count=count+1
	174	endif
	175	if (noms(iq).eq."stormdust_number") then
	176	igcm_stormdust_number=iq
	177	count=count+1
	178	endif
	179	enddo
[2199]	180	endif ! of if (rdstorm)
[2628]	181	if (topflows) then
[2199]	182	do iq=1,nq
	183	if (noms(iq).eq."topdust_mass") then
	184	igcm_topdust_mass=iq
	185	count=count+1
	186	endif
	187	if (noms(iq).eq."topdust_number") then
	188	igcm_topdust_number=iq
	189	count=count+1
	190	endif
	191	enddo
[2628]	192	endif ! of if (topflows)
[38]	193	! 2. find chemistry and water tracers
[1036]	194	do iq=1,nq
[38]	195	if (noms(iq).eq."co2") then
	196	igcm_co2=iq
	197	mmol(igcm_co2)=44.
	198	count=count+1
	199	endif
	200	if (noms(iq).eq."co") then
	201	igcm_co=iq
	202	mmol(igcm_co)=28.
	203	count=count+1
	204	endif
	205	if (noms(iq).eq."o") then
	206	igcm_o=iq
	207	mmol(igcm_o)=16.
	208	count=count+1
	209	endif
	210	if (noms(iq).eq."o1d") then
	211	igcm_o1d=iq
	212	mmol(igcm_o1d)=16.
	213	count=count+1
	214	endif
	215	if (noms(iq).eq."o2") then
	216	igcm_o2=iq
	217	mmol(igcm_o2)=32.
	218	count=count+1
	219	endif
	220	if (noms(iq).eq."o3") then
	221	igcm_o3=iq
	222	mmol(igcm_o3)=48.
	223	count=count+1
	224	endif
	225	if (noms(iq).eq."h") then
	226	igcm_h=iq
	227	mmol(igcm_h)=1.
	228	count=count+1
	229	endif
	230	if (noms(iq).eq."h2") then
	231	igcm_h2=iq
	232	mmol(igcm_h2)=2.
	233	count=count+1
	234	endif
	235	if (noms(iq).eq."oh") then
	236	igcm_oh=iq
	237	mmol(igcm_oh)=17.
	238	count=count+1
	239	endif
	240	if (noms(iq).eq."ho2") then
	241	igcm_ho2=iq
	242	mmol(igcm_ho2)=33.
	243	count=count+1
	244	endif
	245	if (noms(iq).eq."h2o2") then
	246	igcm_h2o2=iq
	247	mmol(igcm_h2o2)=34.
	248	count=count+1
	249	endif
	250	if (noms(iq).eq."n2") then
	251	igcm_n2=iq
	252	mmol(igcm_n2)=28.
	253	count=count+1
	254	endif
[324]	255	if (noms(iq).eq."ch4") then
	256	igcm_ch4=iq
	257	mmol(igcm_ch4)=16.
	258	count=count+1
	259	endif
[38]	260	if (noms(iq).eq."ar") then
	261	igcm_ar=iq
	262	mmol(igcm_ar)=40.
	263	count=count+1
	264	endif
[324]	265	if (noms(iq).eq."n") then
[171]	266	igcm_n=iq
	267	mmol(igcm_n)=14.
	268	count=count+1
	269	endif
	270	if (noms(iq).eq."no") then
	271	igcm_no=iq
	272	mmol(igcm_no)=30.
	273	count=count+1
	274	endif
	275	if (noms(iq).eq."no2") then
	276	igcm_no2=iq
	277	mmol(igcm_no2)=46.
	278	count=count+1
	279	endif
	280	if (noms(iq).eq."n2d") then
	281	igcm_n2d=iq
	282	mmol(igcm_n2d)=28.
	283	count=count+1
	284	endif
[1660]	285	if (noms(iq).eq."he") then
	286	igcm_he=iq
	287	mmol(igcm_he)=4.
	288	count=count+1
	289	endif
[171]	290	if (noms(iq).eq."co2plus") then
	291	igcm_co2plus=iq
	292	mmol(igcm_co2plus)=44.
	293	count=count+1
	294	endif
	295	if (noms(iq).eq."oplus") then
	296	igcm_oplus=iq
	297	mmol(igcm_oplus)=16.
	298	count=count+1
	299	endif
	300	if (noms(iq).eq."o2plus") then
	301	igcm_o2plus=iq
	302	mmol(igcm_o2plus)=32.
	303	count=count+1
	304	endif
	305	if (noms(iq).eq."coplus") then
	306	igcm_coplus=iq
	307	mmol(igcm_coplus)=28.
	308	count=count+1
	309	endif
	310	if (noms(iq).eq."cplus") then
	311	igcm_cplus=iq
	312	mmol(igcm_cplus)=12.
	313	count=count+1
	314	endif
	315	if (noms(iq).eq."nplus") then
	316	igcm_nplus=iq
	317	mmol(igcm_nplus)=14.
	318	count=count+1
	319	endif
	320	if (noms(iq).eq."noplus") then
	321	igcm_noplus=iq
	322	mmol(igcm_noplus)=30.
	323	count=count+1
	324	endif
[324]	325	if (noms(iq).eq."n2plus") then
[171]	326	igcm_n2plus=iq
	327	mmol(igcm_n2plus)=28.
	328	count=count+1
	329	endif
	330	if (noms(iq).eq."hplus") then
	331	igcm_hplus=iq
	332	mmol(igcm_hplus)=1.
	333	count=count+1
	334	endif
[635]	335	if (noms(iq).eq."hco2plus") then
	336	igcm_hco2plus=iq
	337	mmol(igcm_hco2plus)=45.
	338	count=count+1
	339	endif
[2284]	340	if (noms(iq).eq."hcoplus") then
	341	igcm_hcoplus=iq
	342	mmol(igcm_hcoplus)=29.
	343	count=count+1
	344	endif
[2302]	345	if (noms(iq).eq."h2oplus") then
	346	igcm_h2oplus=iq
	347	mmol(igcm_h2oplus)=18.
	348	count=count+1
	349	endif
[2321]	350	if (noms(iq).eq."h3oplus") then
	351	igcm_h3oplus=iq
	352	mmol(igcm_h3oplus)=19.
	353	count=count+1
	354	endif
	355	if (noms(iq).eq."ohplus") then
	356	igcm_ohplus=iq
	357	mmol(igcm_ohplus)=17.
	358	count=count+1
	359	endif
[171]	360	if (noms(iq).eq."elec") then
	361	igcm_elec=iq
	362	mmol(igcm_elec)=1./1822.89
	363	count=count+1
	364	endif
[38]	365	if (noms(iq).eq."h2o_vap") then
	366	igcm_h2o_vap=iq
	367	mmol(igcm_h2o_vap)=18.
	368	count=count+1
	369	endif
[2312]	370	if (noms(iq).eq."hdo_vap") then
	371	igcm_hdo_vap=iq
[2462]	372	mmol(igcm_hdo_vap)=19.
[2312]	373	count=count+1
	374	endif
[2461]	375	if (noms(iq).eq."od") then
	376	igcm_od=iq
	377	mmol(igcm_od)=18.
	378	count=count+1
	379	endif
	380	if (noms(iq).eq."d") then
	381	igcm_d=iq
	382	mmol(igcm_d)=2.
	383	count=count+1
	384	endif
	385	if (noms(iq).eq."hd") then
	386	igcm_hd=iq
	387	mmol(igcm_hd)=3.
	388	count=count+1
	389	endif
	390	if (noms(iq).eq."do2") then
	391	igcm_do2=iq
	392	mmol(igcm_do2)=34.
	393	count=count+1
	394	endif
	395	if (noms(iq).eq."hdo2") then
	396	igcm_hdo2=iq
	397	mmol(igcm_hdo2)=35.
	398	count=count+1
	399	endif
[1617]	400	if (noms(iq).eq."co2_ice") then
	401	igcm_co2_ice=iq
	402	mmol(igcm_co2_ice)=44.
	403	count=count+1
	404	endif
[38]	405	if (noms(iq).eq."h2o_ice") then
	406	igcm_h2o_ice=iq
	407	mmol(igcm_h2o_ice)=18.
	408	count=count+1
	409	endif
[2312]	410	if (noms(iq).eq."hdo_ice") then
	411	igcm_hdo_ice=iq
[2463]	412	mmol(igcm_hdo_ice)=19.
[2312]	413	count=count+1
	414	endif
[38]	415	! Other stuff: e.g. for simulations using co2 + neutral gaz
	416	if (noms(iq).eq."Ar_N2") then
	417	igcm_ar_n2=iq
	418	mmol(igcm_ar_n2)=30.
	419	count=count+1
	420	endif
[1720]	421	if (co2clouds) then
[1617]	422	if (noms(iq).eq."ccnco2_mass") then
	423	igcm_ccnco2_mass=iq
	424	count=count+1
	425	endif
	426	if (noms(iq).eq."ccnco2_number") then
	427	igcm_ccnco2_number=iq
	428	count=count+1
	429	endif
[2589]	430	if (meteo_flux) then
	431	if (noms(iq).eq."ccnco2_meteor_mass") then
	432	igcm_ccnco2_meteor_mass=iq
	433	count=count+1
	434	endif
	435	if (noms(iq).eq."ccnco2_meteor_number") then
	436	igcm_ccnco2_meteor_number=iq
	437	count=count+1
	438	endif
	439	end if
[2562]	440	if (co2useh2o) then
	441	if (noms(iq).eq."ccnco2_h2o_number") then
	442	igcm_ccnco2_h2o_number=iq
	443	count=count+1
	444	endif
	445	if (noms(iq).eq."ccnco2_h2o_mass_ice") then
	446	igcm_ccnco2_h2o_mass_ice=iq
	447	count=count+1
	448	endif
	449	if (noms(iq).eq."ccnco2_h2o_mass_ccn") then
	450	igcm_ccnco2_h2o_mass_ccn=iq
	451	count=count+1
	452	endif
	453	end if
[1617]	454	endif
	455	enddo ! of do iq=1,nq
	456
[38]	457	! check that we identified all tracers:
[1036]	458	if (count.ne.nq) then
[38]	459	write(,) "initracer: found only ",count," tracers"
[1036]	460	write(,) " expected ",nq
[38]	461	do iq=1,count
	462	write(,)' ',iq,' ',trim(noms(iq))
	463	enddo
[2302]	464	call abort_physic("initracer","tracer mismatch",1)
[38]	465	else
	466	write(,) "initracer: found all expected tracers, namely:"
[1036]	467	do iq=1,nq
[38]	468	write(,)' ',iq,' ',trim(noms(iq))
	469	enddo
	470	endif
	471
	472	! if water cycle but iceparty=.false., there will nevertheless be
	473	! water ice at the surface (iceparty is not used anymore, but this
	474	! part is still relevant, as we want to stay compatible with the
	475	! older versions).
	476	if (water.and.(igcm_h2o_ice.eq.0)) then
	477	igcm_h2o_ice=igcm_h2o_vap ! so that qsurf(i_h2o_ice) is identified
	478	! even though there is no q(i_h2o_ice)
	479	else
	480	! surface ice qsurf(i_h2o_ice) was loaded twice by phyetat0,
	481	! as qsurf(i_h2o_vap) & as qsurf(i_h2o_ice), so to be clean:
	482	if (igcm_h2o_vap.ne.0) then
[1036]	483	qsurf(1:ngrid,igcm_h2o_vap)=0
[38]	484	endif
	485	endif
	486
[2312]	487	! Additional test required for HDO
	488	! We need to compute some things for H2O before HDO
	489	if (hdo) then
	490	if (igcm_h2o_vap.gt.igcm_hdo_vap) then
[2398]	491	call abort_physic("initracer",
	492	& "Tracer H2O must be initialized before HDO",1)
[2324]	493	else if ((nqfils(igcm_h2o_ice).lt.1)
	494	& .or. (nqfils(igcm_h2o_vap).lt.1)) then
	495	write(,) "HDO must be transported as a son",
	496	& " of H2O: complete traceur.def"
[2398]	497	call abort_physic("initracer","adapt your tracer.def",1)
[2324]	498	else if ((igcm_hdo_ice.lt.nq-2)
	499	& .or. (igcm_hdo_vap.lt.nq-2)) then
	500	write(,) "The isotopes (HDO) must be placed at",
	501	& " the end of the list in traceur.def"
[2398]	502	call abort_physic("initracer","adapt your tracer.def",1)
[2312]	503	endif
	504	endif
	505
[38]	506	c------------------------------------------------------------
[1036]	507	c Initialize tracers .... (in tracer_mod)
[38]	508	c------------------------------------------------------------
[1005]	509	! start by setting everything to (default) zero
[1036]	510	rho_q(1:nq)=0 ! tracer density (kg.m-3)
	511	radius(1:nq)=0. ! tracer particle radius (m)
	512	alpha_lift(1:nq) =0. ! tracer saltation vertical flux/horiz flux ratio (m-1)
	513	alpha_devil(1:nq)=0. ! tracer lifting coefficient by dust devils
[38]	514
[1005]	515
	516	! some reference values
[38]	517	rho_dust=2500. ! Mars dust density (kg.m-3)
	518	rho_ice=920. ! Water ice density (kg.m-3)
[1685]	519	rho_ice_co2=1650.
[38]	520
	521	if (doubleq) then
	522	c "doubleq" technique
	523	c -------------------
	524	c (transport of mass and number mixing ratio)
	525	c iq=1: Q mass mixing ratio, iq=2: N number mixing ratio
	526
[2312]	527	if( (nq.lt.2).or.(water.and.(nq.lt.4))
	528	* .or.(hdo.and.(nq.lt.6) )) then
[1036]	529	write(,)'initracer: nq is too low : nq=', nq
[38]	530	write(,)'water= ',water,' doubleq= ',doubleq
	531	end if
	532
	533	nueff_lift = 0.5
	534	varian=sqrt(log(1.+nueff_lift))
	535
	536	rho_q(igcm_dust_mass)=rho_dust
	537	rho_q(igcm_dust_number)=rho_dust
	538
	539	c Intermediate calcul for computing geometric mean radius r0
	540	c as a function of mass and number mixing ratio Q and N
	541	c (r0 = (r3n_q * Q/ N)^(1/3))
	542	r3n_q = exp(-4.5varian2)(3./4.)/(pi*rho_dust)
	543
	544	c Intermediate calcul for computing effective radius reff
	545	c from geometric mean radius r0
	546	c (reff = ref_r0 * r0)
	547	ref_r0 = exp(2.5varian*2)
	548
	549	c lifted dust :
	550	c '''''''''''
	551	reff_lift = 3.0e-6 !3.e-6 !Effective radius of lifted dust (m)
	552	alpha_devil(igcm_dust_mass)=9.e-9 ! dust devil lift mass coeff
	553	c alpha_lift(igcm_dust_mass)=3.0e-15 ! Lifted mass coeff
[1455]	554
	555	!! default lifting settings
	556	!! -- GCM: alpha_lift not zero because large-scale lifting by default
	557	!! -- MESOSCALE: alpha_lift zero because no lifting at all in mesoscale by default
	558	#ifdef MESOSCALE
	559	alpha_lift(igcm_dust_mass)=0.0
	560	#else
[38]	561	alpha_lift(igcm_dust_mass)=1.e-6 !1.e-6 !Lifted mass coeff
[1974]	562	IF (dustinjection.ge.1) THEN
	563	reff_lift = 3.0e-6 ! Effective radius of lifted dust (m)
	564	alpha_lift(igcm_dust_mass)=(4/3.)reff_liftrho_dust
	565	& /2.4
	566	ENDIF
[1455]	567	#endif
[38]	568
	569	r0_lift = reff_lift/ref_r0
	570	alpha_devil(igcm_dust_number)=r3n_q*
	571	& alpha_devil(igcm_dust_mass)/r0_lift**3
	572	alpha_lift(igcm_dust_number)=r3n_q*
	573	& alpha_lift(igcm_dust_mass)/r0_lift**3
	574
	575	radius(igcm_dust_mass) = reff_lift
	576	radius(igcm_dust_number) = reff_lift
	577
	578	write(,) "initracer: doubleq_param reff_lift:", reff_lift
	579	write(,) "initracer: doubleq_param nueff_lift:", nueff_lift
	580	write(,) "initracer: doubleq_param alpha_lift:",
	581	& alpha_lift(igcm_dust_mass)
[1974]	582	!c ----------------------------------------------------------------------
	583	!c rocket dust storm scheme
	584	!c lifting tracer stormdust using same distribution than
	585	!c normal dust
	586	if (rdstorm) then
	587	reff_storm=3.e-6 ! reff_lift !3.e-6
	588	r0_storm=reff_storm/ref_r0
	589	rho_q(igcm_stormdust_mass)=rho_dust
	590	rho_q(igcm_stormdust_number)=rho_dust
	591
	592	alpha_devil(igcm_stormdust_mass)=9.e-9 ! dust devil lift mass coeff
	593	alpha_lift(igcm_stormdust_mass)=4./3./2.4reff_stormrho_dust
	594
	595	write(,) 'alpha_lift(rds):',alpha_lift(igcm_stormdust_mass)
	596
	597	alpha_devil(igcm_stormdust_number)=r3n_q*
	598	& alpha_devil(igcm_stormdust_mass)/r0_storm**3
	599	alpha_lift(igcm_stormdust_number)=r3n_q*
	600	& alpha_lift(igcm_stormdust_mass)/r0_storm**3
	601
	602	radius(igcm_stormdust_mass) = reff_storm
	603	radius(igcm_stormdust_number) = reff_storm
	604	end if !(rdstorm)
	605	!c ----------------------------------------------------------------------
[2628]	606	!c mountain top dust flows scheme
[2199]	607	!c you need a radius value for topdust to active its sedimentation
	608	!c we take the same value as for the normal dust
[2628]	609	if (topflows) then
[2199]	610	rho_q(igcm_topdust_mass)=rho_dust
	611	rho_q(igcm_topdust_number)=rho_dust
	612	radius(igcm_topdust_mass) = 3.e-6
	613	radius(igcm_topdust_number) = 3.e-6
[2628]	614	end if !(topflows)
[2199]	615	!c ----------------------------------------------------------------------
[1974]	616
[38]	617	else
	618
[648]	619	! initialize varian, which may be used (e.g. by surfacearea)
	620	! even with conrath dust
	621	nueff_lift = 0.5
	622	varian=sqrt(log(1.+nueff_lift))
	623
[38]	624	if (dustbin.gt.1) then
	625	print*,'initracer: STOP!',
	626	$ ' properties of dust need to be set in initracer !!!'
[2302]	627	call abort_physic("initracer","dustbin properties issue",1)
[38]	628
	629	else if (dustbin.eq.1) then
	630
	631	c This will be used for 1 dust particle size:
	632	c ------------------------------------------
	633	radius(igcm_dustbin(1))=3.e-6
	634	alpha_lift(igcm_dustbin(1))=0.0e-6
	635	alpha_devil(igcm_dustbin(1))=7.65e-9
	636	rho_q(igcm_dustbin(1))=rho_dust
	637
	638	endif
	639	end if ! (doubleq)
	640
[358]	641
	642	c Scavenging of dust particles by H2O clouds:
	643	c ------------------------------------------
	644	c Initialize the two tracers used for the CCNs
	645	if (water.AND.doubleq.AND.scavenging) then
	646	radius(igcm_ccn_mass) = radius(igcm_dust_mass)
	647	alpha_lift(igcm_ccn_mass) = 1e-30
	648	alpha_devil(igcm_ccn_mass) = 1e-30
	649	rho_q(igcm_ccn_mass) = rho_dust
	650
	651	radius(igcm_ccn_number) = radius(igcm_ccn_mass)
	652	alpha_lift(igcm_ccn_number) = alpha_lift(igcm_ccn_mass)
	653	alpha_devil(igcm_ccn_number) = alpha_devil(igcm_ccn_mass)
	654	rho_q(igcm_ccn_number) = rho_q(igcm_ccn_mass)
	655	endif ! of if (water.AND.doubleq.AND.scavenging)
	656
[38]	657	c Submicron dust mode:
	658	c --------------------
	659
	660	if (submicron) then
	661	radius(igcm_dust_submicron)=0.1e-6
	662	rho_q(igcm_dust_submicron)=rho_dust
	663	if (doubleq) then
	664	c If doubleq is also active, we use the population ratio:
	665	alpha_lift(igcm_dust_submicron) =
	666	& alpha_lift(igcm_dust_number)popratio
	667	& rho_q(igcm_dust_submicron)4./3.pi*
	668	& radius(igcm_dust_submicron)**3.
	669	alpha_devil(igcm_dust_submicron)=1.e-30
	670	else
	671	alpha_lift(igcm_dust_submicron)=1e-6
	672	alpha_devil(igcm_dust_submicron)=1.e-30
	673	endif ! (doubleq)
	674	end if ! (submicron)
	675
	676	c Initialization for water vapor
	677	c ------------------------------
	678	if(water) then
	679	radius(igcm_h2o_vap)=0.
	680	alpha_lift(igcm_h2o_vap) =0.
	681	alpha_devil(igcm_h2o_vap)=0.
[1036]	682	if(water.and.(nq.ge.2)) then
[38]	683	radius(igcm_h2o_ice)=3.e-6
	684	rho_q(igcm_h2o_ice)=rho_ice
	685	alpha_lift(igcm_h2o_ice) =0.
	686	alpha_devil(igcm_h2o_ice)=0.
[1036]	687	elseif(water.and.(nq.lt.2)) then
	688	write(,) 'nq is too low : nq=', nq
[38]	689	write(,) 'water= ',water
	690	endif
	691
	692	end if ! (water)
[2312]	693
	694	c Initialization for hdo vapor
	695	c ------------------------------
	696	if (hdo) then
	697	radius(igcm_hdo_vap)=0.
	698	alpha_lift(igcm_hdo_vap) =0.
	699	alpha_devil(igcm_hdo_vap)=0.
	700	if(water.and.(nq.ge.2)) then
	701	radius(igcm_hdo_ice)=3.e-6
	702	rho_q(igcm_hdo_ice)=rho_ice
	703	alpha_lift(igcm_hdo_ice) =0.
	704	alpha_devil(igcm_hdo_ice)=0.
	705	elseif(hdo.and.(nq.lt.6)) then
	706	write(,) 'nq is too low : nq=', nq
	707	write(,) 'hdo= ',hdo
	708	endif
	709
	710	end if ! (hdo)
	711
[1617]	712
[38]	713
	714	c Output for records:
	715	c ~~~~~~~~~~~~~~~~~~
	716	write(,)
	717	Write(,) '******** initracer : dust transport parameters :'
	718	write(,) 'alpha_lift = ', alpha_lift
	719	write(,) 'alpha_devil = ', alpha_devil
	720	write(,) 'radius = ', radius
	721	if(doubleq) then
	722	write(,) 'reff_lift (um) = ', reff_lift
	723	write(,) 'size distribution variance = ', varian
	724	write(,) 'r3n_q , ref_r0 : ', r3n_q , ref_r0
	725	end if
	726
	727	!
	728	! some extra (possibly redundant) sanity checks for tracers:
	729	! ---------------------------------------------------------
	730
	731	if (doubleq) then
	732	! verify that we indeed have dust_mass and dust_number tracers
	733	if (igcm_dust_mass.eq.0) then
	734	write(,) "initracer: error !!"
	735	write(,) " cannot use doubleq option without ",
	736	& "a dust_mass tracer !"
[2302]	737	call abort_physic("initracer","doubleq issue",1)
[38]	738	endif
	739	if (igcm_dust_number.eq.0) then
	740	write(,) "initracer: error !!"
	741	write(,) " cannot use doubleq option without ",
	742	& "a dust_number tracer !"
[2302]	743	call abort_physic("initracer","doubleq issue",1)
[38]	744	endif
	745	endif
	746
	747	if ((.not.doubleq).and.(dustbin.gt.0)) then
	748	! verify that we indeed have 'dustbin' dust tracers
	749	count=0
	750	do iq=1,dustbin
	751	if (igcm_dustbin(iq).ne.0) then
	752	count=count+1
	753	endif
	754	enddo
	755	if (count.ne.dustbin) then
	756	write(,) "initracer: error !!"
[2302]	757	write(,) " dustbin is set to ",dustbin,
[38]	758	& " but we only have the following dust tracers:"
	759	do iq=1,count
	760	write(,)" ",trim(noms(igcm_dustbin(iq)))
	761	enddo
[2302]	762	call abort_physic("initracer","dustbin issue",1)
[38]	763	endif
	764	endif
	765
	766	if (water) then
	767	! verify that we indeed have h2o_vap and h2o_ice tracers
	768	if (igcm_h2o_vap.eq.0) then
	769	write(,) "initracer: error !!"
	770	write(,) " cannot use water option without ",
	771	& "an h2o_vap tracer !"
[2302]	772	call abort_physic("initracer","water cycle issue",1)
[38]	773	endif
	774	if (igcm_h2o_ice.eq.0) then
	775	write(,) "initracer: error !!"
	776	write(,) " cannot use water option without ",
	777	& "an h2o_ice tracer !"
[2302]	778	call abort_physic("initracer","water cycle issue",1)
[38]	779	endif
	780	endif
	781
[2312]	782	if (hdo) then
	783	! verify that we indeed have hdo_vap and hdo_ice tracers
	784	if (igcm_hdo_vap.eq.0) then
	785	write(,) "initracer: error !!"
	786	write(,) " cannot use hdo option without ",
	787	& "an hdo_vap tracer !"
[2398]	788	call abort_physic("initracer","hdo cycle issue",1)
[2312]	789	endif
	790	if (igcm_hdo_ice.eq.0) then
	791	write(,) "initracer: error !!"
	792	write(,) " cannot use hdo option without ",
	793	& "an hdo_ice tracer !"
[2398]	794	call abort_physic("initracer","hdo cycle issue",1)
[2312]	795	endif
	796	endif
	797
	798
[1617]	799	if (co2clouds) then
	800	!verify that we have co2_ice and co2 tracers
	801	if (igcm_co2 .eq. 0) then
	802	write(,) "initracer: error !!"
	803	write(,) " cannot use co2 clouds option without ",
	804	& "a co2 tracer !"
[2302]	805	call abort_physic("initracer","co2 clouds issue",1)
[2589]	806	end if
[1617]	807	if (igcm_co2_ice .eq. 0) then
	808	write(,) "initracer: error !!"
	809	write(,) " cannot use co2 clouds option without ",
	810	& "a co2_ice tracer !"
[2302]	811	call abort_physic("initracer","co2 clouds issue",1)
[2589]	812	end if
[2659]	813	if (igcm_ccnco2_number .eq. 0) then
	814	write(,) "initracer: error !!"
	815	write(,) " cannot use co2 clouds option without ",
	816	& "a ccnco2_number tracer !"
	817	call abort_physic("initracer","co2 clouds issue",1)
	818	end if
	819	if (igcm_ccnco2_mass .eq. 0) then
	820	write(,) "initracer: error !!"
	821	write(,) " cannot use co2 clouds option without ",
	822	& "a ccnco2_mass tracer !"
	823	call abort_physic("initracer","co2 clouds issue",1)
	824	end if
	825	if (co2useh2o) then
	826	if (igcm_ccnco2_h2o_number .eq. 0) then
	827	write(,) "initracer: error !!"
	828	write(,) " cannot use co2 clouds option without ",
	829	& "a ccnco2_h2o_number tracer !"
	830	call abort_physic("initracer","co2 clouds issue",1)
	831	end if
	832	if (igcm_ccnco2_h2o_mass_ice .eq. 0) then
	833	write(,) "initracer: error !!"
	834	write(,) " cannot use co2 clouds option without ",
	835	& "a ccnco2_h2o_mass_ice tracer !"
	836	call abort_physic("initracer","co2 clouds issue",1)
	837	end if
	838	if (igcm_ccnco2_h2o_mass_ccn .eq. 0) then
	839	write(,) "initracer: error !!"
	840	write(,) " cannot use co2 clouds option without ",
	841	& "a ccnco2_h2o_mass_ccn tracer !"
	842	call abort_physic("initracer","co2 clouds issue",1)
	843	end if
	844	end if
	845	if (meteo_flux) then
	846	if (igcm_ccnco2_meteor_number .eq. 0) then
	847	write(,) "initracer: error !!"
	848	write(,) " cannot use co2 clouds option without ",
	849	& "a ccnco2_meteor_number tracer !"
	850	call abort_physic("initracer","co2 clouds issue",1)
	851	end if
	852	if (igcm_ccnco2_meteor_mass .eq. 0) then
	853	write(,) "initracer: error !!"
	854	write(,) " cannot use co2 clouds option without ",
	855	& "a ccnco2_h2o_mass_ice tracer !"
	856	call abort_physic("initracer","co2 clouds issue",1)
	857	end if
	858	if (igcm_ccnco2_meteor_mass .eq. 0) then
	859	write(,) "initracer: error !!"
	860	write(,) " cannot use co2 clouds option without ",
	861	& "a ccnco2_meteor_mass tracer !"
	862	call abort_physic("initracer","co2 clouds issue",1)
	863	end if
	864	end if
[1617]	865	endif
[1974]	866
	867	if (rdstorm) then
	868	! verify that we indeed have stormdust_mass and stormdust_number tracers
	869	if (igcm_stormdust_mass.eq.0) then
	870	write(,) "initracer: error !!"
	871	write(,) " cannot use rdstorm option without ",
	872	& "a stormdust_mass tracer !"
[2302]	873	call abort_physic("initracer","rdstorm issue",1)
[1974]	874	endif
	875	if (igcm_stormdust_number.eq.0) then
	876	write(,) "initracer: error !!"
	877	write(,) " cannot use rdstorm option without ",
	878	& "a stormdust_number tracer !"
[2302]	879	call abort_physic("initracer","rdstorm issue",1)
[1974]	880	endif
	881	endif
[2199]	882
[2628]	883	if (topflows) then
[2199]	884	! verify that we indeed have topdust_mass and topdust_number tracers
	885	if (igcm_topdust_mass.eq.0) then
	886	write(,) "initracer: error !!"
[2628]	887	write(,) " cannot use topflows option without ",
[2199]	888	& "a topdust_mass tracer !"
[2628]	889	call abort_physic("initracer","topflows issue",1)
[2199]	890	endif
	891	if (igcm_topdust_number.eq.0) then
	892	write(,) "initracer: error !!"
[2628]	893	write(,) " cannot use topflows option without ",
[2199]	894	& "a topdust_number tracer !"
[2628]	895	call abort_physic("initracer","topflows issue",1)
[2199]	896	endif
	897	endif
[1974]	898
[1380]	899	if (callnlte) then ! NLTE requirements
	900	if (nltemodel.ge.1) then
	901	! check that co2, co, o and n2 tracers are available
	902	if (igcm_co2.eq.0) then
	903	write(,) "initracer: error !!"
	904	write(,) " with nltemodel>0, we need the co2 tracer!"
[2302]	905	call abort_physic("initracer","missing co2 tracer",1)
[1380]	906	endif
	907	if (igcm_co.eq.0) then
	908	write(,) "initracer: error !!"
	909	write(,) " with nltemodel>0, we need the co tracer!"
[2302]	910	call abort_physic("initracer","missing co tracer",1)
[1380]	911	endif
	912	if (igcm_o.eq.0) then
	913	write(,) "initracer: error !!"
	914	write(,) " with nltemodel>0, we need the o tracer!"
[2302]	915	call abort_physic("initracer","missing o tracer",1)
[1380]	916	endif
	917	if (igcm_n2.eq.0) then
	918	write(,) "initracer: error !!"
	919	write(,) " with nltemodel>0, we need the n2 tracer!"
[2302]	920	call abort_physic("initracer","missing n2 tracer",1)
[1380]	921	endif
	922	endif
	923	endif
	924
[358]	925	if (scavenging) then
	926	! verify that we indeed have ccn_mass and ccn_number tracers
[1617]	927	if (igcm_ccn_mass.eq.0 .and. igcm_ccnco2_mass.eq.0) then
[358]	928	write(,) "initracer: error !!"
	929	write(,) " cannot use scavenging option without ",
[1617]	930	& "a ccn_mass or ccnco2_mass tracer !"
[2302]	931	call abort_physic("initracer","scavenging issue",1)
[358]	932	endif
[1617]	933	if (igcm_ccn_number.eq.0 .and. igcm_ccnco2_number.eq.0 ) then
[358]	934	write(,) "initracer: error !!"
	935	write(,) " cannot use scavenging option without ",
[1617]	936	& "a ccn_number or ccnco2_number tracer !"
[2302]	937	call abort_physic("initracer","scavenging issue",1)
[358]	938	endif
	939	endif ! of if (scavenging)
	940
[38]	941	if (photochem .or. callthermos) then
	942	! verify that we indeed have the chemistry tracers
	943	if (igcm_co2.eq.0) then
	944	write(,) "initracer: error !!"
	945	write(,) " cannot use chemistry option without ",
	946	& "a co2 tracer !"
[2302]	947	call abort_physic("initracer","missing co2 tracer",1)
[38]	948	endif
	949	if (igcm_co.eq.0) then
	950	write(,) "initracer: error !!"
	951	write(,) " cannot use chemistry option without ",
	952	& "a co tracer !"
[2302]	953	call abort_physic("initracer","missing co tracer",1)
[38]	954	endif
	955	if (igcm_o.eq.0) then
	956	write(,) "initracer: error !!"
	957	write(,) " cannot use chemistry option without ",
	958	& "a o tracer !"
[2302]	959	call abort_physic("initracer","missing o tracer",1)
[38]	960	endif
	961	if (igcm_o1d.eq.0) then
	962	write(,) "initracer: error !!"
	963	write(,) " cannot use chemistry option without ",
	964	& "a o1d tracer !"
[2302]	965	call abort_physic("initracer","missing o1d tracer",1)
[38]	966	endif
	967	if (igcm_o2.eq.0) then
	968	write(,) "initracer: error !!"
	969	write(,) " cannot use chemistry option without ",
	970	& "an o2 tracer !"
[2302]	971	call abort_physic("initracer","missing o2 tracer",1)
[38]	972	endif
	973	if (igcm_o3.eq.0) then
	974	write(,) "initracer: error !!"
	975	write(,) " cannot use chemistry option without ",
	976	& "an o3 tracer !"
[2302]	977	call abort_physic("initracer","missing o3 tracer",1)
[38]	978	endif
	979	if (igcm_h.eq.0) then
	980	write(,) "initracer: error !!"
	981	write(,) " cannot use chemistry option without ",
	982	& "a h tracer !"
[2302]	983	call abort_physic("initracer","missing h tracer",1)
[38]	984	endif
	985	if (igcm_h2.eq.0) then
	986	write(,) "initracer: error !!"
	987	write(,) " cannot use chemistry option without ",
	988	& "a h2 tracer !"
[2302]	989	call abort_physic("initracer","missing h2 tracer",1)
[38]	990	endif
	991	if (igcm_oh.eq.0) then
	992	write(,) "initracer: error !!"
	993	write(,) " cannot use chemistry option without ",
	994	& "an oh tracer !"
[2302]	995	call abort_physic("initracer","missing oh tracer",1)
[38]	996	endif
	997	if (igcm_ho2.eq.0) then
	998	write(,) "initracer: error !!"
	999	write(,) " cannot use chemistry option without ",
	1000	& "a ho2 tracer !"
[2302]	1001	call abort_physic("initracer","missing ho2 tracer",1)
[1720]	1002	endif
[38]	1003	if (igcm_h2o2.eq.0) then
	1004	write(,) "initracer: error !!"
	1005	write(,) " cannot use chemistry option without ",
	1006	& "a h2o2 tracer !"
[2302]	1007	call abort_physic("initracer","missing h2o2 tracer",1)
[38]	1008	endif
	1009	if (igcm_n2.eq.0) then
	1010	write(,) "initracer: error !!"
	1011	write(,) " cannot use chemistry option without ",
	1012	& "a n2 tracer !"
[2302]	1013	call abort_physic("initracer","missing n2 tracer",1)
[38]	1014	endif
	1015	if (igcm_ar.eq.0) then
	1016	write(,) "initracer: error !!"
	1017	write(,) " cannot use chemistry option without ",
	1018	& "an ar tracer !"
[2302]	1019	call abort_physic("initracer","missing ar tracer",1)
[38]	1020	endif
	1021	endif ! of if (photochem .or. callthermos)
	1022
[2659]	1023	! Initialisation for CO2 clouds
	1024	if (co2clouds) then
	1025	radius(igcm_ccnco2_mass) = radius(igcm_dust_mass)
	1026	alpha_lift(igcm_ccnco2_mass) = 1e-30
	1027	alpha_devil(igcm_ccnco2_mass) = 1e-30
	1028	rho_q(igcm_ccnco2_mass) = rho_dust
	1029	radius(igcm_ccnco2_number) = radius(igcm_ccnco2_mass)
	1030	alpha_lift(igcm_ccnco2_number) = alpha_lift(igcm_ccnco2_mass)
	1031	alpha_devil(igcm_ccnco2_number) = alpha_devil(igcm_ccnco2_mass)
	1032	rho_q(igcm_ccnco2_number) = rho_q(igcm_ccnco2_mass)
	1033
	1034	radius(igcm_co2)=0.
	1035	alpha_lift(igcm_co2) =0.
	1036	alpha_devil(igcm_co2)=0.
	1037	radius(igcm_co2_ice)=1.e-8
	1038	rho_q(igcm_co2_ice)=rho_ice_co2
	1039	alpha_lift(igcm_co2_ice) =0.
	1040	alpha_devil(igcm_co2_ice)=0.
	1041
	1042	endif
	1043
[38]	1044	end

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

Download in other formats:

Original Format