346
Appendix B
C
C
C
CFA(NI+I)=CFA(NI+I)+VV
ENDIF
CFB(NI+I)=CFB(NI+I)
ENDIF
END DO
DO I=l,NS
IF(MCPN(I,MI) .EQ. MK) THEN
CFA(NI+I)=CFA(NI+I)
1
+((D11*BBJK+D12*CCJK)*(BI-BJ)
2
+ (D12*BBJK+D22*CCJK)* (CI-CJ)
3
+ (D11*BBKK+D12*CCKK)*(BI-BK)
4
+ (D12*BBKK+D22*CCKK)*(CI-CK))/(2.*TRS)
IF (IHC .EQ. 2) THEN
VV=(VX(N)*(BBJK+BBKK)+VY(N)*(CCJK+CCKK))/4.
CFA(NI+I) =CFA(NI+I) +VV
ENDIF
CFB(NI+I)=CFB(NI+I)
ENDIF
END DO
J=ME(l)
ME(1)=ME(2)
ME(2)=ME(3)
ME(3)=J
END DO
END DO
DO 1=1, NP
PSE(I)=PSE(I)/PS(I)
PPR(I)=PPR(I)/PS(I)
PAD(I)=PAD(I)/PS(I)
END DO
RETURN
END
C***************************************************** ***************
C
C
Subroutine WEIT generates adaptive upstream weights.
C
C********************************************************************
C
C
SUBROUTINE WEIT (BI,BJ,BK,CI,CJ,CK,DXX,VX,VY,WI,WJ,WK,PT)
IMPLICIT REAL*8 (A-H, O-Z)
VIJ=CK*VX-BK*VY
VIK=BJ*VY-CJ*VX
VKJ=BI*VY-CI*VX
IF (DXX .GT. 1.E-20) THEN
WI=( VIJ+VIK)/DXX
WJ= (-VIJ -VKJ) /DXX
WK=(-VIK+VKJ)/DXX
ELSE
WI=O.
WJ=O.
WK=O.
ENDIF
WI=1./3.+PT*WI
WJ=1.!3. +PT*WJ
WK=l.!3. +PT*WK
RETURN
END