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 - 5. fejezethez
05.02. EGYENLETÁLLANDÓK SZÁMÍTÁSÁNAK SZUBRUTINJA, FÉLDINAMIKUS MÓDSZERHEZ (EQCONSTP) ln ooai LN 0002 C LN 0003 C LN 00 0<» LN 0005 LN 0006 LN 0007 LN 0006 LN 0009 LN 0010 LN 0011 LN 0012 LN 0013 LN 00H* LN 0015 LN 0016 LN 0017 LN 0018 LN 0019 LN 0020 LN 0021 LN 0022 LN 0023 LN 002A- LN 0025 LN 0026 LN 00 27 LN 0028 LN 0029 LN 0030 LN 0031 LN 00 32 LN 0033 LN 003J» SUBROUTINE EQCONSTP ÍA , B ,Q, V,Z,S2,Bl,32, CC,W,DX ,TX , TB,BS,E,L2) 05.02 EGYENLET ALL ANDOK SZAMITASA COMPUTATION OF EQUTION CONSTANT INTEGER SZ» SZM1, BS , P ,R LOGIOAL BREAK DIMENSION A (L2, 8) , BCSZ ,33) COMMON /10/ ITER,BREAK B (1 ,8)-Q 8(1,9)=V B(SZ,10)=Z SZM1 = SZ-1 DO 97 I=1,SZH1 M=A*(I-13+1 A(M,2>=0. A (M ,6) =0 • A (M» 1,21=0. A (H +i,4) =0. A(M+1,6J=0. A (M , 1) =“1. A (M ,5) =1. A(M,3)=B1*W A(M ,7)=B1*W A(M,<»)=-2*DX A (M,B)=AIM,3J*(B (I, 6 )*B < I *1 ,6) ) FB (I, A)-B (I + l»«,) IF( BREAK. AN D.TX.GE.TB. ANt.I.EQ.BS) A (M, 8) =A (M, 8) - (B (1,29) +6 (I, 30 ) 1) *PX M=M+1 AC= 0 BC= 0 0C=0 FO 0 R=0 P=0 DO 53 K= P,7,7 334
