fix some missing ones

Author: sbryngelson

src/common/m_derived_types.fpp

@@ -188,7 +188,7 @@ module m_derived_types
         !! Identity (id) of the patch with which current patch is to get smoothed
 
         real(wp) :: smooth_coeff !<
-        !! Smoothing coefficient (coeff) adminstrating the size of the stencil of
+        !! Smoothing coefficient (coeff) for the size of the stencil of
         !! cells across which boundaries of the current patch will be smeared out
 
         real(wp), dimension(num_fluids_max) :: alpha_rho

src/common/m_helper.fpp

@@ -229,16 +229,16 @@ contains
         real(wp), intent(in) :: omega, peclet
         real(wp), intent(out) :: Re_trans, Im_trans
 
-        complex :: trans, c1, c2, c3
-        complex :: imag = (0._wp, 1._wp)
+        complex(wp) :: trans, c1, c2, c3
+        complex(wp) :: imag = (0._wp, 1._wp)
         real(wp) :: f_transcoeff
 
         c1 = imag*omega*peclet
         c2 = CSQRT(c1)
         c3 = (CEXP(c2) - CEXP(-c2))/(CEXP(c2) + CEXP(-c2)) ! TANH(c2)
         trans = ((c2/c3 - 1._wp)**(-1) - 3._wp/c1)**(-1) ! transfer function
 
-        Re_trans = dble(trans)
+        Re_trans = trans
         Im_trans = aimag(trans)
 
     end subroutine s_transcoeff
@@ -279,7 +279,7 @@ contains
         ! phi = ln( R0 ) & return R0
         do ir = 1, nb
             phi(ir) = log(R0mn) &
-                      + dble(ir - 1)*log(R0mx/R0mn)/dble(nb - 1)
+                      + (ir - 1._wp)*log(R0mx/R0mn)/(nb - 1._wp)
             R0(ir) = exp(phi(ir))
         end do
         dphi = phi(2) - phi(1)

src/common/m_precision_select.f90

@@ -19,7 +19,7 @@ module m_precision_select
     integer, parameter :: sp = single_precision
     integer, parameter :: dp = double_precision
 
-    ! Set the working precision (wp) to single or double precision
+    ! Set the working precision (wp) to single or double
 #ifdef MFC_SINGLE_PRECISION
     integer, parameter :: wp = single_precision  ! Change to single_precision if needed
 #else

src/simulation/m_fftw.fpp

@@ -51,9 +51,9 @@ module m_fftw
 #if defined(MFC_OpenACC)
     !$acc declare create(real_size, cmplx_size, x_size, batch_size, Nfq)
 
-    real(kind(0d0)), allocatable, target :: data_real_gpu(:)
-    complex(kind(0d0)), allocatable, target :: data_cmplx_gpu(:)
-    complex(kind(0d0)), allocatable, target :: data_fltr_cmplx_gpu(:)
+    real(dp), allocatable, target :: data_real_gpu(:)
+    complex(dp), allocatable, target :: data_cmplx_gpu(:)
+    complex(dp), allocatable, target :: data_fltr_cmplx_gpu(:)
 !$acc declare create(data_real_gpu, data_cmplx_gpu, data_fltr_cmplx_gpu)
 
 #if defined(__PGI)
@@ -146,7 +146,7 @@ contains
         do k = 1, sys_size
             do j = 0, m
                 do l = 1, cmplx_size
-                    data_fltr_cmplx_gpu(l + j*cmplx_size + (k - 1)*cmplx_size*x_size) = (0d0, 0d0)
+                    data_fltr_cmplx_gpu(l + j*cmplx_size + (k - 1)*cmplx_size*x_size) = (0_dp, 0_dp)
                 end do
             end do
         end do
@@ -198,7 +198,7 @@ contains
         do k = 1, sys_size
             do j = 0, m
                 do l = 0, p
-                    data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size) = data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size)/real(real_size, kind(0d0))
+                    data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size) = data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size)/real(real_size, dp)
                     q_cons_vf(k)%sf(j, 0, l) = data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size)
                 end do
             end do
