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

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

@@ -34 +34 @@
       INTEGER, INTENT(IN), DIMENSION(kdep) :: &
          & kkco                 ! Array indicies for interpolation
-      REAL(KIND=wp), INTENT(IN), DIMENSION(kdep) :: &
+      REAL(KIND=dp), INTENT(IN), DIMENSION(kdep) :: &
          & pobsdep              ! Depth of the observation
-      REAL(KIND=wp), INTENT(IN), DIMENSION(kpk) :: &
+      REAL(KIND=dp), INTENT(IN), DIMENSION(kpk) :: &
          & pobsk,  &            ! Model profile at a given (lon,lat)
          & pobs2k, &            ! 2nd derivative of the interpolating function
          & pdep,   &            ! Model depth array
          & pobsmask             ! Vertical mask
-      REAL(KIND=wp), INTENT(OUT), DIMENSION(kdep) :: &
+      REAL(KIND=dp), INTENT(OUT), DIMENSION(kdep) :: &
          & pobs                 ! Model equivalent at observation point
   
       !! * Local declarations
-      REAL(KIND=wp) :: z1dm       ! Distance above and below obs to model grid points
-      REAL(KIND=wp) :: z1dp         
-      REAL(KIND=wp) :: zsum       ! Dummy variables for computation
-      REAL(KIND=wp) :: zsum2
+      REAL(KIND=dp) :: z1dm       ! Distance above and below obs to model grid points
+      REAL(KIND=dp) :: z1dp         
+      REAL(KIND=dp) :: zsum       ! Dummy variables for computation
+      REAL(KIND=dp) :: zsum2
       INTEGER :: jdep             ! Observation depths loop variable
     
@@ -122 +122 @@
       !! * Arguments
       INTEGER, INTENT(IN) :: kpk               ! Number of vertical levels
-      REAL(KIND=wp), INTENT(IN), DIMENSION(kpk) :: &
+      REAL(KIND=dp), INTENT(IN), DIMENSION(kpk) :: &
          & pobsk, &          ! Model profile at a given (lon,lat)
          & pdep,  &          ! Model depth array
          & pobsmask          ! Vertical mask
-      REAL(KIND=wp), INTENT(OUT), DIMENSION(kpk) :: &
+      REAL(KIND=dp), INTENT(OUT), DIMENSION(kpk) :: &
          & pobs2k            ! 2nd derivative of the interpolating function
   
       !! * Local declarations
       INTEGER :: jk
-      REAL(KIND=wp) :: za
-      REAL(KIND=wp) :: zb
-      REAL(KIND=wp) :: zc
-      REAL(KIND=wp) :: zpa
-      REAL(KIND=wp) :: zkm
-      REAL(KIND=wp) :: zkp
-      REAL(KIND=wp) :: zk
-      REAL(KIND=wp), DIMENSION(kpk-1) :: &
+      REAL(KIND=dp) :: za
+      REAL(KIND=dp) :: zb
+      REAL(KIND=dp) :: zc
+      REAL(KIND=dp) :: zpa
+      REAL(KIND=dp) :: zkm
+      REAL(KIND=dp) :: zkp
+      REAL(KIND=dp) :: zk
+      REAL(KIND=dp), DIMENSION(kpk-1) :: &
          & zs, &
          & zp, &