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

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

@@ -74 +74 @@
       CHARACTER(len = 8)    ::   clname    ! Name of the coupling field
       CHARACTER(len = 1)    ::   clgrid    ! Grid type
-      REAL(wp)              ::   nsgn      ! Control of the sign change
+      REAL(dp)              ::   nsgn      ! Control of the sign change
       INTEGER, DIMENSION(nmaxcat,nmaxcpl) ::   nid   ! Id of the field (no more than 9 categories and 9 extrena models)
       INTEGER               ::   nct       ! Number of categories in field
@@ -82 +82 @@
    TYPE(FLD_CPL), DIMENSION(nmaxfld), PUBLIC ::   srcv, ssnd   !: Coupling fields
 
-   REAL(wp), DIMENSION(:,:), ALLOCATABLE ::   exfld   ! Temporary buffer for receiving
+   REAL(dp), DIMENSION(:,:), ALLOCATABLE ::   exfld   ! Temporary buffer for receiving
 
    !!----------------------------------------------------------------------
@@ -89 +89 @@
    !! Software governed by the CeCILL license (see ./LICENSE)
    !!----------------------------------------------------------------------
+#  include "single_precision_substitute.h90"
 CONTAINS
 
@@ -316 +317 @@
       INTEGER                   , INTENT(  out) ::   kinfo     ! OASIS3 info argument
       INTEGER                   , INTENT(in   ) ::   kstep     ! ocean time-step in seconds
-      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   pdata
+      REAL(dp), DIMENSION(:,:,:), INTENT(in   ) ::   pdata
       !!
       INTEGER                                   ::   jc,jm     ! local loop index
@@ -361 +362 @@
       INTEGER                   , INTENT(in   ) ::   kid       ! variable index in the array
       INTEGER                   , INTENT(in   ) ::   kstep     ! ocean time-step in seconds
-      REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   pdata     ! IN to keep the value if nothing is done
-      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   pmask     ! coupling mask
+      REAL(dp), DIMENSION(:,:,:), INTENT(inout) ::   pdata     ! IN to keep the value if nothing is done
+      REAL(dp), DIMENSION(:,:,:), INTENT(in   ) ::   pmask     ! coupling mask
       INTEGER                   , INTENT(  out) ::   kinfo     ! OASIS3 info argument
       !!
@@ -419 +420 @@
          !--- we must call lbc_lnk to fill the halos that where not received.
          IF( .NOT. ll_1st ) THEN
-            CALL lbc_lnk( 'cpl_oasis3', pdata(:,:,jc), srcv(kid)%clgrid, srcv(kid)%nsgn )
+            CALL lbc_lnk( 'cpl_oasis3', pdata(:,:,jc), srcv(kid)%clgrid, CASTSP(srcv(kid)%nsgn) )
          ENDIF
 
@@ -547 +548 @@
 
    SUBROUTINE oasis_put(k1,k2,p1,k3)
-      REAL(wp), DIMENSION(:,:), INTENT(in   ) ::  p1
+      REAL(dp), DIMENSION(:,:), INTENT(in   ) ::  p1
       INTEGER                 , INTENT(in   ) ::  k1,k2
       INTEGER                 , INTENT(  out) ::  k3
@@ -555 +556 @@
 
    SUBROUTINE oasis_get(k1,k2,p1,k3)
-      REAL(wp), DIMENSION(:,:), INTENT(  out) ::  p1
+      REAL(dp), DIMENSION(:,:), INTENT(  out) ::  p1
       INTEGER                 , INTENT(in   ) ::  k1,k2
       INTEGER                 , INTENT(  out) ::  k3