2022年新安江模型程序C++代码.docx

《2022年新安江模型程序C++代码.docx》由会员分享，可在线阅读，更多相关《2022年新安江模型程序C++代码.docx（31页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选学习资料 - - - - - - - - - 新安江模型程序 C+代码以下是类的声明：class XinanjiangModel private: / FORCING double *m_pP; / 降水数据double *m_pEm; / 水面蒸发数据/ long m_nSteps; / 模型要运行的步长 一共 m_nSteps步 long steps; / OUTPUT double *m_pR; / 流域内每一步长的产流量 径流深度 double *m_pRs; / 每一步长的地表径流深 毫米 double *m_pRi; / 每一步长的壤中流深毫米double *m_pRg; /

2、每一步长的地下径流深 毫米 double *m_pE; / 每一步长的蒸发 毫米 double *m_pQrs; / 流域出口地表径流量double *m_pQri; / 流域出口壤中流径流流量double *m_pQrg; / 流域出口地下径流量double *m_pQ; / 流域出口的总流量double m_U; / SOIL double *m_pW; / 流域内土壤湿度double *m_pWu; / 流域内上层土壤湿度double *m_pWl; / 流域内下层土壤适度double *m_pWd; / 流域内深层土壤湿度double m_Wum; / 流域内上层土壤蓄水容量名师归纳总

3、结 double m_Wlm; / 流域内下层土壤蓄水容量第 1 页,共 17 页- - - - - - -精选学习资料 - - - - - - - - - double m_Wdm; / 流域内深层土壤蓄水容量,WDM=WM-WUM-WLM / EVAPORATION double *m_pEu; / 上层土壤蒸发量毫米double *m_pEl; / 下层土壤蒸发量毫米double *m_pEd; / 深层土壤蒸发量毫米/runoff double *RF; / PARAMETER double m_Kc; / 流域蒸散发才能与实测蒸散发值的比double m_IM; / 不透水面积占全流

4、域面积之比double m_B; / 蓄水容量曲线的方次,小流域几平方公里B0.1左右/ 中等面积平方公里以内,较大面积 .30.4 double m_WM; / 流域平均蓄水容量毫米WM=WUM+WLM+WDM double m_C; / 流域内深层土壤蒸发系数,江南潮湿地区：0,/华北半潮湿地区：double m_SM; /自由水蓄水容量double m_EX; /自由水蓄水容量面积分布曲线指数double m_KG; /地下水日出流系数double m_KI; /壤中流日出流系数double m_CG; /地下水消退系数double m_CI; /壤中流消退系数double *m_UH;

5、 / 单元流域上地面径流的单位线double m_WMM; / 流域内最大蓄水容量double m_Area; / 流域面积int m_DeltaT; / 每一步长的小时数int m_PD; / 给定数据,用以判定是否时行河道汇流运算public: XinanjiangModel void; XinanjiangModel void; 名师归纳总结 - - - - - - -第 2 页,共 17 页精选学习资料 - - - - - - - - - / 初始化模型 void InitModel long nSteps, double Area,int DeltaT, int PD, char *

6、ForcingFile; / 设置模型参数 void SetParametersdouble *Params; / 运行新安江模型 void RunModelvoid; / 储存模拟结果到文件 void SaveResultschar *FileName; / 记录出流数据,用以作图分析 void Runoffchar *runoff; private: / 进行汇流运算,将径流深度转换为流域出口的流量 void Routingvoid; ; 以下是类的定义 #include stdafx.h #include xinanjiangmodel.h #include #include #incl

7、ude using namespace std; #include math.h #include stdio.h #include conio.h XinanjiangModel:XinanjiangModel void this-m_pP = NULL; this-m_pEm = NULL; this-m_pE = NULL; this-m_pEd = NULL; 名师归纳总结 - - - - - - -第 3 页,共 17 页精选学习资料 - - - - - - - - - this-m_pEl = NULL; this-m_pEu = NULL; this-m_pW = NULL; t

