56 type(CoefficientsType) :: phi_face
58 class(CentralDifferenceSchemeType),
target :: this
59 integer(I4B),
intent(in) :: n
60 integer(I4B),
intent(in) :: m
61 integer(I4B),
intent(in) :: iposnm
63 real(DP) :: lnm, lmn, omega
67 if (this%dis%con%ihc(this%dis%con%jas(iposnm)) == 0)
then
69 lnm = dhalf * (this%dis%top(n) - this%dis%bot(n))
70 lmn = dhalf * (this%dis%top(m) - this%dis%bot(m))
73 lnm = this%dis%con%cl1(this%dis%con%jas(iposnm))
74 lmn = this%dis%con%cl2(this%dis%con%jas(iposnm))
77 omega = lmn / (lnm + lmn)
80 phi_face%c_m = done - omega