Skip to content

Commit fc3b249

Browse files
committed
added more module precision
1 parent 1fafb74 commit fc3b249

File tree

2 files changed

+4
-5
lines changed

2 files changed

+4
-5
lines changed

src/common/m_variables_conversion.fpp

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1045,10 +1045,10 @@ contains
10451045
#ifdef MFC_POST_PROCESS
10461046
if (hyperelasticity) then
10471047
! to save von Mises stress instead of elastic internal energy
1048-
qK_prim_vf(xiend+1)%sf(j, k, l) = sqrt((3d0/2d0)* (qK_prim_vf(strxb)%sf(j, k, l)**2d0 + &
1049-
2d0*qK_prim_vf(strxb+1)%sf(j, k, l)**2d0 + qK_prim_vf(strxb+2)%sf(j, k, l)**2d0 + &
1050-
2d0*qK_prim_vf(strxb+3)%sf(j, k, l)**2d0 + 2d0*qK_prim_vf(strxb+4)%sf(j, k, l)**2d0 + &
1051-
qK_prim_vf(strxe)%sf(j, k, l)**2d0) )
1048+
qK_prim_vf(xiend+1)%sf(j, k, l) = sqrt((3_wp/2_wp)* (qK_prim_vf(strxb)%sf(j, k, l)**2_wp + &
1049+
2_wp*qK_prim_vf(strxb+1)%sf(j, k, l)**2_wp + qK_prim_vf(strxb+2)%sf(j, k, l)**2_wp + &
1050+
2_wp*qK_prim_vf(strxb+3)%sf(j, k, l)**2_wp + 2_wp*qK_prim_vf(strxb+4)%sf(j, k, l)**2_wp + &
1051+
qK_prim_vf(strxe)%sf(j, k, l)**2_wp) )
10521052
end if
10531053
#endif
10541054

src/pre_process/m_patches.fpp

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1444,7 +1444,6 @@ contains
14441444
! smoothed out, the pseudo volume fraction is set to 1 to make sure
14451445
! that only the current patch contributes to the fluid state in the
14461446
! cells that this patch covers.
1447-
eps = 1.d-32
14481447
eta = 1._wp
14491448
eps = 1.e-32_wp
14501449

0 commit comments

Comments
 (0)