Added feedback for if variables are not the same in non-MPI

Author: danieljvickers

src/common/m_mpi_common.fpp

@@ -492,6 +492,8 @@ contains
                                mpi_p, MPI_SUM, MPI_COMM_WORLD, ierr)
         end if
 
+#else
+        var_glb = var_loc
 #endif
 
     end subroutine s_mpi_allreduce_vectors_sum

src/simulation/m_ibm.fpp

@@ -1017,10 +1017,11 @@ contains
                     ib_idx = ib_markers%sf(i, j, k)
                     if (ib_idx /= 0) then ! only need to compute the gradient for cells inside a IB
                         if (patch_ib(ib_idx)%moving_ibm == 2) then ! make sure that this IB has 2-way coupling enabled
+                            ! get the vector pointing to the grid cell from the IB centroid
                             if (num_dims == 3) then
-                                radial_vector = [x_cc(i), y_cc(j), z_cc(k)] - [patch_ib(ib_idx)%x_centroid, patch_ib(ib_idx)%y_centroid, patch_ib(ib_idx)%z_centroid] ! get the vector pointing to the grid cell
+                                radial_vector = [x_cc(i), y_cc(j), z_cc(k)] - [patch_ib(ib_idx)%x_centroid, patch_ib(ib_idx)%y_centroid, patch_ib(ib_idx)%z_centroid]
                             else
-                                radial_vector = [x_cc(i), y_cc(j), 0._wp] - [patch_ib(ib_idx)%x_centroid, patch_ib(ib_idx)%y_centroid, 0._wp] ! get the vector pointing to the grid cell
+                                radial_vector = [x_cc(i), y_cc(j), 0._wp] - [patch_ib(ib_idx)%x_centroid, patch_ib(ib_idx)%y_centroid, 0._wp]
                             end if
 
                             ! use a finite difference to compute the 2D components of the gradient of the pressure and cell volume
@@ -1040,7 +1041,7 @@ contains
                             $:GPU_ATOMIC(atomic='update')
                             forces(ib_idx, :) = forces(ib_idx, :) - (pressure_divergence*cell_volume)
                             $:GPU_ATOMIC(atomic='update')
-                            torques(ib_idx, :) = torques(ib_idx, :) + (cross_product(radial_vector, pressure_divergence)*cell_volume)
+                            torques(ib_idx, :) = torques(ib_idx, :) - (cross_product(radial_vector, pressure_divergence)*cell_volume)
                         end if
                     end if
                 end do