gaussian_points.f90

Go to the documentation of this file.

subroutine gaussian_points
!-----------------------------------------------------------------------
! calculate the sixteen functions at the sixteen gaussian points
!-----------------------------------------------------------------------
  
  use itm_types
  use mod_corners
  use mod_gauss
  
  implicit none
  
  real(r8) :: r0, s0
  real(r8) :: r, s
  real(r8) :: dumrr, dumss
  integer(itm_i4) :: irs, im, ngr, ngs
  integer(itm_i4) :: mi, mj
  
  do irs = 1, 4
    r0 = rs(irs, 1)
    s0 = rs(irs, 2)
    do im = 1, 4
      mi = ij(im, 1)
      mj = ij(im, 2)
      do ngr = 1, 4
        r = xgauss(ngr)
        do ngs = 1, 4
          s = xgauss(ngs)
          call cubic_polynomials(mi, mj, r0, s0, r, s, &
           h(im, irs, ngr, ngs), hr(im, irs, ngr, ngs), &
           hs(im, irs, ngr, ngs), hrs(im, irs, ngr, ngs), dumrr, &
           dumss)
        end do
      end do
    end do
  end do
  
  return
end subroutine gaussian_points