C
C
C
C
C
00 I=l, NP
IF ( ABS(DWR(I)) .GT. 0.0001)
END 00
RETURN
END IF
END 00
END IF
RETURN
END
B.6. Source Programs
341
DWC(I)=DWC(I)/DWR(I)
C***************************************************** ***************
C
C
Subroutine FDS forms source/sink data for each time step
C
C********************************************************************
C
*
*
C
C
C
C
C
C
SUBROUTINE FDS (TSUM,X,Y,IJM,SRAT,SCON,DSR,DSC,STEPS,
NSZON,NP,NE,NS,NTS,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION X(NP),Y(NP) ,IJM(3,NE),SRAT(NS,NTS) ,SCON(NS,NTS),
DSR(NP) ,DSC(NP) ,NSZON(NE),STEPS(NTS)
00 I=l, NP
DSR(I)=O.
DSC(I)=O.
END 00
IF (NS .GT. 0) THEN
00 K=l, NTS
IF (STEPS(K) .GT. TSUM) THEN
00 J=l, NE
MI=IJM(l, J)
MJ=IJM(2, J)
MK=IJM(3, J)
BI=Y (MJ) -Y (MK)
BJ=Y(MK) -Y(MI)
CI=X(MK)-X(MJ)
CJ=X(MI)-X(MK)
TRS=ABS(CJ*BI-BJ*CI)
R=SRAT(NSZON(J), K)
R=R*TRS/6.
DSR(MI)=DSR(MI) + R
DSR(MJ)=DSR(MJ) + R
DSR(MK)=DSR(MK) + R
IF (IHC .EQ. 2 ) THEN
C=SCON(NSZON(J), K)
DSC(MI)=DSC(MI) + C*R
DSC(MJ)=DSC(MJ) + C*R
DSC(MK)=DSC(MK) + C*R
END IF
END 00
00 I=l, NP
IF ( ABS(DSR(I)) .GT. 0.00001) DSC(I)=DSC(I)/DSR(I)