天然河道计算.docVIP

题目： 某天然河道有四个水文站，如图所示，测得四个断面的水文资料（列于表1）。四个断面处的局部阻力系数分别为0，0，0.3，0.3；四个断面之间的距离为：△L1-2=35000m，△L2-3=33500m，△L3-4=31500m。当河道通过流量Q=3800m3/s时，测得站4处闸前水位Z4=39.27m，求上游断面3、2、1处的水位。 变 量 说 明 M----------------------------------------断面总数 N-----------------------------------------每个断面内已知水文资料总数 Q-----------------------------------------流量 Z4----------------------------------------断面4的已知水位 EPS--------------------------------------允许误差 SN---------------------------------------粗糙系数 G-----------------------------------------重力加速度 SS(I)-------------------------------------局部阻力系数 SL(N)------------------------------------各断面间的距离 A(I)--------------------------------------过水断面面积 B(I)--------------------------------------水面宽度 R(I)----------------------------------------水力半径 Z(I)---------------------------------------水位 Z(M)-------------------------------------断面4的水位 A2(J),B2(J),Z2(J)-------------------------------------中间变量,断面上已知水力要素 A0(I,J)-----------------------------------断面内已知水位的过水断面面积 B0(I,J)-----------------------------------断面内已知水位的水面宽度 Z0(I,J)------------------------------------断面内已知水位 E2------------------------------------------断面2的能量函数 G1------------------------------------------断面1的能量函数 L----------------------------------------迭代次数 Z1-------------------------------------------计算段上游断面水位 X(N),Y(N)---------------------------------拉格朗日插值公式中的已点值 R1,R2,R3-----------------------------------拉格朗日插值公式中的中间变量 Q1,Q2,Q3----------------------------------拉格朗日插值公式中的中间量 源程序 PARAMETER (Q=3800.0, M=4,N=3,Z4=39.27,EPS=0.01,SN=0.025,G=9.8) DIMENSION SL(3),SS(4),A(4),B(4),R(4),Z2(3),A2(3) DIMENSION B2(3),Z(4),Z0(4,3),A0(4,3),B0(4,3) OPEN(10,FILE=HMX.DAT) OPEN(10,FILE=HMX.OUT) DO 100 I=1,4 READ(10,*) (A0(I,J),J=1,3) 100 CONTINUE DO 200 I=1,4 READ(10,*) (B0(I,J),J=1,3) 200 CONTINUE DO 300 I=1,4 READ(10,*) (Z0(I,J),J=1,3) CONTINUE READ(10,*)(SL(I),I=1,3) READ(10,*)(SS(I),I=1,3) WRITE(20,1100) 1100 FORMAT(/5X,计算结果:/) WRITE(20,12