tvdschememodule Module Reference

Data Types
interface	tvdschemetype

Functions/Subroutines
type(tvdschemetype) function	constructor (dis, fmi, ibound)
type(coefficientstype) function, target	compute (this, n, m, iposnm, phi)

Function/Subroutine Documentation

compute()

type(coefficientstype) function, target tvdschememodule::compute ( class(tvdschemetype), target  this, integer(i4b), intent(in)  n, integer(i4b), intent(in)  m, integer(i4b), intent(in)  iposnm, real(dp), dimension(:), intent(in)  phi  )

private

Definition at line 42 of file TVDScheme.f90.

    !-- return
    type(CoefficientsType), target :: phi_face
    ! -- dummy
    class(TVDSchemeType), target :: this
    integer(I4B), intent(in) :: n
    integer(I4B), intent(in) :: m
    integer(I4B), intent(in) :: iposnm
    real(DP), intent(in), dimension(:) :: phi
    ! -- local
    integer(I4B) :: ipos, isympos, iup, idn, i2up, j
    real(DP) :: qnm, qmax, qupj, elupdn, elup2up
    real(DP) :: smooth, cdiff, alimiter
    real(DP), pointer :: coef_up, coef_dn
    !
    ! -- Find upstream node
    isympos = this%dis%con%jas(iposnm)
    qnm = this%fmi%gwfflowja(iposnm)
    if (qnm > dzero) then
      ! -- positive flow into n means m is upstream
      iup = m
      idn = n
 
      coef_up => phi_face%c_m
      coef_dn => phi_face%c_n
    else
      iup = n
      idn = m
 
      coef_up => phi_face%c_n
      coef_dn => phi_face%c_m
    end if
    elupdn = this%dis%con%cl1(isympos) + this%dis%con%cl2(isympos)
    !
    ! -- Add low order terms
    coef_up = done
    !
    ! -- Add high order terms
    !
    ! -- Find second node upstream to iup
    i2up = 0
    qmax = dzero
    do ipos = this%dis%con%ia(iup) + 1, this%dis%con%ia(iup + 1) - 1
      j = this%dis%con%ja(ipos)
      if (this%ibound(j) == 0) cycle
      qupj = this%fmi%gwfflowja(ipos)
      isympos = this%dis%con%jas(ipos)
      if (qupj > qmax) then
        qmax = qupj
        i2up = j
        elup2up = this%dis%con%cl1(isympos) + this%dis%con%cl2(isympos)
      end if
    end do
    !
    ! -- Calculate flux limiting term
    if (i2up > 0) then
      smooth = dzero
      cdiff = abs(phi(idn) - phi(iup))
      if (cdiff > dprec) then
        smooth = (phi(iup) - phi(i2up)) / elup2up * &
                 elupdn / (phi(idn) - phi(iup))
      end if
      if (smooth > dzero) then
        alimiter = dtwo * smooth / (done + smooth)
        phi_face%rhs = -dhalf * alimiter * (phi(idn) - phi(iup))
      end if
    end if
 

constructor()

type(tvdschemetype) function tvdschememodule::constructor ( class(disbasetype), intent(in), pointer  dis, type(tspfmitype), intent(in), pointer  fmi, integer(i4b), dimension(:), intent(in), pointer, contiguous  ibound  )

private

Definition at line 28 of file TVDScheme.f90.

    ! -- return
    type(TVDSchemeType) :: interpolation_scheme
    ! --dummy
    class(DisBaseType), pointer, intent(in) :: dis
    type(TspFmiType), pointer, intent(in) :: fmi
    integer(I4B), dimension(:), pointer, contiguous, intent(in) :: ibound
 
    interpolation_scheme%dis => dis
    interpolation_scheme%fmi => fmi
    interpolation_scheme%ibound => ibound
 

Here is the caller graph for this function: