Context Navigation

readreleases.f90 @ 9484483

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

Last change on this file since 9484483 was 0f20c31, checked in by flexpart <>, 9 years ago
Updated SPECIES files
Property mode set to `100644`
File size: 24.6 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 readreleases
23
24	!*****************************************************************************
25	! *
26	! This routine reads the release point specifications for the current *
27	! model run. Several release points can be used at the same time. *
28	! *
29	! Author: A. Stohl *
30	! *
31	! 18 May 1996 *
32	! *
33	! Update: 29 January 2001 *
34	! Release altitude can be either in magl or masl *
35	! HSO, 12 August 2013
36	! Added optional namelist input
37	! *
38	!*****************************************************************************
39	! *
40	! Variables: *
41	! decay decay constant of species *
42	! dquer [um] mean particle diameters *
43	! dsigma e.g. dsigma=10 or dsigma=0.1 means that 68% of the mass*
44	! are between 0.1dquer and 10dquer *
45	! ireleasestart, ireleaseend [s] starting time and ending time of each *
46	! release *
47	! kindz 1: zpoint is in m agl, 2: zpoint is in m asl, 3: zpoint*
48	! is in hPa *
49	! npart number of particles to be released *
50	! nspec number of species to be released *
51	! density [kg/m3] density of the particles *
52	! rm [s/m] Mesophyll resistance *
53	! species name of species *
54	! xmass total mass of each species *
55	! xpoint1,ypoint1 geograf. coordinates of lower left corner of release *
56	! area *
57	! xpoint2,ypoint2 geograf. coordinates of upper right corner of release *
58	! area *
59	! weta, wetb parameters to determine the below-cloud scavenging *
60	! weta_in, wetb_in parameters to determine the in-cloud scavenging *
61	! wetc_in, wetd_in parameters to determine the in-cloud scavenging *
62	! zpoint1,zpoint2 height range, over which release takes place *
63	! num_min_discrete if less, release cannot be randomized and happens at *
64	! time mid-point of release interval *
65	! lroot true if serial version, or if MPI and root process *
66	! *
67	!*****************************************************************************
68
69	use point_mod
70	use xmass_mod
71	use par_mod
72	use com_mod
73
74	implicit none
75
76	integer :: numpartmax,i,j,id1,it1,id2,it2,idum,stat
77	integer,parameter :: num_min_discrete=100
78	real :: vsh(ni),fracth(ni),schmih(ni),releaserate,xdum,cun
79	real(kind=dp) :: jul1,jul2,julm,juldate
80	character(len=50) :: line
81	logical :: old
82
83	! help variables for namelist reading
84	integer :: numpoints, parts, readerror
85	integer*2 :: zkind
86	integer :: idate1, itime1, idate2, itime2
87	real :: lon1,lon2,lat1,lat2,z1,z2
88	character*40 :: comment
89	integer,parameter :: unitreleasesout=2
90	real,allocatable, dimension (:) :: mass
91	integer,allocatable, dimension (:) :: specnum_rel,specnum_rel2
92
93	! declare namelists
94	namelist /releases_ctrl/ &
95	nspec, &
96	specnum_rel
97
98	namelist /release/ &
99	idate1, itime1, &
100	idate2, itime2, &
101	lon1, lon2, &
102	lat1, lat2, &
103	z1, z2, &
104	zkind, &
105	mass, &
106	parts, &
107	comment
108
109	numpoint=0
110
111	! allocate with maxspec for first input loop
112	allocate(mass(maxspec),stat=stat)
113	if (stat.ne.0) write(,)'ERROR: could not allocate mass'
114	allocate(specnum_rel(maxspec),stat=stat)
115	if (stat.ne.0) write(,)'ERROR: could not allocate specnum_rel'
116
117	! presetting namelist releases_ctrl
118	nspec = -1 ! use negative value to determine failed namelist input
119	specnum_rel = 0
120
121	!sec, read release to find how many releasepoints should be allocated
122	open(unitreleases,file=path(1)(1:length(1))//'RELEASES',status='old',form='formatted',err=999)
123
124	! check if namelist input provided
125	read(unitreleases,releases_ctrl,iostat=readerror)
126
127	! prepare namelist output if requested
128	if (nmlout.and.lroot) then
129	open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',access='append',status='replace',err=1000)
130	endif
131
132	if ((readerror.ne.0).or.(nspec.lt.0)) then
133
134	! no namelist format, close file and allow reopening in old format
135	close(unitreleases)
136	open(unitreleases,file=path(1)(1:length(1))//'RELEASES',status='old',err=999)
137
138	readerror=1 ! indicates old format
139
140	! Check the format of the RELEASES file (either in free format,
141	! or using a formatted mask)
142	! Use of formatted mask is assumed if line 10 contains the word 'DIRECTION'
143	!**************************************************************************
144
145	call skplin(12,unitreleases)
146	read (unitreleases,901) line
147	901 format (a)
148	if (index(line,'Total') .eq. 0) then
149	old = .false.
150	else
151	old = .true.
152	endif
153	rewind(unitreleases)
154
155
156	! Skip first 11 lines (file header)
157	!**********************************
158
159	call skplin(11,unitreleases)
160
161	read(unitreleases,*,err=998) nspec
162	if (old) call skplin(2,unitreleases)
163	do i=1,nspec
164	read(unitreleases,*,err=998) specnum_rel(i)
165	if (old) call skplin(2,unitreleases)
166	end do
167
168	numpoint=0
169	100 numpoint=numpoint+1
170	read(unitreleases,*,end=25)
171	read(unitreleases,*,err=998,end=25) idum,idum
172	if (old) call skplin(2,unitreleases)
173	read(unitreleases,*,err=998) idum,idum
174	if (old) call skplin(2,unitreleases)
175	read(unitreleases,*,err=998) xdum
176	if (old) call skplin(2,unitreleases)
177	read(unitreleases,*,err=998) xdum
178	if (old) call skplin(2,unitreleases)
179	read(unitreleases,*,err=998) xdum
180	if (old) call skplin(2,unitreleases)
181	read(unitreleases,*,err=998) xdum
182	if (old) call skplin(2,unitreleases)
183	read(unitreleases,*,err=998) idum
184	if (old) call skplin(2,unitreleases)
185	read(unitreleases,*,err=998) xdum
186	if (old) call skplin(2,unitreleases)
187	read(unitreleases,*,err=998) xdum
188	if (old) call skplin(2,unitreleases)
189	read(unitreleases,*,err=998) idum
190	if (old) call skplin(2,unitreleases)
191	do i=1,nspec
192	read(unitreleases,*,err=998) xdum
193	if (old) call skplin(2,unitreleases)
194	end do
195	!save compoint only for the first 1000 release points
196	read(unitreleases,'(a40)',err=998) compoint(1)(1:40)
197	if (old) call skplin(1,unitreleases)
198
199	goto 100
200
201	25 numpoint=numpoint-1
202
203	else
204
205	readerror=0
206	do while (readerror.eq.0)
207	idate1=-1
208	read(unitreleases,release,iostat=readerror)
209	if ((idate1.lt.0).or.(readerror.ne.0)) then
210	readerror=1
211	else
212	numpoint=numpoint+1
213	endif
214	end do
215	readerror=0
216	endif ! if namelist input
217
218	rewind(unitreleases)
219
220	! allocate arrays of matching size for number of species (namelist output)
221	deallocate(mass)
222	allocate(mass(nspec),stat=stat)
223	if (stat.ne.0) write(,)'ERROR: could not allocate mass'
224	allocate(specnum_rel2(nspec),stat=stat)
225	if (stat.ne.0) write(,)'ERROR: could not allocate specnum_rel2'
226	specnum_rel2=specnum_rel(1:nspec)
227	deallocate(specnum_rel)
228	allocate(specnum_rel(nspec),stat=stat)
229	if (stat.ne.0) write(,)'ERROR: could not allocate specnum_rel'
230	specnum_rel=specnum_rel2
231	deallocate(specnum_rel2)
232
233	!allocate memory for numpoint releaspoints
234	allocate(ireleasestart(numpoint),stat=stat)
235	if (stat.ne.0) write(,)'ERROR: could not allocate ireleasestart'
236	allocate(ireleaseend(numpoint),stat=stat)
237	if (stat.ne.0) write(,)'ERROR: could not allocate ireleaseend'
238	allocate(xpoint1(numpoint),stat=stat)
239	if (stat.ne.0) write(,)'ERROR: could not allocate xpoint1'
240	allocate(xpoint2(numpoint),stat=stat)
241	if (stat.ne.0) write(,)'ERROR: could not allocate xpoint2'
242	allocate(ypoint1(numpoint),stat=stat)
243	if (stat.ne.0) write(,)'ERROR: could not allocate ypoint1'
244	allocate(ypoint2(numpoint),stat=stat)
245	if (stat.ne.0) write(,)'ERROR: could not allocate ypoint2'
246	allocate(zpoint1(numpoint),stat=stat)
247	if (stat.ne.0) write(,)'ERROR: could not allocate zpoint1'
248	allocate(zpoint2(numpoint),stat=stat)
249	if (stat.ne.0) write(,)'ERROR: could not allocate zpoint2'
250	allocate(kindz(numpoint),stat=stat)
251	if (stat.ne.0) write(,)'ERROR: could not allocate kindz'
252	allocate(xmass(numpoint,maxspec),stat=stat)
253	if (stat.ne.0) write(,)'ERROR: could not allocate xmass'
254	allocate(rho_rel(numpoint),stat=stat)
255	if (stat.ne.0) write(,)'ERROR: could not allocate rho_rel'
256	allocate(npart(numpoint),stat=stat)
257	if (stat.ne.0) write(,)'ERROR: could not allocate npart'
258	allocate(xmasssave(numpoint),stat=stat)
259	if (stat.ne.0) write(,)'ERROR: could not allocate xmasssave'
260
261	if (lroot) write (,) 'Releasepoints : ', numpoint
262
263	do i=1,numpoint
264	xmasssave(i)=0.
265	end do
266
267	!now save the information
268	DEP=.false.
269	DRYDEP=.false.
270	WETDEP=.false.
271	OHREA=.false.
272	do i=1,maxspec
273	DRYDEPSPEC(i)=.false.
274	end do
275
276	if (readerror.ne.0) then
277	! Skip first 11 lines (file header)
278	!**********************************
279
280	call skplin(11,unitreleases)
281
282	! Assign species-specific parameters needed for physical processes
283	!*************************************************************************
284
285	read(unitreleases,*,err=998) nspec
286	if (nspec.gt.maxspec) goto 994
287	if (old) call skplin(2,unitreleases)
288	endif
289
290	! namelist output
291	if (nmlout.and.lroot) then
292	write(unitreleasesout,nml=releases_ctrl)
293	endif
294
295	do i=1,nspec
296	if (readerror.ne.0) then
297	read(unitreleases,*,err=998) specnum_rel(i)
298	if (old) call skplin(2,unitreleases)
299	call readspecies(specnum_rel(i),i)
300	else
301	call readspecies(specnum_rel(i),i)
302	endif
303
304	! For backward runs, only 1 species is allowed
305	!*********************************************
306
307	!if ((ldirect.lt.0).and.(nspec.gt.1)) then
308	!write(,) '#####################################################'
309	!write(,) '#### FLEXPART MODEL SUBROUTINE READRELEASES: ####'
310	!write(,) '#### FOR BACKWARD RUNS, ONLY 1 SPECIES IS ALLOWED####'
311	!write(,) '#####################################################'
312	! stop
313	!endif
314
315	! Molecular weight
316	!*****************
317
318	if (((iout.eq.2).or.(iout.eq.3)).and.(weightmolar(i).lt.0.)) then
319	write(,) 'For mixing ratio output, valid molar weight'
320	write(,) 'must be specified for all simulated species.'
321	write(,) 'Check table SPECIES or choose concentration'
322	write(,) 'output instead if molar weight is not known.'
323	stop
324	endif
325
326	! Radioactive decay
327	!******************
328
329	decay(i)=0.693147/decay(i) !conversion half life to decay constant
330
331
332	! Dry deposition of gases
333	!************************
334
335	if (reldiff(i).gt.0.) rm(i)=1./(henry(i)/3000.+100.*f0(i)) ! mesophyll resistance
336
337	! Dry deposition of particles
338	!****************************
339
340	vsetaver(i)=0.
341	cunningham(i)=0.
342	dquer(i)=dquer(i)*1000000. ! Conversion m to um
343	if (density(i).gt.0.) then ! Additional parameters
344	call part0(dquer(i),dsigma(i),density(i),fracth,schmih,cun,vsh)
345	do j=1,ni
346	fract(i,j)=fracth(j)
347	schmi(i,j)=schmih(j)
348	vset(i,j)=vsh(j)
349	cunningham(i)=cunningham(i)+cun*fract(i,j)
350	vsetaver(i)=vsetaver(i)-vset(i,j)*fract(i,j)
351	end do
352	if (lroot) write(,) 'Average settling velocity: ',i,vsetaver(i)
353	endif
354
355	! Dry deposition for constant deposition velocity
356	!************************************************
357
358	dryvel(i)=dryvel(i)*0.01 ! conversion to m/s
359
360	! Check if wet deposition or OH reaction shall be calculated
361	!***********************************************************
362	! NIK 15.02.2015, with new wet dep scheme either A or B parameters can be assigned a positive value which switches on wet dep
363	if (weta(i).gt.0. .or. wetb(i).gt.0.) then
364	WETDEP=.true.
365	if (lroot) then
366	write (,) 'Below-cloud scavenging: ON'
367	! write (,) 'Below-cloud scavenging coefficients: ',weta(i),i
368	end if
369	else
370	if (lroot) write (,) 'Below-cloud scavenging: OFF'
371	endif
372
373	! NIK 31.01.2013 + 10.12.2013 + 15.02.2015
374	if (weta_in(i).gt.0. .or. wetb_in(i).gt.0.) then
375	WETDEP=.true.
376	if (lroot) then
377	write (,) 'In-cloud scavenging: ON'
378	! write (,) 'In-cloud scavenging coefficients: ',&
379	! &weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i
380	end if
381	else
382	if (lroot) write (,) 'In-cloud scavenging: OFF'
383	endif
384
385	if (ohcconst(i).gt.0.) then
386	OHREA=.true.
387	if (lroot) write (,) 'OHreaction switched on: ',ohcconst(i),i
388	endif
389
390	if ((reldiff(i).gt.0.).or.(density(i).gt.0.).or.(dryvel(i).gt.0.)) then
391	DRYDEP=.true.
392	DRYDEPSPEC(i)=.true.
393	endif
394
395	end do
396
397	if (WETDEP.or.DRYDEP) DEP=.true.
398
399	! Read specifications for each release point
400	!*******************************************
401	numpoints=numpoint
402	numpoint=0
403	numpartmax=0
404	releaserate=0.
405	101 numpoint=numpoint+1
406
407	if (readerror.lt.1) then ! reading namelist format
408
409	if (numpoint.gt.numpoints) goto 250
410	zkind = 1
411	mass = 0
412	parts = 0
413	comment = ' '
414	read(unitreleases,release,iostat=readerror)
415	id1=idate1
416	it1=itime1
417	id2=idate2
418	it2=itime2
419	xpoint1(numpoint)=lon1
420	xpoint2(numpoint)=lon2
421	ypoint1(numpoint)=lat1
422	ypoint2(numpoint)=lat2
423	zpoint1(numpoint)=z1
424	zpoint2(numpoint)=z2
425	kindz(numpoint)=zkind
426	do i=1,nspec
427	xmass(numpoint,i)=mass(i)
428	end do
429	npart(numpoint)=parts
430	compoint(min(1001,numpoint))=comment
431
432	! namelist output
433	if (nmlout.and.lroot) then
434	write(unitreleasesout,nml=release)
435	endif
436
437	else
438
439	read(unitreleases,*,end=250)
440	read(unitreleases,*,err=998,end=250) id1,it1
441	if (old) call skplin(2,unitreleases)
442	read(unitreleases,*,err=998) id2,it2
443	if (old) call skplin(2,unitreleases)
444	read(unitreleases,*,err=998) xpoint1(numpoint)
445	if (old) call skplin(2,unitreleases)
446	read(unitreleases,*,err=998) ypoint1(numpoint)
447	if (old) call skplin(2,unitreleases)
448	read(unitreleases,*,err=998) xpoint2(numpoint)
449	if (old) call skplin(2,unitreleases)
450	read(unitreleases,*,err=998) ypoint2(numpoint)
451	if (old) call skplin(2,unitreleases)
452	read(unitreleases,*,err=998) kindz(numpoint)
453	if (old) call skplin(2,unitreleases)
454	read(unitreleases,*,err=998) zpoint1(numpoint)
455	if (old) call skplin(2,unitreleases)
456	read(unitreleases,*,err=998) zpoint2(numpoint)
457	if (old) call skplin(2,unitreleases)
458	read(unitreleases,*,err=998) npart(numpoint)
459	if (old) call skplin(2,unitreleases)
460	do i=1,nspec
461	read(unitreleases,*,err=998) xmass(numpoint,i)
462	if (old) call skplin(2,unitreleases)
463	mass(i)=xmass(numpoint,i)
464	end do
465	!save compoint only for the first 1000 release points
466	if (numpoint.le.1000) then
467	read(unitreleases,'(a40)',err=998) compoint(numpoint)(1:40)
468	comment=compoint(numpoint)(1:40)
469	else
470	read(unitreleases,'(a40)',err=998) compoint(1001)(1:40)
471	comment=compoint(1001)(1:40)
472	endif
473	if (old) call skplin(1,unitreleases)
474
475	! namelist output
476	if (nmlout.and.lroot) then
477	idate1=id1
478	itime1=it1
479	idate2=id2
480	itime2=it2
481	lon1=xpoint1(numpoint)
482	lon2=xpoint2(numpoint)
483	lat1=ypoint1(numpoint)
484	lat2=ypoint2(numpoint)
485	z1=zpoint1(numpoint)
486	z2=zpoint2(numpoint)
487	zkind=kindz(numpoint)
488	parts=npart(numpoint)
489	write(unitreleasesout,nml=release)
490	endif
491
492	if (numpoint.le.1000) then
493	if((xpoint1(numpoint).eq.0.).and.(ypoint1(numpoint).eq.0.).and. &
494	(xpoint2(numpoint).eq.0.).and.(ypoint2(numpoint).eq.0.).and. &
495	(compoint(numpoint)(1:8).eq.' ')) goto 250
496	else
497	if((xpoint1(numpoint).eq.0.).and.(ypoint1(numpoint).eq.0.).and. &
498	(xpoint2(numpoint).eq.0.).and.(ypoint2(numpoint).eq.0.)) goto 250
499	endif
500
501	endif ! if namelist format
502
503	! If a release point contains no particles, stop and issue error message
504	!***********************************************************************
505
506	if (npart(numpoint).eq.0) then
507	write(,) 'FLEXPART MODEL ERROR'
508	write(,) 'RELEASES file is corrupt.'
509	write(,) 'At least for one release point, there are zero'
510	write(,) 'particles released. Make changes to RELEASES.'
511	stop
512	endif
513
514	! Check whether x coordinates of release point are within model domain
515	!*********************************************************************
516
517	if (xpoint1(numpoint).lt.xlon0) &
518	xpoint1(numpoint)=xpoint1(numpoint)+360.
519	if (xpoint1(numpoint).gt.xlon0+(nxmin1)*dx) &
520	xpoint1(numpoint)=xpoint1(numpoint)-360.
521	if (xpoint2(numpoint).lt.xlon0) &
522	xpoint2(numpoint)=xpoint2(numpoint)+360.
523	if (xpoint2(numpoint).gt.xlon0+(nxmin1)*dx) &
524	xpoint2(numpoint)=xpoint2(numpoint)-360.
525
526	! Determine relative beginning and ending times of particle release
527	!******************************************************************
528
529	jul1=juldate(id1,it1)
530	jul2=juldate(id2,it2)
531	julm=(jul1+jul2)/2.
532	if (jul1.gt.jul2) then
533	write(,) 'FLEXPART MODEL ERROR'
534	write(,) 'Release stops before it begins.'
535	write(,) 'Make changes to file RELEASES.'
536	stop
537	endif
538	if (mdomainfill.eq.0) then ! no domain filling
539	if (ldirect.eq.1) then
540	if ((jul1.lt.bdate).or.(jul2.gt.edate)) then
541	write(,) 'FLEXPART MODEL ERROR'
542	write(,) 'Release starts before simulation begins or ends'
543	write(,) 'after simulation stops.'
544	write(,) 'Make files COMMAND and RELEASES consistent.'
545	stop
546	endif
547	if (npart(numpoint).gt.num_min_discrete) then
548	ireleasestart(numpoint)=int((jul1-bdate)*86400.)
549	ireleaseend(numpoint)=int((jul2-bdate)*86400.)
550	else
551	ireleasestart(numpoint)=int((julm-bdate)*86400.)
552	ireleaseend(numpoint)=int((julm-bdate)*86400.)
553	endif
554	else if (ldirect.eq.-1) then
555	if ((jul1.lt.edate).or.(jul2.gt.bdate)) then
556	write(,) 'FLEXPART MODEL ERROR'
557	write(,) 'Release starts before simulation begins or ends'
558	write(,) 'after simulation stops.'
559	write(,) 'Make files COMMAND and RELEASES consistent.'
560	stop
561	endif
562	if (npart(numpoint).gt.num_min_discrete) then
563	ireleasestart(numpoint)=int((jul1-bdate)*86400.)
564	ireleaseend(numpoint)=int((jul2-bdate)*86400.)
565	else
566	ireleasestart(numpoint)=int((julm-bdate)*86400.)
567	ireleaseend(numpoint)=int((julm-bdate)*86400.)
568	endif
569	endif
570	endif
571
572	! Determine the release rate (particles per second) and total number
573	! of particles released during the simulation
574	!*******************************************************************
575
576	if (ireleasestart(numpoint).ne.ireleaseend(numpoint)) then
577	releaserate=releaserate+real(npart(numpoint))/ &
578	real(ireleaseend(numpoint)-ireleasestart(numpoint))
579	else
580	releaserate=99999999
581	endif
582	numpartmax=numpartmax+npart(numpoint)
583	goto 101
584
585	250 close(unitreleases)
586
587	if (nmlout.and.lroot) then
588	close(unitreleasesout)
589	endif
590
591	if (lroot) write (,) 'Particles allocated (maxpart) : ',maxpart
592	if (lroot) write (,) 'Particles released (numpartmax): ',numpartmax
593	numpoint=numpoint-1
594
595	if (ioutputforeachrelease.eq.1) then
596	maxpointspec_act=numpoint
597	else
598	maxpointspec_act=1
599	endif
600
601	! Check, whether the total number of particles may exceed totally allowed
602	! number of particles at some time during the simulation
603	!************************************************************************
604
605	if (releaserate.gt. &
606	0.99*real(maxpart)/real(lage(nageclass))) then
607	if (numpartmax.gt.maxpart.and.lroot) then
608	write(,) '#####################################################'
609	write(,) '#### FLEXPART MODEL SUBROUTINE READRELEASES: ####'
610	write(,) '#### ####'
611	write(,) '####WARNING - TOTAL NUMBER OF PARTICLES SPECIFIED####'
612	write(,) '#### IN FILE "RELEASES" MAY AT SOME POINT DURING ####'
613	write(,) '#### THE SIMULATION EXCEED THE MAXIMUM ALLOWED ####'
614	write(,) '#### NUMBER (MAXPART).IF RELEASES DO NOT OVERLAP,####'
615	write(,) '#### FLEXPART CAN POSSIBLY COMPLETE SUCCESSFULLY.####'
616	write(,) '#### HOWEVER, FLEXPART MAY HAVE TO STOP ####'
617	write(,) '#### AT SOME TIME DURING THE SIMULATION. PLEASE ####'
618	write(,) '#### MAKE SURE THAT YOUR SETTINGS ARE CORRECT. ####'
619	write(,) '#####################################################'
620	write(,) 'Maximum release rate may be: ',releaserate, &
621	' particles per second'
622	write(,) 'Maximum allowed release rate is: ', &
623	real(maxpart)/real(lage(nageclass)),' particles per second'
624	write(,) &
625	'Total number of particles released during the simulation is: ', &
626	numpartmax
627	write(,) 'Maximum allowed number of particles is: ',maxpart
628	endif
629	endif
630
631	return
632
633	994 write(,) '#####################################################'
634	write(,) '#### FLEXPART MODEL SUBROUTINE READRELEASES: ####'
635	write(,) '#### ####'
636	write(,) '#### ERROR - MAXIMUM NUMBER OF EMITTED SPECIES IS####'
637	write(,) '#### TOO LARGE. PLEASE REDUCE NUMBER OF SPECIES. ####'
638	write(,) '#####################################################'
639	stop
640
641	998 write(,) '#####################################################'
642	write(,) '#### FLEXPART MODEL SUBROUTINE READRELEASES: ####'
643	write(,) '#### ####'
644	write(,) '#### FATAL ERROR - FILE "RELEASES" IS ####'
645	write(,) '#### CORRUPT. PLEASE CHECK YOUR INPUTS FOR ####'
646	write(,) '#### MISTAKES OR GET A NEW "RELEASES"- ####'
647	write(,) '#### FILE ... ####'
648	write(,) '#####################################################'
649	stop
650
651
652	999 write(,) '#####################################################'
653	write(,) ' FLEXPART MODEL SUBROUTINE READRELEASES: '
654	write(,)
655	write(,) 'FATAL ERROR - FILE CONTAINING PARTICLE RELEASE POINTS'
656	write(,) 'POINTS IS NOT AVAILABLE OR YOU ARE NOT'
657	write(,) 'PERMITTED FOR ANY ACCESS'
658	write(,) '#####################################################'
659	stop
660
661	1000 write(,) ' #### FLEXPART MODEL ERROR! FILE "RELEASES" #### '
662	write(,) ' #### CANNOT BE OPENED IN THE DIRECTORY #### '
663	write(*,'(a)') path(2)(1:length(2))
664	stop
665
666	end subroutine readreleases

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

Download in other formats:

Original Format