Fix 3D IBM Infinite CFL Number on GPUs (#519)

Co-authored-by: Samuel Briney <[email protected]>

Author: Sam-Briney

src/simulation/m_start_up.fpp

@@ -1417,6 +1417,10 @@ contains
             !$acc update device(palpha_eps, ptgalpha_eps)
         end if
 
+        if (ib) then
+            !$acc update device(ib_markers%sf)
+        end if
+
     end subroutine s_initialize_gpu_vars
 
     subroutine s_finalize_modules

src/simulation/m_time_steppers.fpp

@@ -303,14 +303,6 @@ contains
 
         call s_compute_rhs(q_cons_ts(1)%vf, q_prim_vf, rhs_vf, pb_ts(1)%sf, rhs_pb, mv_ts(1)%sf, rhs_mv, t_step, time_avg)
 
-        if (ib .and. t_step == 1) then
-            if (qbmm .and. .not. polytropic) then
-                call s_ibm_correct_state(q_cons_ts(1)%vf, q_prim_vf, pb_ts(1)%sf, mv_ts(1)%sf)
-            else
-                call s_ibm_correct_state(q_cons_ts(1)%vf, q_prim_vf)
-            end if
-        end if
-
 #ifdef DEBUG
         print *, 'got rhs'
 #endif
@@ -420,14 +412,6 @@ contains
 
         call s_compute_rhs(q_cons_ts(1)%vf, q_prim_vf, rhs_vf, pb_ts(1)%sf, rhs_pb, mv_ts(1)%sf, rhs_mv, t_step, time_avg)
 
-        if (ib .and. t_step == 1) then
-            if (qbmm .and. .not. polytropic) then
-                call s_ibm_correct_state(q_cons_ts(1)%vf, q_prim_vf, pb_ts(1)%sf, mv_ts(1)%sf)
-            else
-                call s_ibm_correct_state(q_cons_ts(1)%vf, q_prim_vf)
-            end if
-        end if
-
         if (run_time_info) then
             call s_write_run_time_information(q_prim_vf, t_step)
         end if
@@ -609,14 +593,6 @@ contains
 
         call s_compute_rhs(q_cons_ts(1)%vf, q_prim_vf, rhs_vf, pb_ts(1)%sf, rhs_pb, mv_ts(1)%sf, rhs_mv, t_step, time_avg)
 
-        if (ib .and. t_step == 1) then
-            if (qbmm .and. .not. polytropic) then
-                call s_ibm_correct_state(q_cons_ts(1)%vf, q_prim_vf, pb_ts(1)%sf, mv_ts(1)%sf)
-            else
-                call s_ibm_correct_state(q_cons_ts(1)%vf, q_prim_vf)
-            end if
-        end if
-
         if (run_time_info) then
             call s_write_run_time_information(q_prim_vf, t_step)
         end if