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 - 6. fejezethez
06.07. A KEZDETI FELTÉTELT SZÁMÍTÓ PROGRAM ÖSSZETETT PRIZMA. TIKUS ÉS SZIMMETRIKUS HULLÁMTERŰ MEDREKRE L N 0001 C KEZDETI FEL TETEL, PRIZMATIKUS, 0E3SZETETT MEDERRE L N 3002 C COMPUTATION OF INITIAL CONDITION IN COMPOLNO CROSSECTION LN 0003 INTEGER SEC ,0,1 ,N,L,F '.N 03 09 REAL sl.zo,dx,kf»kh,bf,bh,dz»ozh,zl,QERR,ES,K, knf,KNH.KKF.KKH, I.N 00 05 DIMENSION A <27, 7) LN 0006 READ (6 0, 101 )SEC,D,0X,SL ,Z0,KKF,KKH, 3F, BH , ZL ,QERR , A (1,1) ,A (1,2) LN 0007 1A (1,3) ,A (1,9) LN 00 08 101 F0RMAT(I3»I3,F8.1,F9.5,F7.2,F6.1,F6.i»F6.1»F6,l,F6, 2 .F8.9/9F7 • L N 00 09 RO = 0. 0 L N 0010 ES=DX*SL LN 0011 OO 1 1=2,SEC LN 0012 A (I ,1) =A(I-1,1)+DX LN 0013 A (1,2) =A (1-1,2) +-ES L N 0019 A (I , 3) = A (1-1,3) +ES LN 3015 1 A (I ,9) =A (1-1,9) LN 0016 A (1,6) =ZD LN 0017 A (1,7)=0.0 L N 0018 0Z=ZD-A(1,2) LN 0019 OZH = ZD-A (1, 3) LN 0020 KF=BF»DZ*KKF»(8F*DZ/ (BF + 2*0Z) )** (2/3) LN 0021 KH= BH*DZH*KKH» (8H»DZH/ (8H+2*0ZH) ) ** (2/3) LN 0022 K=( KF + KH) /2 LN 0023 C SZAKASZOK SZÁMIT ASA LN 3029 C COM PUTION OF STRECHER LN 0025 DO 2 N= 2,SEC LN 0026 F=0 LN 0027 L =0 LN 0028 A (N ,6) =A (N-1,6) +ZL LN 0029 3 DZ= A (N ,6 ) -A (N, 2 ) L N 0030 DZH = A (N, 6) -A (N, 3) LN 00 31 IF( OZH.L E.O .0) GO TO 9 LN 0032 KNF = BF»OZ*KKF«(BF*0Z/(8F +2»DZ))**(2/3) L N 00 33 KNH = BH*0ZH*KKH* (8H*DZH/( BH(-2*DZN) ) ** (2/3) LN 00 39 SO= (K* (KNF + KNH) /2) »SORT ( (A(N,6) - A (N-l,6) )/DX) LN 0 0 35 IF(ABS(SQ-A (N,9) )/A(N,9) .LT.QERR) GO TO 5 L N 0 0 36 IF(SQ-A(N,9))6,7,7 L N 00 37 6 IF(F.EQ. 1) ZL-ZL/2 LN 0038 A (N ,6) =A (N, 6) +ZL LN 0039 L =1 LN C3N0 GO TO 3 LN 00 LI 7 IF(L.EQ.l) ZL=ZL/2 LN 0092 A (N, 6) =A (N, 6) -ZL LN 0093 F=1 LN 0099 GO TO 3 LN 0095 5 ES=SQRT( (A(N,6)-A(N-1*6) >/OX) LN 0096 A (N,9)=(KF«-KNF)/2*ES L N 0097 A (N,5) = (KHtKNH)/2*ES LN 0098 A( N,7>= 0.0 LN 0099 KF= KNF LN 0050 KH=KNH LN 0051 K = (KF*KH)/2 LN 0052 2 ZL= A (N ,6) -A (N-1,6) LN 0053 9 A (1,9) =A(2,9) LN 0059 A (1, 5)=A (2, 5) LN 0055 WRITE (61,102) BF.BH.RO, KKF.KKH LN 005b 10 2 F0RMATÍ1H1///36H HYDRAULIC PARAMETER OF COMPOUND BED//2X, 2H0F, LN 0057 12H8H.9X, 2HR0.9X, 2HKF.9X, 2HKH/F5. 1,F8. 1.3F6.1) LN 0058 WRITE (61,10 3) LN 0059 103 FORMAT (1 HO, 3X,9HX(M) ,5X, 6HZ0F (M),9X, 6HZ0H(M),6X,7HQ<M3/S) , LN 0060 19X, 7HV(M3/S),5X,9HZ(M) ,5X,2HQL) LN 0061 HRITE(61,109) ( (A(I,N),N=1,0) ,1=1,SEC) LN 0062 109 FORMAT (1H0, F8.1,2F10.3,2F11.3,FI3. 9,F6.1 > L N 0063 END 06.07 5.6- 1X 5.6- 2X 5.6- 1X 5.6- 2X b.2-11 6.2-10 5.6- 1 5.6- 2 USASI FORTRAN DIAGNOSTIC RESULTS FOR FTN. MAIN NO ERRORS
