2024年5月7日发(作者:)

%%%新安江模型程序核心源代码

function Qr=XAJ_JUN(DAREA,DT,EM,WwU,WwL,WwD,P,S0, FR0, Qrs0, Qrss0,

Qrg0,parameter,Qm)

% XAJ是新安江的运行程序,用于单纯形和遗传算法调用,也用于新安江模型的预报

Imp1= ; %流域不透水面积比: 次洪

Kc= ; %流域蒸散发折算系数:多年总径流量决定

WMU= ; %流域上层蓄水容量

WML= ; %流域中层蓄水容量

WMD = ; %流域下层蓄水容量

B = parameter.B ; %流域蓄水容量分布曲线指数

C = parameter.C ; %流域深层蒸发系数

Ex = ; %流域自由水分布曲线指数

SM = ; %流域自由水平均蓄水容量

Ki = ; %自由水箱壤中流出流系数

Kg = ; %自由水箱地下水出流系数

Cs = ; %地面水线性水库汇流系数

Ci = ; %壤中流线性水库汇流系数

Cg = ; %地下水线性水库汇流系数

Ke = ; %马斯京根法河段传播时间

Xe = ; %马斯京根法流量比重系数

L = parameter.L ; %滞后演算法参数

%次洪决定:WM,B,Imp

%WwU(0)=WwU;WwL(0)=WwL;WwD(0)=WwD;

%由于日模型与次洪模型的计算时段长不同,参数值不能全部通用,但K、WM、WUM、

WLM、B、IMP、EX、C与时段长无关,可以直接引用,

%Kc SM、KG、KSS、CS、CI、CG与时段长相关,不能直接引用,需要另外率定

%junjunzhu-XAJ-MODEL

U=DAREA/(DT*3.6); %单位转换