We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 624c36d commit 7134cc0Copy full SHA for 7134cc0
src/common/m_variables_conversion.fpp
@@ -942,7 +942,7 @@ contains
942
if (model_eqns /= 4) then
943
qK_prim_vf(i)%sf(j, k, l) = qK_cons_vf(i)%sf(j, k, l) &
944
/rho_K
945
- dyn_pres_K = dyn_pres_K + 5e-1_wp*qK_cons_vf(i)%sf(j, k, l) &
+ dyn_pres_K = dyn_pres_K + 0.5_wp*qK_cons_vf(i)%sf(j, k, l) &
946
*qK_prim_vf(i)%sf(j, k, l)
947
else
948
@@ -1341,7 +1341,7 @@ contains
1341
1342
! Computing the energy from the pressure
1343
E_K = gamma_K*pres_K + pi_inf_K &
1344
- + 5e-1_wp*rho_K*vel_K_sum + qv_K
+ + 0.5_wp*rho_K*vel_K_sum + qv_K
1345
1346
! mass flux, this should be \alpha_i \rho_i u_i
1347
!$acc loop seq
@@ -1466,7 +1466,7 @@ contains
1466
(rho*(1._wp - adv(num_fluids)))
1467
end if
1468
1469
- c = ((H - 5e-1*vel_sum)/gamma)
+ c = ((H - 0.5_wp*vel_sum)/gamma)
1470
1471
1472
if (mixture_err .and. c < 0._wp) then
0 commit comments