AverageSides_Vector3D_t Subroutine

public subroutine AverageSides_Vector3D_t(this)

Arguments

TypeIntentOptionalAttributesName
class(Vector3D_t), intent(inout) :: this

Contents


Source Code

  subroutine AverageSides_Vector3D_t(this)
    implicit none
    class(Vector3D_t),intent(inout) :: this
    ! Local
    integer :: iel
    integer :: iside
    integer :: ivar
    integer :: i,j
    integer :: idir

    do concurrent(i=1:this%N+1,j=1:this%N+1, &
                  iside=1:6,iel=1:this%nelem,ivar=1:this%nvar, &
                  idir=1:3)
      this%boundary(i,j,iside,iel,ivar,idir) = 0.5_prec*( &
                                               this%boundary(i,j,iside,iel,ivar,idir)+ &
                                               this%extBoundary(i,j,iside,iel,ivar,idir))
    enddo

  endsubroutine AverageSides_Vector3D_t