fix

Author: sbryngelson

src/common/m_helper.fpp

@@ -503,12 +503,12 @@ contains
         real(wp), intent(in) :: x, phi
         real(wp) :: Y, prefactor, pi
 
-        pi = acos(-1_wp)
+        pi = acos(-1._wp)
         prefactor = sqrt((2*l + 1)/(4*pi)*factorial(l - m)/factorial(l + m)); 
         if (m == 0) then
             Y = prefactor*associated_legendre(x, l, m); 
         elseif (m > 0) then
-            Y = (-1_wp)**m*sqrt(2_wp)*prefactor*associated_legendre(x, l, m)*cos(m*phi); 
+            Y = (-1._wp)**m*sqrt(2._wp)*prefactor*associated_legendre(x, l, m)*cos(m*phi); 
         end if
 
     end function spherical_harmonic_func
@@ -530,7 +530,7 @@ contains
         elseif (l == 1 .and. m <= 0) then
             P = x; 
         elseif (l == 1 .and. m == 1) then
-            P = -(1 - x**2)**(1/2); 
+            P = -(1 - x**2)**(1._wp/2._wp); 
         elseif (m == l) then
             P = (-1)**l*double_factorial(2*l - 1)*(1 - x**2)**(l/2); 
         elseif (m == l - 1) then

src/common/m_variables_conversion.fpp

@@ -1416,6 +1416,7 @@ contains
 
         real(wp) :: blkmod1, blkmod2
         real(wp) :: Tolerance
+
         integer :: q
 
         if (chemistry) then

src/post_process/m_data_output.fpp

@@ -1139,6 +1139,8 @@ contains
         dV = 0_wp
         pres_av = 0_wp
         pres = 0_wp
+        c = 0._wp
+
         do k = 0, p
             do j = 0, n
                 do i = 0, m
@@ -1168,10 +1170,10 @@ contains
 
                     call s_compute_speed_of_sound(pres, rho, &
                                                   gamma, pi_inf, &
-                                                  H, adv, 0_wp, 0_wp, c)
+                                                  H, adv, 0._wp, 0._wp, c)
 
                     Ma = maxvel/c
-                    if (Ma > MaxMa .and. adv(1) > 1.0_wp - 1.0d-10) then
+                    if (Ma > MaxMa .and. (adv(1) > (1.0_wp - 1.0e-10_wp))) then
                         MaxMa = Ma
                     end if
                     Vl = Vl + adv(1)*dV

src/simulation/m_riemann_solvers.fpp

@@ -572,7 +572,7 @@ contains
                             !        G_R = G_R + alpha_R(i)*Gs(i)
                             !    end do
                             !    ! Elastic contribution to energy if G large enough
-                            !    if ((G_L > 1d-3) .and. (G_R > 1d-3)) then
+                            !    if ((G_L > 1e-3_wp) .and. (G_R > 1e-3_wp)) then
                             !    E_L = E_L + G_L*qL_prim_rs${XYZ}$_vf(j, k, l, xiend + 1)
                             !    E_R = E_R + G_R*qR_prim_rs${XYZ}$_vf(j + 1, k, l, xiend + 1)
                             !    !$acc loop seq
@@ -1299,13 +1299,13 @@ contains
 
                                 ! goes with numerical star velocity in x/y/z directions
                                 ! xi_P/M = 0.5 +/m sgn(0.5,s_star)
-                                xi_M = (5d-1 + sign(5d-1, s_S))
-                                xi_P = (5d-1 - sign(5d-1, s_S))
+                                xi_M = (5e-1_wp + sign(0.5_wp, s_S))
+                                xi_P = (5e-1_wp - sign(0.5_wp, s_S))
 
                                 ! goes with the numerical velocity in x/y/z directions
                                 ! xi_P/M (pressure) = min/max(0. sgn(1,sL/sR))
-                                xi_MP = -min(0_wp, sign(1_wp, s_L))
-                                xi_PP = max(0_wp, sign(1_wp, s_R))
+                                xi_MP = -min(0._wp, sign(1._wp, s_L))
+                                xi_PP = max(0._wp, sign(1._wp, s_R))
 
                                 E_star = xi_M*(E_L + xi_MP*(xi_L*(E_L + (s_S - vel_L(dir_idx(1)))* &
                                                                   (rho_L*s_S + pres_L/(s_L - vel_L(dir_idx(1))))) - E_L)) + &