mpp_sum.inc

! -*-f90-*-
 
!***********************************************************************
!*                   GNU Lesser General Public License
!*
!* This file is part of the GFDL Flexible Modeling System (FMS).
!*
!* FMS is free software: you can redistribute it and/or modify it under
!* the terms of the GNU Lesser General Public License as published by
!* the Free Software Foundation, either version 3 of the License, or (at
!* your option) any later version.
!*
!* FMS is distributed in the hope that it will be useful, but WITHOUT
!* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
!* FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
!* for more details.
!*
!* You should have received a copy of the GNU Lesser General Public
!* License along with FMS.  If not, see <http://www.gnu.org/licenses/>.
!***********************************************************************
 
!> @file
!> @brief Summation routines for @ref mpp_mod
 
!> @addtogroup mpp_mod
!> @{
!#######################################################################
 
    !> Sums array a when only first element is passed: this routine just converts to a call to MPP_SUM_
    subroutine mpp_sum_scalar_( a, pelist )
      mpp_type_, intent(inout) :: a
      integer, intent(in), optional :: pelist(:)
      mpp_type_ :: b(1)
 
      b(1) = a
      if( debug )call mpp_error( note, 'MPP_SUM_SCALAR_: calling MPP_SUM_ ...' )
      call mpp_sum_( b, 1, pelist )
      a = b(1)
      return
    end subroutine mpp_sum_scalar_
 
!#######################################################################
    !> Sums 2d array across pes
    subroutine mpp_sum_2d_( a, length, pelist )
      mpp_type_, intent(inout) :: a(:,:) !< 2d array to sum
      integer, intent(in) :: length !< amount of indices in given 2d array
      integer, intent(in), optional :: pelist(:) !< pelist to calculate sum across
      mpp_type_ :: a1d(length)
 
      pointer( ptr, a1d )
      ptr = loc(a)
      call mpp_sum( a1d, length, pelist )
 
      return
    end subroutine mpp_sum_2d_
 
!#######################################################################
    !> Sums 3d array across pes
    subroutine mpp_sum_3d_( a, length, pelist )
      mpp_type_, intent(inout) :: a(:,:,:) !< 3d array to sum
      integer, intent(in) :: length !< amount of indices in given 3d array
      integer, intent(in), optional :: pelist(:) !< pelist to calculate sum across
      mpp_type_ :: a1d(length)
 
      pointer( ptr, a1d )
      ptr = loc(a)
      call mpp_sum( a1d, length, pelist )
 
      return
    end subroutine mpp_sum_3d_
 
!#######################################################################
    !> Sums 4d array across pes
    subroutine mpp_sum_4d_( a, length, pelist )
      mpp_type_, intent(inout) :: a(:,:,:,:) !< 4d array to sum
      integer, intent(in) :: length !< amount of indices in given 4d array
      integer, intent(in), optional :: pelist(:) !< pelist to calculate sum across
      mpp_type_ :: a1d(length)
 
      pointer( ptr, a1d )
      ptr = loc(a)
      call mpp_sum( a1d, length, pelist )
 
      return
    end subroutine mpp_sum_4d_
 
!#######################################################################
    !> Sums 5d array across pes
    subroutine mpp_sum_5d_( a, length, pelist )
      mpp_type_, intent(inout) :: a(:,:,:,:,:) !< 5d array to sum
      integer, intent(in) :: length !< amount of indices in given 5d array
      integer, intent(in), optional :: pelist(:) !< pelist to calculate sum across
      mpp_type_ :: a1d(length)
 
      pointer( ptr, a1d )
      ptr = loc(a)
      call mpp_sum( a1d, length, pelist )
 
      return
    end subroutine mpp_sum_5d_
!> @}