Made it through 66% of simulation

Author: danieljvickers

src/common/m_variables_conversion.fpp

@@ -741,7 +741,7 @@ contains
 
         type(scalar_field), dimension(sys_size), intent(in) :: qK_cons_vf
 
-        real(wp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(inout) :: mv
+        real(stp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(inout) :: mv
 
         integer :: i, j, k, l
         real(wp) :: mu, sig, nbub_sc
@@ -773,8 +773,8 @@ contains
     pure subroutine s_initialize_pb(qK_cons_vf, mv, pb)
         type(scalar_field), dimension(sys_size), intent(in) :: qK_cons_vf
 
-        real(wp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(in) :: mv
-        real(wp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(inout) :: pb
+        real(stp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(in) :: mv
+        real(stp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(inout) :: pb
 
         integer :: i, j, k, l
         real(wp) :: mu, sig, nbub_sc

src/pre_process/m_patches.fpp

@@ -2182,7 +2182,7 @@ contains
 
                             ! Correct the sign of the levelset
                             if (patch_id_fp(i, j, k) > 0) then
-                                STL_levelset%sf(i, j, k, patch_id) = -abs(STL_levelset%sf(i, j, k, patch_id))
+                                STL_levelset%sf(i, j, k, patch_id) = -real(abs(real(STL_levelset%sf(i, j, k, patch_id), kind=wp)), kind=stp) ! TODO :: Make sure this isn't slow. It feels bad.
                             end if
 
                             ! Correct the sign of the levelset_norm
@@ -2208,7 +2208,7 @@ contains
 
                             ! Correct the sign of the levelset
                             if (patch_id_fp(i, j, k) > 0) then
-                                STL_levelset%sf(i, j, 0, patch_id) = -abs(STL_levelset%sf(i, j, 0, patch_id))
+                                STL_levelset%sf(i, j, 0, patch_id) = -real(abs(real(STL_levelset%sf(i, j, 0, patch_id), kind=wp)), kind=stp)
                             end if
 
                             ! Get the boundary normals

src/pre_process/m_perturbation.fpp

@@ -98,7 +98,7 @@ contains
         do q = 1, elliptic_smoothing_iters
 
             ! Communication of buffer regions and apply boundary conditions
-            call s_populate_variables_buffers(bc_type, q_prim_vf, pb%sf, mv%sf)
+            call s_populate_variables_buffers(bc_type, q_prim_vf, real(pb%sf, kind=wp), real(mv%sf, kind=wp))
 
             ! Perform smoothing and store in temp array
             if (n == 0) then