pcm_functions Module Reference

Detailed Description

Piecewise constant reconstruction functions.

Functions/Subroutines
subroutine, public	pcm_reconstruction (N, u, ppoly_E, ppoly_coef)
	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. More...

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_coef
	)

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.

It is assumed that the dimension of 'u' is equal to the number of cells defining 'grid' and 'ppoly'. No consistency check is performed.

Parameters

[in]	n	Number of cells
[in]	u	cell averages
[in,out]	ppoly_e	Edge value of polynomial, with the same units as u.
[in,out]	ppoly_coef	Coefficients of polynomial, with the same units as u.

Definition at line 19 of file PCM_functions.F90.

  integer,              intent(in)    :: N !< Number of cells
  real, dimension(:),   intent(in)    :: u !< cell averages
  real, dimension(:,:), intent(inout) :: ppoly_E    !< Edge value of polynomial,
                                           !! with the same units as u.
  real, dimension(:,:), intent(inout) :: ppoly_coef !< Coefficients of polynomial,
                                           !! with the same units as u.
 
  ! Local variables
  integer :: k
 
  ! The coefficients of the piecewise constant polynomial are simply
  ! the cell averages.
  ppoly_coef(:,1) = u(:)
 
  ! The edge values are equal to the cell average
  do k = 1,n
    ppoly_e(k,:) = u(k)
  enddo
 

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

Here is the caller graph for this function: