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 - 8. fejezethez

.02. A MAGÁNYOS HULLÁM LEVONULÁSÁT IMPLICIT MÓDSZERREL SZÁMÍTÓ FŐPROGRAM A 8.1. PÉLDÁHOZ

I 0001 I 0002 I 0003 I 0004 ^ 0005 4 0006 4 00 07 ^ 0008 Y 00 09 N 0010 4 0011 N 0012 N 0013 N 0014 N 0015 N 0016 N 0017 N 0018 N 0019 N 0020 N 0021 N 0022 N 0023 N 0024

C NEMPERMANENS VIZMOZGAS FOLYÓKBAN IMPLICIT MOOSZER

C UNSTEADY FLOH IN RIVERS IMPLICIT METHOD

INTEGER SEC ,P»R,TD»REK,OAY REAL K1

DIMENSION Bt51,28),A(100,5) ,C (51,12 > ,0 <2 ,4)

LOGICAL ITER EXTERNAL RESPRINT COMMON 710/ ITER

CALL DATA (8 ,SEC, OAY, Bl ,K 1,0H, OX, OT, RO,ZI,REK.TM , QERR,Z0i, Kll, ILII ,QL,CC,H ,E,TO,Rl, JM.a.KEK, OZ)

CALL CINPAR (B,ZO,SEC)

IF(REK.LT«5) GO TO 31 CALL OATAPRIN(B,SEC)

31 CONTINUE

C DIRECTION OF PROGRAM BY TIME DO 170 J = 1,JM TX= (J-l) *TO ITER=.FALSE.

KC= 0

CALL BOUNOC ON (Q,QR,TX,TM,Z, ZD)

33 CONTINUE

CALL SECPARIM(B,SEC,0Z)

8 (1» 15 )*QR 8 (1,26) =QR

.N 0025 .N 0026 .N 0027 .N 0028 .N 0029 .N 0030 -N 0031 lN 0032 .N 0033 LN 0034 LN 0035 LN 0036 LN 0037 LN 0038 LN 0039 LN 0040 LN 0041 LN 0042 LN 0043 LN 0044 LN 0045 LN 0046 LN 0047 LN 0048 LN 0049 LN 0050

301 B (S EO,14)=Z O L2=2*(SEC-1>

CALL EQCONI MP(A,B,SEC,OX,DT ,E,CC,L2)

CALL MATCORR4(B»A,SEC,L2)

CALL EQS0L4 (A,L2)

,C DISPOSITION OF RESULT

BCl ,14» = A (1 ,5)

304 B(SEC,26)=A(L2,5)

Nl= SEC-1

DO 161 1=2, N1

M=2*(I-1)

B (I , 26) = A (M ,5 )

161 B(I,14) = A(MH,5)

C HIBAELLENORZES

C CONTROL! OF ERROR

DO 162 1 = 1, N1

0Q=ABS((8(I,26)-B(I,15) ) / B( 1,151 )

IF(OQ-OERR) 162, 162,164 164 KC=KC*1

ITER=.TRUE.

GO TO 33

162 CONTINUE

CALL RESORD (B ,C, O, KEK.TX ,KC ,K 11, L11, SEC, RESPRIN T) 170 CONTINUE STOP ENO

0 6.92

08.01

08.04

08.03

04.01

04.02 0 8.07 04*03

08.08

USASI FORTRAN DIAGNOSTIC RESULTS FOR FTN.MAIN

NO ERRORS

375