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

《新安江模型程序C++代码.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、新安江模型程序C++代码以下是类的声明：classXinanjiangModel{private://FORCINGdouble*m_pP;//降水数据double*m_pEm;//水面蒸发数据//longm_nSteps;//模型要运行的步长(一共m_nSteps步)longsteps;//OUTPUTdouble*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;//流域出口的总流量doublem_U;//for24h.U=A(km^2)/3.6/delta_t//SOILdouble*m_pW;//流域内土壤湿度double*m_pWu;//流域内上层土壤湿度double*m_pWl;//流域内下层土壤适度doubl

3、e*m_pWd;//流域内深层土壤湿度doublem_Wum;//流域内上层土壤蓄水容量doublem_Wlm;//流域内下层土壤蓄水容量doublem_Wdm;//流域内深层土壤蓄水容量，WDM=WM-WUM-WLM//EVAPORATIONdouble*m_pEu;//上层土壤蒸发量（毫米）double*m_pEl;//下层土壤蒸发量（毫米）double*m_pEd;//深层土壤蒸发量（毫米）//runoffdouble*RF;//PARAMETERdoublem_Kc;//流域蒸散发能力与实测蒸散发值的

4、比doublem_IM;//不透水面积占全流域面积之比doublem_B;//蓄水容量曲线的方次，小流域（几平方公里）B0.1左右//中等面积（平方公里以内）.2~0.3，较大面积.3~0.4doublem_WM;//流域平均蓄水容量（毫米）(WM=WUM+WLM+WDM)doublem_C;//流域内深层土壤蒸发系数，江南湿润地区：0.15-0.2，//华北半湿润地区：.09-0.12doublem_SM;//自由水蓄水容量doublem_EX;//自由水蓄水容量～面积分布曲线指数doublem_KG;//

5、地下水日出流系数doublem_KI;//壤中流日出流系数doublem_CG;//地下水消退系数doublem_CI;//壤中流消退系数double*m_UH;//单元流域上地面径流的单位线doublem_WMM;//流域内最大蓄水容量doublem_Area;//流域面积intm_DeltaT;//每一步长的小时数intm_PD;//给定数据，用以判断是否时行河道汇流计算public:XinanjiangModel(void);~XinanjiangModel(void);//初始化模型voidInitM

6、odel(longnSteps,doubleArea,intDeltaT,intPD,char*ForcingFile);//设置模型参数voidSetParameters(double*Params);//运行新安江模型voidRunModel(void);//保存模拟结果到文件voidSaveResults(char*FileName);//记录出流数据，用以作图分析voidRunoff(char*runoff);private://进行汇流计算，将径流深度转换为流域出口的流量voidRouting(vo

7、id);};以下是类的定义#include"stdafx.h"#include"xinanjiangmodel.h"#include#include#includeusingnamespacestd;#include"math.h"#include"stdio.h"#include"conio.h"XinanjiangModel::XinanjiangModel(void){this->m_pP=NULL;this->m_pEm=NULL;this->m

8、_pE=NULL;this->m_pEd=NULL;this->m_pEl=NULL;this->m_pEu=NULL;this->m_pW=NULL;this->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