Context Navigation

initracer.F90 @ 3557

Last change on this file since 3557 was 3557, checked in by debatzbr, 5 days ago
Miscellaneous cleans + Set-up the physics for the implementation of the microphysical model.
File size: 23.7 KB

Rev	Line
[3184]	1	SUBROUTINE initracer(ngrid,nq)
	2
	3	use surfdat_h, ONLY: dryness
	4	USE tracer_h
[3557]	5	USE callkeys_mod, only: aerohaze,nb_monomer,haze,fractal,fasthaze,rad_haze,callmufi
[3184]	6	USE recombin_corrk_mod, ONLY: ini_recombin
	7	USE mod_phys_lmdz_para, only: is_master, bcast
	8	use generic_cloud_common_h
[3195]	9	use aerosol_mod, only: iaero_haze,i_haze
[3184]	10	IMPLICIT NONE
	11	!=======================================================================
	12	! subject:
	13	! --------
[3195]	14	! Initialization related to tracer
[3184]	15	! (transported dust, water, chemical species, ice...)
	16	!
	17	! Name of the tracer
	18	!
	19	! Test of dimension :
	20	! Initialize COMMON tracer in tracer.h, using tracer names provided
	21	! by the argument nametrac
	22	!
	23	! author: F.Forget
	24	! ------
[3557]	25	! Ehouarn Millour (oct. 2008): identify tracers by their names
	26	! Y. Jaziri & J. Vatant d'Ollone (2020) : modern traceur.def
	27	! B. de Batz de Trenquelléon (2024): specific microphysical tracers
[3184]	28	!=======================================================================
	29
	30	integer,intent(in) :: ngrid,nq
	31
	32	character(len=500) :: tracline ! to read traceur.def lines
	33	integer :: blank !to store the index of 1st blank when reading tracers names
	34	integer iq,ig,count,ierr
[3195]	35	real r0_lift , reff_lift, rho_haze
	36	integer nqhaze(nq) ! to store haze tracers
[3557]	37	integer i, ia, block, j
[3195]	38	character(len=20) :: txt ! to store some text
[3557]	39	character(LEN=20) :: tracername ! to temporarily store text
	40	character(LEN=20) :: str
[3184]	41
	42	!-----------------------------------------------------------------------
	43	! radius(nq) ! aerosol particle radius (m)
	44	! rho_q(nq) ! tracer densities (kg.m-3)
	45	! qext(nq) ! Single Scat. Extinction coeff at 0.67 um
	46	! alpha_lift(nq) ! saltation vertical flux/horiz flux ratio (m-1)
	47	! alpha_devil(nq) ! lifting coeeficient by dust devil
	48	! rho_dust ! Mars dust density
	49	! rho_ice ! Water ice density
	50	! doubleq ! if method with mass (iq=1) and number(iq=2) mixing ratio
	51	! varian ! Characteristic variance of log-normal distribution
	52	!-----------------------------------------------------------------------
	53
	54	if (is_master) then ! only the master proc/thread needs do this
	55
	56	moderntracdef=.false. ! For modern traceur.def (default false, old type)
	57
	58	open(407, form = 'formatted', status = 'old', &
	59	file = 'traceur.def', iostat=ierr)
	60	if (ierr /=0) then
	61	! call abort_physic('initracer',&
	62	! 'Problem in opening traceur.def',1)
	63	print*,'Problem in opening traceur.def'
	64	stop
	65	end if
	66	!! - Modif. by JVO and YJ for modern planetary traceur.def ---------------
	67	READ(407,'(A)') tracline
	68	IF (trim(tracline).ne.'#ModernTrac-v1') THEN ! Test modern traceur.def
	69	READ(tracline,*) nqtot ! Try standard traceur.def
	70	ELSE
	71	moderntracdef = .true.
	72	DO
	73	READ(407,'(A)',iostat=ierr) tracline
	74	IF (ierr==0) THEN
	75	IF (index(tracline,'#').ne.1) THEN ! Allows arbitary number of comments lines in the header
	76	READ(tracline,*) nqtot
	77	! Temporary not implemented solution
	78	if (nqtot/=nq) then
	79	! call abort_physic('initracer','Different number of '// &
	80	! 'tracers in dynamics and physics not managed yet',1)
	81	print*,'!= nbr oftracers in dynamics and physics not managed yet'
	82	stop
	83	endif
	84	EXIT
	85	ENDIF
	86	ELSE ! If pb, or if reached EOF without having found nqtot
	87	! call abort_physic('initracer','Unable to read numbers '// &
	88	! 'of tracers in traceur.def',1)
	89	print*,"unable to read numbers of tracer in tracer.def"
	90	stop
	91	ENDIF
	92	ENDDO
	93	ENDIF ! if modern or standard traceur.def
