Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/lib_fortran.F90

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/lib_fortran.F90 (modified) (6 diffs)

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/lib_fortran.F90

@@ -91 +91 @@
    FUNCTION local_sum_2d( ptab )
       !!----------------------------------------------------------------------
-      REAL(wp),  INTENT(in   ) ::   ptab(:,:) ! array on which operation is applied
+      REAL(dp),  INTENT(in   ) ::   ptab(:,:) ! array on which operation is applied
       COMPLEX(dp)              ::  local_sum_2d
       !
@@ -97 +97 @@
       !
       COMPLEX(dp)::   ctmp
-      REAL(wp)   ::   ztmp
+      REAL(dp)   ::   ztmp
       INTEGER    ::   ji, jj    ! dummy loop indices
       INTEGER    ::   ipi, ipj  ! dimensions
@@ -126 +126 @@
       !
       COMPLEX(dp)::   ctmp
-      REAL(wp)   ::   ztmp
+      REAL(dp)   ::   ztmp
       INTEGER    ::   ji, jj, jk   ! dummy loop indices
       INTEGER    ::   ipi, ipj, ipk    ! dimensions
@@ -158 +158 @@
       !! ** Purpose : sum over 3x3 boxes
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION (:,:), INTENT(inout) ::   p2d
+      REAL(dp), DIMENSION (:,:), INTENT(inout) ::   p2d
       !
       INTEGER ::   ji, ji2, jj, jj2     ! dummy loop indices
@@ -210 +210 @@
       !! ** Purpose : sum over 3x3 boxes
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION (:,:,:), INTENT(inout) ::   p3d
+      REAL(dp), DIMENSION (:,:,:), INTENT(inout) ::   p3d
       !
       INTEGER ::   ji, ji2, jj, jj2, jn     ! dummy loop indices
@@ -292 +292 @@
       !
       ! The result is t1 + t2, after normalization.
-      yddb = CMPLX( zt1 + zt2, zt2 - ((zt1 + zt2) - zt1), wp )
+      yddb = CMPLX( zt1 + zt2, zt2 - ((zt1 + zt2) - zt1), dp )
       !
    END SUBROUTINE DDPDD