蒸发模型

//以下计算蒸散发，按照水文教材p156流程图
EU[i] = k * em[i]; // 上层按蒸散发能力蒸发
EL[i] = 0;
ED[i] = 0;
if (ua[i] + PE - R[i] <UM)// 如果上层土壤含水量没有饱和
{
ua[i + 1] = ua[i] + PE - R[i];
la[i + 1] = la[i];
da[i + 1] = da[i];
} else {//如果上层土壤含水量已经饱和
if(ua[i] + la[i] + PE- R[i] - UM>LM)// 如果第二层土壤含水量已经饱和
{
ua[i + 1] = UM;
la[i + 1] = LM;
da[i + 1] = wmt + PE - R[i] - ua[i + 1] - la[i + 1];
} else// 如果第二层土壤含水量没有饱和
{
ua[i + 1] = UM;
la[i+ 1] = ua[i] + la[i] + PE - R[i] - UM;
da[i + 1] = da[i];
}
}
}
// 以下计算产流量PE<0
elseif (PE< 0)
{//？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？
R[i] = 0;
if (ua[i] + PE> 0)// 上层没有被蒸发完
{
EU[i] = EP;// 上层按蒸散发能力蒸发
EL[i] = 0;
ED[i] = 0;
ua[i + 1] = ua[i] + PE;
la[i + 1] = la[i];
da[i + 1] = da[i];
} else// 上层已经被蒸发完
{
EU[i] = ua[i] + P[i];
ua[i + 1] = 0;
if (la[i] >c * LM) {
EL[i] = (k * em[i] - EU[i]) * la[i] / LM;
la[i + 1] = la[i] - EL[i];
da[i + 1] = da[i];
ED[i] = 0;
} else {
if (la[i] >c * (k * em[i] - EU[i])) {
EL[i] = c * (k * em[i] - EU[i]);
ED[i] = 0;
la[i + 1] = la[i] - EL[i];
da[i + 1] = da[i];
} else {
EL[i] = la[i];
la[i + 1] = 0;
ED[i] = c * (k * em[i] - EU[i]) - EL[i];
da[i + 1] = da[i] - ED[i];
}
}
}
}
if (ua[i + 1] < 0)
ua[i + 1] = 0;
if (la[i + 1] < 0)
la[i + 1] = 0;
if (da[i + 1] < 0)
da[i + 1] = 0;
if (ua[i + 1] >UM)
ua[i + 1] = UM;
if (la[i + 1] >LM)
la[i + 1] = LM;
if (da[i + 1] >DM)
da[i + 1] = DM;
}
}
// writeFilewF = new writeFile(R,i_periods);
for (int i = 0; i<i_periods; i++){
System.out.println(R[i]+" "+i);
}
// for(int i=0;i<i_periods;i++){
// System.out.println(EU[i]+" "+EL[i]+" "+ED[i]+" "+i+" "+"陈");
// }
return R;
}
if (PE<=0)
{
ED[i] = 0;
la[i + 1] = la[i] - EL[i];
da[i + 1] = da[i]; }。