MOM_checksums.F90 File Reference

#include "version_variable.h"

Include dependency graph for MOM_checksums.F90:

Go to the source code of this file.

Data Types
interface	mom_checksums::hchksum_pair
interface	mom_checksums::uvchksum
interface	mom_checksums::uchksum
interface	mom_checksums::vchksum
interface	mom_checksums::bchksum_pair
interface	mom_checksums::hchksum
interface	mom_checksums::bchksum
interface	mom_checksums::qchksum
interface	mom_checksums::chksum
interface	mom_checksums::chk_sum_msg
interface	mom_checksums::is_nan

Modules
module	mom_checksums

Functions/Subroutines
subroutine	mom_checksums::chksum_pair_h_2d (mesg, arrayA, arrayB, HI, haloshift, omit_corners, scale)
subroutine	mom_checksums::chksum_pair_h_3d (mesg, arrayA, arrayB, HI, haloshift, omit_corners, scale)
subroutine	mom_checksums::chksum_h_2d (array, mesg, HI, haloshift, omit_corners, scale)
	chksum_h_2d performs checksums on a 2d array staggered at tracer points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, scale)
subroutine	mom_checksums::chksum_pair_b_2d (mesg, arrayA, arrayB, HI, haloshift, symmetric, omit_corners, scale)
subroutine	mom_checksums::chksum_pair_b_3d (mesg, arrayA, arrayB, HI, haloshift, symmetric, omit_corners, scale)
subroutine	mom_checksums::chksum_b_2d (array, mesg, HI, haloshift, symmetric, omit_corners, scale)
	chksum_B_2d performs checksums on a 2d array staggered at corner points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, sym_stats, scale)
subroutine	mom_checksums::chksum_uv_2d (mesg, arrayU, arrayV, HI, haloshift, symmetric, omit_corners, scale)
subroutine	mom_checksums::chksum_uv_3d (mesg, arrayU, arrayV, HI, haloshift, symmetric, omit_corners, scale)
subroutine	mom_checksums::chksum_u_2d (array, mesg, HI, haloshift, symmetric, omit_corners, scale)
	chksum_u_2d performs checksums on a 2d array staggered at C-grid u points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, sym_stats, scale)
subroutine	mom_checksums::chksum_v_2d (array, mesg, HI, haloshift, symmetric, omit_corners, scale)
	chksum_v_2d performs checksums on a 2d array staggered at C-grid v points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, sym_stats, scale)
subroutine	mom_checksums::chksum_h_3d (array, mesg, HI, haloshift, omit_corners, scale)
	chksum_h_3d performs checksums on a 3d array staggered at tracer points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, scale)
subroutine	mom_checksums::chksum_b_3d (array, mesg, HI, haloshift, symmetric, omit_corners, scale)
	chksum_B_3d performs checksums on a 3d array staggered at corner points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, sym_stats, scale)
subroutine	mom_checksums::chksum_u_3d (array, mesg, HI, haloshift, symmetric, omit_corners, scale)
	chksum_u_3d performs checksums on a 3d array staggered at C-grid u points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, sym_stats, scale)
subroutine	mom_checksums::chksum_v_3d (array, mesg, HI, haloshift, symmetric, omit_corners, scale)
	chksum_v_3d performs checksums on a 3d array staggered at C-grid v points. More...
integer function	subchk (array, HI, di, dj, scale)
subroutine	substats (HI, array, mesg, sym_stats, scale)
subroutine	mom_checksums::chksum1d (array, mesg, start_i, end_i, compare_PEs)
	chksum1d does a checksum of a 1-dimensional array. More...
subroutine	mom_checksums::chksum2d (array, mesg)
	chksum2d does a checksum of all data in a 2-d array. More...
subroutine	mom_checksums::chksum3d (array, mesg)
	chksum3d does a checksum of all data in a 2-d array. More...
logical function	mom_checksums::is_nan_0d (x)
	This function returns .true. if x is a NaN, and .false. otherwise. More...
logical function	mom_checksums::is_nan_1d (x, skip_mpp)
	This function returns .true. if any element of x is a NaN, and .false. otherwise. More...
logical function	mom_checksums::is_nan_2d (x)
	This function returns .true. if any element of x is a NaN, and .false. otherwise. More...
logical function	mom_checksums::is_nan_3d (x)
	This function returns .true. if any element of x is a NaN, and .false. otherwise. More...
