需求

计算饱和水汽压差VPD

代码

%处理计算vpd数据，参考世界气象组织仪器和观测方法委员会指南转换方程
%By Yang，2019/5/13
load rh02_05d.mat          %加载湿度数据
load tem02_05d.mat         %加载温度数据
%前四行为站点号、年月日、第五列为数据
%% 计算饱和水汽压差
for i = 1:length(tem02_05d)
     temp = tem02_05d{1,i};
     temp2 = rh02_05d{1,i};
     Bh = 6112*0.1*exp(22.46*temp(:,5)./(272.62+temp(:,5)));   %Bh饱和蒸气压,单位帕pa
     Sj = Bh.*temp2(:,5);                                      %Sj实际蒸气压     
     vpd = Bh - Sj;                                            %vpd饱和水汽压差
     vpd02_05d{1,i}(:,1:4) = temp(:,1:4);
     vpd02_05d{1,i}(:,5) = vpd;                                %一个站点一个元胞，排为一行
end

VPD计算参考

世界气象组织仪器和观测方法委员会指南转换方程