Context Navigation

initracer.F @ 524

Last change on this file since 524 was 459, checked in by emillour, 13 years ago

Mars GCM: Include changes and updates for photochemistry by FL:

aeronomars/calchim.F : change in units of surface density.
aeronomars/surfacearea.F : new routine to compute ice and dust surface area

(m2/m3) available for heterogeneous reactions.

phymars/initracer.F : bug correction: initialize igcm_ch4 and change loop

bounds (when guessing tracer names/properties with
old input files).

phymars/watercloud.F : cleanup.
phymars/physiq.F : add call to surfacearea; photochemistry is now called

after sedimentation to take into acount updated rdust
and rice.

File size: 25.2 KB

Rev	Line
[38]	1	SUBROUTINE initracer(qsurf,co2ice)
	2
	3	IMPLICIT NONE
	4	c=======================================================================
	5	c subject:
	6	c --------
	7	c Initialization related to tracer
	8	c (transported dust, water, chemical species, ice...)
	9	c
	10	c Name of the tracer
	11	c
	12	c Test of dimension :
	13	c Initialize COMMON tracer in tracer.h, using tracer names provided
	14	c by the dynamics in "advtrac.h"
	15	c
	16	c Old conventions: (not used any more)
	17	c
	18	c If water=T : q(iq=nqmx) is the water mass mixing ratio
	19	c and q(iq=nqmx-1) is the ice mass mixing ratio
	20
	21	c If there is transported dust, it uses iq=1 to iq=dustbin
	22	c If there is no transported dust : dustbin=0
	23	c If doubleq=T : q(iq=1) is the dust mass mixing ratio
	24	c q(iq=2) is the dust number mixing ratio
	25
	26	c If (photochem.or.thermochem) there is "ncomp" chemical species (ncomp
	27	c is set in aeronomars/chimiedata.h) using the ncomp iq values starting at
	28	c iq=nqchem_min = dustbin+1 (nqchem_min is defined in inifis.F)
	29	c
	30	c
	31	c author: F.Forget
	32	c ------
	33	c Modifs: Franck Montmessin, Sebastien Lebonnois (june 2003)
	34	c Ehouarn Millour (oct. 2008) identify tracers by their names
	35	c=======================================================================
	36
	37
	38	#include "dimensions.h"
	39	#include "dimphys.h"
	40	#include "comcstfi.h"
	41	#include "callkeys.h"
	42	#include "tracer.h"
	43	#include "advtrac.h"
	44	#include "comgeomfi.h"
	45	#include "chimiedata.h"
	46
[91]	47	#include "surfdat.h"
[38]	48
	49	real qsurf(ngridmx,nqmx) ! tracer on surface (e.g. kg.m-2)
	50	real co2ice(ngridmx) ! co2 ice mass on surface (e.g. kg.m-2)
