FORTRAN在罗氏法洪水计算中应用.docVIP

FORTRAN在罗氏法洪水计算中应用 　　摘要：以吉林省某流域为例，采用FORTRAN语言对罗氏法洪水计算进行程序化设计，并与手工试算成果进行对比分析，验证了该罗氏法洪水计算程序的正确性。与传统的手工试算方法相比，FORTRAN程序计算更加高效、精确。 　　（Taking a water basin in Jilin province as an example, applying FORTRAN program in Roche Method for Flood Calculation, then making comparative analysis on the results of program calculation and manual calculation. Compared to traditional manual calculation method, FORTRAN program calculation method is more effective and accurate. ） 　　关键词：罗氏法；洪水计算；FORTRAN；程序设计 　　（Roche Method; Flood Calculation; FORTRAN; Program Design） 　　中图分类号： S611 文献标识码： A 　　国内小流域设计洪水计算方法主要采用水科所法[1~2]和罗斯托莫夫法[3~5]。罗氏法在我省应用较广，但其在实际工作中往往采用手工试算的方法，需进行多次试算，比较麻烦。本文采用FORTRAN[6]对罗氏法洪水计算进行程序化设计，以提高运算效率。 　　1 基本公式 　　罗氏法是依据全面汇流理论推导得出，但没有考虑部分汇流的产峰情况，因而不宜用于过大的集水面积，适用于流域面积小于200km2的洪水计算。其洪峰流量值与洪水调查值作过验证，一般可采用此法计算。 　　 　　;;; 　　;; 　　;; 　　式中: 为设计流量, m3/s; 为径流系数; 为暴雨不均匀系数; 为面积, km2; 为频率对应的雨力, mm/min; 为频率对应的24h雨量, mm; 为暴雨强度, mm/min; 为暴雨递减指数; 总为汇流历时, s; 为河道汇流历时, s; 坡面汇流历时, s; 为土壤系数; 为坡面粗糙系数; 为河道长度, km; 为河道坡度, ‰; 为坡面长度, km; 为坡面坡度, %. 　　洪水计算方法不同，河长和河道坡度的计算方法也不相同。其中罗氏法的河道长度是从河源点开始计算，如果河源点不清楚，可采用量至分水岭的总长度扣除平均坡面长度作为河长。⑴当L500m时，按500m扣除；⑵当L10km时按1km扣除。 　　2 程序代码 　　罗氏法洪水计算程序可分为三个部分：⑴初始数据录入；⑵迭代计算；⑶结果输出。 　　2.1 数据录入部分 　　新建一个文本文件，将原始数据录入文本文件，并将数据读入。具体代码如下： 　　OPEN (5,FILE=DATA.DAT,STATUS=OLD) 　　READ (5,*) (P(K),F(K),L(K),J(K),LS(K),JS(K),XI(K),PHI(K),H24(K),N(K),K=1,6) 　　2.2 迭代计算部分 　　迭代原理：⑴先调整假定暴雨不均匀系数β0，使β0=β；⑵再调整假定汇流历时T0，使T0=T，当T0=T时，停止迭代。迭代计算代码如下： 　　DO K=1,6 　　 BETA0(K)=0.85 　　 T0(K)=180. 　　 SP(K)=H24(K)/1440**(1-N(K)) 　　 B(K)=SQRT(LS(K)/JS(K)**0.6/XI(K))/PHI(K) 　　100 A(K)=XI(K)*F(K)*BETA0(K) 　　 C(K)=L(K)/0.222/A(K)**0.2/J(K)**(0.24+1.6*J(K)) 　　 I(K)=SP(K)/T0(K)**N(K) 　　 ALPHA(K)=XI(K)*(I(K)+0.1)**0.345*T0(K)**0.15 　　 D(K)=1/I(K)/(I(K)+0.1)**0.345/T0(K)**0.15 　　 TR(K)=C(K)*D(K)**0.2 　　 TS(K)=B(K)*D(K)**0.5 　　 T(K)=TR(K)+TS(K) 　　 BETA(K)=1-0.28*(F(K)/T0(K))**0.4 　　 IF ((BETA0(K).NE.BETA(K)).AND.(T(K).NE.T0(K))) THEN 　　 BETA0(K)=BETA(K) 　　 GOTO 100 　　 ENDIF 　　 IF (