fix

Author: Archith Iyer

misc/m_silo_proxy.f90

@@ -189,9 +189,9 @@ function DBPUTQM(dbid, name, lname, xname, lxname, yname, lyname, & !! -
         integer, intent(IN) :: lyname
         character(LEN=*), intent(IN) :: zname
         integer, intent(IN) :: lzname
-        real(kind(0d0)), dimension(:), intent(IN) :: x
-        real(kind(0d0)), dimension(:), intent(IN) :: y
-        real(kind(0d0)), dimension(:), intent(IN) :: z
+        real(wp), dimension(:), intent(IN) :: x
+        real(wp), dimension(:), intent(IN) :: y
+        real(wp), dimension(:), intent(IN) :: z
         integer, dimension(:), intent(IN) :: dims
         integer, intent(IN) :: ndims
         integer, intent(IN) :: datatype
@@ -215,8 +215,8 @@ function DBPUTCURVE(dbid, curvename, lcurvename, xvals, yvals, & !! ----
         integer, intent(IN) :: dbid
         character(LEN=*), intent(IN) :: curvename
         integer, intent(IN) :: lcurvename
-        real(kind(0d0)), dimension(:), intent(IN) :: xvals
-        real(kind(0d0)), dimension(:), intent(IN) :: yvals
+        real(wp), dimension(:), intent(IN) :: xvals
+        real(wp), dimension(:), intent(IN) :: yvals
         integer, intent(IN) :: datatype
         integer, intent(IN) :: npoints
         integer, intent(IN) :: optlist_id
@@ -264,7 +264,7 @@ function DBPUTQV1(dbid, name, lname, meshname, lmeshname, var, & !! ----
         integer, intent(IN) :: lname
         character(LEN=*), intent(IN) :: meshname
         integer, intent(IN) :: lmeshname
-        real(kind(0d0)), dimension(:, :, :), intent(IN) :: var
+        real(wp), dimension(:, :, :), intent(IN) :: var
         integer, dimension(:), intent(IN) :: dims
         integer, intent(IN) :: ndims
         integer, intent(IN) :: mixvar

src/common/m_constants.fpp

@@ -35,10 +35,10 @@ module m_constants
     integer, parameter :: Ifactor_3D = 5 !< Multiple factor of the ratio (edge to cell width) for interpolation along edges for 3D models
     integer, parameter :: Ifactor_bary_3D = 20 !< Multiple factor of the ratio (triangle area to cell face area) for interpolation on triangle facets for 3D models
     integer, parameter :: num_ray = 20 !< Default number of rays traced per cell
-    real(kind(0d0)), parameter :: ray_tracing_threshold = 0.9d0 !< Threshold above which the cell is marked as the model patch
-    real(kind(0d0)), parameter :: threshold_vector_zero = 1d-10 !< Threshold to treat the component of a vector to be zero
-    real(kind(0d0)), parameter :: threshold_edge_zero = 1d-10 !< Threshold to treat two edges to be overlapped
-    real(kind(0d0)), parameter :: threshold_bary = 1d-1 !< Threshold to interpolate a barycentric facet
-    real(kind(0d0)), parameter :: initial_distance_buffer = 1d12 !< Initialized levelset distance for the shortest path pair algorithm
+    real(wp), parameter :: ray_tracing_threshold = 0.9d0 !< Threshold above which the cell is marked as the model patch
+    real(wp), parameter :: threshold_vector_zero = 1d-10 !< Threshold to treat the component of a vector to be zero
+    real(wp), parameter :: threshold_edge_zero = 1d-10 !< Threshold to treat two edges to be overlapped
+    real(wp), parameter :: threshold_bary = 1d-1 !< Threshold to interpolate a barycentric facet
+    real(wp), parameter :: initial_distance_buffer = 1d12 !< Initialized levelset distance for the shortest path pair algorithm
 
 end module m_constants

src/common/m_derived_types.fpp

@@ -168,7 +168,6 @@ module m_derived_types
         !! patch geometries. It is specified through its x-, y-, and z-components
         !! respectively.
 
-
         real(wp) :: epsilon, beta !<
         !! The spherical harmonics eccentricity parameters.
 
@@ -239,7 +238,7 @@ module m_derived_types
         integer :: model_spc !<
         !! Number of samples per cell to use when discretizing the STL object.
 
-        real(kind(0d0)) :: model_threshold !<
+        real(wp) :: model_threshold !<
         !! Threshold to turn on smoothen STL patch.
 
     end type ic_patch_parameters
@@ -277,7 +276,7 @@ module m_derived_types
         integer :: model_spc !<
         !! Number of samples per cell to use when discretizing the STL object.
 
-        real(kind(0d0)) :: model_threshold !<
+        real(wp) :: model_threshold !<
         !! Threshold to turn on smoothen STL patch.
     end type ib_patch_parameters
 

src/common/m_helper_basic.f90

@@ -25,6 +25,7 @@ module m_helper_basic
     !! @param tol_input Relative error (default = 1e-6_wp).
     !! @return Result of the comparison.
     logical function f_approx_equal(a, b, tol_input) result(res)
+        !$acc routine seq
         real(wp), intent(in) :: a, b
         real(wp), optional, intent(in) :: tol_input
         real(wp) :: tol
@@ -47,6 +48,7 @@ end function f_approx_equal
     !> Checks if a real(wp) variable is of default value.
     !! @param var Variable to check.
     logical function f_is_default(var) result(res)
+        !$acc routine seq
         real(wp), intent(in) :: var
 
         res = f_approx_equal(var, dflt_real)

src/pre_process/m_compute_levelset.fpp

@@ -360,9 +360,9 @@ contains
         type(levelset_norm_field), intent(INOUT) :: levelset_norm
 
         integer, intent(IN) :: ib_patch_id
-        real(kind(0d0)) :: Right, Left, Bottom, Top, Front, Back
-        real(kind(0d0)) :: x, y, z, min_dist
-        real(kind(0d0)) :: side_dists(6)
+        real(wp) :: Right, Left, Bottom, Top, Front, Back
+        real(wp) :: x, y, z, min_dist
+        real(wp) :: side_dists(6)
 
         integer :: i, j, k !< Loop index variables