Context Navigation

readcommand.f90 @ 36

Last change on this file since 36 was 36, checked in by pesei, 9 years ago
Implement switch for incremental deposition, see ticket:113 and many small changes, see changelog.txt
File size: 23.3 KB

Rev	Line
[4]	1	!**********************************************************************
[36]	2	! Copyright 1998-2015 *
[4]	3	! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, *
	4	! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann *
	5	! *
	6	! This file is part of FLEXPART. *
	7	! *
	8	! FLEXPART is free software: you can redistribute it and/or modify *
	9	! it under the terms of the GNU General Public License as published by*
	10	! the Free Software Foundation, either version 3 of the License, or *
	11	! (at your option) any later version. *
	12	! *
	13	! FLEXPART is distributed in the hope that it will be useful, *
	14	! but WITHOUT ANY WARRANTY; without even the implied warranty of *
	15	! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
	16	! GNU General Public License for more details. *
	17	! *
	18	! You should have received a copy of the GNU General Public License *
	19	! along with FLEXPART. If not, see <http://www.gnu.org/licenses/>. *
	20	!**********************************************************************
	21
	22	subroutine readcommand
	23
	24	!*****************************************************************************
	25	! *
	26	! This routine reads the user specifications for the current model run. *
	27	! *
	28	! Author: A. Stohl *
	29	! *
	30	! 18 May 1996 *
[36]	31	!
	32	! HSO, 1 July 2014: Added optional namelist input
	33	! PS 2/2015: add ldep_incr as optional command input
	34	! make new parameters (last 3) optional for bwd compatibility
[4]	35	! *
	36	!*****************************************************************************
	37	! *
	38	! Variables: *
	39	! bdate beginning date as Julian date *
	40	! ctl factor by which time step must be smaller than *
	41	! Lagrangian time scale *
	42	! ibdate,ibtime beginnning date and time (YYYYMMDD, HHMISS) *
	43	! ideltas [s] modelling period *
	44	! iedate,ietime ending date and time (YYYYMMDD, HHMISS) *
	45	! ifine reduction factor for vertical wind time step *
	46	! outputforeachrel for forward runs it is possible either to create *
	47	! one outputfield or several for each releasepoint *
	48	! iflux switch to turn on (1)/off (0) flux calculations *
	49	! iout 1 for conc. (residence time for backward runs) output,*
	50	! 2 for mixing ratio output, 3 both, 4 for plume *
	51	! trajectory output, 5 = options 1 and 4 *
	52	! ipin 1 continue simulation with dumped particle data, 0 no *
	53	! ipout 0 no particle dump, 1 every output time, 3 only at end*
	54	! itsplit [s] time constant for particle splitting *
	55	! loutaver [s] concentration output is an average over loutaver *
	56	! seconds *
	57	! loutsample [s] average is computed from samples taken every [s] *
	58	! seconds *
	59	! loutstep [s] time interval of concentration output *
	60	! lsynctime [s] synchronisation time interval for all particles *
	61	! lagespectra switch to turn on (1)/off (0) calculation of age *
	62	! spectra *
	63	! lconvection value of either 0 and 1 indicating mixing by *
	64	! convection *
	65	! = 0 .. no convection *
	66	! + 1 .. parameterisation of mixing by subgrid-scale *
	67	! convection = on *
	68	! lsubgrid switch to turn on (1)/off (0) subgrid topography *
	69	! parameterization *
	70	! method method used to compute the particle pseudovelocities *
	71	! mdomainfill 1 use domain-filling option, 0 not, 2 use strat. O3 *
	72	! *
[36]	73	! ldep_incr .true. for incremental,
	74	! .false. for accumulated deposition output (DEFAULT)
[4]	75	! Constants: *
	76	! unitcommand unit connected to file COMMAND *
	77	! *
	78	!*****************************************************************************
	79
	80	use par_mod
	81	use com_mod
	82
	83	implicit none
	84
