Checked the 0_wp issues

Author: XZTian64

src/common/m_helper.fpp

@@ -488,7 +488,7 @@ contains
         real(wp) :: P
 
         if (l == 0) then
-            P = 1_wp
+            P = 1._wp
         else if (l == 1) then
             P = x
         else

src/post_process/m_data_output.fpp

@@ -1202,7 +1202,7 @@ contains
         allocate (x_d1(m*n))
         allocate (y_d1(m*n))
         counter = 0
-        maxalph_loc = 0_wp
+        maxalph_loc = 0._wp
         do k = 0, p
             do j = 0, n
                 do i = 0, m
@@ -1274,7 +1274,7 @@ contains
                         x_td(i), y_td(i), size(x_td)
                 else
                     write (211, '(F12.9,1X,F12.9,1X,F3.1)') &
-                        x_td(i), y_td(i), 0_wp
+                        x_td(i), y_td(i), 0._wp
                 end if
             end do
         end if
@@ -1289,29 +1289,29 @@ contains
         real(wp), dimension(num_fluids) :: adv
         integer :: i, j, k, l, s !looping indices
 
-        Egk = 0_wp
-        Elp = 0_wp
-        Egint = 0_wp
-        Vb = 0_wp
-        maxvel = 0_wp
-        MaxMa = 0_wp
-        Vl = 0_wp
-        Elk = 0_wp
-        Et = 0_wp
-        Vb = 0_wp
-        dV = 0_wp
-        pres_av = 0_wp
-        pres = 0_wp
+        Egk = 0._wp
+        Elp = 0._wp
+        Egint = 0._wp
+        Vb = 0._wp
+        maxvel = 0._wp
+        MaxMa = 0._wp
+        Vl = 0._wp
+        Elk = 0._wp
+        Et = 0._wp
+        Vb = 0._wp
+        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
-                    pres = 0_wp
+                    pres = 0._wp
                     dV = dx(i)*dy(j)*dz(k)
-                    rho = 0_wp
-                    gamma = 0_wp
-                    pi_inf = 0_wp
+                    rho = 0._wp
+                    gamma = 0._wp
+                    pi_inf = 0._wp
                     pres = q_prim_vf(E_idx)%sf(i, j, k)
                     Egint = Egint + q_prim_vf(E_idx + 2)%sf(i, j, k)*(fluid_pp(2)%gamma*pres)*dV
                     do s = 1, num_vels

src/pre_process/m_assign_variables.fpp

@@ -502,7 +502,7 @@ contains
                 theta = atan2(y_cc(k), x_cc(j))
                 phi = atan2(sqrt(x_cc(j)**2 + y_cc(k)**2), z_cc(l))
                 !spherical coord, assuming Rmax=1
-                xi_sph = (rcoord**3 - R0ref**3 + 1_wp)**(1._wp/3._wp)
+                xi_sph = (rcoord**3 - R0ref**3 + 1._wp)**(1._wp/3._wp)
                 xi_cart(1) = xi_sph*sin(phi)*cos(theta)
                 xi_cart(2) = xi_sph*sin(phi)*sin(theta)
                 xi_cart(3) = xi_sph*cos(phi)
@@ -515,7 +515,7 @@ contains
             ! assigning the reference map to the q_prim vector field
             do i = 1, num_dims
                 q_prim_vf(i + xibeg - 1)%sf(j, k, l) = eta*xi_cart(i) + &
-                                                       (1_wp - eta)*orig_prim_vf(i + xibeg - 1)
+                                                       (1._wp - eta)*orig_prim_vf(i + xibeg - 1)
             end do
         end if
 

src/pre_process/m_patches.fpp

@@ -1648,7 +1648,7 @@ contains
 
                     if (non_axis_sym) then
                         phi = atan(((y_cc(j) - y_centroid) + eps)/((x_cc(i) - x_centroid) + eps))
-                        r = sqrt((x_cc(i) - x_centroid)**2_wp + (y_cc(j) - y_centroid)**2_wp) + eps
+                        r = sqrt((x_cc(i) - x_centroid)**2._wp + (y_cc(j) - y_centroid)**2._wp) + eps
                         x_p = (eps)/r
                         Ps(2) = spherical_harmonic_func(x_p, phi, 2, 2)
                         Ps(3) = spherical_harmonic_func(x_p, phi, 3, 3)
@@ -1659,7 +1659,7 @@ contains
                         Ps(8) = spherical_harmonic_func(x_p, phi, 8, 8)
                         Ps(9) = spherical_harmonic_func(x_p, phi, 9, 9)
                     else
-                        r = sqrt((x_cc(i) - x_centroid)**2_wp + (y_cc(j) - y_centroid)**2_wp) + eps
+                        r = sqrt((x_cc(i) - x_centroid)**2._wp + (y_cc(j) - y_centroid)**2._wp) + eps
                         x_p = abs(x_cc(i) - x_centroid + eps)/r
                         Ps(2) = unassociated_legendre(x_p, 2)
                         Ps(3) = unassociated_legendre(x_p, 3)

src/simulation/m_hyperelastic.fpp

@@ -126,7 +126,7 @@ contains
                     if (G > verysmall) then
                         !$acc loop seq
                         do i = 1, tensor_size
-                            tensora(i) = 0_wp
+                            tensora(i) = 0._wp
                         end do
                         ! STEP 1: computing the grad_xi tensor using finite differences
                         ! grad_xi definition / organization
@@ -173,7 +173,7 @@ contains
                             end do
 
                             ! STEP 2d: computing the J = det(F) = 1/det(\grad{\xi})
-                            tensorb(tensor_size) = 1_wp/tensorb(tensor_size)
+                            tensorb(tensor_size) = 1._wp/tensorb(tensor_size)
 
                             ! STEP 3: computing F transpose F
                             tensorb(1) = tensora(1)**2 + tensora(2)**2 + tensora(3)**2

src/simulation/m_hypoelastic.fpp

@@ -165,8 +165,8 @@ contains
                     do q = 0, p
                         do l = 0, n
                             do k = 0, m
-                                du_dz(k, l, q) = 0_wp; dv_dz(k, l, q) = 0_wp; dw_dx(k, l, q) = 0_wp; 
-                                dw_dy(k, l, q) = 0_wp; dw_dz(k, l, q) = 0_wp; 
+                                du_dz(k, l, q) = 0._wp; dv_dz(k, l, q) = 0._wp; dw_dx(k, l, q) = 0._wp; 
+                                dw_dy(k, l, q) = 0._wp; dw_dz(k, l, q) = 0._wp; 
                             end do
                         end do
                     end do