新安江模型程序c代碼-資料下載頁(yè)
2025-07-07 14:43本頁(yè)面
【正文】 thism_pE[i],thism_pEu[i],thism_pEl[i],thism_pEd[i], thism_pW[i],thism_pWu[i],thism_pWl[i],thism_pWd[i], thism_pR[i],thism_pRs[i],thism_pRi[i],thism_pRg[i], thism_pQ[i],thism_pQrs[i],thism_pQri[i],thism_pQrg[i])。 } fclose(fp)。}// 進(jìn)行匯流計(jì)算，將徑流深度轉(zhuǎn)換為流域出口的流量void XinanjiangModel::Routing(void){ ///////////// 地面徑流匯流計(jì)算：?jiǎn)挝痪€法 /////////////////////// int i,j。 double B[10000] = {}。 if (thism_PD == 1) { double UH[] ={,,, , , , ,}。 for(i = 0。ithism_nSteps。i++ ) { for(j = 0。j18。j++ ) { B[i + j] += thism_pRs[i] * UH[j]/。 } } } else { double UH[] ={ ,,,,,}。 for(i = 0。ithism_nSteps。i++ ) { for(j = 0。j18。j++ ) { B[i + j] += thism_pRs[i] * UH[j]/。 } } } for(i = 0。ithissteps。i++ ) thism_pQrs[i] = B[i]。 ///// 壤中流匯流計(jì)算:線性水庫(kù) for(i = 1。ithissteps。i++ ) {thism_pQri[i] = thism_CI * thism_pQri[i 1] + ( thism_CI) * thism_pRi[i] * thism_U。 } ///// 地下徑流匯流計(jì)算：線性水庫(kù) for(i = 1。ithissteps。i++ ) {thism_pQrg[i] = thism_pQrg[i 1] * thism_CG + thism_pRg[i] * ( thism_CG) * thism_U。 } //////單元面積總?cè)肓饔?jì)算 for(i = 0。ithissteps。i++ ) { thism_pQ[i] = thism_pQrs[i] + thism_pQri[i] + thism_pQrg[i]。 }}void XinanjiangModel::Runoff(char * runoff){ int i。 ofstream outfile。 (runoff)。 if(()){ for(i = 0。ithissteps。i++ ) { outfilesetprecision(3)setiosflags(ios::fixed)thism_pQ[i]endl。 }} ()。}以下是main()函數(shù)語(yǔ)句int _tmain(int argc, _TCHAR * argv[]){ long nSteps = 942。 int DeltaT = 24。 double Area1 = 1603。 XinanjiangModel Model1。 (nSteps, Area1,DeltaT,1,)。//模型參數(shù)/*Kc,IM,B,m_Wum,Wlm,Wdm,C,SM,EX,KG,KI,CG,CI double Params1[] = {,10,60,40,32,,}。 (Params1)。 ()。 ()。 ()。 Model1.~XinanjiangModel()。 double Area2 = 2991。 XinanjiangModel Model2。 (nSteps, Area2,DeltaT,0,)。//模型參數(shù)/*Kc,IM,B,m_Wum,Wlm,Wdm,C,SM,EX,KG,KI,CG,CI double Params2[] = {,10,60,40,27,,}。 (Params2)。 ()。 ()。 ()。 Model2.~XinanjiangModel()。 FILE *fp1,*fp2。 double Q1[1000],Q2[1000],Q[1000]={}。 ofstream outfile。 ()。 if((fp1=fopen(,r))==NULL) { printf(Can not open the file!\n)。 return 0。} if((fp2=fopen(,r))==NULL) { printf(Can not open the file!\n)。 return 0。 } if(()) { for(int i=0。i960。i++) { fscanf(fp1,%lf,amp。Q1[i])。 fscanf(fp2,%lf,amp。Q2[i])。 Q[i]=Q1[i]+Q2[i]。 outfilesetprecision(3)setiosflags(ios::fixed)Q[i]endl。 } fclose(fp1)。 fclose(fp2)。 } () 。 return 0。}
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1