[36]	85	real (kind=dp) :: juldate
	86	character (len=50) :: line
[4]	87	logical :: old
[36]	88	integer :: ireaderror, icmdstat
[4]	89
[20]	90	namelist /command/ &
[24]	91	ldirect, &
	92	ibdate,ibtime, &
	93	iedate,ietime, &
	94	loutstep, &
	95	loutaver, &
	96	loutsample, &
	97	itsplit, &
	98	lsynctime, &
	99	ctl, &
	100	ifine, &
	101	iout, &
	102	ipout, &
	103	lsubgrid, &
	104	lconvection, &
	105	lagespectra, &
	106	ipin, &
	107	ioutputforeachrelease, &
	108	iflux, &
	109	mdomainfill, &
	110	ind_source, &
	111	ind_receptor, &
	112	mquasilag, &
	113	nested_output, &
	114	linit_cond, &
[30]	115	lnetcdfout, &
[36]	116	surf_only, &
	117	ldep_incr
[4]	118
[36]	119	! Set default values for namelist command
	120	ldirect = 0
	121	ibdate = 20000101
	122	ibtime = 0
	123	iedate = 20000102
	124	ietime = 0
	125	loutstep = 10800
	126	loutaver = 10800
	127	loutsample = 900
	128	itsplit = 999999999
	129	lsynctime = 900
	130	ctl = -5.0
	131	ifine = 4
	132	iout = 3
	133	ipout = 0
	134	lsubgrid = 1
	135	lconvection = 1
	136	lagespectra = 0
	137	ipin = 1
	138	ioutputforeachrelease = 1
	139	iflux = 1
	140	mdomainfill = 0
	141	ind_source = 1
	142	ind_receptor = 1
	143	mquasilag = 0
	144	nested_output = 0
	145	linit_cond = 0
	146	lnetcdfout = 0
	147	surf_only = 0
	148	ldep_incr = .false.
[20]	149
[4]	150	! Open the command file and read user options
[20]	151	! Namelist input first: try to read as namelist file
	152	!**************************************************************************
[27]	153	open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old',form='formatted',err=999)
[4]	154
[27]	155	! try namelist input (default)
[36]	156	read(unitcommand,command,iostat=ireaderror)
[27]	157	close(unitcommand)
[4]	158
[27]	159	! distinguish namelist from fixed text input
[36]	160	if ((ireaderror .ne. 0) .or. (ldirect .eq. 0)) then
	161
	162	! parse as text file format
[24]	163
[27]	164	open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old', err=999)
[24]	165
[20]	166	! Check the format of the COMMAND file (either in free format,
	167	! or using formatted mask)
	168	! Use of formatted mask is assumed if line 10 contains the word 'DIRECTION'
	169	!**************************************************************************
[4]	170
[20]	171	call skplin(9,unitcommand)
	172	read (unitcommand,901) line
	173	901 format (a)
	174	if (index(line,'LDIRECT') .eq. 0) then
	175	old = .false.
[27]	176	write(,) 'COMMAND in old short format, please update to namelist format'
[20]	177	else
	178	old = .true.
[27]	179	write(,) 'COMMAND in old long format, please update to namelist format'
[20]	180	endif
	181	rewind(unitcommand)
[4]	182
[24]	183
[20]	184	! Read parameters
	185	!****************