[283]	51	integer iq,ig,count
[38]	52	real r0_lift , reff_lift, nueff_lift
	53	c Ratio of small over large dust particles (used when both
	54	c doubleq and the submicron mode are active); In Montmessin
	55	c et al. (2002), a value of 25 has been deduced;
	56	real, parameter :: popratio = 25.
	57	logical :: oldnames ! =.true. if old tracer naming convention (q01,...)
	58	character(len=20) :: txt ! to store some text
	59
	60	c-----------------------------------------------------------------------
	61	c radius(nqmx) ! aerosol particle radius (m)
	62	c rho_q(nqmx) ! tracer densities (kg.m-3)
	63	c alpha_lift(nqmx) ! saltation vertical flux/horiz flux ratio (m-1)
	64	c alpha_devil(nqmx) ! lifting coeeficient by dust devil
	65	c rho_dust ! Mars dust density
	66	c rho_ice ! Water ice density
	67	c nuice_ref ! Effective variance nueff of the
	68	c ! water-ice size distributions
	69	c doubleq ! if method with mass (iq=1) and number(iq=2) mixing ratio
	70	c varian ! Characteristic variance of log-normal distribution
	71	c-----------------------------------------------------------------------
	72
	73	integer :: nqchem_start
	74
	75	! Initialization: get tracer names from the dynamics and check if we are
	76	! using 'old' tracer convention ('q01',q02',...)
	77	! or new convention (full tracer names)
	78	! check if tracers have 'old' names
	79	oldnames=.false.
	80	count=0
	81	do iq=1,nqmx
	82	txt=" "
	83	write(txt,'(a1,i2.2)') 'q',iq
	84	if (txt.eq.tnom(iq)) then
	85	count=count+1
	86	endif
	87	enddo ! of do iq=1,nqmx
	88
	89	if (count.eq.nqmx) then
	90	write(,) "initracer: tracers seem to follow old naming ",
	91	& "convention (q01,q02,...)"
	92	write(,) " => will work for now ... "
	93	write(,) " but you should run newstart to rename them"
	94	oldnames=.true.
	95	endif
	96
	97	! copy/set tracer names
	98	if (oldnames) then ! old names (derived from iq & option values)
	99	! 1. dust:
	100	if (dustbin.ne.0) then ! transported dust
	101	do iq=1,dustbin
	102	txt=" "
	103	write(txt,'(a4,i2.2)') 'dust',iq
	104	noms(iq)=txt
	105	mmol(iq)=100.
	106	enddo
	107	! special case if doubleq
	108	if (doubleq) then
	109	if (dustbin.ne.2) then
	110	write(,) 'initracer: error expected dustbin=2'
	111	else
[164]	112	! noms(1)='dust_mass' ! dust mass mixing ratio
	113	! noms(2)='dust_number' ! dust number mixing ratio
	114	! same as above, but avoid explict possible out-of-bounds on array
	115	noms(1)='dust_mass' ! dust mass mixing ratio
	116	do iq=2,2
	117	noms(iq)='dust_number' ! dust number mixing ratio
	118	enddo
[38]	119	endif
	120	endif
	121	endif
	122	! 2. water & ice
	123	if (water) then
	124	noms(nqmx)='h2o_vap'
	125	mmol(nqmx)=18.
[164]	126	! noms(nqmx-1)='h2o_ice'
	127	! mmol(nqmx-1)=18.
	128	!"loop" to avoid potential out-of-bounds in array
	129	do iq=nqmx-1,nqmx-1
	130	noms(iq)='h2o_ice'
	131	mmol(iq)=18.
	132	enddo
[38]	133	endif
	134	! 3. Chemistry
	135	if (photochem .or. callthermos) then
	136	if (doubleq) then
	137	nqchem_start=3
	138	else
	139	nqchem_start=dustbin+1
	140	end if
[459]	141	do iq=nqchem_start,nqchem_start+ncomp-2
[38]	142	noms(iq)=nomchem(iq-nqchem_start+1)
	143	mmol(iq)=mmolchem(iq-nqchem_start+1)
	144	enddo
	145	endif ! of if (photochem .or. callthermos)
	146	! 4. Other tracers ????
	147	if ((dustbin.eq.0).and.(.not.water)) then
	148	noms(1)='co2'
	149	mmol(1)=44
	150	if (nqmx.eq.2) then
	151	noms(nqmx)='Ar_N2'
	152	mmol(nqmx)=30
	153	endif
	154	endif
	155	else
	156	! copy tracer names from dynamics
	157	do iq=1,nqmx
	158	noms(iq)=tnom(iq)
	159	enddo
	160	! mmol(:) array is initialized later (see below)
	161	endif ! of if (oldnames)
	162
	163	! special modification when using 'old' tracers:
	164	! qsurf(nqmx) was h2o ice whereas q(nqmx) was water vapour
	165	! and (if iceparty) q(nqmx-1) was null whereas q(nqmx-1) was water ice
	166	if (oldnames.and.water) then
	167	write(,)'initracer: moving surface water ice to index ',nqmx-1
[164]	168	! "loop" to avoid potential out-of-bounds on arrays
	169	do iq=nqmx-1,nqmx-1
	170	qsurf(1:ngridmx,iq)=qsurf(1:ngridmx,iq+1)
	171	enddo
[38]	172	qsurf(1:ngridmx,nqmx)=0
	173	endif
	174
	175	c------------------------------------------------------------
	176	c Test Dimensions tracers
	177	c------------------------------------------------------------
	178
	179	! Ehouarn: testing number of tracers is obsolete...
	180	! if(photochem.or.thermochem) then
	181	! if (water) then
	182	! if ((dustbin+ncomp+2).ne.nqmx) then
	183	! print*,'initracer: tracer dimension problem:'
	184	! print*,"(dustbin+ncomp+2).ne.nqmx"
	185	! print*,"ncomp: ",ncomp
	186	! print*,"dustbin: ",dustbin
	187	! print*,"nqmx: ",nqmx
	188	! print*,'Change ncomp in chimiedata.h'
	189	! endif
	190	! else
	191	! if ((dustbin+ncomp+1).ne.nqmx) then
	192	! print*,'initracer: tracer dimension problem:'
	193	! print*,"(dustbin+ncomp+1).ne.nqmx"
	194	! print*,"ncomp: ",ncomp
	195	! print*,"dustbin: ",dustbin
	196	! print*,"nqmx: ",nqmx
	197	! print*,'Change ncomp in chimiedata.h'
	198	! STOP
	199	! endif
	200	! endif
	201	! endif
	202
	203	c------------------------------------------------------------
	204	c NAME and molar mass of the tracer
	205	c------------------------------------------------------------
	206
	207
	208	! Identify tracers by their names: (and set corresponding values of mmol)
	209	! 0. initialize tracer indexes to zero:
	210	do iq=1,nqmx
	211	igcm_dustbin(iq)=0
	212	enddo
	213	igcm_dust_mass=0
	214	igcm_dust_number=0
[358]	215	igcm_ccn_mass=0
	216	igcm_ccn_number=0
[38]	217	igcm_dust_submicron=0
	218	igcm_h2o_vap=0
	219	igcm_h2o_ice=0
	220	igcm_co2=0
	221	igcm_co=0
	222	igcm_o=0
	223	igcm_o1d=0
	224	igcm_o2=0
	225	igcm_o3=0
	226	igcm_h=0
	227	igcm_h2=0
	228	igcm_oh=0
	229	igcm_ho2=0
	230	igcm_h2o2=0
[459]	231	igcm_ch4=0
[38]	232	igcm_n2=0
	233	igcm_ar=0
	234	igcm_ar_n2=0
[171]	235	igcm_n=0
	236	igcm_no=0
	237	igcm_no2=0
	238	igcm_n2d=0
	239	igcm_co2plus=0
	240	igcm_oplus=0
	241	igcm_o2plus=0
	242	igcm_coplus=0
	243	igcm_cplus=0
	244	igcm_nplus=0
	245	igcm_noplus=0
	246	igcm_n2plus=0
	247	igcm_hplus=0
	248	igcm_elec=0
[38]	249
[171]	250
[38]	251	! 1. find dust tracers
	252	count=0
	253	if (dustbin.gt.0) then
	254	do iq=1,nqmx
	255	txt=" "
	256	write(txt,'(a4,i2.2)')'dust',count+1
	257	if (noms(iq).eq.txt) then
	258	count=count+1
	259	igcm_dustbin(count)=iq
	260	mmol(iq)=100.
	261	endif
	262	enddo !do iq=1,nqmx
	263	endif ! of if (dustbin.gt.0)
	264	if (doubleq) then
	265	do iq=1,nqmx
	266	if (noms(iq).eq."dust_mass") then
	267	igcm_dust_mass=iq
	268	count=count+1
	269	endif
	270	if (noms(iq).eq."dust_number") then
	271	igcm_dust_number=iq
	272	count=count+1
	273	endif
	274	enddo
	275	endif ! of if (doubleq)
[358]	276	if (scavenging) then
	277	do iq=1,nqmx
	278	if (noms(iq).eq."ccn_mass") then
	279	igcm_ccn_mass=iq
	280	count=count+1
	281	endif
	282	if (noms(iq).eq."ccn_number") then
	283	igcm_ccn_number=iq
	284	count=count+1
	285	endif
	286	enddo
	287	endif ! of if (scavenging)
[38]	288	if (submicron) then
	289	do iq=1,nqmx
	290	if (noms(iq).eq."dust_submicron") then
	291	igcm_dust_submicron=iq
	292	mmol(iq)=100.
	293	count=count+1
	294	endif
	295	enddo
	296	endif ! of if (submicron)
	297	! 2. find chemistry and water tracers
	298	do iq=1,nqmx
	299	if (noms(iq).eq."co2") then
	300	igcm_co2=iq
	301	mmol(igcm_co2)=44.
	302	count=count+1
	303	endif
	304	if (noms(iq).eq."co") then
	305	igcm_co=iq
	306	mmol(igcm_co)=28.
	307	count=count+1
	308	endif
	309	if (noms(iq).eq."o") then
	310	igcm_o=iq
	311	mmol(igcm_o)=16.
	312	count=count+1
	313	endif
	314	if (noms(iq).eq."o1d") then
	315	igcm_o1d=iq
	316	mmol(igcm_o1d)=16.
	317	count=count+1
	318	endif
	319	if (noms(iq).eq."o2") then
	320	igcm_o2=iq
	321	mmol(igcm_o2)=32.
	322	count=count+1
	323	endif
	324	if (noms(iq).eq."o3") then
	325	igcm_o3=iq
	326	mmol(igcm_o3)=48.
	327	count=count+1
	328	endif
	329	if (noms(iq).eq."h") then
	330	igcm_h=iq
	331	mmol(igcm_h)=1.
	332	count=count+1
	333	endif
	334	if (noms(iq).eq."h2") then
	335	igcm_h2=iq
	336	mmol(igcm_h2)=2.
	337	count=count+1
	338	endif
	339	if (noms(iq).eq."oh") then
	340	igcm_oh=iq
	341	mmol(igcm_oh)=17.
	342	count=count+1
	343	endif
	344	if (noms(iq).eq."ho2") then
	345	igcm_ho2=iq
	346	mmol(igcm_ho2)=33.
	347	count=count+1
	348	endif
	349	if (noms(iq).eq."h2o2") then
	350	igcm_h2o2=iq
	351	mmol(igcm_h2o2)=34.
	352	count=count+1
	353	endif
	354	if (noms(iq).eq."n2") then
	355	igcm_n2=iq
	356	mmol(igcm_n2)=28.
	357	count=count+1
	358	endif
[324]	359	if (noms(iq).eq."ch4") then
	360	igcm_ch4=iq
	361	mmol(igcm_ch4)=16.
	362	count=count+1
	363	endif
[38]	364	if (noms(iq).eq."ar") then
	365	igcm_ar=iq
	366	mmol(igcm_ar)=40.
	367	count=count+1
	368	endif
[324]	369	if (noms(iq).eq."n") then
[171]	370	igcm_n=iq
	371	mmol(igcm_n)=14.
	372	count=count+1
	373	endif
	374	if (noms(iq).eq."no") then
	375	igcm_no=iq
	376	mmol(igcm_no)=30.
	377	count=count+1
	378	endif
	379	if (noms(iq).eq."no2") then
	380	igcm_no2=iq
	381	mmol(igcm_no2)=46.
	382	count=count+1
	383	endif
	384	if (noms(iq).eq."n2d") then
	385	igcm_n2d=iq
	386	mmol(igcm_n2d)=28.
	387	count=count+1
	388	endif
	389	if (noms(iq).eq."co2plus") then
	390	igcm_co2plus=iq
	391	mmol(igcm_co2plus)=44.
	392	count=count+1
	393	endif
	394	if (noms(iq).eq."oplus") then
	395	igcm_oplus=iq
	396	mmol(igcm_oplus)=16.
	397	count=count+1
	398	endif
	399	if (noms(iq).eq."o2plus") then
	400	igcm_o2plus=iq
	401	mmol(igcm_o2plus)=32.
	402	count=count+1
	403	endif
	404	if (noms(iq).eq."coplus") then
	405	igcm_coplus=iq
	406	mmol(igcm_coplus)=28.
	407	count=count+1
	408	endif
	409	if (noms(iq).eq."cplus") then
	410	igcm_cplus=iq
	411	mmol(igcm_cplus)=12.
	412	count=count+1
	413	endif
	414	if (noms(iq).eq."nplus") then
	415	igcm_nplus=iq
	416	mmol(igcm_nplus)=14.
	417	count=count+1
	418	endif
	419	if (noms(iq).eq."noplus") then
	420	igcm_noplus=iq
	421	mmol(igcm_noplus)=30.
	422	count=count+1
	423	endif
[324]	424	if (noms(iq).eq."n2plus") then
[171]	425	igcm_n2plus=iq
	426	mmol(igcm_n2plus)=28.
	427	count=count+1
	428	endif
	429	if (noms(iq).eq."hplus") then
	430	igcm_hplus=iq
	431	mmol(igcm_hplus)=1.
	432	count=count+1
	433	endif
	434	if (noms(iq).eq."elec") then
	435	igcm_elec=iq
	436	mmol(igcm_elec)=1./1822.89
	437	count=count+1
	438	endif
[38]	439	if (noms(iq).eq."h2o_vap") then
	440	igcm_h2o_vap=iq
	441	mmol(igcm_h2o_vap)=18.
	442	count=count+1
	443	endif
	444	if (noms(iq).eq."h2o_ice") then
	445	igcm_h2o_ice=iq
	446	mmol(igcm_h2o_ice)=18.
	447	count=count+1
	448	endif
	449	! Other stuff: e.g. for simulations using co2 + neutral gaz
	450	if (noms(iq).eq."Ar_N2") then
	451	igcm_ar_n2=iq
	452	mmol(igcm_ar_n2)=30.
	453	count=count+1
	454	endif
[171]	455
[324]	456
[38]	457	enddo ! of do iq=1,nqmx
	458	! count=count+nbqchem
	459
	460	! check that we identified all tracers:
	461	if (count.ne.nqmx) then
	462	write(,) "initracer: found only ",count," tracers"
	463	write(,) " expected ",nqmx
	464	do iq=1,count
	465	write(,)' ',iq,' ',trim(noms(iq))
	466	enddo
	467	stop
	468	else
	469	write(,) "initracer: found all expected tracers, namely:"
	470	do iq=1,nqmx
	471	write(,)' ',iq,' ',trim(noms(iq))
	472	enddo
	473	endif
	474
	475	! if water cycle but iceparty=.false., there will nevertheless be
	476	! water ice at the surface (iceparty is not used anymore, but this
	477	! part is still relevant, as we want to stay compatible with the
	478	! older versions).
	479	if (water.and.(igcm_h2o_ice.eq.0)) then
	480	igcm_h2o_ice=igcm_h2o_vap ! so that qsurf(i_h2o_ice) is identified
	481	! even though there is no q(i_h2o_ice)
	482	else
	483	! surface ice qsurf(i_h2o_ice) was loaded twice by phyetat0,
	484	! as qsurf(i_h2o_vap) & as qsurf(i_h2o_ice), so to be clean:
	485	if (igcm_h2o_vap.ne.0) then
	486	qsurf(1:ngridmx,igcm_h2o_vap)=0
	487	endif
	488	endif
	489
	490	c------------------------------------------------------------
	491	c Initialisation tracers ....
	492	c------------------------------------------------------------
	493	call zerophys(nqmx,rho_q)
	494
	495	rho_dust=2500. ! Mars dust density (kg.m-3)
	496	rho_ice=920. ! Water ice density (kg.m-3)
	497	nuice_ref=0.1 ! Effective variance nueff of the
[358]	498	! water-ice size distribution
[455]	499	!!!nuice_sed=0.45 ! Sedimentation effective variance
[358]	500	! of the water-ice size distribution
[38]	501
	502	if (doubleq) then
	503	c "doubleq" technique
	504	c -------------------
	505	c (transport of mass and number mixing ratio)
	506	c iq=1: Q mass mixing ratio, iq=2: N number mixing ratio
	507
	508	if( (nqmx.lt.2).or.(water.and.(nqmx.lt.4)) ) then
	509	write(,)'initracer: nqmx is too low : nqmx=', nqmx
	510	write(,)'water= ',water,' doubleq= ',doubleq
	511	end if
	512
	513	nueff_lift = 0.5
	514	varian=sqrt(log(1.+nueff_lift))
	515
	516	rho_q(igcm_dust_mass)=rho_dust
	517	rho_q(igcm_dust_number)=rho_dust
	518
	519	c Intermediate calcul for computing geometric mean radius r0
	520	c as a function of mass and number mixing ratio Q and N
	521	c (r0 = (r3n_q * Q/ N)^(1/3))
	522	r3n_q = exp(-4.5varian2)(3./4.)/(pi*rho_dust)
	523
	524	c Intermediate calcul for computing effective radius reff
	525	c from geometric mean radius r0
	526	c (reff = ref_r0 * r0)
	527	ref_r0 = exp(2.5varian*2)
	528
	529	c lifted dust :
	530	c '''''''''''
	531	reff_lift = 3.0e-6 !3.e-6 !Effective radius of lifted dust (m)
	532	alpha_devil(igcm_dust_mass)=9.e-9 ! dust devil lift mass coeff
	533	c alpha_lift(igcm_dust_mass)=3.0e-15 ! Lifted mass coeff
	534	alpha_lift(igcm_dust_mass)=1.e-6 !1.e-6 !Lifted mass coeff
	535
	536	r0_lift = reff_lift/ref_r0
	537	alpha_devil(igcm_dust_number)=r3n_q*
	538	& alpha_devil(igcm_dust_mass)/r0_lift**3
	539	alpha_lift(igcm_dust_number)=r3n_q*
	540	& alpha_lift(igcm_dust_mass)/r0_lift**3
	541
	542	radius(igcm_dust_mass) = reff_lift
	543	radius(igcm_dust_number) = reff_lift
	544
	545	write(,) "initracer: doubleq_param reff_lift:", reff_lift
	546	write(,) "initracer: doubleq_param nueff_lift:", nueff_lift
	547	write(,) "initracer: doubleq_param alpha_lift:",
	548	& alpha_lift(igcm_dust_mass)
	549	else
	550
	551	if (dustbin.gt.1) then
	552	print*,'initracer: STOP!',
	553	$ ' properties of dust need to be set in initracer !!!'
	554	stop
	555
	556	else if (dustbin.eq.1) then
	557
	558	c This will be used for 1 dust particle size:
	559	c ------------------------------------------
	560	radius(igcm_dustbin(1))=3.e-6
	561	alpha_lift(igcm_dustbin(1))=0.0e-6
	562	alpha_devil(igcm_dustbin(1))=7.65e-9
	563	rho_q(igcm_dustbin(1))=rho_dust
	564
	565	endif
	566	end if ! (doubleq)
	567
[358]	568
	569	c Scavenging of dust particles by H2O clouds:
	570	c ------------------------------------------
	571	c Initialize the two tracers used for the CCNs
	572	if (water.AND.doubleq.AND.scavenging) then
	573	radius(igcm_ccn_mass) = radius(igcm_dust_mass)
	574	alpha_lift(igcm_ccn_mass) = 1e-30
	575	alpha_devil(igcm_ccn_mass) = 1e-30
	576	rho_q(igcm_ccn_mass) = rho_dust
	577
	578	radius(igcm_ccn_number) = radius(igcm_ccn_mass)
	579	alpha_lift(igcm_ccn_number) = alpha_lift(igcm_ccn_mass)
	580	alpha_devil(igcm_ccn_number) = alpha_devil(igcm_ccn_mass)
	581	rho_q(igcm_ccn_number) = rho_q(igcm_ccn_mass)
	582	endif ! of if (water.AND.doubleq.AND.scavenging)
	583
[38]	584	c Submicron dust mode:
	585	c --------------------
	586
	587	if (submicron) then
	588	radius(igcm_dust_submicron)=0.1e-6
	589	rho_q(igcm_dust_submicron)=rho_dust
	590	if (doubleq) then
	591	c If doubleq is also active, we use the population ratio:
	592	alpha_lift(igcm_dust_submicron) =
	593	& alpha_lift(igcm_dust_number)popratio
	594	& rho_q(igcm_dust_submicron)4./3.pi*
	595	& radius(igcm_dust_submicron)**3.
	596	alpha_devil(igcm_dust_submicron)=1.e-30
	597	else
	598	alpha_lift(igcm_dust_submicron)=1e-6
	599	alpha_devil(igcm_dust_submicron)=1.e-30
	600	endif ! (doubleq)
	601	end if ! (submicron)
	602
	603	c Initialization for photochemistry:
	604	c ---------------------------------
	605	if (photochem) then
	606	! initialize chemistry+water (water will be correctly initialized below)
	607	! by initializing everything which is not dust ...
	608	do iq=1,nqmx
	609	txt=noms(iq)
	610	if (txt(1:4).ne."dust") then
	611	radius(iq)=0.
	612	alpha_lift(iq) =0.
	613	alpha_devil(iq)=0.
	614	endif
	615	enddo ! do iq=1,nqmx
	616	endif
	617
	618	c Initialization for water vapor
	619	c ------------------------------
	620	if(water) then
	621	radius(igcm_h2o_vap)=0.
	622	alpha_lift(igcm_h2o_vap) =0.
	623	alpha_devil(igcm_h2o_vap)=0.
	624	if(water.and.(nqmx.ge.2)) then
	625	radius(igcm_h2o_ice)=3.e-6
	626	rho_q(igcm_h2o_ice)=rho_ice
	627	alpha_lift(igcm_h2o_ice) =0.
	628	alpha_devil(igcm_h2o_ice)=0.
	629	elseif(water.and.(nqmx.lt.2)) then
	630	write(,) 'nqmx is too low : nqmx=', nqmx
	631	write(,) 'water= ',water
	632	endif
	633
	634	end if ! (water)
	635
	636	c Output for records:
	637	c ~~~~~~~~~~~~~~~~~~
	638	write(,)
	639	Write(,) '******** initracer : dust transport parameters :'
	640	write(,) 'alpha_lift = ', alpha_lift
	641	write(,) 'alpha_devil = ', alpha_devil
	642	write(,) 'radius = ', radius
	643	if(doubleq) then
	644	write(,) 'reff_lift (um) = ', reff_lift
	645	write(,) 'size distribution variance = ', varian
	646	write(,) 'r3n_q , ref_r0 : ', r3n_q , ref_r0
	647	end if
	648
	649	!
	650	! some extra (possibly redundant) sanity checks for tracers:
	651	! ---------------------------------------------------------
	652
	653	if (doubleq) then
	654	! verify that we indeed have dust_mass and dust_number tracers
	655	if (igcm_dust_mass.eq.0) then
	656	write(,) "initracer: error !!"
	657	write(,) " cannot use doubleq option without ",
	658	& "a dust_mass tracer !"
	659	stop
	660	endif
	661	if (igcm_dust_number.eq.0) then
	662	write(,) "initracer: error !!"
	663	write(,) " cannot use doubleq option without ",
	664	& "a dust_number tracer !"
	665	stop
	666	endif
	667	endif
	668
	669	if ((.not.doubleq).and.(dustbin.gt.0)) then
	670	! verify that we indeed have 'dustbin' dust tracers
	671	count=0
	672	do iq=1,dustbin
	673	if (igcm_dustbin(iq).ne.0) then
	674	count=count+1
	675	endif
	676	enddo
	677	if (count.ne.dustbin) then
	678	write(,) "initracer: error !!"
	679	write(,) " dusbin is set to ",dustbin,
	680	& " but we only have the following dust tracers:"
	681	do iq=1,count
	682	write(,)" ",trim(noms(igcm_dustbin(iq)))
	683	enddo
	684	stop
	685	endif
	686	endif
	687
	688	if (water) then
	689	! verify that we indeed have h2o_vap and h2o_ice tracers
	690	if (igcm_h2o_vap.eq.0) then
	691	write(,) "initracer: error !!"
	692	write(,) " cannot use water option without ",
	693	& "an h2o_vap tracer !"
	694	stop
	695	endif
	696	if (igcm_h2o_ice.eq.0) then
	697	write(,) "initracer: error !!"
	698	write(,) " cannot use water option without ",
	699	& "an h2o_ice tracer !"
	700	stop
	701	endif
	702	endif
	703
[358]	704	if (scavenging) then
	705	! verify that we indeed have ccn_mass and ccn_number tracers
	706	if (igcm_ccn_mass.eq.0) then
	707	write(,) "initracer: error !!"
	708	write(,) " cannot use scavenging option without ",
	709	& "a ccn_mass tracer !"
	710	stop
	711	endif
	712	if (igcm_ccn_number.eq.0) then
	713	write(,) "initracer: error !!"
	714	write(,) " cannot use scavenging option without ",
	715	& "a ccn_number tracer !"
	716	stop
	717	endif
	718	endif ! of if (scavenging)
	719
[38]	720	if (photochem .or. callthermos) then
	721	! verify that we indeed have the chemistry tracers
	722	if (igcm_co2.eq.0) then
	723	write(,) "initracer: error !!"
	724	write(,) " cannot use chemistry option without ",
	725	& "a co2 tracer !"
	726	stop
	727	endif
	728	if (igcm_co.eq.0) then
	729	write(,) "initracer: error !!"
	730	write(,) " cannot use chemistry option without ",
	731	& "a co tracer !"
	732	stop
	733	endif
	734	if (igcm_o.eq.0) then
	735	write(,) "initracer: error !!"
	736	write(,) " cannot use chemistry option without ",
	737	& "a o tracer !"
	738	stop
	739	endif
	740	if (igcm_o1d.eq.0) then
	741	write(,) "initracer: error !!"
	742	write(,) " cannot use chemistry option without ",
	743	& "a o1d tracer !"
	744	stop
	745	endif
	746	if (igcm_o2.eq.0) then
	747	write(,) "initracer: error !!"
	748	write(,) " cannot use chemistry option without ",
	749	& "an o2 tracer !"
	750	stop
	751	endif
	752	if (igcm_o3.eq.0) then
	753	write(,) "initracer: error !!"
	754	write(,) " cannot use chemistry option without ",
	755	& "an o3 tracer !"
	756	stop
	757	endif
	758	if (igcm_h.eq.0) then
	759	write(,) "initracer: error !!"
	760	write(,) " cannot use chemistry option without ",
	761	& "a h tracer !"
	762	stop
	763	endif
	764	if (igcm_h2.eq.0) then
	765	write(,) "initracer: error !!"
	766	write(,) " cannot use chemistry option without ",
	767	& "a h2 tracer !"
	768	stop
	769	endif
	770	if (igcm_oh.eq.0) then
	771	write(,) "initracer: error !!"
	772	write(,) " cannot use chemistry option without ",
	773	& "an oh tracer !"
	774	stop
	775	endif
	776	if (igcm_ho2.eq.0) then
	777	write(,) "initracer: error !!"
	778	write(,) " cannot use chemistry option without ",
	779	& "a ho2 tracer !"
	780	stop
	781	endif
	782	if (igcm_h2o2.eq.0) then
	783	write(,) "initracer: error !!"
	784	write(,) " cannot use chemistry option without ",
	785	& "a h2o2 tracer !"
	786	stop
	787	endif
	788	if (igcm_n2.eq.0) then
	789	write(,) "initracer: error !!"
	790	write(,) " cannot use chemistry option without ",
	791	& "a n2 tracer !"
	792	stop
	793	endif
	794	if (igcm_ar.eq.0) then
	795	write(,) "initracer: error !!"
	796	write(,) " cannot use chemistry option without ",
	797	& "an ar tracer !"
	798	stop
	799	endif
	800	endif ! of if (photochem .or. callthermos)
	801
	802	end

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

Download in other formats:

Original Format