新安江模型程序c代碼(參考版)
2025-07-10 14:43本頁面
【正文】 }。 } () 。 } fclose(fp1)。 Q[i]=Q1[i]+Q2[i]。 fscanf(fp2,%lf,amp。i++) { fscanf(fp1,%lf,amp。 } if(()) { for(int i=0。} if((fp2=fopen(,r))==NULL) { printf(Can not open the file!\n)。 if((fp1=fopen(,r))==NULL) { printf(Can not open the file!\n)。 ofstream outfile。 FILE *fp1,*fp2。 ()。 ()。//模型參數(shù)/*Kc,IM,B,m_Wum,Wlm,Wdm,C,SM,EX,KG,KI,CG,CI double Params2[] = {,10,60,40,27,}。 XinanjiangModel Model2。 Model1.~XinanjiangModel()。 ()。 (Params1)。 (nSteps, Area1,DeltaT,1,)。 double Area1 = 1603。}以下是main()函數(shù)語句int _tmain(int argc, _TCHAR * argv[]){ long nSteps = 942。i++ ) { outfilesetprecision(3)setiosflags(ios::fixed)thism_pQ[i]endl。 if(()){ for(i = 0。 ofstream outfile。i++ ) { thism_pQ[i] = thism_pQrs[i] + thism_pQri[i] + thism_pQrg[i]。 } //////單元面積總?cè)肓饔?jì)算 for(i = 0。ithissteps。i++ ) {thism_pQri[i] = thism_CI * thism_pQri[i 1] + ( thism_CI) * thism_pRi[i] * thism_U。 ///// 壤中流匯流計(jì)算:線性水庫 for(i = 1。ithissteps。j++ ) { B[i + j] += thism_pRs[i] * UH[j]/。i++ ) { for(j = 0。 for(i = 0。j++ ) { B[i + j] += thism_pRs[i] * UH[j]/。i++ ) { for(j = 0。 for(i = 0。 double B[10000] = {}。 } fclose(fp)。ithissteps。 fprintf(fp, \n)。 } fprintf(fp, \n)。 if((fp = fopen(FileName,w)) == NULL) { printf(Can not create output file!\n)。}// 保存模擬結(jié)果到文件void XinanjiangModel::SaveResults(char* FileName){ int i。 // 當(dāng)前步長(zhǎng)流域地表徑流徑流深度 /* 12 */ thism_pR[i] = R。 // 當(dāng)前步長(zhǎng)流域地下徑流深度 /* 10 */ thism_pRi[i] = RI。 // 當(dāng)前步長(zhǎng)流域下層土壤含水量 /* 8 */ thism_pWd[i] = WD。 // 當(dāng)前步長(zhǎng)流域平均土壤含水量 /* 6 */ thism_pWu[i] = WU。 // 當(dāng)前步長(zhǎng)下層土壤蒸發(fā) /* 4 */ thism_pEd[i] = ED。 // 當(dāng)前步長(zhǎng)的蒸發(fā)（模型重要輸出） /* 2 */ thism_pEu[i] = EU。 FRo = FR。 RS = RB。 RG = 。 FR = 1 pow((1 – S / MS) , m_EX )。 FRo = FR。 R = RS + RI + RG。 RG = m_KG * S * FR。 S = So * Fro / FR + ( R – RS ) / FR。 if(PE + AUMS) RS = FR * (PE + So * FRo/FR m_SM + m_SM * pow((1 (PE + AU) / MS),m_EX + 1))。 ////三水源劃分匯流計(jì)算 if(PE0) { FR = (R RB) / PE。 else WD = m_Wdm。 } else { WL = m_Wlm。 if(WL EL + ( WU + P EU R m_Wum ) = m_Wlm ) { WL = WL – EL + ( WU + P EU R m_Wum )。 WD = WD – ED。 /* 不透水面上產(chǎn)生的徑流*/ } } //三層蓄水量的計(jì)
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1