8、his-m_pWd = NULL; this-m_pWl = NULL; this-m_pWu = NULL; this-m_pR = NULL; this-m_pRg = NULL; this-m_pRi = NULL; this-m_pRs = NULL; this-m_pQ = NULL; this-m_pQrg = NULL; this-m_pQri = NULL; this-m_pQrs = NULL; XinanjiangModel:XinanjiangModel void delete this-m_pP; delete this-m_pEm; delete this-m_pE;

9、 delete this-m_pEd; delete this-m_pEl; delete this-m_pEu; delete this-m_pW; delete this-m_pWd; delete this-m_pWl; delete this-m_pWu; delete this-m_pR; delete this-m_pRg; delete this-m_pRi; 名师归纳总结 - - - - - - -第 4 页,共 17 页精选学习资料 - - - - - - - - - delete this-m_pRs; delete this-m_pQ; delete this-m_pQr

10、g; delete this-m_pQrs; delete this-m_pQri; / 初始化模型 void XinanjiangModel:InitModel long nSteps, double Area, int DeltaT,int PD, char * ForcingFile FILE * fp; int i; this-m_nSteps = nSteps; this-steps = this-m_nSteps + 18; / 驱动数据 this-m_pP = new doublethis-steps; this-m_pEm = new doublethis-steps; / 模

11、型输出,蒸散发项 this-m_pE = new doublethis-steps; this-m_pEd = new doublethis-steps; this-m_pEl = new doublethis-steps; this-m_pEu = new doublethis-steps; / 模型输出,出流项,经过汇流的产流 this-m_pQrg = new doublethis-steps; this-m_pQrs = new doublethis-steps; this-m_pQri = new doublethis-steps; this-m_pQ = new doublethi

12、s-steps; / 模型输出,产流项 this-m_pR = new doublethis-steps; this-m_pRg= new doublethis-steps; 名师归纳总结 - - - - - - -第 5 页,共 17 页精选学习资料 - - - - - - - - - this-m_pRi= new doublethis-steps; this-m_pRs = new doublethis-steps; / 模型状态量,土壤湿度 this-m_pW = new doublethis-steps; this-m_pWd = new doublethis-steps; this

13、-m_pWl = new doublethis-steps; this-m_pWu = new doublethis-steps; /runoff 值 this-RF = new doublethis-steps; fori = 0;isteps;i+ / 驱动数据 this-m_pP i = 0.00; this-m_pEm i = 0.00; / 模型输出,蒸散发项 this-m_pE i = 0.00; this-m_pEd i = 0.00; this-m_pEl i = 0.00; this-m_pEu i = 0.00; / 模型输出,出流项,经过汇流的产流 this-m_pQrg

14、i = 0.00; this-m_pQrsi = 0.00; this-m_pQrii = 0.00; this-m_pQi = 0.00; / 模型输出,产流项 this-m_pR i = 0.00; this-m_pRg i = 0.00; this-m_pRi i = 0.00; this-m_pRs i = 0.00; 名师归纳总结 - - - - - - -第 6 页,共 17 页精选学习资料 - - - - - - - - - / 模型状态量,土壤湿度 this-m_pW i = 0.00; this-m_pWdi = 0.00; this-m_pWli = 0.00; this-

15、m_pWui = 0.00; this-m_Area = Area; this-m_DeltaT = DeltaT; this-m_PD = PD; this-m_U = this-m_Area/3.6 * this-m_DeltaT; / Forcing文件格式： 第一列：降水单位毫米 空格其次列水面蒸发 毫米iffp = fopenForcingFile, r = NULL printf Can not open forcing file.n;return; fori = 0;im_nSteps;i+ fscanffp,%lf%lf ,& this-m_pPi,& this-m_pEmi;

