水文分析与计算——年均流量及年极值洪量

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
//年均流量及年极值洪量.h
void NianJunZhi()
{
using namespace std;
ofstream outfile;
outfile.open("outfile_YearQ_W.txt");
outfile<<setw(10)<<"年份"<<setw(10)<<"天数"<<setw(10)<<"年均流量";
for(int k = 0; k < J; k++)
outfile<<setw(7)<<"最大"<<NumJ[k]<<"日";
outfile<<endl;
ifstream infile;
infile.open("infile_DayQ.txt");
for(int i = 0; i <
Y; i++)
{
Year = StartYear + i;
SumDay[i] =
Sum[12];//初始化
if(Year == StartYear)
SumDay[i] -= (Sum[StartMonth -1] + StartDate - 1);
if(Year == EndYear)
SumDay[i] = Sum[EndMonth - 1]+ EndDate;
if(Year%400==0||(Year%4==0&&Year%100!=0))
{
if (Year == StartYear&&StartMonth< 2) SumDay[i]++;
if(Year == EndYear&&EndMonth > 2) SumDay[i]++;
if(Year!=StartYear&&Year!=EndYear) //&&SumDay[i] == Sum[12]
SumDay[i]++;//闰年且含2月份加一天
}
YearQ[i] = 0;//初始化
for(int j = 0; j < SumDay[i]; j++)
{
infile>>Q[j];
YearQ[i] += Q[j];
}
YearQ[i] /= SumDay[i];
outfile<<setw(10)<<Year<<setw(10)<<SumDay[i]<<setw(10)<<YearQ[i];
for(int k = 0; k<J; k++)
{//求年极值洪量
MaxYearW[i][k] = 0;//初始化
for(int j = 0; j < SumDay[i] - NumJ[k] + 1; j++)
{
temp_W = 0;//初始化
for(int n = 0; n < NumJ[k]; n++)
temp_W += Q[j + n];
if(temp_W > MaxYearW[i][k])
MaxYearW[i][k] = temp_W;//(m3/s*月)
}
outfile<<setw(10)<<MaxYearW[i][k];
}
outfile<<endl;
}
infile.close();
outfile.close();
}