Context Navigation

readcommand.f90 @ 4c64400

10.4.1_peseiGFS_025bugfixes+enhancementsdevrelease-10release-10.4.1scaling-bugunivie

Last change on this file since 4c64400 was 4c64400, checked in by Espen Sollum ATMOS <eso@…>, 7 years ago
Bugfix for double precision dry deposition calculation. Added (hardcoded) option to not use output kernel. Version/date string updated.
Property mode set to `100644`
File size: 23.5 KB

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

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

Download in other formats:

Original Format