Changes in src/readreleases.f90 [242571d:b4d29ce] in flexpart.git

File:

• src/readreleases.f90 (modified) (6 diffs)

src/readreleases.f90

@@ -126 +126 @@
   ! prepare namelist output if requested
   if (nmlout.eqv..true.) then
-    !open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',access='append',status='new',err=1000)
-    open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',err=1000)
+    open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',access='append',status='new',err=1000)
   endif
 
@@ -259 +258 @@
   if (stat.ne.0) write(*,*)'ERROR: could not allocate xmasssave'
 
-  write (*,*) 'readreleases> Releasepoints : ', numpoint
+  write (*,*) 'Releasepoints : ', numpoint
 
   do i=1,numpoint
@@ -294 +293 @@
 
   do i=1,nspec
-    if (verbosity.gt.0) then
-      print*, 'readreleases> call readspecies', i
-    endif
- 
     if (readerror.ne.0) then
       read(unitreleases,*,err=998) specnum_rel(i)
@@ -499 +494 @@
   endif ! if namelist format
 
-
-  if (verbosity.gt.1 .and. numpoint.eq.1) then ! verbosity 2 or larger
-    write(*,*) 'numpoint=', numpoint
-    print*,  id1,it1
-    print*,  id2,it2
-    print*,  xpoint1(numpoint)
-    print*,  ypoint1(numpoint)
-    print*,  xpoint2(numpoint)
-    print*,  ypoint2(numpoint)
-    print*,  'kindz=' , kindz(numpoint)
-    print*,  zpoint1(numpoint)
-    print*,  zpoint2(numpoint)
-    print*,  npart(numpoint)
-    do i=1,nspec
-      !mass(i)=
-      print*, 'xmass=', xmass(numpoint,i)
-    end do
-    print*, compoint(numpoint) 
-  endif
-
-
   ! If a release point contains no particles, stop and issue error message
   !***********************************************************************
@@ -562 +536 @@
         write(*,*) 'after simulation stops.'
         write(*,*) 'Make files COMMAND and RELEASES consistent.'
-        write(*,*) jul1, ' < ' , bdate 
-        write(*,*) ' .or. '
-        write(*,*) jul2 , ' > ', edate
-        
         stop
       endif
@@ -593 +563 @@
   endif
 
-  if (verbosity.gt.1 .and. numpoint.eq.1) then ! verbosity 2 or larger
-    print*, 'ireleasestart(',numpoint,')', ireleasestart(numpoint) 
-    print*, 'ireleaseend(',numpoint,')', ireleaseend(numpoint) 
-  endif
-
   ! Determine the release rate (particles per second) and total number
   ! of particles released during the simulation