按土壤垂向分布的不均匀性将土层分为三层,用三层蒸散发模型计算蒸散发量。参数有流域平均张力水容量WM(mm),上层张力水容量WUM(mm),下层张力水容量WLM(mm),深层张力水容量WDM(mm),深层蒸散发扩散系数C。
具体计算为
若P+WU>=EP,则EU=EP,EL=0,ED=0; 若P+WU
若WL>C*LM,则WL=(EP-EU)WL/LM,ED=0;
若WL=C*(EP-EU),则EL=C*(EP-EU),ED=0; 若WL
>> WUM=20;WLM=60;WDM=40;C=1/6;
P=[0,0,0.8,0,0,0.2,0,0,0,0,0,0.3,0,0,0,7.4,0.1,0];
Ep=[7.9,7.4,5.9,6.1,6.2,5.8,5.0,5.2,5.4,6.9,6.7,4.1,5.8,4.0,4.3,5.9,4.2,6.3];
WU=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1.5,0];
n=length(Ep);
WL(1)=14.9;EL(1)=2;
for(i=2:n)
WL(i)=WL(i-1)-EL(i-1)
if(WU[i]+P[i])>=Ep[i];
EU[i]=Ep[i];EL[i]=o;ED[i]=o;
elseif (WL[i]>=C*WLM);
EU[i]=WU[i]+P[i];EL[i]=(Ep[i]-(WU[i]+P[i]))*(WL[i]/WLM);ED[i]=o;
elseif(C*(Ep[i]-(WU[i]+P[i]))<=WL[i])
EU[i]=WU[i]+P[i];EL[i]=(Ep[i]-(WU[i]+P[i]))*C;ED[i]=o;
else
EU[i]=WU[i]+P[i];EL[i]=WL[i];ED[i]=(Ep[i]-(WU[i]+P[i]))*C-WL[i];
end
end
end
end
??? if(WU[i]+P[i])>=Ep[i];
|
Error: Unbalanced or misused parentheses or brackets.
============================================================
实在不知道怎么写。。。初学者求教大神