From ecb6f0de41b46fffb53cb24b7d686bf511217eb8 Mon Sep 17 00:00:00 2001
From: Juan ESCOBAR <juan.escobar@aero.obs-mip.fr>
Date: Tue, 14 Nov 2023 19:10:57 +0100
Subject: [PATCH] Juan 14/11/2023:discretisation.f90 , for real*4 compilation ,
replace DSQRT -> SQRT
---
.../discretisation.f90 | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/src/LIB/tensorproductmultigrid_Source/discretisation.f90 b/src/LIB/tensorproductmultigrid_Source/discretisation.f90
index 06352e12b..0952a99a3 100644
--- a/src/LIB/tensorproductmultigrid_Source/discretisation.f90
+++ b/src/LIB/tensorproductmultigrid_Source/discretisation.f90
@@ -503,27 +503,27 @@ contains
sigma_j = 2.0_rl*(1.0_rl*iy-0.5_rl)/grid_param%n-1.0_rl
! alpha_{i+1,j}
if (ix == grid_param%n) then
- alpha_T(1) = 2.0_rl*DSQRT((1.0_rl+(rho_i+0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
+ alpha_T(1) = 2.0_rl*SQRT((1.0_rl+(rho_i+0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
else
- alpha_T(1) = DSQRT((1.0_rl+(rho_i+0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
+ alpha_T(1) = SQRT((1.0_rl+(rho_i+0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
end if
! alpha_{i-1,j}
if (ix == 1) then
- alpha_T(2) = 2.0_rl*DSQRT((1.0_rl+(rho_i-0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
+ alpha_T(2) = 2.0_rl*SQRT((1.0_rl+(rho_i-0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
else
- alpha_T(2) = DSQRT((1.0_rl+(rho_i-0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
+ alpha_T(2) = SQRT((1.0_rl+(rho_i-0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
end if
! alpha_{i,j+1}
if (iy == grid_param%n) then
- alpha_T(3) = 2.0_rl*DSQRT((1.0_rl+(sigma_j+0.25_rl*h)**2)/(1.0_rl+rho_i**2))
+ alpha_T(3) = 2.0_rl*SQRT((1.0_rl+(sigma_j+0.25_rl*h)**2)/(1.0_rl+rho_i**2))
else
- alpha_T(3) = DSQRT((1.0_rl+(sigma_j+0.5_rl*h)**2)/(1.0_rl+rho_i**2))
+ alpha_T(3) = SQRT((1.0_rl+(sigma_j+0.5_rl*h)**2)/(1.0_rl+rho_i**2))
end if
! alpha_{i,j-1}
if (iy == 1) then
- alpha_T(4) = 2.0_rl*DSQRT((1.0_rl+(sigma_j-0.25_rl*h)**2)/(1.0_rl+rho_i**2))
+ alpha_T(4) = 2.0_rl*SQRT((1.0_rl+(sigma_j-0.25_rl*h)**2)/(1.0_rl+rho_i**2))
else
- alpha_T(4) = DSQRT((1.0_rl+(sigma_j-0.5_rl*h)**2)/(1.0_rl+rho_i**2))
+ alpha_T(4) = SQRT((1.0_rl+(sigma_j-0.5_rl*h)**2)/(1.0_rl+rho_i**2))
end if
#endif
alpha_T(5) = alpha_T(1) + alpha_T(2) + alpha_T(3) + alpha_T(4)
--
GitLab