mus_intp_error Subroutine

private subroutine mus_intp_error(pdf, level, tree, general, fac, scheme, layout, intp, ind, varSys, derVarPos, eType)

Determine the numerical error of the interpolated quantities to the given initial conditions

Arguments

TypeIntentOptionalAttributesName
type(pdf_data_type), intent(in) :: pdf

global fluid parameters

integer, intent(in) :: level

level

type(treelmesh_type), intent(in) :: tree

treelm tree

type(tem_general_type), intent(in) :: general
type(mus_convertFac_type), intent(in) :: fac
type(mus_scheme_type), intent(in) :: scheme

the current scheme

type(mus_scheme_layout_type), intent(in) :: layout

the layout used

type(mus_interpolation_type), intent(in) :: intp

interpolation method info

type(grw_intarray_type), intent(in) :: ind

indirection list

type(tem_varSys_type), intent(in) :: varSys

global variable system

type(mus_derVarPos_type), intent(in) :: derVarPos(:)

required variable system which maps to global varsys

integer, intent(in) :: eType

element type, eT_ghostFromFiner, eT_ghostFromCoarser


Calls

proc~~mus_intp_error~~CallsGraph proc~mus_intp_error mus_intp_error proc~tem_baryofid tem_BaryOfId proc~mus_intp_error->proc~tem_baryofid proc~mus_open_config mus_open_config proc~mus_intp_error->proc~mus_open_config proc~newunit newunit proc~mus_intp_error->proc~newunit mpi_reduce mpi_reduce proc~mus_intp_error->mpi_reduce proc~set_momentindices set_momentIndices proc~mus_intp_error->proc~set_momentindices interface~tem_spatial_for tem_spatial_for proc~mus_intp_error->interface~tem_spatial_for proc~tem_open_distconf_array tem_open_distconf_array proc~mus_open_config->proc~tem_open_distconf_array proc~mus_create_funcstr mus_create_funcStr proc~mus_open_config->proc~mus_create_funcstr proc~open_config_chunk open_config_chunk proc~mus_open_config->proc~open_config_chunk

Called by

proc~~mus_intp_error~~CalledByGraph proc~mus_intp_error mus_intp_error proc~mus_testinterpolation mus_testInterpolation proc~mus_testinterpolation->proc~mus_intp_error proc~mus_init_flow mus_init_flow proc~mus_init_flow->proc~mus_testinterpolation proc~mus_initialize mus_initialize proc~mus_initialize->proc~mus_init_flow program~musubi musubi program~musubi->proc~mus_initialize

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: targetElem
integer, private :: iVal
integer, private :: iPos
integer, private :: iElem
integer, private :: nElems
integer, private :: globnElems
integer, private :: nStresses
integer, private :: indElem
integer, private :: fUnit
real(kind=rk), private :: targetPdf(layout%fStencil%QQ)
real(kind=rk), private :: targetMom(layout%fStencil%QQ)
integer, private :: t_target
integer, private :: iErr
integer, private :: nChunkELems
integer, private :: nDims
integer, private :: iVelMin
integer, private :: iVelMax
integer, private :: iStressMin
integer, private :: iStressMax
integer, private :: iPress
integer(kind=long_k), private :: targetID
integer(kind=long_k), private, allocatable:: minErrorID(:)
integer(kind=long_k), private, allocatable:: maxErrorID(:)
real(kind=rk), private :: targetVel(3)
real(kind=rk), private :: refVel(3)
real(kind=rk), private :: targetPress(1)
real(kind=rk), private :: refPress
real(kind=rk), private :: targetStrain(6)
type(flu_state), private, allocatable:: conf(:)
real(kind=rk), private, allocatable:: xc(:,:)
real(kind=rk), private, allocatable:: ux(:)
real(kind=rk), private, allocatable:: uy(:)
real(kind=rk), private, allocatable:: uz(:)
real(kind=rk), private, allocatable:: rho(:)
real(kind=rk), private, allocatable:: Sxx(:,:)
real(kind=rk), private, allocatable:: maxError(:)
real(kind=rk), private, allocatable:: minError(:)
real(kind=rk), private, allocatable:: error(:)
real(kind=rk), private, allocatable:: errorNorm(:)
real(kind=rk), private, allocatable:: globMaxError(:)
real(kind=rk), private, allocatable:: globminError(:)
real(kind=rk), private, allocatable:: globErrorNorm(:)
character(len=labelLen), private :: intpString
character(len=PathLen), private :: buffer
character(len=PathLen), private :: filename
character(len=PathLen), private :: header
logical, private :: file_exists
integer, private :: bound(2)
integer, private :: nSize
integer, private :: QQ