Hidrológiai Közlöny 1987 (67. évfolyam)

2-3. szám - Gauzer Balázs: Néhány FORTRAN szubrutin hidrológiai adatok elsődleges feldolgozásához

GAUZER B.: Néhány FORTRAN szubrutin 155 $STORAGE:2 SUBROUTINE SPLINE(N,RC) DIMENSION RC(1) , AA(33), A(33), AB(33) DATA A(1) /2./» AB(l) /0./ C IF(N-2) 1, 7, 2 1 STOP 'STOP S/P SPLINE' 2 NB=N-1 1 = 4 DO 3 K=2,

NB DA=RC(I)-RC( 1-3) DB=RC(1+3)-RC(I) DD=DA+DB AA(K)=DA/DD A (K)=2. AB(K)=DB/DD RC( I + 2 ) = 6. *( (RC(I+4)-RC(I+l))/DB(RC(I+l)-RC(I-2))/DA)/DD 3 1=1+3 A (N)=2. AA(N)=0. RC(3)=0. RC(1+2)=0. 1 = 1 DO 4 K=2,N 1 = 1 + 3 DD=AA(K)/A(K-1) A(K)=A(K)-DD*AB(K-1) 4 RC( I

4 )=RC(I + 2)-DD*RC(

I-1) RC(I '2)=RC(1+2)/A(N) K = N 5 K=K-1 IF(K) 100, 100, 6 6 1=1-3 RC(1+2 ) = (RC(1 + 2)-AB(K)+RC(1 + 5))/A(K, GOTO 5 7 RC(3)=0. RC(6)=0. 100 RETURN END $STORAGE:2 FUNCTION FSPLN(X,N,RC) DIMENSION RC(l) C 1 = 1 DO 1 K = 2,N 1 = 1 + 3 IF(X-RCU)) 2, 2, 1 1 CONTINUE 2 D =RC(I)-RC(1-3) DA=RC(

I)-X DB=D-DA DD=D*D FSPLN=((RC(I-1)*(DA*(DA*DA-DD))+RC(1+2)* (DB*(DB*DB-DD)))/6.+ / PC(I-2)*DA+RC(I+1)*DB)/D C RETURN END 10. ábra. A SPLINE és az FSPLN szubrutin A SPLINE szubrutin A SPLINE függvény paramétereit számítja (10.ábra). Paraméterek: N — Az adatsor pontjainak a száma (max. 34) —INP— RC — A SPLINE függvény paramétereit tartalmazó 3 dimenziós tömb. Tartalma: RC(1, I) — a független változó az I-edik pontban —INP— RC(2, I) — a függő változó az I-edik pontban —INP— RC(3, I) — a függvény görbülete az I-edik pontban — OUT— Hívott szubrutin: NINCS. Az FSPLN függvény Számítja a SPLINE függvény értékét. Paraméterek: X — A függvény független változója —INP— N — Az adatsor pontjainak a száma (max. 34) RC — A SPLINE függvény paramétereit tartalmazó kétdimenziós tömb. —INP— Hívott szubrutin: NINCS. Irodalom 1. Futaki K., Vargay Z. 1984. Primary processing of hydrological data oil microcomputers. In: Microprocessors in Operational Hydrology, pp. 113—127. Reidel Publishing, Dordrecht. 2. Ganzer B. 1985. FORTRAN nyelvű programcsomag vízminőségi adatok elsődleges feldolgozása, Kézirat VITUKI, Budapest. 3. VITUKI 1985. Flood innundation forecasting HOMS sequence, Kézirat, Budapest. Kézirat beérkezett: 1986. szeptember 26. Átdolgozás beérkezett: 1987. március 15. Közlésre elfogadva: 1987. március 18.