Context Navigation

initracer.F @ 2628

Last change on this file since 2628 was 2628, checked in by abierjon, 3 years ago

Mars GCM:
Big changes on mountain top dust flows for GCM6:

the scheme now activates only in grid meshes that contain a mountain among a hard-written list, instead of every meshes. This is done to prevent strong artificial reinjections of dust in places that don't present huge converging slopes enabling the concentration of dust (ex: Valles Marineris, Hellas). Topdust is now also detrained as soon as it leaves the column it originated from.
the list of the 19 allowed mountains is used by the subroutine topmons_setup in module topmons_mod, to compute a logical array contains_mons(ngrid). alpha_hmons and hsummit are also set up once and for all by this subroutine, which is called in physiq_mod's firstcall.
contains_mons, alpha_hmons and hsummit are now saved variables of the module surfdat_h, and are called as such and not as arguments in the subroutines using them.
the logical flag "slpwind" and the comments in the code have also been updated to the new terminology "mountain top dust flows", accordingly to ticket #71. The new flag read in callphys.def is "topflows".

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

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

Download in other formats:

Original Format