[3195]	94
[3184]	95	endif ! of if (is_master)
[3195]	96
[3184]	97	! share the information with other procs/threads (if any)
	98	CALL bcast(nqtot)
	99	CALL bcast(moderntracdef)
[3195]	100
[3184]	101	!! -----------------------------------------------------------------------
	102	!! For the moment number of tracers in dynamics and physics are equal
	103	nqtot=nq
	104	!! we allocate once for all arrays in common in tracer_h.F90
	105	!! (supposedly those are not used before call to initracer)
	106	IF (.NOT.ALLOCATED(noms)) ALLOCATE(noms(nq))
	107	IF (.NOT.ALLOCATED(mmol)) ALLOCATE(mmol(nq))
	108	IF (.NOT.ALLOCATED(aki)) ALLOCATE(aki(nqtot))
	109	IF (.NOT.ALLOCATED(cpi)) ALLOCATE(cpi(nqtot))
	110	IF (.NOT.ALLOCATED(radius)) ALLOCATE(radius(nq))
	111	IF (.NOT.ALLOCATED(rho_q)) ALLOCATE(rho_q(nq))
	112	IF (.NOT.ALLOCATED(qext)) ALLOCATE(qext(nq))
	113	IF (.NOT.ALLOCATED(alpha_lift)) ALLOCATE(alpha_lift(nq))
	114	IF (.NOT.ALLOCATED(alpha_devil)) ALLOCATE(alpha_devil(nq))
	115	IF (.NOT.ALLOCATED(qextrhor)) ALLOCATE(qextrhor(nq))
[3195]	116	! IF (.NOT.ALLOCATED(igcm_dustbin)) ALLOCATE(igcm_dustbin(nq))
[3184]	117	IF (.NOT.ALLOCATED(is_chim)) ALLOCATE(is_chim(nqtot))
	118	IF (.NOT.ALLOCATED(is_rad)) ALLOCATE(is_rad(nqtot))
	119	IF (.NOT.ALLOCATED(is_recomb)) ALLOCATE(is_recomb(nqtot))
	120	IF (.NOT.ALLOCATED(is_recomb_qset)) THEN
	121	ALLOCATE(is_recomb_qset(nqtot))
	122	ENDIF
	123	IF (.NOT.ALLOCATED(is_recomb_qotf)) THEN
	124	ALLOCATE(is_recomb_qotf(nqtot))
	125	ENDIF
	126	IF (.NOT. allocated(is_condensable)) allocate(is_condensable(nq)) !LT
	127	IF (.NOT. allocated(is_rgcs)) allocate(is_rgcs(nq)) !LT
	128	IF (.NOT. allocated(constants_mass)) allocate(constants_mass(nq))
[3275]	129	IF (.NOT. allocated(constants_delta_gasH)) allocate(constants_delta_gasH(nq))
[3184]	130	IF (.NOT. allocated(constants_Tref)) allocate(constants_Tref(nq))
	131	IF (.NOT. allocated(constants_Pref)) allocate(constants_Pref(nq))
	132	IF (.NOT. allocated(constants_epsi_generic)) allocate(constants_epsi_generic(nq))
	133	IF (.NOT. allocated(constants_RLVTT_generic)) allocate(constants_RLVTT_generic(nq))
	134	IF (.NOT. allocated(constants_metallicity_coeff)) allocate(constants_metallicity_coeff(nq))
	135	IF (.NOT. allocated(constants_RCPV_generic)) allocate(constants_RCPV_generic(nq))
	136
	137	!! initialization
	138	alpha_lift(:) = 0.
	139	alpha_devil(:) = 0.
	140	mmol(:) = 0.
	141	aki(:) = 0.
	142	cpi(:) = 0.
	143	is_chim(:) = 0
	144	is_rad(:) = 0
	145	is_recomb(:) = 0
	146	is_recomb_qset(:) = 0
	147	is_recomb_qotf(:) = 0
