Changeset 12032 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src

Timestamp:

2019-12-02T19:08:44+01:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_nemo: improve isf debug with a basic checksum

File:

• NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfutils.F90 (modified) (4 diffs)

NEMO/branches/2019/ENHANCE-02_ISF_nemo/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
@@ -60 +62 @@
       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)
+      !
+      ! print out
+      IF (lwp) WRITE(numout,*) TRIM(cdtxt),' (min, max, sum, tag) : ',zmin, zmax, zsum, itmps
       !
       CALL FLUSH(numout)
@@ -84 +104 @@
       !!--------------------------------------------------------------------
       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)
+      !
+      ! print out
+      IF (lwp) WRITE(numout,*) TRIM(cdtxt),' (min, max, sum, tag) : ',zmin, zmax, zsum, itmps
       !
       CALL FLUSH(numout)
@@ -97 +137 @@
 
 END MODULE isfutils
-