bug fixes

Author: wilfonba

src/pre_process/m_global_parameters.fpp

@@ -825,6 +825,8 @@ contains
             sys_size = eqn_idx%species%end
         end if
 
+        if (bubbles_lagrange) fd_number = max(1, fd_order/2)
+
         call s_configure_coordinate_bounds(recon_type, weno_polyn, muscl_polyn, igr_order, buff_size, idwint, idwbuff, viscous, &
                                            & bubbles_lagrange, m, n, p, num_dims, igr, ib, fd_number)
 

src/simulation/m_bubbles_EL.fpp

@@ -703,7 +703,7 @@ contains
                         case (1)
                             mtn_dposdt(k, l, stage) = f_interpolate_velocity(myPos(l), cell, l, q_prim_vf)
                             mtn_dveldt(k, l, stage) = 0._wp
-                            mtn_vel(k, l, stage) = mtn_dposdt(k, l, stage)
+                            mtn_vel(k, l, 1) = mtn_dposdt(k, l, stage)
                         case (2)
                             mtn_dposdt(k, l, stage) = myVel(l)
                             f_b = f_get_bubble_force(myPos(l), myR, myV, myVel(l), myMass_n, myMass_v, Re(1), myRho, cell, l, &

toolchain/mfc/test/case.py

@@ -446,7 +446,7 @@ def create_input_lagrange(path_test):
         os.mkdir(folder_path_lagrange)
 
     with open(file_path_lagrange, "w") as file:
-        file.write("0.5\t0.5\t0.5\t0.0\t0.0\t0.0\t8.0e-03\t0.0")
+        file.write("0.5\t0.5\t0.5\t0.0\t0.0\t0.0\t8.0e-03\t0.0\n")
 
 
 def copy_input_lagrange(path_example_input, path_test):