Enable Restart Probe Files

This allows data to be appended in the existing probe files instead of overwriting them when the simulation restarts.
MFlowCode · Nov 5, 2024 · aafd5bc · aafd5bc
1 parent 8aaf0d6
commit aafd5bc
Showing 1 changed file with 14 additions and 3 deletions.
diff --git a/src/simulation/m_data_output.fpp b/src/simulation/m_data_output.fpp
@@ -175,6 +175,7 @@ contains
             !! Relative path to the probe data file in the case directory
 
         integer :: i !< Generic loop iterator
+        logical :: file_exist
 
         do i = 1, num_probes
             ! Generating the relative path to the data file
@@ -183,9 +184,19 @@ contains
 
             ! Creating the formatted data file and setting up its
             ! structure
-            open (i + 30, FILE=trim(file_path), &
-                  FORM='formatted', &
-                  STATUS='unknown')
+            inquire(file=trim(file_path), exist=file_exist)
+
+            if (file_exist) then
+                open (i + 30, FILE=trim(file_path), &
+                            FORM='formatted', &
+                            STATUS='old', &
+                            POSITION='append')
+            else
+                open (i + 30, FILE=trim(file_path), &
+                                FORM='formatted', &
+                                STATUS='unknown')
+            endif
+
             ! POSITION = 'append', &
             !WRITE(i+30,'(A,I0,A)') 'Probe ',i, ' located at:'
             !WRITE(i+30,'(A,F10.6)') 'x = ',probe(i)%x