[4]	186
[20]	187	call skplin(7,unitcommand)
	188	if (old) call skplin(1,unitcommand)
	189	read(unitcommand,*) ldirect
	190	if (old) call skplin(3,unitcommand)
	191	read(unitcommand,*) ibdate,ibtime
	192	if (old) call skplin(3,unitcommand)
	193	read(unitcommand,*) iedate,ietime
	194	if (old) call skplin(3,unitcommand)
	195	read(unitcommand,*) loutstep
	196	if (old) call skplin(3,unitcommand)
	197	read(unitcommand,*) loutaver
	198	if (old) call skplin(3,unitcommand)
	199	read(unitcommand,*) loutsample
	200	if (old) call skplin(3,unitcommand)
	201	read(unitcommand,*) itsplit
	202	if (old) call skplin(3,unitcommand)
	203	read(unitcommand,*) lsynctime
	204	if (old) call skplin(3,unitcommand)
	205	read(unitcommand,*) ctl
	206	if (old) call skplin(3,unitcommand)
	207	read(unitcommand,*) ifine
	208	if (old) call skplin(3,unitcommand)
	209	read(unitcommand,*) iout
	210	if (old) call skplin(3,unitcommand)
	211	read(unitcommand,*) ipout
	212	if (old) call skplin(3,unitcommand)
	213	read(unitcommand,*) lsubgrid
	214	if (old) call skplin(3,unitcommand)
	215	read(unitcommand,*) lconvection
	216	if (old) call skplin(3,unitcommand)
	217	read(unitcommand,*) lagespectra
	218	if (old) call skplin(3,unitcommand)
	219	read(unitcommand,*) ipin
	220	if (old) call skplin(3,unitcommand)
	221	read(unitcommand,*) ioutputforeachrelease
	222	if (old) call skplin(3,unitcommand)
	223	read(unitcommand,*) iflux
	224	if (old) call skplin(3,unitcommand)
	225	read(unitcommand,*) mdomainfill
	226	if (old) call skplin(3,unitcommand)
	227	read(unitcommand,*) ind_source
	228	if (old) call skplin(3,unitcommand)
	229	read(unitcommand,*) ind_receptor
	230	if (old) call skplin(3,unitcommand)
	231	read(unitcommand,*) mquasilag
	232	if (old) call skplin(3,unitcommand)
	233	read(unitcommand,*) nested_output
	234	if (old) call skplin(3,unitcommand)
[36]	235	read(unitcommand,*,iostat=icmdstat) linit_cond
	236	if (icmdstat .gt. 0) &
	237	print*, 'readcommand: linit_cond not read properly',icmdstat,linit_cond
[24]	238	if (old) call skplin(3,unitcommand)
[36]	239	read(unitcommand,*,iostat=icmdstat) surf_only
	240	if (icmdstat .gt. 0) &
	241	print*, 'readcommand: linit_cond not read properly',icmdstat,surf_only
	242	if (old) call skplin(3,unitcommand)
	243	read(unitcommand,*,iostat=icmdstat) ldep_incr
	244	if (icmdstat .gt. 0) &
	245	print*, 'readcommand: linit_cond not read properly',icmdstat, ldep_incr
[20]	246	close(unitcommand)
	247
	248	endif ! input format
	249
	250	! write command file in namelist format to output directory if requested
