@@ -1379,7 +1379,7 @@ contains
13791379
13801380 ! ENERGY ADJUSTMENTS FOR HYPOELASTIC/HYPERELASTIC ENERGY
13811381 if (hypoelasticity .or. hyperelasticity) then
1382- G_L = 0_wp ; G_R = 0_wp
1382+ G_L = 0._wp ; G_R = 0._wp
13831383 !$acc loop seq
13841384 do i = 1, num_fluids
13851385 G_L = G_L + alpha_L(i)*Gs(i)
@@ -1395,12 +1395,12 @@ contains
13951395 do i = 1, strxe - strxb + 1
13961396 ! Elastic contribution to energy if G large enough
13971397 if ((G_L > verysmall) .and. (G_R > verysmall)) then
1398- E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4_wp *G_L)
1399- E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4_wp *G_R)
1398+ E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4._wp *G_L)
1399+ E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4._wp *G_R)
14001400 ! Additional terms in 2D and 3D
14011401 if ((i == 2) .or. (i == 4) .or. (i == 5)) then
1402- E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4_wp *G_L)
1403- E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4_wp *G_R)
1402+ E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4._wp *G_L)
1403+ E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4._wp *G_R)
14041404 end if
14051405 end if
14061406 end do
@@ -1410,7 +1410,7 @@ contains
14101410 xi_field_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, xibeg - 1 + i)
14111411 xi_field_R(i) = qR_prim_rs${XYZ}$_vf(j + 1, k, l, xibeg - 1 + i)
14121412 end do
1413- G_L = 0_wp ; G_R = 0_wp;
1413+ G_L = 0._wp ; G_R = 0._wp;
14141414 !$acc loop seq
14151415 do i = 1, num_fluids
14161416 ! Mixture left and right shear modulus
@@ -2498,7 +2498,7 @@ contains
24982498
24992499 ! ENERGY ADJUSTMENTS FOR HYPOELASTIC/HYPERELASTIC ENERGY
25002500 if (hypoelasticity .or. hyperelasticity) then
2501- G_L = 0_wp ; G_R = 0_wp
2501+ G_L = 0._wp ; G_R = 0._wp
25022502 !$acc loop seq
25032503 do i = 1, num_fluids
25042504 G_L = G_L + alpha_L(i)*Gs(i)
@@ -2514,12 +2514,12 @@ contains
25142514 do i = 1, strxe - strxb + 1
25152515 ! Elastic contribution to energy if G large enough
25162516 if ((G_L > verysmall) .and. (G_R > verysmall)) then
2517- E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4_wp *G_L)
2518- E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4_wp *G_R)
2517+ E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4._wp *G_L)
2518+ E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4._wp *G_R)
25192519 ! Additional terms in 2D and 3D
25202520 if ((i == 2) .or. (i == 4) .or. (i == 5)) then
2521- E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4_wp *G_L)
2522- E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4_wp *G_R)
2521+ E_L = E_L + (tau_e_L(i)*tau_e_L(i))/(4._wp *G_L)
2522+ E_R = E_R + (tau_e_R(i)*tau_e_R(i))/(4._wp *G_R)
25232523 end if
25242524 end if
25252525 end do
@@ -2529,7 +2529,7 @@ contains
25292529 xi_field_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, xibeg - 1 + i)
25302530 xi_field_R(i) = qR_prim_rs${XYZ}$_vf(j + 1, k, l, xibeg - 1 + i)
25312531 end do
2532- G_L = 0_wp ; G_R = 0_wp;
2532+ G_L = 0._wp ; G_R = 0._wp;
25332533 !$acc loop seq
25342534 do i = 1, num_fluids
25352535 ! Mixture left and right shear modulus
0 commit comments