爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: sunflower

[分享资料] 做了很久沙氏指数si的计算,最后一个问题请教大神

[复制链接]

新浪微博达人勋

发表于 2016-10-31 08:28:39 | 显示全部楼层
海蓝蓝 发表于 2016-10-29 23:12
建议检查一下 抬升高度处假相当位温算的对不对。再就是是绝热过程的迭代算的对不对。

都检查过了,确实没有问题,怀疑这个方法不是很适用,问一下,你也是用过这个方法的不?。。。另外对比了探空和fnl的资料,他们的相对湿度有的时候差的较大。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-31 19:11:05 | 显示全部楼层
往昔皇后 发表于 2016-10-31 08:28
都检查过了,确实没有问题,怀疑这个方法不是很适用,问一下,你也是用过这个方法的不?。。。另外对比了 ...

我是按照文章的方法算的,没问题。你再检查一下抬升凝结高度处的气温
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-11-1 08:53:20 | 显示全部楼层
海蓝蓝 发表于 2016-10-31 19:11
我是按照文章的方法算的,没问题。你再检查一下抬升凝结高度处的气温

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

新浪微博达人勋

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

如果计算SI,应该从850hPa开始上升到达抬升凝结高度。然后计算这个抬升凝结高度处的假相当位温。然后在根据它的守恒做后续计算。
而不是计算850hPa的假相当位温
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-11-1 14:51:43 | 显示全部楼层
海蓝蓝 发表于 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))

这个思路不对吗?


密码修改失败请联系微信:mofangbao

新浪微博达人勋

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

讲真,我验证了一下,不论是计算抬升凝结高度处还是850hpa处的假相当位搵,最后得出的SI的值几乎没有差别
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-11-2 00:23:56 | 显示全部楼层
先搞清楚si的定义。必须由850沿着干绝热线上升到lcl
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-11-2 08:45:23 | 显示全部楼层
海蓝蓝 发表于 2016-11-2 00:23
先搞清楚si的定义。必须由850沿着干绝热线上升到lcl

这个我晓得啊,必须由850hpa沿干绝热上升到抬升凝结高度再沿湿绝热上升。
另外假相当位搵在干湿绝热中守恒啊,既然是守恒的话,计算850hpa还是LCL的theta没什么差别。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-21 15:15:01 | 显示全部楼层
'define t=etd' 这行不应该是 'define e = etd'吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-7-10 11:30:39 | 显示全部楼层
好好好好哦啊好
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表