@@ -210,7 +210,7 @@ contains
             do k = 1, sys_size
                 do j = 0, m
                     do l = 1, cmplx_size
-                        data_fltr_cmplx_gpu(l + j*cmplx_size + (k - 1)*cmplx_size*x_size) = (0d0, 0d0)
+                        data_fltr_cmplx_gpu(l + j*cmplx_size + (k - 1)*cmplx_size*x_size) = (0_dp, 0_dp)
                     end do
                 end do
             end do
@@ -233,7 +233,7 @@ contains
 #endif
             !$acc end host_data
 
-            Nfq = min(floor(2d0*real(i, kind(0d0))*pi), cmplx_size)
+            Nfq = min(floor(2_dp*real(i, dp)*pi), cmplx_size)
             !$acc update device(Nfq)
 
             !$acc parallel loop collapse(3) gang vector default(present)
@@ -258,7 +258,7 @@ contains
             do k = 1, sys_size
                 do j = 0, m
                     do l = 0, p
-                        data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size) = data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size)/real(real_size, kind(0d0))
+                        data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size) = data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size)/real(real_size, dp)
                         q_cons_vf(k)%sf(j, i, l) = data_real_gpu(l + j*real_size + 1 + (k - 1)*real_size*x_size)
                     end do
                 end do
@@ -270,27 +270,27 @@ contains
         Nfq = 3
         do j = 0, m
             do k = 1, sys_size
-                data_fltr_cmplx(:) = (0d0, 0d0)
+                data_fltr_cmplx(:) = (0_dp, 0_dp)
                 data_real(1:p + 1) = q_cons_vf(k)%sf(j, 0, 0:p)
                 call fftw_execute_dft_r2c(fwd_plan, data_real, data_cmplx)
                 data_fltr_cmplx(1:Nfq) = data_cmplx(1:Nfq)
                 call fftw_execute_dft_c2r(bwd_plan, data_fltr_cmplx, data_real)
-                data_real(:) = data_real(:)/real(real_size, kind(0d0))
+                data_real(:) = data_real(:)/real(real_size, dp)
                 q_cons_vf(k)%sf(j, 0, 0:p) = data_real(1:p + 1)
             end do
         end do
 
         ! Apply Fourier filter to additional rings
         do i = 1, fourier_rings
-            Nfq = min(floor(2d0*real(i, kind(0d0))*pi), cmplx_size)
+            Nfq = min(floor(2_dp*real(i, dp)*pi), cmplx_size)
             do j = 0, m
                 do k = 1, sys_size
-                    data_fltr_cmplx(:) = (0d0, 0d0)
+                    data_fltr_cmplx(:) = (0_dp, 0_dp)
                     data_real(1:p + 1) = q_cons_vf(k)%sf(j, i, 0:p)
                     call fftw_execute_dft_r2c(fwd_plan, data_real, data_cmplx)
                     data_fltr_cmplx(1:Nfq) = data_cmplx(1:Nfq)
                     call fftw_execute_dft_c2r(bwd_plan, data_fltr_cmplx, data_real)
-                    data_real(:) = data_real(:)/real(real_size, kind(0d0))
+                    data_real(:) = data_real(:)/real(real_size, dp)
                     q_cons_vf(k)%sf(j, i, 0:p) = data_real(1:p + 1)
                 end do
             end do

src/simulation/m_rhs.fpp

@@ -899,7 +899,7 @@ contains
 
         if (chemistry) then
             !$acc kernels
-            rhs_vf(T_idx)%sf(:, :, :) = 0d0
+            rhs_vf(T_idx)%sf(:, :, :) = 0.0_wp
             !$acc end kernels
 
             if (chem_params%reactions) then

src/syscheck/syscheck.fpp

@@ -53,7 +53,7 @@ program syscheck
 
     @:ACC(integer(acc_device_kind) :: devtype)
     @:ACC(integer :: i, num_devices)
-    @:ACC(real(kind(0d0)), allocatable, dimension(:) :: arr)
+    @:ACC(real(8), allocatable, dimension(:) :: arr)
     @:ACC(integer, parameter :: N = 100)
 
     @:MPIC(call mpi_init(ierr))