mom_coms::reproducing_sum real function real function mom_coms::reproducing_sum::reproducing_sum_2d (array, isr, ier, jsr, jer, EFP_sum, reproducing, overflow_check, err, only_on_PE) reproducing_sum_2d array array isr isr ier ier jsr jsr jer jer EFP_sum EFP_sum reproducing reproducing overflow_check overflow_check err err only_on_PE only_on_PE This subroutine uses a conversion to an integer representation of real numbers to give an order-invariant sum of distributed 2-D arrays that reproduces across domain decomposition. This technique is described in Hallberg & Adcroft, 2014, Parallel Computing, doi:10.1016/j.parco.2014.04.007. array The array to be summed isr The starting i-index of the sum, noting that the array indices starts at 1 ier The ending i-index of the sum, noting that the array indices starts at 1 jsr The starting j-index of the sum, noting that the array indices starts at 1 jer The ending j-index of the sum, noting that the array indices starts at 1 efp_sum The result in extended fixed point format reproducing If present and false, do the sum using the naive non-reproducing approach overflow_check If present and false, disable checking for overflows in incremental results. This can speed up calculations if the number of values being summed is small enough err If present, return an error code instead of triggering any fatal errors directly from this routine. only_on_pe If present and true, do not do the sum across processors, only reporting the local sum Result real function real function mom_coms::reproducing_sum::reproducing_sum_3d (array, isr, ier, jsr, jer, sums, EFP_sum, EFP_lay_sums, err, only_on_PE) reproducing_sum_3d array array isr isr ier ier jsr jsr jer jer sums sums EFP_sum EFP_sum EFP_lay_sums EFP_lay_sums err err only_on_PE only_on_PE This subroutine uses a conversion to an integer representation of real numbers to give an order-invariant sum of distributed 3-D arrays that reproduces across domain decomposition. This technique is described in Hallberg & Adcroft, 2014, Parallel Computing, doi:10.1016/j.parco.2014.04.007. array The array to be summed isr The starting i-index of the sum, noting that the array indices starts at 1 ier The ending i-index of the sum, noting that the array indices starts at 1 jsr The starting j-index of the sum, noting that the array indices starts at 1 jer The ending j-index of the sum, noting that the array indices starts at 1 sums The sums by vertical layer efp_sum The result in extended fixed point format efp_lay_sums The sums by vertical layer in EFP format err If present, return an error code instead of triggering any fatal errors directly from this routine. only_on_pe If present and true, do not do the sum across processors, only reporting the local sum Result Find an accurate and order-invariant sum of a distributed 2d or 3d field. mom_coms::reproducing_sumreproducing_sum_2d mom_coms::reproducing_sumreproducing_sum_3d