Kozák Miklós: A szabadfelszínű nempermanens vízmozgások számítása digitális számítógépek felhasználásával (Akadémiai Kiadó, Budapest, 1977)
Mellékletek. Programok, szubrutinok, nyomtatási sémák 313 - 4. fejezethez
04.02, EGYENLETÁLLANDÓK SZÁMÍTÁSÁNAK SZUBRUTINJA, IMPLICIT MÓDSZERHEZ (EQCONIMP) LN 0001 SU3 ROUTINE E OCONIH? (A » B ,SE C, DX, OT, E ,CC, L 2) LN 0002 INTEGER SEC«SZM 1,P LN 0003 DIMENSION A (L2,5) ,B(SEC, 28) LN 0004 SZM 1=SEC-1 LN 0005 OO 97 1=1,SZM1 LN 0006 M=2* (1-1) *1 LN 0007 DX= 8 (1*1,1) -B (I, 1) L N 0008 W=OX/OT LN 0009 A(M,1)=-1 LN 0010 A (N,2) =9*8(1,17) LN 0011 A(M , 3)=1 LN 0012 A (H , 4) =W*B( 1*1, 17) LN 0013 A(H,5)=*(B(IH,9)-8tI,9))«’(B(I,ll)»8(I,8)»B(IH,ll)*0(I*l, LN 0014 1( B( I , 10) *B( I, 16) ) *OX L N 0015 M=M * 1 LN 0016 C1=A9S(B(I,9))/B(I,13) LN 0017 C2=A9S(3(I*1,9))/B(I*1«13) L N 0018 C3=A8S (8(1, 15) )/B(I, 19) L N 0019 C4=A8S(B(I* 1,15) )/B(I*l, 19) LN 0020 AC= 0 • LN 0021 8 C= 0 • LN 0022 DC=0. LN 0023 R=0 LN 0024 P=0 LN 0025 OO 53 K= P, 6 »6 LN 0026 M1=9*K LN 0027 M2= 12*K LN 0028 AC= AC*B(I,M1)/B(I,M2)**2*B( I«-1, Ml) / 8 (I *1 , F2> **2 LN 0029 BC=BC*1. /B( I, M2) U./B(I*1,H2) LN 0030 5 3 0C=0C*8(I,M1)**2/B(I,M2) **3 *B (1*1, M 1) **2/8( I »1, M2) ** 3 L N OOSi A (M, 1) =C3*DX/2.*E*BC-CC*AC LN 0032 A (M, 2)=-l+CC*B (1,11)*OC L N 0033 A (M, 3) =C4*0X/2.*E*8C*Ce*AC L N 0034 A (M, 4) =1-CC *8(1*1,11 >*OC L N 0035 A (M ,5) =- (8(1*1,8)-B( I, 8) > *CC*OC* (2. * (B (I * 1,12) -8(1,12) > ♦ ( B (I LN 0036 18(I,11)-B(I*1,8)*B(I*1,11)) )*E*BC»(8(I,9) *8 (1*1,9) ) - (C1»B (I, L N 0037 2C2*B (I *1,9) )*OX/2.-CC*AC»(8 (I *1, 9) -B (1,9 ) ) LN 0038 97 CONTINUE LN 00 39 RETURN LN 0040 ENO USASI FORTRAN DIAGNOSTIC RESULTS FOR ECCONIMP NO ERRORS 329
