41 type(CoefficientsType) :: phi_face
43 class(CentralDifferenceSchemeType),
target :: this
44 integer(I4B),
intent(in) :: n
45 integer(I4B),
intent(in) :: m
46 integer(I4B),
intent(in) :: iposnm
47 real(DP),
intent(in),
dimension(:) :: phi
49 real(DP) :: lnm, lmn, omega
53 if (this%dis%con%ihc(this%dis%con%jas(iposnm)) == 0)
then
55 lnm = dhalf * (this%dis%top(n) - this%dis%bot(n))
56 lmn = dhalf * (this%dis%top(m) - this%dis%bot(m))
59 lnm = this%dis%con%cl1(this%dis%con%jas(iposnm))
60 lmn = this%dis%con%cl2(this%dis%con%jas(iposnm))
63 omega = lmn / (lnm + lmn)
66 phi_face%c_m = done - omega