[3195]	148
[3184]	149	! Added by JVO 2017 : these arrays are handled later
	150	! -> initialization is the least we can do, please !!!
	151	radius(:)=0.
	152	qext(:)=0.
	153
	154	! For condensable tracers, by Lucas Teinturier and Noé Clément (2022)
	155
	156	is_condensable(:)= 0
	157	is_rgcs(:) = 0
	158	constants_mass(:)=0
[3275]	159	constants_delta_gasH(:)=0
[3184]	160	constants_Tref(:)=0
	161	constants_Pref(:)=0
	162	constants_epsi_generic(:)=0
	163	constants_RLVTT_generic(:)=0
	164	constants_metallicity_coeff(:)=0
	165	constants_RCPV_generic(:)=0
	166
	167	rho_q(:) = 0. !need to be init here if we want to read it from modern traceur with get_tracdat
	168
	169
	170	! Initialization: Read tracers names from traceur.def
[3195]	171	do iq=1,nq
[3184]	172	if (is_master) read(407,'(A)') tracline
	173	call bcast(tracline)
	174	blank = index(tracline,' ') ! Find position of 1st blank in tracline
	175	noms(iq) = tracline(1:blank) !ensure that in Modern-trac case, noms is actually the name and not all properties
	176	enddo
	177
	178	! Identify tracers by their names: (and set corresponding values of mmol)
	179	! 0. initialize tracer indexes to zero:
	180	! NB: igcm_* indexes are commons in 'tracer.h'
	181	igcm_n2=0
[3195]	182	igcm_ch4_gas=0
	183	igcm_ch4_ice=0
	184	igcm_prec_haze=0
	185	igcm_co_gas=0
	186	igcm_co_ice=0
[3184]	187
[3195]	188	nqhaze(:)=0
	189	i=0
	190	DO iq=1,nq
	191	txt=noms(iq)
	192	IF (txt(1:4).eq."haze") THEN
	193	i=i+1
	194	nqhaze(i)=iq
	195	ENDIF
	196	ENDDO
	197	if ((haze.or.fasthaze).and.i==0) then
	198	print*, 'Haze active but no haze tracer in traceur.def'
	199	stop
	200	endif
[3184]	201
[3195]	202	igcm_haze=0
	203	igcm_haze10=0
	204	igcm_haze30=0
	205	igcm_haze50=0
	206	igcm_haze100=0
	207
	208	! Eddy diffusion tracers
	209	igcm_eddy1e6=0
	210	igcm_eddy1e7=0
	211	igcm_eddy5e7=0
	212	igcm_eddy1e8=0
	213	igcm_eddy5e8=0
	214	write(,) 'initracer: noms() ', noms
	215	rho_n2=1000 ! n2 ice
	216	rho_ch4_ice=520. ! rho ch4 ice
	217	rho_co_ice=520. ! rho ch4 ice
	218	rho_haze=800. ! haze
	219
[3184]	220	! 1. find dust tracers
	221	count=0
	222
	223	! 2. find chemistry and water tracers
	224	do iq=1,nq
	225	if (noms(iq).eq."n2") then
	226	igcm_n2=iq
	227	mmol(igcm_n2)=28.
	228	count=count+1
[3195]	229	write(,) 'Tracer ',count,' = n2'
[3184]	230	endif
[3275]	231	if (noms(iq).eq."ch4_gas") then
[3195]	232	igcm_ch4_gas=iq
	233	mmol(igcm_ch4_gas)=16.
[3184]	234	count=count+1
[3195]	235	write(,) 'Tracer ',count,' = ch4 gas'
[3184]	236	endif
[3195]	237	if (noms(iq).eq."ch4_ice") then
	238	igcm_ch4_ice=iq
	239	mmol(igcm_ch4_ice)=16.
	240	radius(iq)=3.e-6
	241	rho_q(iq)=rho_ch4_ice
[3184]	242	count=count+1
[3195]	243	write(,) 'Tracer ',count,' = ch4 ice'
[3184]	244	endif
[3275]	245	if (noms(iq).eq."co_gas") then
[3195]	246	igcm_co_gas=iq
	247	mmol(igcm_co_gas)=28.
[3184]	248	count=count+1
[3195]	249	write(,) 'Tracer ',count,' = co gas'
[3184]	250	endif
[3195]	251	if (noms(iq).eq."co_ice") then
	252	igcm_co_ice=iq
	253	mmol(igcm_co_ice)=28.
	254	radius(iq)=3.e-6
	255	rho_q(iq)=rho_co_ice