subroutine	mom_checksums::chk_sum_msg1 (fmsg, bc0, mesg)
subroutine	mom_checksums::chk_sum_msg5 (fmsg, bc0, bcSW, bcSE, bcNW, bcNE, mesg)
subroutine	mom_checksums::chk_sum_msg_nsew (fmsg, bc0, bcN, bcS, bcE, bcW, mesg)
subroutine	mom_checksums::chk_sum_msg_s (fmsg, bc0, bcS, mesg)
subroutine	mom_checksums::chk_sum_msg_w (fmsg, bc0, bcW, mesg)
subroutine	mom_checksums::chk_sum_msg2 (fmsg, bc0, bcSW, mesg)
subroutine	mom_checksums::chk_sum_msg3 (fmsg, aMean, aMin, aMax, mesg)
subroutine, public	mom_checksums::mom_checksums_init (param_file)
	MOM_checksums_init initializes the MOM_checksums module. As it happens, the only thing that it does is to log the version of this module. More...
subroutine	mom_checksums::chksum_error (signal, message)

Variables
integer, parameter	mom_checksums::default_shift =0
logical	mom_checksums::calculatestatistics =.true.
logical	mom_checksums::writechksums =.true.
logical	mom_checksums::checkfornans =.true.

Function/Subroutine Documentation

subchk() [1/8]

integer function chksum_h_2d::subchk ( real, dimension(hi%isd:,hi%jsd:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 193 of file MOM_checksums.F90.

Referenced by mom_checksums::chksum_b_2d(), mom_checksums::chksum_b_3d(), mom_checksums::chksum_h_2d(), mom_checksums::chksum_h_3d(), mom_checksums::chksum_u_2d(), mom_checksums::chksum_u_3d(), mom_checksums::chksum_v_2d(), mom_checksums::chksum_v_3d(), and subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%jsd:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, bc
    subchk = 0
    do j=hi%jsc+dj,hi%jec+dj; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j)))
      subchk = subchk + bc
    enddo; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the caller graph for this function:

subchk() [2/8]

integer function chksum_b_2d::subchk ( real, dimension(hi%isdb:,hi%jsdb:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 362 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%JsdB:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, bc
    subchk = 0
    ! This line deliberately uses the h-point computational domain.
    do j=hi%jsc+dj,hi%jec+dj; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j)))
      subchk = subchk + bc
    enddo; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

subchk() [3/8]

integer function chksum_u_2d::subchk ( real, dimension(hi%isdb:,hi%jsd:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 531 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%jsd:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, bc
    subchk = 0
    ! This line deliberately uses the h-point computational domain.
    do j=hi%jsc+dj,hi%jec+dj; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j)))
      subchk = subchk + bc
    enddo; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

subchk() [4/8]

integer function chksum_v_2d::subchk ( real, dimension(hi%isd:,hi%jsdb:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 659 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%JsdB:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, bc
    subchk = 0
    ! This line deliberately uses the h-point computational domain.
    do j=hi%jsc+dj,hi%jec+dj; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j)))
      subchk = subchk + bc
    enddo; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

subchk() [5/8]

integer function chksum_h_3d::subchk ( real, dimension(hi%isd:,hi%jsd:,:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 770 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%jsd:,:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, k, bc
    subchk = 0
    do k=lbound(array,3),ubound(array,3) ; do j=hi%jsc+dj,hi%jec+dj ; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j,k)))
      subchk = subchk + bc
    enddo ; enddo ; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

subchk() [6/8]

integer function chksum_b_3d::subchk ( real, dimension(hi%isdb:,hi%jsdb:,:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 894 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%JsdB:,:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, k, bc
    subchk = 0
    ! This line deliberately uses the h-point computational domain.
    do k=lbound(array,3),ubound(array,3) ; do j=hi%jsc+dj,hi%jec+dj ; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j,k)))
      subchk = subchk + bc
    enddo ; enddo ; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

subchk() [7/8]

integer function chksum_u_3d::subchk ( real, dimension(hi%isdb:,hi%jsd:,:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 1022 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%jsd:,:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, k, bc
    subchk = 0
    ! This line deliberately uses the h-point computational domain.
    do k=lbound(array,3),ubound(array,3) ; do j=hi%jsc+dj,hi%jec+dj ; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j,k)))
      subchk = subchk + bc
    enddo ; enddo ; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

subchk() [8/8]

integer function chksum_v_3d::subchk ( real, dimension(hi%isd:,hi%jsdb:,:), intent(in)  array, type(hor_index_type), intent(in)  HI, integer, intent(in)  di, integer, intent(in)  dj, real, intent(in)  scale  )

private

Definition at line 1150 of file MOM_checksums.F90.

