比湿计算结果偏大是什么原因

zhangyuxin0430

求助！！！本人现有2m高度比湿数据，但是想要一个2m地表的比湿梯度，因此想根据地表温度计算出地表比湿，但是最后的结果偏大，是什么原因呢?附上程序及2m比湿结果及自己求出的地表结果（带有黄色），求帮忙。。。

ca=6.112;
cb=273.16;
c70=21.874;
c71=7.66;
c72=17.269;
c73=35.86;

if TS1<=273.15     %0°以下冰面     
   eog=ca*exp(21.874*(TS1-cb)./(TS1-7.66));           %surface (sat) vapor pressure
   else
   eog=ca*exp(17.269*(TS1-cb)./(TS1-35.86));           %surface (sat) vapor pressure
   end

qg=0.622*eog./(PS1/100-0.378*eog);

delq=Q1-qg;