[3184]	256	count=count+1
[3195]	257	write(,) 'Tracer ',count,' = co ice'
[3184]	258	endif
[3195]	259	if (noms(iq).eq."prec_haze") then
	260	igcm_prec_haze=iq
[3184]	261	count=count+1
[3195]	262	write(,) 'Tracer ',count,' = prec haze'
[3184]	263	endif
[3195]	264	if (noms(iq).eq."haze") then
	265	igcm_haze=iq
[3184]	266	count=count+1
[3195]	267	radius(iq)=rad_haze
	268	rho_q(iq)=rho_haze
	269	write(,) 'Tracer ',count,' = haze'
[3184]	270	endif
[3195]	271	if (noms(iq).eq."haze10") then
	272	igcm_haze10=iq
[3184]	273	count=count+1
[3195]	274	radius(iq)=10.e-9
	275	rho_q(iq)=rho_haze
	276	write(,) 'Tracer ',count,' = haze10'
[3184]	277	endif
[3195]	278	if (noms(iq).eq."haze30") then
	279	igcm_haze30=iq
[3184]	280	count=count+1
[3195]	281	radius(iq)=30.e-9
	282	rho_q(iq)=rho_haze
	283	write(,) 'Tracer ',count,' = haze30'
[3184]	284	endif
[3195]	285	if (noms(iq).eq."haze50") then
	286	igcm_haze50=iq
[3184]	287	count=count+1
[3195]	288	radius(iq)=50.e-9
	289	rho_q(iq)=rho_haze
	290	write(,) 'Tracer ',count,' = haze50'
[3184]	291	endif
[3195]	292	if (noms(iq).eq."haze100") then
	293	igcm_haze100=iq
[3184]	294	count=count+1
[3195]	295	radius(iq)=100.e-9
	296	rho_q(iq)=rho_haze
	297	write(,) 'Tracer ',count,' = haze100'
[3184]	298	endif
[3195]	299	! Eddy diffusion tracers
	300	if (noms(iq).eq."eddy1e6") then
	301	igcm_eddy1e6=iq
[3184]	302	count=count+1
[3195]	303	write(,) 'Tracer ',count,' = eddy1e6'
[3184]	304	endif
[3195]	305	if (noms(iq).eq."eddy1e7") then
	306	igcm_eddy1e7=iq
[3184]	307	count=count+1
[3195]	308	write(,) 'Tracer ',count,' = eddy1e7'
[3184]	309	endif
[3195]	310	if (noms(iq).eq."eddy5e7") then
	311	igcm_eddy5e7=iq
[3184]	312	count=count+1
[3195]	313	write(,) 'Tracer ',count,' = eddy5e7'
[3184]	314	endif
[3195]	315	if (noms(iq).eq."eddy1e8") then
	316	igcm_eddy1e8=iq
[3184]	317	count=count+1
[3195]	318	write(,) 'Tracer ',count,' = eddy1e8'
[3184]	319	endif
[3195]	320	if (noms(iq).eq."eddy5e8") then
	321	igcm_eddy5e8=iq
[3184]	322	count=count+1
[3195]	323	write(,) 'Tracer ',count,' = eddy5e8'
[3184]	324	endif
	325	enddo ! of do iq=1,nq
	326
[3241]	327	! ! 3. find condensable traceurs different from h2o and n2
	328	! do iq=1,nq
	329	! if ((index(noms(iq),"vap") .ne. 0) .and. (index(noms(iq),"n2") .eq. 0)) then
	330	! count=count+1
	331	! endif
	332	! if ((index(noms(iq),"ice") .ne. 0) .and. (index(noms(iq),"n2") .eq. 0)) then
	333	! count=count+1
	334	! endif
[3184]	335
[3241]	336	! enddo ! of do iq=1,nq
[3195]	337
[3184]	338	! check that we identified all tracers:
	339	if (count.ne.nq) then
	340	write(,) "initracer: found only ",count," tracers"
	341	write(,) " expected ",nq
	342	do iq=1,count
	343	write(,)' ',iq,' ',trim(noms(iq))
	344	enddo
	345	! stop
	346	else
	347	write(,) "initracer: found all expected tracers, namely:"
	348	do iq=1,nq
	349	write(,)' ',iq,' ',trim(noms(iq))
	350	enddo
	351	endif
