水文预报——竹溪坡流域洪水预报

最新推荐文章于 2021-12-12 21:46:52 发布

iSuwen

最新推荐文章于 2021-12-12 21:46:52 发布

阅读量3.4k

点赞数

分类专栏：爱进步 C++程序文章标签： c c++

本文链接：https://blog.csdn.net/Superwen_go/article/details/7629194

版权

这是一个C++程序，用于实现竹溪坡流域的洪水预报。程序中包含三层蒸发模式下的蓄流产流模型，根据不同的降雨和蒸发数据进行计算，并通过单位线数据进行汇流分析。程序涉及读取和写入数据文件，如'infile_P_Ep.txt'和'infile_q.txt'，并输出结果到'outfile_PE_R.txt'，'outfile_R_SIG.txt'和'outfile_Q_SIG.txt'。

摘要由CSDN通过智能技术生成

// 竹溪坡流域水文预报C++程序
//本流域中，不同次洪需要修改M值，其为次洪降雨起止间时段数
//本流域中，不同次洪需要在"infile_P_Ep.txt"文件中修改为对应的数据
//若单位线需要调整，在"infile_q.txt"文件中修改


#include "stdafx.h"
#include <fstream>
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;

 void main()
 {//5
	 const int M =44; //每个次洪的值都不一样，其为次洪降雨起止间时段数
	 const double FE = 0.8;
	 ifstream infile;
	 ofstream outfile;

	 
//**********三层蒸发模式下的蓄流产流模型(开始)************	
	 const double 	WUM = 30,	WLM = 80,	WDM = 40, WM = 150, C = 0.15, b = 1.5;
	 double WMM = WM*(1 + b); 
	 double P[M], Ep[M], EU[M], EL[M], ED[M], E[M], PE[M], WU[M + 1], WL[M + 1], WD[M + 1], W[M], a[M], R[M];
	 double af[M];//af[i]指产流面积比率(%),三水源划分中需要的数据
	 WU[0] = FE*WUM;
	 WL[0] = FE*WLM;
	 WD[0] = FE*WDM;

	 infile.open("infile_P_Ep.txt");//读人文件，次洪降雨、蒸发能力数据
	 for (int i = 0; i < M; i++)   
	 {
		 infile>>P[i]>>Ep[i];
	 }
	 infile.close();
 
	 outfile.open("outfile_PE_R.txt");//写出文件，次洪蓄满产流数据
	 outfile<<setw(10)<<"P"<<setw(10)<<"Ep"
			  <<setw(10)<<"EU"<<setw(10)<<"EL"<<setw(10)<<"ED"<<setw(10)<<"E"<<setw(10)<<"PE"
			 <<setw(10)<<"WU"<<setw(10)<<"WL"<<setw(10)<<"WD"<<setw(10)<<"W"
			 <<setw(10)<<"a" <<setw(10)<<"