16、 fclosefp; / 设置模型参数 void XinanjiangModel:SetParametersdouble* Params this-m_Kc = Params0; this-m_IM = Params1; this-m_B = Params2; / 1 流域蒸散发才能与实测水面蒸发之比 / 2 流域不透水面积占全流域面积之比 / 3 蓄水容量曲线的方次this-m_Wum = Params3; / 4 上层蓄水容量this-m_Wlm = Params4; / 5 下层蓄水容量this-m_Wdm = Params5; / 6 深层蓄水容量this-m_C = Params6;

17、 / 7 深层蒸散发系数this-m_SM = Params7; / 8自由水蓄水容量this-m_EX = Params8; / 9自由水蓄水容量面积分布曲线指数名师归纳总结 - - - - - - -第 7 页,共 17 页精选学习资料 - - - - - - - - - this-m_KG = Params9; this-m_KI = Params10; this-m_CG = Params11; / 10地下水日出流系数 / 11壤中流日出流系数 / 12地下水消退系数this-m_CI = Params12; / 13壤中流消退系数 this-m_WM = this-m_Wum +

18、this-m_Wlm + this-m_Wdm; this-m_WMM = this-m_WM * 1.0 + this-m_B/1.0 - this-m_IM; / 运行新安江模型 void XinanjiangModel:RunModel void long i; / 模型的状态变量double PE; / 0 时为净雨量 ; 0 为蒸发不足量 mmdouble Ep; /m_Kc * m_pEmi double P; double R; / 产流深度,包括地表径流、壤中流和地下径流mmdouble RB; / 不透水面上产生的径流深度mmdouble RG; / 地下径流深度 mmdou

19、ble RI; / 壤中流深度 mmdouble RS; / 地表径流深 mmdouble A; /土壤湿度为 W时土壤含水量折算成的径流深度mmdouble E = 0.0; / 蒸散发 mm double EU = 0.0; / 上层土壤蒸散发量 mmdouble EL = 0.0; / 下层土壤蒸散发量 mmdouble ED =0.0; / 深层土壤蒸散发量 mmdouble S; double FRo; double FR; double MS; double AU; 名师归纳总结 - - - - - - -第 8 页,共 17 页精选学习资料 - - - - - - - - - d

20、ouble WU = 5.0; / 流域内上层土壤湿度 double WL = 55.0; / 流域内下层土壤适度 double WD = 40.0; / 流域内深层土壤湿度 double W = 100.0; double So = 5.0; MS = m_SM * 1 + m_EX; FRo = 1 - pow1 - So/MS,m_EX; fori = 0;im_nSteps;i+ / 蒸散发运算 / RB = m_pPi * m_IM; / RB是降在不透水面的降雨量P = m_pPi * 1 - m_IM; Ep = m_Kc * m_pEmi; if WU + P= Ep EU =

21、 Ep; EL = 0; ED = 0; elseif WU + P= m_C * m_Wlm EL = Ep - EU * WL/m_Wlm; ED = 0; elseif m_C * Ep - EU= WL&WLm_C * m_Wlm EL = m_C * Ep - EU; ED = 0; elseif WLm_C * Ep - EU EL = WL; ED = m_C * Ep - EU - EL; E = EU + EL + ED; PE = P - E; /* 蒸散发运算终止 */ / 子流域产流量运算 / if PE= 0 名师归纳总结 - - - - - - -第 9 页,共 1

22、7 页精选学习资料 - - - - - - - - - R = 0.00; W = W + PE; else A = m_WMM * 1 - pow 1.0 - W/m_WM, 1.0/1 + m_B ; / 土壤湿度折算净雨量 +降水后蒸发剩余雨量 流域内最大含水容量 if A + PEm_WMM / 流域内的产流深度运算 R = PE /* 降水蒸发后的剩余量 */ + W /* 流域内土壤湿度 */ + m_WM * pow1 - PE + A/m_WMM,1 + m_B - m_WM /* 减去流域平均蓄水容量 m_WM: 参数*/ + RB; /* 不透水面上产生的径流 */ / 土