[3557]	352
	353	! Compute number of microphysics tracers:
	354	! By convention they all have the prefix "mu_" (case sensitive !)
	355	nmicro = 0
	356	IF (callmufi) THEN
	357	DO iq=1,nq
	358	str = noms(iq)
	359	IF (str(1:3) == "mu_") THEN
	360	nmicro = nmicro+1
	361	count = count+1
	362	ENDIF
	363	ENDDO
[3184]	364
[3557]	365	! Checking the expected number of tracers:
	366	! Microphysics moment model: nmicro = 4
	367	IF (nmicro < 4) THEN
	368	WRITE(,) "initracer:error:"," Inconsistent number of microphysical tracers"
	369	WRITE(,) "expected at least 4 tracers,", nmicro, " given"
	370	CALL abort
	371	ELSE IF (nmicro > 4) THEN
	372	WRITE(,) "!!! WARNING !!! initracer: I was expecting only four tracers, you gave me more."
	373	CALL abort
	374	ENDIF
	375
	376	! microphysics indexes share the same values than original tracname.
	377	IF (.NOT.ALLOCATED(micro_indx)) ALLOCATE(micro_indx(nmicro))
	378	j = 1
	379	DO i=1,nq
	380	str = noms(i)
	381	IF (str(1:3) == "mu_") THEN
	382	micro_indx(j) = i
	383	j=j+1
	384	ENDIF
	385	ENDDO
	386
	387	ELSE
	388	IF (.NOT.ALLOCATED(micro_indx)) ALLOCATE(micro_indx(nmicro))
	389
	390	ENDIF ! end of callmufi
	391
[3184]	392	! Get data of tracers. Need to rewind traceur.def first
	393	if (is_master) then
	394	rewind(407)
[3195]	395	do
[3184]	396	read(407,'(A)') tracline
[3195]	397	if (index(tracline,"#") .ne. 1) then
[3184]	398	exit
	399	endif
	400	enddo
	401	endif
	402	do iq=1,nqtot
	403	if (is_master) read(407,'(A)') tracline
	404	call bcast(tracline)
	405	call get_tracdat(iq, tracline)
	406	enddo
	407
	408	if (is_master) close(407)
	409
	410	! Get specific data of condensable tracers
	411	do iq=1,nq
	412	if((is_condensable(iq)==1)) then
	413	write(,) "There is a specie which is condensable, for generic condensation : ", noms(iq)
	414	write(,) 'looking specie parameters for : ',noms(iq)(1:len(trim(noms(iq)))-4)
[3405]	415	! call specie_parameters_table(noms(iq)(1:len(trim(noms(iq)))-4))
	416	! constants_mass(iq)=m
	417	! constants_delta_gasH(iq)=delta_gasH
	418	! constants_Tref(iq)=Tref
	419	! constants_Pref(iq)=Pref
	420	! constants_epsi_generic(iq)=epsi_generic
	421	! constants_RLVTT_generic(iq)=RLVTT_generic
	422	! constants_metallicity_coeff(iq)=metallicity_coeff
	423	! constants_RCPV_generic(iq)=RCPV_generic
[3184]	424	else
	425	write(,) "This tracer is not condensable, for generic condensation : : ", noms(iq)
	426	write(,) "We keep condensable constants at zero"
	427	endif !(is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0))
	428	enddo ! iq=1,nq
	429
	430	! Calculate number of species in the chemistry
	431	nesp = sum(is_chim)
	432	write(,) 'Number of species in the chemistry nesp = ',nesp
	433
	434	! Calculate number of generic tracers
	435	ngt = sum(is_condensable)
	436	write(,) 'Number of generic tracer is ngt = ',ngt
	437
	438	! Calculate number of radiative generic condensable species
	439	n_rgcs = sum(is_rgcs)
	440	write(,)'Number of Radiative Generic Condensable Species is n_rgcs = ',n_rgcs
[3195]	441	if (n_rgcs> ngt/2) then
[3184]	442	write(,) 'You have more Radiative Generic Condensable Species than Generic Condensable Species'
	443	write(,)'This is not possible: check your Modern traceur.def'
	444	call abort_physic("initracer, issue with # of RGCS and GCS")
	445	endif
	446
[3557]	447	! Calculate number of microphysical tracer
	448	write(,) 'Number of microphysical tracer nmicro = ',nmicro
	449	IF (callmufi) THEN
	450	call dumptracers(micro_indx)
	451	ENDIF
	452
