蒸发模型
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
//以下计算蒸散发，按照水文教材p156流程图
EU[i] = k * em[i]; // 上层按蒸散发能力蒸发
EL[i] = 0;
ED[i] = 0;
if (ua[i] + PE - R[i] { 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 System.out.println(R[i]+" "+i); } // for(int i=0;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]; }