342
Appendix B
END DO
C
RETURN
C
END IF
C
END DO
C
END IF
C
RETURN
END
C
C********************************************************************
C
C
Subroutine FVD generates veloeities and dispersion eoeffieients.
C
C********************************************************************
C
C
SUBROUTINE FVD (X,Y,IJM,TK,PR,AL,AT,AD,DFU,DXX,DXY,DYY,VX,VY,
HO,NP,NE,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION X(NP) ,Y(NP) ,IJM(3,NE) ,TK(NE) ,PR(NE),
1
AL(NE) ,AT(NE) ,AD(NE) ,DXX(NE) ,DXY(NE) ,DYY(NE) ,
2
VX(NE) ,VY(NE) ,HO(NP)
C For flow in an anisotropie aquifer
C
C
IF (IHC .EQ. 1) THEN
DO N=l, NE
DXX(N)=AL(N)
DYY(N)=AT(N)
DXY(N)=AD(N)
END DO
ELSE
DO N=l, NE
MI=IJM(l,N)
MJ=IJM(2,N)
MK=IJM(3,N)
BI=Y(MJ) -Y(MK)
BJ=Y(MK) -Y(MI)
BK=Y(MI) -Y(MJ)
CI=X(MK)-X(MJ)
CJ=X(MI) -X(MK)
CK=X (MJ) -x (MI)
TRS=ABS(CJ*BI-BJ*CI)
C Caleulating average values of parameters in eaeh element
C
A=-TK(N)/PR(N)/TRS
RET=l. + (l.-PR(N»*AD(N)/PR(N)
C
C Calculating the veloeity distribution
C
C
VX(N)=A*(BI*HO(MI)+BJ*HO(MJ)+BK*HO(MK»/RET
VY(N)=A*(CI*HO(MI)+CJ*HO(MJ)+CK*HO(MK»/RET
V=SQRT(VX(N)*VX(N)+VY(N)*VY(N»
IF (V .LT. 1.E-10) V=1.E-10
C Caleulating dispersion eoeffieients