attempt fix for double precision with mpi support

Author: arciyer123

src/common/m_phase_change.fpp

@@ -247,7 +247,7 @@ contains
 
                     ! entropy
                     sk(1:num_fluids) = cvs(1:num_fluids)*log((TS**gs_min(1:num_fluids)) &
-                                                              /((pS + ps_inf(1:num_fluids))**(gs_min(1:num_fluids) - 1.0_wp))) + qvps(1:num_fluids)
+                                                             /((pS + ps_inf(1:num_fluids))**(gs_min(1:num_fluids) - 1.0_wp))) + qvps(1:num_fluids)
 
                     ! enthalpy
                     hk(1:num_fluids) = gs_min(1:num_fluids)*cvs(1:num_fluids)*TS &

src/common/m_precision_select.f90

@@ -1,5 +1,4 @@
-!>
-!! @file m_precision_select.f90
+!> @file m_precision_select.f90
 !! @brief Contains module m_precision_select
 
 !> @brief This file contains the definition of floating point used in MFC
@@ -10,14 +9,18 @@ module m_precision_select
 
     implicit none
 
+    ! Define the available precision types
     integer, parameter :: single_precision = selected_real_kind(6, 37)
     integer, parameter :: double_precision = selected_real_kind(15, 307)
 
-    integer, parameter :: wp = double_precision
+    ! Set the working precision (wp) to single or double precision
+    integer, parameter :: wp = double_precision  ! Change to single_precision if needed
+
 #ifdef MFC_MPI
-    integer, parameter :: mpi_p = MPI_DOUBLE_PRECISION
+    ! Set mpi_p based on wp using the merge intrinsic function
+    integer, parameter :: mpi_p = merge(MPI_DOUBLE_PRECISION, MPI_FLOAT, wp == double_precision)
 #else
-    integer, parameter :: mpi_p = -100
+    integer, parameter :: mpi_p = -100  ! Default value when MPI is not used
 #endif
 
 end module m_precision_select

src/pre_process/m_patches.fpp

@@ -663,7 +663,7 @@ contains
         do j = 0, n
             do i = 0, m
                 myr = sqrt((x_cc(i) - x_centroid)**2 &
-                            + (y_cc(j) - y_centroid)**2)
+                           + (y_cc(j) - y_centroid)**2)
 
                 if (myr <= radius + thickness/2._wp .and. &
                     myr >= radius - thickness/2._wp .and. &
@@ -726,7 +726,7 @@ contains
             do j = 0, n
                 do i = 0, m
                     myr = sqrt((x_cc(i) - x_centroid)**2 &
-                                + (y_cc(j) - y_centroid)**2)
+                               + (y_cc(j) - y_centroid)**2)
 
                     if (myr <= radius + thickness/2._wp .and. &
                         myr >= radius - thickness/2._wp .and. &
@@ -1054,8 +1054,8 @@ contains
 
                 if (patch_icpp(patch_id)%smoothen) then
                     eta = 5e-1_wp + 5e-1_wp*tanh(smooth_coeff/min(dx, dy) &
-                                           *(a*x_cc(i) + b*y_cc(j) + c) &
-                                           /sqrt(a**2 + b**2))
+                                                 *(a*x_cc(i) + b*y_cc(j) + c) &
+                                                 /sqrt(a**2 + b**2))
                 end if
 
                 if ((a*x_cc(i) + b*y_cc(j) + c >= 0._wp &
@@ -1949,10 +1949,10 @@ contains
 
                     if (patch_icpp(patch_id)%smoothen) then
                         eta = 5e-1_wp + 5e-1_wp*tanh(smooth_coeff/min(dx, dy, dz) &
-                                               *(a*x_cc(i) + &
-                                                 b*cart_y + &
-                                                 c*cart_z + d) &
-                                               /sqrt(a**2 + b**2 + c**2))
+                                                     *(a*x_cc(i) + &
+                                                       b*cart_y + &
+                                                       c*cart_z + d) &
+                                                     /sqrt(a**2 + b**2 + c**2))
                     end if
 
                     if ((a*x_cc(i) + b*cart_y + c*cart_z + d >= 0._wp &

src/simulation/m_bubbles.fpp

@@ -140,7 +140,7 @@ contains
                             divu%sf(j, k, l) = 0._wp
                             divu%sf(j, k, l) = &
                                 5e-1_wp/dx(j)*(q_prim_vf(contxe + idir)%sf(j + 1, k, l) - &
-                                            q_prim_vf(contxe + idir)%sf(j - 1, k, l))
+                                               q_prim_vf(contxe + idir)%sf(j - 1, k, l))
 
                         end do
                     end do
@@ -155,7 +155,7 @@ contains
                     do j = 0, m
                         divu%sf(j, k, l) = divu%sf(j, k, l) + &
                                            5e-1_wp/dy(k)*(q_prim_vf(contxe + idir)%sf(j, k + 1, l) - &
-                                                       q_prim_vf(contxe + idir)%sf(j, k - 1, l))
+                                                          q_prim_vf(contxe + idir)%sf(j, k - 1, l))
 
                     end do
                 end do
@@ -169,7 +169,7 @@ contains
                     do j = 0, m
                         divu%sf(j, k, l) = divu%sf(j, k, l) + &
                                            5e-1_wp/dz(l)*(q_prim_vf(contxe + idir)%sf(j, k, l + 1) - &
-                                                       q_prim_vf(contxe + idir)%sf(j, k, l - 1))
+                                                          q_prim_vf(contxe + idir)%sf(j, k, l - 1))
 
                     end do
                 end do