Seite 1 Institut für Hydraulik und landeskulturelle Wasserwirtschaft ...

HYDRODYNAMIK SEITE 60PROGRAM FLOODC DATE 88-12-08C COMPUTATION OF A FLOOD WAVE ALONG A CHANNELC PARAMETER LIST:C C C ACCORDING TO HCC HC DEPTH FOR C-VALUESC CC INTERPOLATED C-VALUESC DX LENGTHSTEPSC DT TIMESTEPCH(I,J) DEPTH IN EACH GRIDPOINTC HH INTERPOLATED DEPTH - BOUNDARYC HT DEPTH - BOUNDARYC TH TIME - BOUNDARYC IT OUTPUT FOR SELECTED TIMESTEPSC K NUMBER OF VALUES FOR BOUNDARYC L NUMBER OF C-VALUES FOR INTERPOLATIONC M NUMBER OF TIMESTEPSC N NUMBER OF DXC R R=CC*DT/DXC RL RIVERLENGTHCRMAX MAXIMUM R - VALUECRMIN MINIMUM R - VALUEC RR RIVERRIDGECTIME TIMERANGE CONSIDEREDCTSTEP TIMESTEPSCPARAMETER (IMAX=100,ITAB=50)REALHH,H(0:IMAX,0:ITAB),C(ITAB),HC(ITAB),HT(ITAB),TH(ITAB),RLREAL CC,DX,DT,TIME,R,RR(0:IMAX),TSTEP(0:IMAX),RMAX,RMINCHARACTER LINE*6INTEGER I,J,K,L,N,M,ITOPEN(10,FILE='FLODAT')OPEN(11,FILE='FLORES')C -------------------------- READ DATA ------------------READ(10,*) N,M,RL,TIME,ITREAD(10,*) KREAD(10,*) (TH(I),I=1,K)READ(10,*) (HT(I),I=1,K)READ(10,*) LREAD(10,*) ( C(I),I=1,L)READ(10,*) (HC(I),I=1,L)C -------------------------- INITIALIZE -----------------LINE = '------'DX = RL/NDT = TIME/MRR(0) = 0.H(0,0)= HT(1)TSTEP(0) = 0.RMAX = 0.RMIN = 1.C -------------------------- BOUNDARY VALUES RIVERRIDGE ---DO 20 I=1,NRR(I)=RR(I-1)+DXH(0,I)=HT(1)Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.0