Changeset 242571d in flexpart.git for src

Timestamp:

Mar 1, 2015, 3:57:22 PM (9 years ago)

Author:

Ignacio Pisso <Ignacio.Pisso@…>

Branches:

master, 10.4.1_pesei, FPv9.3.1, FPv9.3.1b_testing, FPv9.3.2, GFS_025, bugfixes+enhancements, dev, fp9.3.1-20161214-nc4, grib2nc4_repair, release-10, release-10.4.1, scaling-bug, univie

Children:

18c1336

Parents:

43c8684

Message:

add verbose messages to readreleases.f90. clobbering issue yet to be finished

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',access='append',status='new',err=1000)
+    open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',err=1000)
   endif
 
@@ -258 +259 @@
   if (stat.ne.0) write(*,*)'ERROR: could not allocate xmasssave'
 
-  write (*,*) 'Releasepoints : ', numpoint
+  write (*,*) 'readreleases> Releasepoints : ', numpoint
 
   do i=1,numpoint
@@ -293 +294 @@
 
   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)
@@ -494 +499 @@
   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
   !***********************************************************************
@@ -536 +562 @@
         write(*,*) 'after simulation stops.'
         write(*,*) 'Make files COMMAND and RELEASES consistent.'
+        write(*,*) jul1, ' < ' , bdate 
+        write(*,*) ' .or. '
+        write(*,*) jul2 , ' > ', edate
+        
         stop
       endif
@@ -563 +593 @@
   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