MATLAB经典蒸发模型代码

经典蒸发模型是指Penman-Monteith方法，其MATLAB代码如下：

% 输入参数
T = 25; % 空气温度，单位：摄氏度
RH = 60; % 相对湿度，单位：百分比
U = 2; % 风速，单位：m/s
P = 101325; % 大气压强，单位：Pa
Rn = 200; % 净辐射，单位：W/m2
G = 0; % 地面热通量，单位：W/m2
LAI = 2; % 叶面积指数
z = 2; % 高度，单位：m

% 常数
sigma = 5.67e-8; % 斯特藩-玻尔兹曼常数，单位：W/m2/K4
cp = 1005; % 干空气比热，单位：J/kg/K
lambda = 2450000; % 水的潜热，单位：J/kg

% 饱和水汽压力，单位：Pa
es = exp(77.345 + 0.0057 * T - 7235 / T) / T^8.2;
% 实际水汽压力，单位：Pa
ea = RH / 100 * es;
% 空气密度，单位：kg/m3
rhoa = P / (287.05 * (T + 273.15));
% 热舒适度参数，单位：Pa/K
delta = 4098 * es / (T + 273.15)^2;
% 水汽压斜率曲线，单位：Pa/K
gamma = 0.665e-3 * P / lambda;
% 叶面温度，单位：摄氏度
Ts = T + 0.5 * (LAI * sigma * (T + 273.15)^4 + Rn - G) / (rhoa * cp * U) - 273.15;
% 潜在蒸散量，单位：mm/s
LE = delta * (Rn - G) / (delta + gamma * (1 + 0.34 * U)) / lambda * 1000;
% 实际蒸散量，单位：mm/s
E = 1.26 * (ea - 0.7 * es) / (delta + gamma * (1 + 0.34 * U)) / lambda * 1000;

% 输出结果
disp(['Ts = ' num2str(Ts) ' °C']);
disp(['LE = ' num2str(LE) ' mm/s']);
disp(['E = ' num2str(E) ' mm/s']);