水文分析与计算——年均流量及年极值洪量
//年均流量及年极值洪量.hvoid 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();}