[3184]	453	! Processing modern traceur options
	454	if(moderntracdef) then
	455	call ini_recombin
	456	endif
	457
	458	!------------------------------------------------------------
	459	! Initialisation tracers ....
	460	!------------------------------------------------------------
	461	! rho_q(1:nq)=0
	462
	463	rho_n2=1000. ! N2 ice density (kg.m-3)
	464
[3193]	465	lw_co=274000.
	466	lw_ch4=586700.
	467	lw_n2=2.5e5
[3237]	468	write(,) "lw_n2 = ", lw_n2
[3195]	469
[3193]	470	if (haze) then
	471	! the sedimentation radius remains radius(igcm_haze)
	472	if (fractal) then
	473	nmono=nb_monomer
	474	else
	475	nmono=1
[3557]	476	endif ! end fractal
[3193]	477
	478	ia=0
	479	if (aerohaze) then
	480	ia=ia+1
	481	iaero_haze=ia
	482	write(,) '--- number of haze aerosol = ', iaero_haze
	483
	484	block=0 ! Only one type of haze is active : the first one set in traceur.def
[3195]	485	do iq=1,nq
[3193]	486	tracername=noms(iq)
[3195]	487	write(,) "--> tracername ",iq,'/',nq,' = ',tracername
[3193]	488	if (tracername(1:4).eq."haze".and.block.eq.0) then
	489	i_haze=iq
	490	block=1
	491	write(,) "i_haze=",i_haze
[3557]	492	write(,) "Careful: if you set many haze traceurs in &
	493	traceur.def,only ",tracername," will be radiatively active &
	494	(first one in traceur.def)"
[3193]	495	endif
	496	enddo
[3557]	497	endif ! end aerohaze
	498	endif ! end haze
[3193]	499
[3184]	500	! Output for records:
	501	! ~~~~~~~~~~~~~~~~~~
	502	write(,)
	503	Write(,) '******** initracer : dust transport parameters :'
	504	write(,) 'alpha_lift = ', alpha_lift
	505	write(,) 'alpha_devil = ', alpha_devil
	506	write(,) 'radius = ', radius
[3195]	507	write(,) 'Qext = ', qext
[3184]	508	write(,)
	509
	510	contains
	511
	512	subroutine get_tracdat(iq,tracline)
	513	!-------------------ADDING NEW OPTIONS-------------------
	514	! Duplicate if sentence for adding new options
	515	! Do not forget to add the new variables in tracer_h.F90
	516	! Do not forget to allocate and initialize the new variables above
	517	! Please update list of options in "LMDZ.GENERIC/deftank/traceur.def.modern"
	518	!-------------------------------------------------------
	519	use tracer_h
	520	implicit none
	521	integer, intent(in) :: iq ! tracer index
	522	character(len=500),intent(in) :: tracline ! traceur.def lines with parameters
[3195]	523
[3184]	524	read(tracline,*) noms(iq)
	525	! JVO 20 : We should add a sanity check aborting when duplicates in names !
	526	write(,)"initracer: iq=",iq,"noms(iq)=",trim(noms(iq))
	527	! option mmol
	528	if (index(tracline,'mmol=' ) /= 0) then
	529	read(tracline(index(tracline,'mmol=')+len('mmol='):),*)&
	530	mmol(iq)
	531	write(,) ' Parameter value (traceur.def) : mmol=', &
	532	mmol(iq)
	533	else
	534	write(,) ' Parameter value (default) : mmol=', &
	535	mmol(iq)
	536	end if
[3195]	537	! option aki
[3184]	538	if (index(tracline,'aki=' ) /= 0) then
	539	read(tracline(index(tracline,'aki=')+len('aki='):),*) &
	540	aki(iq)
	541	write(,) ' Parameter value (traceur.def) : aki=', &
	542	aki(iq)
	543	else
	544	write(,) ' Parameter value (default) : aki=', &
	545	aki(iq)
	546	end if
[3195]	547	! option cpi
[3184]	548	if (index(tracline,'cpi=' ) /= 0) then
	549	read(tracline(index(tracline,'cpi=')+len('cpi='):),*) &
	550	cpi(iq)
	551	write(,) ' Parameter value (traceur.def) : cpi=', &
	552	cpi(iq)
	553	else
	554	write(,) ' Parameter value (default) : cpi=', &
	555	cpi(iq)
	556	end if
