做了很久沙氏指数si的计算，最后一个问题请教大神

往昔皇后

海蓝蓝 发表于 2016-10-29 23:12
建议检查一下 抬升高度处假相当位温算的对不对。再就是是绝热过程的迭代算的对不对。

都检查过了，确实没有问题，怀疑这个方法不是很适用，问一下，你也是用过这个方法的不？。。。另外对比了探空和fnl的资料，他们的相对湿度有的时候差的较大。

海蓝蓝

往昔皇后 发表于 2016-10-31 08:28
都检查过了，确实没有问题，怀疑这个方法不是很适用，问一下，你也是用过这个方法的不？。。。另外对比了 ...

我是按照文章的方法算的，没问题。你再检查一下抬升凝结高度处的气温

往昔皇后

海蓝蓝 发表于 2016-10-31 19:11
我是按照文章的方法算的，没问题。你再检查一下抬升凝结高度处的气温

又检查了下抬升凝结高度的气温，文献上的抬升凝结高度气温有点小错误，我直接按照刘建文李耀东那本书上的正确公式来的。再就是你说检查抬升凝结高度处的假相当位搵，我并不是计算的这个位搵，而是计算的起始高度850hpa的假相当位搵。方便留下联系方式不，真的想请教一下把这个问题解决了。

海蓝蓝

往昔皇后 发表于 2016-11-1 08:53
又检查了下抬升凝结高度的气温，文献上的抬升凝结高度气温有点小错误，我直接按照刘建文李耀东那本书上的 ...

如果计算SI，应该从850hPa开始上升到达抬升凝结高度。然后计算这个抬升凝结高度处的假相当位温。然后在根据它的守恒做后续计算。
而不是计算850hPa的假相当位温

往昔皇后

海蓝蓝 发表于 2016-11-1 13:41
如果计算SI，应该从850hPa开始上升到达抬升凝结高度。然后计算这个抬升凝结高度处的假相当位温。然后在根 ...

我是计算SI，但是我理解的是假相当位搵是干湿绝热都守恒啊，不论是850还是抬升凝结高度处的假相当位搵都可以吧。
文献中假相当位搵是求850和500的假相当位搵差，我是按照下面这样计算的
!计算850hpa的假相当位温 K     tc是抬升凝结温度                                                       
eqsita1=(t850*((1000/850)**(0.2854*(1-0.00028*r1))))*exp((3.376/tc-0.00254)*r1*(1+0.00081*r1))
!求迭代温度tk1处的假相当位温 K           
eqsita2=(tk1*((1000/500)**(0.2854*(1-0.00028*r))))*exp((3.376/tc-0.00254)*r*(1+0.00081*r))

这个思路不对吗？

往昔皇后

海蓝蓝 发表于 2016-11-1 13:41
如果计算SI，应该从850hPa开始上升到达抬升凝结高度。然后计算这个抬升凝结高度处的假相当位温。然后在根 ...

讲真，我验证了一下，不论是计算抬升凝结高度处还是850hpa处的假相当位搵，最后得出的SI的值几乎没有差别

海蓝蓝

先搞清楚si的定义。必须由850沿着干绝热线上升到lcl

往昔皇后

海蓝蓝 发表于 2016-11-2 00:23
先搞清楚si的定义。必须由850沿着干绝热线上升到lcl

这个我晓得啊，必须由850hpa沿干绝热上升到抬升凝结高度再沿湿绝热上升。
另外假相当位搵在干湿绝热中守恒啊，既然是守恒的话，计算850hpa还是LCL的theta没什么差别。

joeal

'define t=etd' 这行不应该是 'define e = etd'吗

chenxuejing

好好好好哦啊好