Changeset 14219 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DOM/domutl.F90

Timestamp:

2020-12-18T18:52:57+01:00 (4 years ago)

Author:

mcastril

Message:

Add Mixed Precision support by Oriol Tintó

File:

• NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DOM/domutl.F90 (modified) (3 diffs)

NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DOM/domutl.F90

@@ -22 +22 @@
 
    INTERFACE is_tile
-      MODULE PROCEDURE is_tile_2d, is_tile_3d, is_tile_4d
+      MODULE PROCEDURE is_tile_2d_sp, is_tile_3d_sp, is_tile_4d_sp
+      MODULE PROCEDURE is_tile_2d_dp, is_tile_3d_dp, is_tile_4d_dp
    END INTERFACE is_tile
 
@@ -108 +109 @@
       !
       puniq(:,:) = ztstref(:,:)                    ! default definition
-      CALL lbc_lnk( 'domwri', puniq, cdgrd, 1. )   ! apply boundary conditions
+      CALL lbc_lnk( 'domwri', puniq, cdgrd, 1._wp )   ! apply boundary conditions
       lluniq(:,:,1) = puniq(:,:) == ztstref(:,:)   ! check which values have not been changed
       !
@@ -116 +117 @@
 
 
-   FUNCTION is_tile_2d( pt )
+   FUNCTION is_tile_2d_sp( pt )
       !!
-      REAL(wp), DIMENSION(:,:), INTENT(in) ::   pt
-      INTEGER :: is_tile_2d
+      REAL(sp), DIMENSION(:,:), INTENT(in) ::   pt
+      INTEGER :: is_tile_2d_sp
       !!
       IF( ln_tile .AND. (SIZE(pt, 1) < jpi .OR. SIZE(pt, 2) < jpj) ) THEN
-         is_tile_2d = 1
+         is_tile_2d_sp = 1
       ELSE
-         is_tile_2d = 0
+         is_tile_2d_sp = 0
       ENDIF
-   END FUNCTION is_tile_2d
+   END FUNCTION is_tile_2d_sp
 
 
-   FUNCTION is_tile_3d( pt )
+   FUNCTION is_tile_3d_sp( pt )
       !!
-      REAL(wp), DIMENSION(:,:,:), INTENT(in) ::   pt
-      INTEGER :: is_tile_3d
+      REAL(sp), DIMENSION(:,:,:), INTENT(in) ::   pt
+      INTEGER :: is_tile_3d_sp
       !!
       IF( ln_tile .AND. (SIZE(pt, 1) < jpi .OR. SIZE(pt, 2) < jpj) ) THEN
-         is_tile_3d = 1
+         is_tile_3d_sp = 1
       ELSE
-         is_tile_3d = 0
+         is_tile_3d_sp = 0
       ENDIF
-   END FUNCTION is_tile_3d
+   END FUNCTION is_tile_3d_sp
 
 
-   FUNCTION is_tile_4d( pt )
+   FUNCTION is_tile_4d_sp( pt )
       !!
-      REAL(wp), DIMENSION(:,:,:,:), INTENT(in) ::   pt
-      INTEGER :: is_tile_4d
+      REAL(sp), DIMENSION(:,:,:,:), INTENT(in) ::   pt
+      INTEGER :: is_tile_4d_sp
       !!
       IF( ln_tile .AND. (SIZE(pt, 1) < jpi .OR. SIZE(pt, 2) < jpj) ) THEN
-         is_tile_4d = 1
+         is_tile_4d_sp = 1
       ELSE
-         is_tile_4d = 0
+         is_tile_4d_sp = 0
       ENDIF
-   END FUNCTION is_tile_4d
+   END FUNCTION is_tile_4d_sp
+
+   FUNCTION is_tile_2d_dp( pt )
+      !!
+      REAL(dp), DIMENSION(:,:), INTENT(in) ::   pt
+      INTEGER :: is_tile_2d_dp
+      !!
+      IF( ln_tile .AND. (SIZE(pt, 1) < jpi .OR. SIZE(pt, 2) < jpj) ) THEN
+         is_tile_2d_dp = 1
+      ELSE
+         is_tile_2d_dp = 0
+      ENDIF
+   END FUNCTION is_tile_2d_dp
+
+
+   FUNCTION is_tile_3d_dp( pt )
+      !!
+      REAL(dp), DIMENSION(:,:,:), INTENT(in) ::   pt
+      INTEGER :: is_tile_3d_dp
+      !!
+      IF( ln_tile .AND. (SIZE(pt, 1) < jpi .OR. SIZE(pt, 2) < jpj) ) THEN
+         is_tile_3d_dp = 1
+      ELSE
+         is_tile_3d_dp = 0
+      ENDIF
+   END FUNCTION is_tile_3d_dp
+
+
+   FUNCTION is_tile_4d_dp( pt )
+      !!
+      REAL(dp), DIMENSION(:,:,:,:), INTENT(in) ::   pt
+      INTEGER :: is_tile_4d_dp
+      !!
+      IF( ln_tile .AND. (SIZE(pt, 1) < jpi .OR. SIZE(pt, 2) < jpj) ) THEN
+         is_tile_4d_dp = 1
+      ELSE
+         is_tile_4d_dp = 0
+      ENDIF
+   END FUNCTION is_tile_4d_dp
+
 
    !!======================================================================