23、壤湿度折算净雨量 +降水后蒸发剩余雨量 流域内最大含水容量 else / 流域内的产流深度运算R = PE /* 降水蒸发后的剩余量*/ + W /* 流域内土壤湿度 */ - m_WM /* 减去流域平均蓄水容量+ RB; /* 不透水面上产生的径流 */ /三层蓄水量的运算 : WU, WL, WD if WU + P - EU - R = m_Wum WU = WU + P - EU - R; else WU = m_Wum; WL = WL - EL; WD = WD ED; 名师归纳总结 - - - - - - -第 10 页,共 17 页精选学习资料 - - - - - - - -

24、 - if WL - EL + WU + P - EU - R - m_Wum = m_Wlm WL = WL EL + WU + P - EU - R - m_Wum ; WD = WD - ED; else WL = m_Wlm; if WD - ED + WL - EL + WU + P - EU - R - m_Wum - m_Wlm 0 FR = R - RB / PE; AU = MS * 1 - pow1 - So * FRo/FR/m_SM,1/1 + m_EX; if PE + AU= MS RS = FR * PE + So * Fro / FR - m_SM ; S =

25、So * Fro / FR + R RS / FR; RI = m_KI * S * FR; RG = m_KG * S * FR; 名师归纳总结 - - - - - - -第 11 页,共 17 页精选学习资料 - - - - - - - - - RS += RB; R = RS + RI + RG; So = S * 1 - m_KI - m_KG ; FRo = FR; else S = So; FR = 1 - pow1 S / MS , m_EX ; RI = 0.00; RG = 0.00; So = S * 1 - m_KI - m_KG ; RS = RB; R = RS +

26、RI + RG; FRo = FR; /三水源划分运算终止/* 以下部分是状态量：总蒸发量、上、下和深层土壤的蒸发的储存 */ /* 1 */ this-m_pEi = E; / 当前步长的蒸发模型重要输出/* 2 */ this-m_pEui = EU; / 当前步长上层土壤蒸发/* 3 */ this-m_pEli = EL; / 当前步长下层土壤蒸发/* 4 */ this-m_pEdi = ED; / 当前步长深层土壤蒸发/* 5 */ this-m_pWi = W; / 当前步长流域平均土壤含水量/* 6 */ this-m_pWui = WU; / 当前步长流域上层土壤含水量/*

27、7 */ this-m_pWli = WL; / 当前步长流域下层土壤含水量/* 8 */ this-m_pWdi = WD; / 当前步长流域深层土壤含水量/* 9 */ this-m_pRgi = RG; / 当前步长流域地下径流深度/* 10 */ this-m_pRii = RI; /* 11 */ this-m_pRsi = RS; / 当前步长流域壤中流深度 / 当前步长流域地表径流径流深度名师归纳总结 - - - - - - -第 12 页,共 17 页精选学习资料 - - - - - - - - - /* 12 */ this-m_pRi = R; / 当前步长的总产流径流深度

28、this-Routing; / 储存模拟结果到文件 void XinanjiangModel:SaveResultschar* FileName int i; FILE * fp; iffp = fopenFileName,w = NULL printfCan not create output file.n; return; fprintffp, - - - - - - - - - - - - - - - - - - - - - - - - - - n; fprintffp, Emm EUmm ELmm EDmm Wmm WUmm WLmm WDmm Rmm RSmm RImm RGmm Q

29、m3/d QSm3/d QIm3/d QGm3/dn; fprintffp, - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -n ; fori = 0;isteps;i+ fprintffp, %9.3lf%9.3lf%9.3lf%9.3lf%9.3lf%9.3lf%9.3lf%9.3lf%9.3lf%9.3lf% 9.3lf%9.3lf%9.3lf%9.3lf%9.3lf%9.3lfn , this-m_pEi, this-m_pEui,this-m_pEli, this-m_pEdi, this-m_pWi, this-

30、m_pWui, this-m_pWli, this-m_pWdi, this-m_pRi, this-m_pRsi,this-m_pRii, this-m_pRgi, this-m_pQi, this-m_pQrsi, this-m_pQrii, this-m_pQrgi; 名师归纳总结 - - - - - - -第 13 页,共 17 页精选学习资料 - - - - - - - - - fclosefp; / 进行汇流运算,将径流深度转换为流域出口的流量 void XinanjiangModel:Routing void / 地面径流汇流运算：单位线法/ int i,j; double B1

