Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/LBC/lib_mpp.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/LBC/lib_mpp.F90 (modified) (9 diffs)

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

@@ -197 +197 @@
    CHARACTER(len=3),  DIMENSION(nbdelay), PUBLIC ::   c_delaycpnt = (/ 'ICE'   , 'OCE' /)
    TYPE, PUBLIC ::   DELAYARR
-      REAL(   wp), POINTER, DIMENSION(:) ::  z1d => NULL()
+      REAL(   dp), POINTER, DIMENSION(:) ::  z1d => NULL()
       COMPLEX(dp), POINTER, DIMENSION(:) ::  y1d => NULL()
    END TYPE DELAYARR
@@ -207 +207 @@
    REAL(dp)              , PUBLIC ::  compute_time = 0._dp, elapsed_time = 0._dp
 
-   REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE ::   tampon   ! buffer in case of bsend
+   REAL(dp), DIMENSION(:), ALLOCATABLE, SAVE ::   tampon   ! buffer in case of bsend
 
    LOGICAL, PUBLIC ::   ln_nnogather                !: namelist control of northfold comms
@@ -290 +290 @@
       !!
       !!----------------------------------------------------------------------
-      REAL(wp), INTENT(inout) ::   pmess(*)   ! array of real
+      REAL(dp), INTENT(inout) ::   pmess(*)   ! array of real
       INTEGER , INTENT(in   ) ::   kbytes     ! size of the array pmess
       INTEGER , INTENT(in   ) ::   kdest      ! receive process number
@@ -365 +365 @@
       !!
       !!----------------------------------------------------------------------
-      REAL(wp), INTENT(inout) ::   pmess(*)   ! array of real
+      REAL(dp), INTENT(inout) ::   pmess(*)   ! array of real
       INTEGER , INTENT(in   ) ::   kbytes     ! suze of the array pmess
       INTEGER , INTENT(in   ) ::   ktyp       ! Tag of the recevied message
@@ -458 +458 @@
       !!
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION(jpi,jpj)      , INTENT(in   ) ::   ptab   ! subdomain input array
+      REAL(dp), DIMENSION(jpi,jpj)      , INTENT(in   ) ::   ptab   ! subdomain input array
       INTEGER                           , INTENT(in   ) ::   kp     ! record length
-      REAL(wp), DIMENSION(jpi,jpj,jpnij), INTENT(  out) ::   pio    ! subdomain input array
+      REAL(dp), DIMENSION(jpi,jpj,jpnij), INTENT(  out) ::   pio    ! subdomain input array
       !!
       INTEGER :: itaille, ierror   ! temporary integer
@@ -484 +484 @@
       !!
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION(jpi,jpj,jpnij)  ::   pio    ! output array
+      REAL(dp), DIMENSION(jpi,jpj,jpnij)  ::   pio    ! output array
       INTEGER                             ::   kp     ! Tag (not used with MPI
-      REAL(wp), DIMENSION(jpi,jpj)        ::   ptab   ! subdomain array input
+      REAL(dp), DIMENSION(jpi,jpj)        ::   ptab   ! subdomain array input
       !!
       INTEGER :: itaille, ierror   ! temporary integer
@@ -513 +513 @@
       CHARACTER(len=*), INTENT(in   )               ::   cdelay  ! name (used as id) of the delayed operation
       COMPLEX(dp),      INTENT(in   ), DIMENSION(:) ::   y_in
-      REAL(wp),         INTENT(  out), DIMENSION(:) ::   pout
+      REAL(dp),         INTENT(  out), DIMENSION(:) ::   pout
       LOGICAL,          INTENT(in   )               ::   ldlast  ! true if this is the last time we call this routine
       INTEGER,          INTENT(in   ), OPTIONAL     ::   kcom
@@ -586 +586 @@
       CHARACTER(len=*), INTENT(in   )                 ::   cdname  ! name of the calling subroutine
       CHARACTER(len=*), INTENT(in   )                 ::   cdelay  ! name (used as id) of the delayed operation
-      REAL(wp),         INTENT(in   ), DIMENSION(:)   ::   p_in    !
-      REAL(wp),         INTENT(  out), DIMENSION(:)   ::   pout    !
+      REAL(dp),         INTENT(in   ), DIMENSION(:)   ::   p_in    !
+      REAL(dp),         INTENT(  out), DIMENSION(:)   ::   pout    !
       LOGICAL,          INTENT(in   )                 ::   ldlast  ! true if this is the last time we call this routine
       INTEGER,          INTENT(in   ), OPTIONAL       ::   kcom
@@ -1234 +1234 @@
 
          ! The result is zt1 + zt2, after normalization.
-         yddb(ji) = cmplx ( zt1 + zt2, zt2 - ((zt1 + zt2) - zt1),wp )
+         yddb(ji) = cmplx ( zt1 + zt2, zt2 - ((zt1 + zt2) - zt1), dp )
       END DO
       !