function?Height=Delte（TakTskRHu）
??z1=6;??????%測量參考高度??單位：m
??z0=0.00015;?????%空氣動力學中表面粗糙度參數??單位：m
??e=inline（‘6.105*exp（25.22*（T-273.2）/T-5.31*log（T/273.2））‘‘T‘）;
??Height=[];
??
??Rib=369*z1*（Tak-Tsk）/（u^2*Tak）;??%總體理查森數
??Gamma=PieceGamma（Rib）;?????
??Lmo=（10*z1*Gamma）/Rib;?????%Monin-Obukhov長度?
??ea=RH/100*e（Tak）;?????????????%海洋表面周圍空氣水汽壓
??es=e（Tsk）;????????????????????%海洋表面海水水汽壓
??
??Na=（77.6/Tak）*（1000+4810*ea/Tak）;????%海洋表面空氣折射率
??Ns=?（77.6/Tsk）*（1000+4810*es/Tsk）;????%海洋表面海水折射率
??Np=（Na-Ns）;

??if?0%???????Height=Np/（-0.125*（log（z1/z0）+5.2*（z1-z0）/Lmo）-5.2*Np/Lmo）;
%???????if?Height>0||Height/Lmo>1
%??????????Height=（6.2*Np+0.65*z1）/（-0.125*log（z1/z0））;
%???????end
??????syms?z;
??????Height=solve（Np*（1+5.2*z/Lmo）/（z*（log（z1/z0）+5.2*（z1-z0）/Lmo））+0.125）;
??????
?elseif?Rib<0
??????B=log（z1/z0）-PiecePsi（z1Lmo）;
??????A=-0.125*B/Np;
??????Height=（A^4-18*A^3/Lmo）^（-1/4）;
??????
??elseif?Rib==0
??????Height?=?Np/（-0.125*log（z1/z0））;
??end
??if?Height<0
??????Height=[];
??end
??????
??%Height=Height;
?
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4389??2014-06-18?19:16??Main_H.m

?????文件????????237??2014-04-27?13:37??PieceGamma.m

?????文件????????390??2014-05-19?22:09??PiecePsi.m

?????文件???????1152??2014-06-18?18:41??Delte.m

-----------?---------??----------?-----??----

?????????????????6168????????????????????4