-
W. Spencer Smith authoredW. Spencer Smith authored
Code owners
Assign users and groups as approvers for specific file changes. Learn more.
pde_solver_constants.f90.svn-base 1.12 KiB
! ------------------------------------------------------------------------
!> \brief Module defining PDE Solver Constants
! ------------------------------------------------------------------------
MODULE pde_solver_constants
USE system_constants !> Global system constants
IMPLICIT NONE
INTEGER, PARAMETER :: NGAUSS_ELEM = 1 !> number of Gaussian integration points per body element
INTEGER, PARAMETER :: NGAUSS_BOUND = 1 !> number of Gaussian integration points per traction element
!> Gaussian integration data for body elements
DOUBLE PRECISION, PARAMETER, DIMENSION(NGAUSS_ELEM,2) :: GAUSS_PT_ELEM = RESHAPE( (/ ONE_THIRD, ONE_THIRD /), SHAPE(GAUSS_PT_ELEM) )
DOUBLE PRECISION, PARAMETER, DIMENSION(NGAUSS_ELEM) :: GAUSS_WT_ELEM = (/ 1.d0 /)
!> Gaussian integration data for traction elements
DOUBLE PRECISION, PARAMETER, DIMENSION(NGAUSS_BOUND) :: GAUSS_PT_BOUND = (/ 0.5d0 /)
DOUBLE PRECISION, PARAMETER, DIMENSION(NGAUSS_BOUND) :: GAUSS_WT_BOUND = (/ 1.d0 /)
!> parameters for Newmark time-stepping
DOUBLE PRECISION, PARAMETER :: GAMA = 0.5d0
DOUBLE PRECISION, PARAMETER :: BETA = 0.25d0
END MODULE pde_solver_constants