Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/OBS/grt_cir_dis.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/grt_cir_dis.h90 (modified) (2 diffs)

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

@@ -5 +5 @@
    !!----------------------------------------------------------------------
 
-   REAL(KIND=wp) FUNCTION grt_cir_dis( pa1, pa2, pb1, pb2, pc1, pc2 )
+   REAL(KIND=dp) FUNCTION grt_cir_dis( pa1, pa2, pb1, pb2, pc1, pc2 )
       !!----------------------------------------------------------------------
       !!                     *** FUNCTION grt_cir_dis ***
@@ -21 +21 @@
       
       !! * Arguments
-      REAL(KIND=wp) :: pa1   !  sin(lat1)
-      REAL(KIND=wp) :: pa2   !  sin(lat2)
-      REAL(KIND=wp) :: pb1   !  cos(lat1) * cos(lon1)
-      REAL(KIND=wp) :: pb2   !  cos(lat2) * cos(lon2)
-      REAL(KIND=wp) :: pc1   !  cos(lat1) * sin(lon1)
-      REAL(KIND=wp) :: pc2   !  cos(lat2) * sin(lon2)
+      REAL(KIND=dp) :: pa1   !  sin(lat1)
+      REAL(KIND=dp) :: pa2   !  sin(lat2)
+      REAL(KIND=dp) :: pb1   !  cos(lat1) * cos(lon1)
+      REAL(KIND=dp) :: pb2   !  cos(lat2) * cos(lon2)
+      REAL(KIND=dp) :: pc1   !  cos(lat1) * sin(lon1)
+      REAL(KIND=dp) :: pc2   !  cos(lat2) * sin(lon2)
 
-      REAL(KIND=wp) :: cosdist ! cosine of great circle distance
+      REAL(KIND=dp) :: cosdist ! cosine of great circle distance
 
       ! Compute cosine of great circle distance, constraining it to be between