Changeset 478e9e6 in flexpart.git for src

Timestamp:

Sep 1, 2015, 9:31:33 AM (9 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

f55fdce

Parents:

cda0919

Message:

dates file is now created at program start, instead of appending to it if it already exists

Location:

src

Files:

• concoutput.f90 (modified) (2 diffs)
• concoutput_mpi.f90 (modified) (2 diffs)

src/concoutput.f90

@@ -97 +97 @@
   real,parameter :: weightair=28.97
   logical :: sp_zer
+  LOGICAL,save :: init=.true.
   character :: adate*8,atime*6
   character(len=3) :: anspec
   integer :: mind
 ! mind        eso:added to ensure identical results between 2&3-fields versions
-
+  CHARACTER(LEN=8),save :: file_stat='REPLACE'
 
   ! Determine current calendar date, needed for the file name
@@ -110 +111 @@
   write(adate,'(i8.8)') jjjjmmdd
   write(atime,'(i6.6)') ihmmss
-  open(unitdates,file=path(2)(1:length(2))//'dates', ACCESS='APPEND')
+  OPEN(unitdates,file=path(2)(1:length(2))//'dates', ACCESS='APPEND', STATUS=file_stat)
   write(unitdates,'(a)') adate//atime
   close(unitdates)  
+
+  ! Overwrite existing dates file on first call, later append to it
+  ! This fixes a bug where the dates file kept growing across multiple runs
+  ! TODO check if the 'always APPEND'-behaviour is useful in other scenarioes
+  ! e.g. (restart?)
+  IF (init) THEN
+    file_stat='OLD'
+    init=.false.
+  END IF
+
+
 
   ! For forward simulations, output fields have dimension MAXSPEC,

src/concoutput_mpi.f90

@@ -102 +102 @@
   real,parameter :: weightair=28.97
   logical :: sp_zer
+  LOGICAL,save :: init=.true.
   character :: adate*8,atime*6
   character(len=3) :: anspec
   integer :: mind
 ! mind        eso:added to ensure identical results between 2&3-fields versions
+  CHARACTER(LEN=8),save :: file_stat='REPLACE'
 
 ! Measure execution time
@@ -118 +120 @@
   write(adate,'(i8.8)') jjjjmmdd
   write(atime,'(i6.6)') ihmmss
-  open(unitdates,file=path(2)(1:length(2))//'dates', ACCESS='APPEND')
+  OPEN(unitdates,file=path(2)(1:length(2))//'dates', ACCESS='APPEND', STATUS=file_stat)
   write(unitdates,'(a)') adate//atime
   close(unitdates)  
+
+  ! Overwrite existing dates file on first call, later append to it
+  ! This fixes a bug where the dates file kept growing across multiple runs
+  ! TODO check if the 'always APPEND'-behaviour is useful in other scenarioes
+  ! e.g. (restart?)
+  IF (init) THEN
+    file_stat='OLD'
+    init=.false.
+  END IF
+
 
 ! For forward simulations, output fields have dimension MAXSPEC,