31、0000 = 0.00; if this-m_PD = 1 double UH =3.71,12.99,38.96,94.63,131.74,154.00,166.99,176.27,178.12, 4,18.55, 9.27, 3.71,0.00; fori = 0;im_nSteps;i+ forj = 0;jm_pRsi * UHj/10.0; else double UH = 7.18,23.38,63.20,143.10,221.75,365.18,447.40,491.29, 506.93,504.82,468.46,388.56,309.91,166.49,84.26,40.37

32、,17.56,3.46; fori = 0;im_nSteps;i+ forj = 0;jm_pRsi * UHj/10.0; fori = 0;isteps;i+ this-m_pQrsi = Bi; 名师归纳总结 - - - - - - -第 14 页,共 17 页精选学习资料 - - - - - - - - - / 壤中流汇流运算 :线性水库1 fori = 1;isteps;i+ this-m_pQrii = this-m_CI * this-m_pQrii - + 1.0 - this-m_CI * this-m_pRii * this-m_U; / 地下径流汇流运算：线性水库for

33、i = 1;isteps;i+ this-m_pQrgi = this-m_pQrgi - 1 * this-m_CG + this-m_pRgi * 1.0 - this-m_CG * this-m_U; /单元面积总入流运算 fori = 0;isteps;i+ this-m_pQi = this-m_pQrsi + this-m_pQrii + this-m_pQrgi; void XinanjiangModel:Runoff char * runoff int i; ofstream outfile; outfile.openrunoff; ifoutfile.is_open fori

34、 = 0;isteps;i+ outfilesetprecision3setiosflagsios:fixedm_pQiendl; outfile.close; 以下是 main函数语句 int _tmainint argc, _TCHAR * argv long nSteps = 942; int DeltaT = 24; double Area1 = 1603; XinanjiangModel Model1; 名师归纳总结 - - - - - - -第 15 页,共 17 页精选学习资料 - - - - - - - - - Model1.InitModelnSteps, Area1,Del

35、taT,1,LFForcingfile.txt ; /模型参数 /*Kc,IM,B,m_Wum,Wlm,Wdm,C,SM,EX,KG,KI,CG,CI double Params1 = 0.50,0.01,0.30,10,60,40,0.18,32,1.2,0.075,0.072,0.94,0.7; Model1.SetParametersParams1; Model1.RunModel; Model1.SaveResults来凤站日模型运算结果 .txt ; Model1.Runoff LF_Q.txt ; Model1.XinanjiangModel; doubleArea2 = 2991

36、; XinanjiangModel Model2; Model2.InitModelnSteps, Area2,DeltaT,0,YCForcingfile.txt ; /模型参数 /*Kc,IM,B,m_Wum,Wlm,Wdm,C,SM,EX,KG,KI,CG,CI double Params2 = 0.75,0.01,0.32,10,60,40,0.18,27,1.2,0.065,0.067,0.96,0.8; Model2.SetParametersParams2; Model2.RunModel; Model2.SaveResultsfile.txt ; Model2.Runoff Y

37、C_Q.txt ; Model2.XinanjiangModel; FILE *fp1,*fp2; double Q11000,Q21000,Q1000=0.00; ofstream outfile; outfile.openQ.txt ; iffp1=fopenLF_Q.txt ,r =NULL printf Can not open the file.n; return 0; iffp2=fopenYC_Q.txt ,r=NULL printf Can not open the file.n; return 0; ifoutfile.is_open forint i=0;i960;i+ 名师归纳总结 - - - - - - -第 16 页,共 17 页精选学习资料 - - - - - - - - - fscanffp1,%lf ,&Q1i; fscanffp2,%lf ,&Q2i; Qi=Q1i+Q2i; outfilesetprecision3setio