[3195]	557	! option is_chim
[3184]	558	if (index(tracline,'is_chim=' ) /= 0) then
	559	read(tracline(index(tracline,'is_chim=')+len('is_chim='):),*) &
	560	is_chim(iq)
	561	write(,) ' Parameter value (traceur.def) : is_chim=', &
	562	is_chim(iq)
	563	else
	564	write(,) ' Parameter value (default) : is_chim=', &
	565	is_chim(iq)
	566	end if
[3195]	567	! option is_rad
[3184]	568	if (index(tracline,'is_rad=') /= 0) then
	569	read(tracline(index(tracline,'is_rad=') &
	570	+len('is_rad='):),*) is_rad(iq)
	571	write(,) ' Parameter value (traceur.def) : is_rad=', &
	572	is_rad(iq)
	573	else
	574	write(,) ' Parameter value (default) : is_rad=', &
	575	is_rad(iq)
	576	end if
	577	! option is_recomb
	578	if (index(tracline,'is_recomb=') /= 0) then
	579	read(tracline(index(tracline,'is_recomb=') &
	580	+len('is_recomb='):),*) is_recomb(iq)
	581	write(,) ' Parameter value (traceur.def) : is_recomb=', &
	582	is_recomb(iq)
	583	else
	584	write(,) ' Parameter value (default) : is_recomb=', &
	585	is_recomb(iq)
	586	end if
	587	! option is_recomb_qset
	588	if (index(tracline,'is_recomb_qset=') /= 0) then
	589	read(tracline(index(tracline,'is_recomb_qset=') &
	590	+len('is_recomb_qset='):),*) is_recomb_qset(iq)
	591	write(,) ' Parameter value (traceur.def) :'// &
	592	' is_recomb_qset=', &
	593	is_recomb_qset(iq)
	594	else
	595	write(,) ' Parameter value (default) :'// &
	596	' is_recomb_qset=', &
	597	is_recomb_qset(iq)
	598	endif
	599	! option is_recomb_qotf
	600	if (index(tracline,'is_recomb_qotf=') /= 0) then
	601	read(tracline(index(tracline,'is_recomb_qotf=') &
	602	+len('is_recomb_qotf='):),*) is_recomb_qotf(iq)
	603	write(,) ' Parameter value (traceur.def) :'// &
	604	' is_recomb_qotf=', &
	605	is_recomb_qotf(iq)
	606	else
	607	write(,) ' Parameter value (default) :'// &
	608	' is_recomb_qotf=', &
	609	is_recomb_qotf(iq)
	610	end if
	611	!option is_condensable (LT)
[3195]	612	if (index(tracline,'is_condensable=') /=0) then
[3184]	613	read(tracline(index(tracline,'is_condensable=') &
	614	+len('is_condensable='):),*) is_condensable(iq)
	615	write(,) ' Parameter value (traceur.def) :'// &
[3195]	616	' is_condensable=', is_condensable(iq)
[3184]	617	else
	618	write(,) ' Parameter value (default) :'// &
[3195]	619	' is_condensable=', is_condensable(iq)
[3184]	620	endif
	621	!option radius
[3195]	622	if (index(tracline,'radius=') .ne. 0) then
[3184]	623	read(tracline(index(tracline,'radius=') &
	624	+len('radius='):),*) radius(iq)
	625	write(,)'Parameter value (traceur.def) :'// &
	626	"radius=",radius(iq), "m"
	627	else
	628	write(,) ' Parameter value (default) :'// &
[3195]	629	' radius=', radius(iq)," m"
[3184]	630	endif
	631	!option rho
[3195]	632	if (index(tracline,'rho=') .ne. 0) then
[3184]	633	read(tracline(index(tracline,'rho=') &
	634	+len('rho='):),*) rho_q(iq)
	635	write(,)'Parameter value (traceur.def) :'// &
	636	"rho=",rho_q(iq)
	637	else
	638	write(,) ' Parameter value (default) :'// &
[3195]	639	' rho=', rho_q(iq)
[3184]	640	endif
[3195]	641	!option is_rgcs
	642	if (index(tracline,'is_rgcs') .ne. 0) then
[3184]	643	read(tracline(index(tracline,'is_rgcs=') &
	644	+len('is_rgcs='):),*) is_rgcs(iq)
	645	write(,)'Parameter value (traceur.def) :'// &
	646	'is_rgcs=',is_rgcs(iq)
	647	else
	648	write(,)'Parameter value (default) : '// &
	649	'is_rgcs = ',is_rgcs(iq)
	650	endif
	651	end subroutine get_tracdat
	652
	653	end
	654

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

Download in other formats:

Original Format