Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/OBS/obs_conv_functions.h90

Timestamp:

2021-11-26T12:27:56+01:00 (3 years ago)

Author:

sparonuz

Message:

Mixed precision version, tested up to 30 years on ORCA2.

File:

• NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/OBS/obs_conv_functions.h90 (modified) (14 diffs)

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/OBS/obs_conv_functions.h90

@@ -5 +5 @@
    !!----------------------------------------------------------------------
 
-  REAL(KIND=wp) FUNCTION potemp( ps, pt, pp, ppr )
+  REAL(KIND=dp) FUNCTION potemp( ps, pt, pp, ppr )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION potemp  ***
@@ -32 +32 @@
       REAL(KIND=wp), INTENT(IN) :: ps
       REAL(KIND=wp), INTENT(IN) :: pt
-      REAL(KIND=wp), INTENT(IN) :: pp
+      REAL(KIND=dp), INTENT(IN) :: pp
       REAL(KIND=wp), INTENT(IN) :: ppr
 
       !! * Local declarations
-      REAL(KIND=wp) :: zpol
-      REAL(KIND=wp), PARAMETER :: a1 =  1.067610e-05
-      REAL(KIND=wp), PARAMETER :: a2 = -1.434297e-06
-      REAL(KIND=wp), PARAMETER :: a3 = -7.566349e-09
-      REAL(KIND=wp), PARAMETER :: a4 = -8.535585e-06
-      REAL(KIND=wp), PARAMETER :: a5 =  3.074672e-08
-      REAL(KIND=wp), PARAMETER :: a6 =  1.918639e-08
-      REAL(KIND=wp), PARAMETER :: a7 =  1.788718e-10
+      REAL(KIND=dp) :: zpol
+      REAL(KIND=dp), PARAMETER :: a1 =  1.067610e-05
+      REAL(KIND=dp), PARAMETER :: a2 = -1.434297e-06
+      REAL(KIND=dp), PARAMETER :: a3 = -7.566349e-09
+      REAL(KIND=dp), PARAMETER :: a4 = -8.535585e-06
+      REAL(KIND=dp), PARAMETER :: a5 =  3.074672e-08
+      REAL(KIND=dp), PARAMETER :: a6 =  1.918639e-08
+      REAL(KIND=dp), PARAMETER :: a7 =  1.788718e-10
 
       zpol = a1 + a2 * ps + a3 * ( pp + ppr ) + a4 * pt &
@@ -52 +52 @@
    END FUNCTION potemp
 
-   REAL(KIND=wp) FUNCTION fspott( pft, pfs, pfp )
+   REAL(KIND=dp) FUNCTION fspott( pft, pfs, pfp )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION fspott  ***
@@ -71 +71 @@
 
       !! * Arguments
-      REAL(KIND=wp) :: pft   ! in situ temperature in degrees Celsius
-      REAL(KIND=wp) :: pfs   ! salinity in psu
-      REAL(KIND=wp) :: pfp   ! pressure in bars
+      REAL(KIND=dp) :: pft   ! in situ temperature in degrees Celsius
+      REAL(KIND=dp) :: pfs   ! salinity in psu
+      REAL(KIND=dp) :: pfp   ! pressure in bars
       
       fspott = &
@@ -91 +91 @@
    END FUNCTION fspott
 
-   REAL(KIND=wp) FUNCTION atg( p_s, p_t, p_p )
+   REAL(KIND=dp) FUNCTION atg( p_s, p_t, p_p )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION atg  ***
@@ -115 +115 @@
 
       REAL(KIND=wp), INTENT(IN) :: p_s    ! Salinity in PSU
-      REAL(KIND=wp), INTENT(IN) :: p_t    ! Temperature in centigrades
-      REAL(KIND=wp), INTENT(IN) :: p_p    ! Pressure in decibars.
-
-      !! * Local declarations
-      
-      REAL(KIND=wp) :: z_ds
+      REAL(KIND=dp), INTENT(IN) :: p_t    ! Temperature in centigrades
+      REAL(KIND=dp), INTENT(IN) :: p_p    ! Pressure in decibars.
+
+      !! * Local declarations
+      
+      REAL(KIND=dp) :: z_ds
 
       z_ds = p_s - 35.0
@@ -131 +131 @@
    END FUNCTION atg
       
-   REAL(KIND=wp) FUNCTION theta( p_s, p_t0, p_p0, p_pr )
+   REAL(KIND=dp) FUNCTION theta( p_s, p_t0, p_p0, p_pr )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION theta  ***
@@ -159 +159 @@
 
       !! * Local declarations
-      REAL(KIND=wp) :: z_p
-      REAL(KIND=wp) :: z_t
-      REAL(KIND=wp) :: z_h
-      REAL(KIND=wp) :: z_xk
-      REAL(KIND=wp) :: z_q
+      REAL(KIND=dp) :: z_p
+      REAL(KIND=dp) :: z_t
+      REAL(KIND=dp) :: z_h
+      REAL(KIND=dp) :: z_xk
+      REAL(KIND=dp) :: z_q
 
       z_p = p_p0
@@ -184 +184 @@
    END FUNCTION theta
 
-   REAL(KIND=wp) FUNCTION depth( p_p, p_lat )
+   REAL(KIND=dp) FUNCTION depth( p_p, p_lat )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION depth  ***
@@ -210 +210 @@
 
       !! * Local declarations
-      REAL(KIND=wp) :: z_x
-      REAL(KIND=wp) :: z_gr
+      REAL(KIND=dp) :: z_x
+      REAL(KIND=dp) :: z_gr
       
       z_x = SIN( p_lat / 57.29578 )
@@ -221 +221 @@
    END FUNCTION depth
 
-   REAL(KIND=wp) FUNCTION p_to_dep( p_p, p_lat )
+   REAL(KIND=dp) FUNCTION p_to_dep( p_p, p_lat )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION p_to_dep  ***
@@ -245 +245 @@
 
       !! * Local declarations
-      REAL(KIND=wp) :: z_x
-      REAL(KIND=wp) :: z_c1
-      REAL(KIND=wp) :: z_c2
+      REAL(KIND=dp) :: z_x
+      REAL(KIND=dp) :: z_c1
+      REAL(KIND=dp) :: z_c2
 
       z_x = SIN( p_lat / 57.29578 )
@@ -257 +257 @@
    END FUNCTION p_to_dep
 
-   REAL(KIND=wp) FUNCTION dep_to_p( p_dep, p_lat )
+   REAL(KIND=dp) FUNCTION dep_to_p( p_dep, p_lat )
       !!----------------------------------------------------------------------
       !!                    ***  FUNCTION dep_to_p  ***
@@ -280 +280 @@
 
       !! * Local declarations
-      REAL(KIND=wp) :: z_x
-      REAL(KIND=wp) :: z_c1
-      REAL(KIND=wp) :: z_c2
-      REAL(KIND=wp) :: z_d
+      REAL(KIND=dp) :: z_x
+      REAL(KIND=dp) :: z_c1
+      REAL(KIND=dp) :: z_c2
+      REAL(KIND=dp) :: z_d
 
       z_x = SIN( p_lat / 57.29578 )