pcm_functions Module Reference

Functions/Subroutines
subroutine, public	pcm_reconstruction (N, u, ppoly_E, ppoly_coefficients)

Function/Subroutine Documentation

pcm_reconstruction()

subroutine, public pcm_functions::pcm_reconstruction	(	integer, intent(in)	N,
		real, dimension(:), intent(in)	u,
		real, dimension(:,:), intent(inout)	ppoly_E,
		real, dimension(:,:), intent(inout)	ppoly_coefficients
	)

Definition at line 24 of file PCM_functions.F90.

Referenced by mom_remapping::build_reconstructions_1d(), and mom_remapping::remapping_unit_tests().

!------------------------------------------------------------------------------
! Reconstruction by constant polynomials within each cell. There is nothing to
! do but this routine is provided to ensure a homogeneous interface
! throughout the regridding toolbox.
!
! N:     number of cells in grid
! h:     thicknesses of grid cells
! u:     cell averages to use in constructing piecewise polynomials
! ppoly_E : edge values of piecewise polynomials
! ppoly_coefficients : coefficients of piecewise polynomials
!
! It is assumed that the dimension of 'u' is equal to the number of cells
! defining 'grid' and 'ppoly'. No consistency check is performed.
!------------------------------------------------------------------------------

  ! Arguments
  integer,              intent(in)    :: n ! Number of cells
  real, dimension(:),   intent(in)    :: u ! cell averages
  real, dimension(:,:), intent(inout) :: ppoly_e            !Edge value of polynomial
  real, dimension(:,:), intent(inout) :: ppoly_coefficients !Coefficients of polynomial

  ! Local variables
  integer :: k

  ! The coefficients of the piecewise constant polynomial are simply
  ! the cell averages.
  ppoly_coefficients(:,1) = u(:)

  ! The edge values are equal to the cell average
  do k = 1,n
    ppoly_e(k,:) = u(k)
  end do

Here is the caller graph for this function: