Changeset 12077 for NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/ISF/isfutils.F90

Timestamp:

2019-12-05T18:41:39+01:00 (4 years ago)

Author:

mathiot

Message:

include ENHANCE-02_ISF_nemo in UKMO merge branch

File:

• NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/ISF/isfutils.F90 (modified) (3 diffs)

NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/ISF/isfutils.F90

@@ -15 +15 @@
    USE lib_fortran   , ONLY: glob_sum, glob_min, glob_max             ! compute global value
    USE par_oce       , ONLY: jpi,jpj,jpk                              ! domain size
-   USE in_out_manager, ONLY: wp, lwp, numout                          ! miscelenious
+   USE dom_oce       , ONLY: nldi, nlei, nldj, nlej                   ! local domain
+   USE in_out_manager, ONLY: i8, wp, lwp, numout                          ! miscelenious
+   USE lib_mpp
 
    IMPLICIT NONE
@@ -54 +56 @@
       !!                  ***  ROUTINE isf_debug2d  ***
       !!
-      !! ** Purpose : add debug print
+      !! ** Purpose : add debug print for 2d variables
       !!
       !!-------------------------- IN  -------------------------------------
-      CHARACTER(LEN=256)          , INTENT(in   ) :: cdtxt
+      CHARACTER(LEN=*)            , INTENT(in   ) :: cdtxt
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pvar
       !!--------------------------------------------------------------------
-      REAL(wp) :: zmin, zmax, zsum
+      REAL(wp)    :: zmin, zmax, zsum
+      INTEGER(i8) :: imodd, ip
+      INTEGER     :: itmps,imods, ji, jj, jk
       !!--------------------------------------------------------------------
       !
+      ! global min/max/sum to check data range and NaN
       zsum = glob_sum( 'debug', pvar(:,:) )
       zmin = glob_min( 'debug', pvar(:,:) )
       zmax = glob_max( 'debug', pvar(:,:) )
       !
-      IF (lwp) WRITE(numout,*) TRIM(cdtxt),' (min, max, sum) : ',zmin, zmax, zsum
+      ! basic check sum to check reproducibility
+      ! TRANSFER function find out the integer corresponding to pvar(i,j) bit pattern
+      ! MOD allow us to keep only the latest digits during the sum
+      ! imod is not choosen to be very large as at the end there is a classic mpp_sum
+      imodd=65521 ! highest prime number < 2**16 with i8 type
+      imods=65521 ! highest prime number < 2**16 with default integer for mpp_sum subroutine
+      itmps=0
+      DO jj=nldj,nlej
+         DO ji=nldi,nlei
+            itmps = MOD(itmps + MOD(TRANSFER(pvar(ji,jj), ip),imodd), imods)
+         END DO
+      END DO
+      CALL mpp_sum('debug',itmps)
       !
-      CALL FLUSH(numout)
+      ! print out
+      IF (lwp) THEN
+         WRITE(numout,*) TRIM(cdtxt),' (min, max, sum, tag) : ',zmin, zmax, zsum, itmps
+         CALL FLUSH(numout)
+      END IF
       !
    END SUBROUTINE debug2d
@@ -77 +98 @@
       !!                  ***  ROUTINE isf_debug3d  ***
       !!
-      !! ** Purpose : add debug print
+      !! ** Purpose : add debug print for 3d variables
       !!
       !!-------------------------- IN  -------------------------------------
-      CHARACTER(LEN=256)              , INTENT(in   ) :: cdtxt
+      CHARACTER(LEN=*)                , INTENT(in   ) :: cdtxt
       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in   ) :: pvar
       !!--------------------------------------------------------------------
       REAL(wp) :: zmin, zmax, zsum
+      INTEGER(i8) :: imodd, ip
+      INTEGER     :: itmps,imods, ji, jj, jk
       !!--------------------------------------------------------------------
       !
+      ! global min/max/sum to check data range and NaN
       zsum = glob_sum( 'debug', pvar(:,:,:) )
       zmin = glob_min( 'debug', pvar(:,:,:) )
       zmax = glob_max( 'debug', pvar(:,:,:) )
       !
-      IF (lwp) WRITE(numout,*) TRIM(cdtxt),' (min, max, sum) : ',zmin, zmax, zsum
+      ! basic check sum to check reproducibility
+      ! TRANSFER function find out the integer corresponding to pvar(i,j) bit pattern
+      ! MOD allow us to keep only the latest digits during the sum
+      ! imod is not choosen to be very large as at the end there is a classic mpp_sum
+      imodd=65521 ! highest prime number < 2**16 with i8 type
+      imods=65521 ! highest prime number < 2**16 with default integer for mpp_sum subroutine
+      itmps=0
+      DO jk=1,jpk
+         DO jj=nldj,nlej
+            DO ji=nldi,nlei
+               itmps = MOD(itmps + MOD(TRANSFER(pvar(ji,jj,jk), ip),imodd), imods)
+            END DO
+         END DO
+      END DO
+      CALL mpp_sum('debug',itmps)
       !
-      CALL FLUSH(numout)
+      ! print out
+      IF (lwp) THEN
+         WRITE(numout,*) TRIM(cdtxt),' (min, max, sum, tag) : ',zmin, zmax, zsum, itmps
+         CALL FLUSH(numout)
+      END IF
       !
    END SUBROUTINE debug3d
 
 END MODULE isfutils
-