Context Navigation

initracer.F @ 2090

Last change on this file since 2090 was 1974, checked in by mvals, 6 years ago
Mars GCM: Integration of the detached dust layer parametrizations (rocket dust storm, slope wind lifting, CW, and dust injection scheme, DB). Still experimental, default behaviour (rdstorm=.false., dustinjection=0) identical to previous revision. NB: Updated newstart requires an updated "surface.nc" containing the "hmons" field. EM+MV
File size: 25.7 KB

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

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

Download in other formats:

Original Format