References subchk().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%JsdB:,:), intent(in) :: array
    integer, intent(in) :: di, dj
    real, intent(in) :: scale
    integer :: bitcount, i, j, k, bc
    subchk = 0
    ! This line deliberately uses the h-point computational domain.
    do k=lbound(array,3),ubound(array,3) ; do j=hi%jsc+dj,hi%jec+dj ; do i=hi%isc+di,hi%iec+di
      bc = bitcount(abs(scale*array(i,j,k)))
      subchk = subchk + bc
    enddo ; enddo ; enddo
    call sum_across_pes(subchk)
    subchk=mod(subchk,1000000000)

Here is the call graph for this function:

substats() [1/8]

subroutine chksum_h_2d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isd:,hi%jsd:), intent(in)  array, character(len=*), intent(in)  mesg, real, intent(in)  scale  )

private

Definition at line 208 of file MOM_checksums.F90.

References mom_error_handler::is_root_pe().

Referenced by mom_checksums::chksum_b_2d(), mom_checksums::chksum_b_3d(), mom_checksums::chksum_h_2d(), mom_checksums::chksum_h_3d(), mom_checksums::chksum_u_2d(), mom_checksums::chksum_u_3d(), mom_checksums::chksum_v_2d(), and mom_checksums::chksum_v_3d().

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%jsd:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    real, intent(in) :: scale
    integer :: i, j, n
    real :: amean, amin, amax

    amin = array(hi%isc,hi%jsc)
    amax = array(hi%isc,hi%jsc)
    n = 0
    do j=hi%jsc,hi%jec ; do i=hi%isc,hi%iec
      amin = min(amin, array(i,j))
      amax = max(amax, array(i,j))
      n = n + 1
    enddo ; enddo
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec))
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("h-point:",amean*scale,amin*scale,amax*scale,mesg)

Here is the call graph for this function:

Here is the caller graph for this function:

substats() [2/8]

subroutine chksum_b_2d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isdb:,hi%jsdb:), intent(in)  array, character(len=*), intent(in)  mesg, logical, intent(in)  sym_stats, real, intent(in)  scale  )

private

Definition at line 378 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%JsdB:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    logical, intent(in) :: sym_stats
    real, intent(in) :: scale
    integer :: i, j, n, isb, jsb
    real :: amean, amin, amax

    isb = hi%isc ; if (sym_stats) isb = hi%isc-1
    jsb = hi%jsc ; if (sym_stats) jsb = hi%jsc-1

    amin = array(hi%isc,hi%jsc) ; amax = amin
    do j=jsb,hi%JecB ; do i=isb,hi%IecB
      amin = min(amin, array(i,j))
      amax = max(amax, array(i,j))
    enddo ; enddo
    ! This line deliberately uses the h-point computational domain.
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec))
    n = (1 + hi%jec - hi%jsc) * (1 + hi%iec - hi%isc)
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("B-point:",amean*scale,amin*scale,amax*scale,mesg)

substats() [3/8]

subroutine chksum_u_2d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isdb:,hi%jsd:), intent(in)  array, character(len=*), intent(in)  mesg, logical, intent(in)  sym_stats, real, intent(in)  scale  )

private

Definition at line 547 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%jsd:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    logical, intent(in) :: sym_stats
    real, intent(in) :: scale
    integer :: i, j, n, isb
    real :: amean, amin, amax

    isb = hi%isc ; if (sym_stats) isb = hi%isc-1

    amin = array(hi%isc,hi%jsc) ; amax = amin
    do j=hi%jsc,hi%jec ; do i=isb,hi%IecB
      amin = min(amin, array(i,j))
      amax = max(amax, array(i,j))
    enddo ; enddo
    ! This line deliberately uses the h-point computational domain.
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec))
    n = (1 + hi%jec - hi%jsc) * (1 + hi%iec - hi%isc)
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("u-point:",amean*scale,amin*scale,amax*scale,mesg)

substats() [4/8]

subroutine chksum_v_2d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isd:,hi%jsdb:), intent(in)  array, character(len=*), intent(in)  mesg, logical, intent(in)  sym_stats, real, intent(in)  scale  )

private

Definition at line 675 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%JsdB:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    logical, intent(in) :: sym_stats
    real, intent(in) :: scale
    integer :: i, j, n, jsb
    real :: amean, amin, amax

    jsb = hi%jsc ; if (sym_stats) jsb = hi%jsc-1

    amin = array(hi%isc,hi%jsc) ; amax = amin
    do j=jsb,hi%JecB ; do i=hi%isc,hi%iec
      amin = min(amin, array(i,j))
      amax = max(amax, array(i,j))
    enddo ; enddo
    ! This line deliberately uses the h-computational domain.
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec))
    n = (1 + hi%jec - hi%jsc) * (1 + hi%iec - hi%isc)
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("v-point:",amean*scale,amin*scale,amax*scale,mesg)

substats() [5/8]

