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)

Második rész. A számítási eljárások gyakorlati alkalmazása - 8. Implicit eljárások alkalmazása

ERRCONT azonosítójú szubrutin 15. és 16. sorát). Iteráció esetén a JEL = 1 és az ITER azonosítójú logikai változó TRUE értéket kap (17 —18. sorok) és ezzel a szubrutin lejár.

Most térjünk vissza a főprogram 38. sorára. Ha az iterációk száma eléri a 100-at akkor, mivel IND — 1, az adott időciklus lejárt. Ha viszont az iterácó szükséges (JEL = 1) akkor a végrehajtás visszacsatolódik a főprogram 33-as jelű címkéjére, ahonnan a számítás a határfeltételek számításától kezdve megismétlődik.

Ha egy adott időcikluson belül a számítás pontossága kielégítő (IND = 0 és JEL = 0), akkor a program végrehajtása a főprogram 40. sorától folytatódik, ahol a 08.27. mellékleten megadott és a számítások eredményeit kinyomtató PRINTRES azonosítójú szubrutin kerül aktivizálásra. A szubrutin 6 — 10. sorában, valamennyi szelvényben a főmeder, a hullámtéri és az egész kereszt- szelvényre vonatkozó középsebesség kerül kiszámításra minden egyes aktivizáláskor.

A PRINTRES szubrutin 11 — 25. sora az eredmények rendezését szervezi és a nyomtatási részt (16., 25. sorok) csak akkor aktivizálja, ha a KEK változó értéke 0, vagyis 4-időciklus eredménye gyűlt össze a megelőző nyomtatás óta. Mivel DT = 900, minden kerek óra számítási eredménye kerül kinyomtatásra, minden második szelvényben.

A 08.28. mellékleten a PRINTRES azonosítójú szubrutin nyomtatási sémáját adtuk meg 61 db szelvényre vonatkozó részeredményekkel, a T = 1 órára vonatkozóan. Az egyes oszlopok jelentése sorrendben a következő:

X a szelvény koordinátája (felülről lefelé növekedve),

QE a főmeder vízhozama,

QH a kétoldali hullámtér eggyüttes vízhozama

Z a vízszint,

HF a főmeder vízmélysége,

HH a hullámtér vízmélysége,

QL a lineáris vízhozam,

VF a főmeder,

VH a hullámtér és

VK a teljes szelvény középsebességei.

8.4.4. A SZÁMÍTÁSI EREDMÉNYEK FELDOLGOZÁSA

A számítások eredményeiből a 8.4—1. ábrán felraktuk a h = h(t)x árhullámképeket néhány (x — 0, 20, 40, 60 km) szelvényben. Ugyancsak felraktuk a felső határfeltételi szelvény Q = Q(t)x=0 vízhozamának időmenti változását is.

Emlékeztetünk arra, hogy ennél a változatnál a hullámtér simasági tényezője kh — 5 volt. Az erre vonatkozó görbéket 8.4—1. ábrán szaggatott vonalakkal ábrázoltuk.

Eredeti feladatunkra visszatérve (lásd 8.4.1. alatti részt) el kell végezni a számítást kh = 20-as hullámtéri simasági tényező esetére is. Ez a változás mindössze azt jelenti, hogy a 08.19. mellékleten megadott főprogram 10. sorában levő beolvasási utasításban K2 értékére most 20.0-at kell beolvastatni, és a programot újból lefuttatni. Ezt is elvégeztük, és a számítások eredményeiből a h — h(t)x görbéket a 8.4—1. ábrán szemléltetjük (folytonos vonallal).

290