[36]	251	if (nmlout .eqv. .true.) then
[20]	252	open(unitcommand,file=path(2)(1:length(2))//'COMMAND.namelist',err=1000)
	253	write(unitcommand,nml=command)
	254	close(unitcommand)
	255	endif
	256
[4]	257	ifine=max(ifine,1)
	258
	259	! Determine how Markov chain is formulated (for w or for w/sigw)
	260	!***************************************************************
	261
[36]	262	if (ctl .ge. 0.1) then
[4]	263	turbswitch=.true.
	264	else
	265	turbswitch=.false.
	266	ifine=1
	267	endif
	268	fine=1./real(ifine)
	269	ctl=1./ctl
	270
	271	! Set the switches required for the various options for input/output units
	272	!*************************************************************************
	273	!AF Set the switches IND_REL and IND_SAMP for the release and sampling
	274	!Af switches for the releasefile:
	275	!Af IND_REL = 1 : xmass * rho
	276	!Af IND_REL = 0 : xmass * 1
	277
	278	!Af switches for the conccalcfile:
	279	!AF IND_SAMP = 0 : xmass * 1
	280	!Af IND_SAMP = -1 : xmass / rho
	281
	282	!AF IND_SOURCE switches between different units for concentrations at the source
	283	!Af NOTE that in backward simulations the release of computational particles
	284	!Af takes place at the "receptor" and the sampling of p[articles at the "source".
	285	!Af 1 = mass units
	286	!Af 2 = mass mixing ratio units
	287	!Af IND_RECEPTOR switches between different units for concentrations at the receptor
	288	!Af 1 = mass units
	289	!Af 2 = mass mixing ratio units
	290
	291	if ( ldirect .eq. 1 ) then ! FWD-Run
	292	!Af set release-switch
	293	if (ind_source .eq. 1 ) then !mass
	294	ind_rel = 0
	295	else ! mass mix
	296	ind_rel = 1
	297	endif
	298	!Af set sampling switch
	299	if (ind_receptor .eq. 1) then !mass
	300	ind_samp = 0
	301	else ! mass mix
	302	ind_samp = -1
	303	endif
	304	elseif (ldirect .eq. -1 ) then !BWD-Run
	305	!Af set sampling switch
	306	if (ind_source .eq. 1 ) then !mass
	307	ind_samp = -1
	308	else ! mass mix
	309	ind_samp = 0
	310	endif
	311	!Af set release-switch
	312	if (ind_receptor .eq. 1) then !mass
	313	ind_rel = 1
	314	else ! mass mix
	315	ind_rel = 0
	316	endif
	317	endif
	318
	319	!*************************************************************
	320	! Check whether valid options have been chosen in file COMMAND
	321	!*************************************************************
	322
	323	! Check options for initial condition output: Switch off for forward runs
	324	!************************************************************************
	325
[36]	326	if (ldirect .eq. 1) linit_cond=0
	327	if ((linit_cond .lt. 0) .or. (linit_cond .gt. 2)) then
[4]	328	write(,) ' #### FLEXPART MODEL ERROR! INVALID OPTION #### '
	329	write(,) ' #### FOR LINIT_COND IN FILE "COMMAND". #### '
	330	stop
	331	endif
	332
	333	! Check input dates
	334	!******************
	335
[36]	336	if (iedate .lt. ibdate) then
[4]	337	write(,) ' #### FLEXPART MODEL ERROR! BEGINNING DATE #### '
	338	write(,) ' #### IS LARGER THAN ENDING DATE. CHANGE #### '
	339	write(,) ' #### EITHER POINT 2 OR POINT 3 IN FILE #### '
	340	write(,) ' #### "COMMAND". #### '
	341	stop
[36]	342	else if (iedate .eq. ibdate) then
	343	if (ietime .lt. ibtime) then
[4]	344	write(,) ' #### FLEXPART MODEL ERROR! BEGINNING TIME #### '
	345	write(,) ' #### IS LARGER THAN ENDING TIME. CHANGE #### '
	346	write(,) ' #### EITHER POINT 2 OR POINT 3 IN FILE #### '
	347	write(,) ' #### "COMMAND". #### '
	348	stop
	349	endif
	350	endif
	351
	352
	353	! Determine kind of dispersion method
	354	!************************************
	355
[36]	356	if (ctl .gt. 0.) then
[4]	357	method=1
	358	mintime=minstep
	359	else
	360	method=0
	361	mintime=lsynctime
	362	endif
	363
[30]	364	! check for netcdf output switch (use for non-namelist input only!)
[36]	365	if (iout .ge. 8) then
[30]	366	lnetcdfout = 1
	367	iout = iout - 8
	368	#ifndef NETCDF_OUTPUT
	369	print*,'ERROR: netcdf output not activated during compile time but used in COMMAND file!'
	370	print*,'Please recompile with netcdf library or use standard output format.'
	371	stop
	372	#endif
	373	endif
	374
[4]	375	! Check whether a valid option for gridded model output has been chosen
	376	!**********************************************************************
	377
[36]	378	if ((iout .lt. 1) .or. (iout .gt. 5)) then
[4]	379	write(,) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### '
[30]	380	write(,) ' #### IOUT MUST BE 1, 2, 3, 4, OR 5 FOR #### '
	381	write(,) ' #### STANDARD FLEXPART OUTPUT OR 9 - 13 #### '
	382	write(,) ' #### FOR NETCDF OUTPUT #### '
[4]	383	stop
	384	endif
	385
	386	!AF check consistency between units and volume mixing ratio
[36]	387	if ( ((iout .eq. 2) .or. (iout .eq. 3)) .and. &
	388	(ind_source .gt. 1 .or. ind_receptor .gt. 1) ) then
[4]	389	write(,) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### '
	390	write(,) ' #### VOLUME MIXING RATIO ONLY SUPPORTED #### '
	391	write(,) ' #### FOR MASS UNITS (at the moment) #### '
	392	stop
	393	endif
	394
	395
	396	! For quasilag output for each release is forbidden
	397	!*****************************************************************************
	398
[36]	399	if ((ioutputforeachrelease .eq. 1) .and. (mquasilag .eq. 1)) then
[4]	400	write(,) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####'
	401	write(,) '#### OUTPUTFOREACHRELEASE AND QUASILAGRANGIAN####'
	402	write(,) '#### MODE IS NOT POSSIBLE ! ####'
	403	stop
	404	endif
	405
	406
	407	! For quasilag backward is forbidden
	408	!*****************************************************************************
	409
[36]	410	if ((ldirect .lt. 0) .and. (mquasilag .eq. 1)) then
[4]	411	write(,) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####'
	412	write(,) '#### FOR BACKWARD RUNS, QUASILAGRANGIAN MODE ####'
	413	write(,) '#### IS NOT POSSIBLE ! ####'
	414	stop
	415	endif
	416
	417
	418	! For backward runs one releasefield for all releases makes no sense,
	419	! For quasilag and domainfill ioutputforechrelease is forbidden
	420	!*****************************************************************************
	421
[36]	422	if ((ldirect .lt. 0) .and. (ioutputforeachrelease .eq. 0)) then
[4]	423	write(,) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####'
	424	write(,) '#### FOR BACKWARD RUNS, IOUTPUTFOREACHRLEASE ####'
	425	write(,) '#### MUST BE SET TO ONE! ####'
	426	stop
	427	endif
	428
	429
	430	! For backward runs one releasefield for all releases makes no sense,
	431	! and is "forbidden"
	432	!*****************************************************************************
	433
[36]	434	if ((mdomainfill .eq. 1) .and. (ioutputforeachrelease .eq. 1)) then
[4]	435	write(,) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####'
	436	write(,) '#### FOR DOMAIN FILLING RUNS OUTPUT FOR ####'
	437	write(,) '#### EACH RELEASE IS FORBIDDEN ! ####'
	438	stop
	439	endif
	440
	441
	442	! For domain-filling trajectories, a plume centroid trajectory makes no sense,
	443	! For backward runs, only residence time output (iout=1) or plume trajectories (iout=4),
	444	! or both (iout=5) makes sense; other output options are "forbidden"
	445	!*****************************************************************************
	446
[36]	447	if (ldirect .lt. 0) then
	448	if ((iout .eq. 2) .or. (iout .eq. 3)) then
[4]	449	write(,) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####'
	450	write(,) '#### FOR BACKWARD RUNS, IOUT MUST BE 1,4,OR 5####'
	451	stop
	452	endif
	453	endif
	454
	455
	456	! For domain-filling trajectories, a plume centroid trajectory makes no sense,
	457	! and is "forbidden"
	458	!*****************************************************************************
	459
[36]	460	if (mdomainfill .ge. 1) then
	461	if ((iout .eq. 4) .or. (iout .eq. 5)) then
[4]	462	write(,) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####'
	463	write(,) '#### FOR DOMAIN-FILLING TRAJECTORY OPTION, ####'
	464	write(,) '#### IOUT MUST NOT BE SET TO 4 OR 5. ####'
	465	stop
	466	endif
	467	endif
	468
	469
	470
	471	! Check whether a valid options for particle dump has been chosen
	472	!****************************************************************
	473
[36]	474	if ((ipout .ne. 0) .and. (ipout .ne. 1) .and. (ipout .ne. 2)) then
[4]	475	write(,) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### '
	476	write(,) ' #### IPOUT MUST BE 1, 2 OR 3! #### '
	477	stop
	478	endif
	479
[36]	480	if(lsubgrid .ne. 1 .and. verbosity .eq. 0) then
[4]	481	write(,) ' ---------------- '
	482	write(,) ' INFORMATION: SUBGRIDSCALE TERRAIN EFFECT IS'
	483	write(,) ' NOT PARAMETERIZED DURING THIS SIMULATION. '
	484	write(,) ' ---------------- '
	485	endif
	486
	487
	488	! Check whether convection scheme is either turned on or off
	489	!***********************************************************
	490
[36]	491	if ((lconvection .ne. 0) .and. (lconvection .ne. 1)) then
[4]	492	write(,) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### '
	493	write(,) ' #### LCONVECTION MUST BE SET TO EITHER 1 OR 0#### '
	494	stop
	495	endif
	496
	497
	498	! Check whether synchronisation interval is sufficiently short
	499	!*************************************************************
	500
[36]	501	if (lsynctime .gt. (idiffnorm/2)) then
[4]	502	write(,) ' #### FLEXPART MODEL ERROR! SYNCHRONISATION #### '
	503	write(,) ' #### TIME IS TOO LONG. MAKE IT SHORTER. #### '
	504	write(,) ' #### MINIMUM HAS TO BE: ', idiffnorm/2
	505	stop
	506	endif
	507
	508
	509	! Check consistency of the intervals, sampling periods, etc., for model output
	510	!*****************************************************************************
	511
[36]	512	if (loutaver .eq. 0) then
[4]	513	write(,) ' #### FLEXPART MODEL ERROR! TIME AVERAGE OF #### '
	514	write(,) ' #### CONCENTRATION FIELD OUTPUT MUST NOT BE #### '
	515	write(,) ' #### ZERO. #### '
	516	write(,) ' #### CHANGE INPUT IN FILE COMMAND. #### '
	517	stop
	518	endif
	519
[36]	520	if (loutaver .gt. loutstep) then
[4]	521	write(,) ' #### FLEXPART MODEL ERROR! TIME AVERAGE OF #### '
	522	write(,) ' #### CONCENTRATION FIELD OUTPUT MUST NOT BE #### '
	523	write(,) ' #### GREATER THAN INTERVAL OF OUTPUT. #### '
	524	write(,) ' #### CHANGE INPUT IN FILE COMMAND. #### '
	525	stop
	526	endif
	527
[36]	528	if (loutsample .gt. loutaver) then
[4]	529	write(,) ' #### FLEXPART MODEL ERROR! SAMPLING TIME OF #### '
	530	write(,) ' #### CONCENTRATION FIELD OUTPUT MUST NOT BE #### '
	531	write(,) ' #### GREATER THAN TIME AVERAGE OF OUTPUT. #### '
	532	write(,) ' #### CHANGE INPUT IN FILE COMMAND. #### '
	533	stop
	534	endif
	535
[36]	536	if (mod(loutaver,lsynctime) .ne. 0) then
[4]	537	write(,) ' #### FLEXPART MODEL ERROR! AVERAGING TIME OF #### '
	538	write(,) ' #### CONCENTRATION FIELD MUST BE A MULTIPLE #### '
	539	write(,) ' #### OF THE SYNCHRONISATION INTERVAL #### '
	540	stop
	541	endif
	542
[36]	543	if ((loutaver/lsynctime) .lt. 2) then
[4]	544	write(,) ' #### FLEXPART MODEL ERROR! AVERAGING TIME OF #### '
	545	write(,) ' #### CONCENTRATION FIELD MUST BE AT LEAST #### '
	546	write(,) ' #### TWICE THE SYNCHRONISATION INTERVAL #### '
	547	stop
	548	endif
	549
[36]	550	if (mod(loutstep,lsynctime) .ne. 0) then
[4]	551	write(,) ' #### FLEXPART MODEL ERROR! INTERVAL BETWEEN #### '
	552	write(,) ' #### CONCENTRATION FIELDS MUST BE A MULTIPLE #### '
	553	write(,) ' #### OF THE SYNCHRONISATION INTERVAL #### '
	554	stop
	555	endif
	556
[36]	557	if ((loutstep/lsynctime) .lt. 2) then
[4]	558	write(,) ' #### FLEXPART MODEL ERROR! INTERVAL BETWEEN #### '
	559	write(,) ' #### CONCENTRATION FIELDS MUST BE AT LEAST #### '
	560	write(,) ' #### TWICE THE SYNCHRONISATION INTERVAL #### '
	561	stop
	562	endif
	563
[36]	564	if (mod(loutsample,lsynctime) .ne. 0) then
[4]	565	write(,) ' #### FLEXPART MODEL ERROR! SAMPLING TIME OF #### '
	566	write(,) ' #### CONCENTRATION FIELD MUST BE A MULTIPLE #### '
	567	write(,) ' #### OF THE SYNCHRONISATION INTERVAL #### '
	568	stop
	569	endif
	570
[36]	571	if (itsplit .lt. loutaver) then
[4]	572	write(,) ' #### FLEXPART MODEL ERROR! SPLITTING TIME FOR#### '
	573	write(,) ' #### PARTICLES IS TOO SHORT. PLEASE INCREASE #### '
	574	write(,) ' #### SPLITTING TIME CONSTANT. #### '
	575	stop
	576	endif
	577
[36]	578	if ((mquasilag .eq. 1) .and. (iout .ge. 4)) then
[4]	579	write(,) ' #### FLEXPART MODEL ERROR! CONFLICTING #### '
	580	write(,) ' #### OPTIONS: IF MQUASILAG=1, PLUME #### '
	581	write(,) ' #### TRAJECTORY OUTPUT IS IMPOSSIBLE. #### '
	582	stop
	583	endif
	584
	585	! Compute modeling time in seconds and beginning date in Julian date
	586	!*******************************************************************
	587
	588	outstep=real(abs(loutstep))
[36]	589	if (ldirect .eq. 1) then
[4]	590	bdate=juldate(ibdate,ibtime)
	591	edate=juldate(iedate,ietime)
	592	ideltas=nint((edate-bdate)*86400.)
[36]	593	else if (ldirect .eq. -1) then
[4]	594	loutaver=-1*loutaver
	595	loutstep=-1*loutstep
	596	loutsample=-1*loutsample
	597	lsynctime=-1*lsynctime
	598	bdate=juldate(iedate,ietime)
	599	edate=juldate(ibdate,ibtime)
	600	ideltas=nint((edate-bdate)*86400.)
	601	else
	602	write(,) ' #### FLEXPART MODEL ERROR! DIRECTION IN #### '
	603	write(,) ' #### FILE "COMMAND" MUST BE EITHER -1 OR 1. #### '
	604	stop
	605	endif
	606
	607	return
	608
	609	999 write(,) ' #### FLEXPART MODEL ERROR! FILE "COMMAND" #### '
	610	write(,) ' #### CANNOT BE OPENED IN THE DIRECTORY #### '
	611	write(*,'(a)') path(1)(1:length(1))
	612	stop
	613
[20]	614	1000 write(,) ' #### FLEXPART MODEL ERROR! FILE "COMMAND" #### '
[27]	615	write(,) ' #### CANNOT BE OPENED IN THE DIRECTORY #### '
	616	write(*,'(a)') path(2)(1:length(1))
	617	stop
[4]	618	end subroutine readcommand

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

Download in other formats:

Original Format