计算沙氏指数跳不出循环，不知哪里出了问题，求帮助

junzhi23

%绝对温度转换为摄氏温度设为气块初始温度
T3=TL(3)-273.15;
%用迭代法来得到850hpa气块先沿干绝热线上升，到达抬升凝结高度后再沿湿绝热线上升至500hpa时所具有的气块温度（摄氏度）
n=0;
while 1
    n=n+1;
    T3=T3-0.0000001*n;
     Tk=T3+273.15;
    est=6.112*exp((17.67*T3)/(T3+243.5));
    etd=6.112*exp((17.67*Td(6)/(Td(6)+243.5)));
    e=etd;
    U1=(etd/est)*100;
    TL1=1/(1/(Tk-55)-log(U1/100)/2840)+55;
    r1=622*(e/(level(6)-e));
    sita2=(Tk.*((1000/level(6))^(0.2854*(1-0.00028*r1)))).*exp((3.376/TL1-0.00254)*r1.*(1+0.00081*r1));%计算500hpa相当位温
    sita3=sita2-sita(3);
    if abs(sita3)<0.002%根据相对位温守恒，当两处的相当位温差小于0.001℃时迭代结束
        sita2=sita2;
        TT=T3; est为饱和水汽压，etd为实际水气压
        break
    end
     
end
si=T(2)-TT; %计算沙瓦特指数

说明：TL（3）为850hPa抬升凝结高度，等于216.74。

chenxuejing

调出来了吗

junzhi23

chenxuejing 发表于 2017-7-10 23:12
调出来了吗

没，好久没做了。。。最近又才开始做。。

chenxuejing

junzhi23 发表于 2017-7-25 13:48
没，好久没做了。。。最近又才开始做。。

我也调不出来

南南偶卖糕Oo

我的也跳不出循环，请问楼主解决了吗