subroutine chksum_h_3d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isd:,hi%jsd:,:), intent(in)  array, character(len=*), intent(in)  mesg, real, intent(in)  scale  )

private

Definition at line 785 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%jsd:,:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    real, intent(in) :: scale
    integer :: i, j, k, n
    real :: amean, amin, amax

    amin = array(hi%isc,hi%jsc,1)
    amax = array(hi%isc,hi%jsc,1)
    n = 0
    do k=lbound(array,3),ubound(array,3) ; do j=hi%jsc,hi%jec ; do i=hi%isc,hi%iec
      amin = min(amin, array(i,j,k))
      amax = max(amax, array(i,j,k))
      n = n + 1
    enddo ; enddo ; enddo
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec,:))
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("h-point:",amean*scale,amin*scale,amax*scale,mesg)

substats() [6/8]

subroutine chksum_b_3d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isdb:,hi%jsdb:,:), intent(in)  array, character(len=*), intent(in)  mesg, logical, intent(in)  sym_stats, real, intent(in)  scale  )

private

Definition at line 910 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%JsdB:,:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    logical, intent(in) :: sym_stats
    real, intent(in) :: scale
    integer :: i, j, k, n, isb, jsb
    real :: amean, amin, amax

    isb = hi%isc ; if (sym_stats) isb = hi%isc-1
    jsb = hi%jsc ; if (sym_stats) jsb = hi%jsc-1

    amin = array(hi%isc,hi%jsc,1) ; amax = amin
    do k=lbound(array,3),ubound(array,3) ; do j=jsb,hi%JecB ; do i=isb,hi%IecB
      amin = min(amin, array(i,j,k))
      amax = max(amax, array(i,j,k))
    enddo ; enddo ; enddo
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec,:))
    n = (1 + hi%jec - hi%jsc) * (1 + hi%iec - hi%isc) * size(array,3)
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("B-point:",amean*scale,amin*scale,amax*scale,mesg)

substats() [7/8]

subroutine chksum_u_3d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isdb:,hi%jsd:,:), intent(in)  array, character(len=*), intent(in)  mesg, logical, intent(in)  sym_stats, real, intent(in)  scale  )

private

Definition at line 1038 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%IsdB:,HI%jsd:,:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    logical, intent(in) :: sym_stats
    real, intent(in) :: scale
    integer :: i, j, k, n, isb
    real :: amean, amin, amax

    isb = hi%isc ; if (sym_stats) isb = hi%isc-1

    amin = array(hi%isc,hi%jsc,1) ; amax = amin
    do k=lbound(array,3),ubound(array,3) ; do j=hi%jsc,hi%jec ; do i=isb,hi%IecB
      amin = min(amin, array(i,j,k))
      amax = max(amax, array(i,j,k))
    enddo ; enddo ; enddo
    ! This line deliberately uses the h-point computational domain.
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec,:))
    n = (1 + hi%jec - hi%jsc) * (1 + hi%iec - hi%isc) * size(array,3)
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("u-point:",amean*scale,amin*scale,amax*scale,mesg)

substats() [8/8]

subroutine chksum_v_3d::substats ( type(hor_index_type), intent(in)  HI, real, dimension(hi%isd:,hi%jsdb:,:), intent(in)  array, character(len=*), intent(in)  mesg, logical, intent(in)  sym_stats, real, intent(in)  scale  )

private

Definition at line 1166 of file MOM_checksums.F90.

    type(hor_index_type), intent(in) :: hi
    real, dimension(HI%isd:,HI%JsdB:,:), intent(in) :: array
    character(len=*), intent(in) :: mesg
    logical, intent(in) :: sym_stats
    real, intent(in) :: scale
    integer :: i, j, k, n, jsb
    real :: amean, amin, amax

    jsb = hi%jsc ; if (sym_stats) jsb = hi%jsc-1

    amin = array(hi%isc,hi%jsc,1) ; amax = amin
    do k=lbound(array,3),ubound(array,3) ; do j=jsb,hi%JecB ; do i=hi%isc,hi%iec
      amin = min(amin, array(i,j,k))
      amax = max(amax, array(i,j,k))
    enddo ; enddo ; enddo
    ! This line deliberately uses the h-point computational domain.
    amean = reproducing_sum(array(hi%isc:hi%iec,hi%jsc:hi%jec,:))
    n = (1 + hi%jec - hi%jsc) * (1 + hi%iec - hi%isc) * size(array,3)
    call sum_across_pes(n)
    call min_across_pes(amin)
    call max_across_pes(amax)
    amean = amean / real(n)
    if (is_root_pe()) call chk_sum_msg("v-point:",amean*scale,amin*scale,amax*scale,mesg)