Fixed IBM not being read

anshgupta1234 · Dec 21, 2023 · 8e13f60 · 8e13f60
1 parent 412f601
commit 8e13f60
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 4 deletions.
diff --git a/examples/2D_ibm/case.py b/examples/2D_ibm/case.py
@@ -121,7 +121,7 @@
     # Uniform medium density, centroid is at the center of the domain
     'patch_ib(1)%y_centroid'     : 3.0E-03,
     'patch_ib(1)%radius'         : 0.2E-03,
-    'patch_ib(1)%slip'           : 'T',
+    'patch_ib(1)%slip'           : 'F',
     # 'patch_ib(1)%length_x'         : 0.2E-03,
     # 'patch_ib(1)%length_y'         : 0.4E-03,
     # # ==========================================================================

diff --git a/src/pre_process/m_data_output.fpp b/src/pre_process/m_data_output.fpp
@@ -500,8 +500,12 @@ contains
             call s_mpi_barrier()
             call DelayFileAccess (proc_rank)
             ! Initialize MPI data I/O
-            call s_initialize_mpi_data(q_cons_vf)
-
+            if (ib) then
+                call s_initialize_mpi_data(q_cons_vf, ib_markers)
+            else
+                call s_initialize_mpi_data(q_cons_vf)
+            end if
+
             ! Open the file to write all flow variables
             write (file_loc, '(I0,A,i7.7,A)') t_step_start, '_', proc_rank, '.dat'
             file_loc = trim(restart_dir)//'/lustre_0'//trim(mpiiofs)//trim(file_loc)
@@ -556,7 +560,11 @@ contains
 
         else
             ! Initialize MPI data I/O
-            call s_initialize_mpi_data(q_cons_vf)
+            if (ib) then
+                call s_initialize_mpi_data(q_cons_vf, ib_markers)
+            else
+                call s_initialize_mpi_data(q_cons_vf)
+            end if
 
             ! Open the file to write all flow variables
             write (file_loc, '(I0,A)') t_step_start, '.dat'