爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: sunflower

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

[复制链接]

新浪微博达人勋

发表于 2016-4-15 16:08:25 | 显示全部楼层
在程序循环迭代部分,求饱和水汽压用的是迭代温度,求水汽压用的是500hpa温度,我在迭代过程中发现在还未满足迭代条件时,相对湿度会出现大于100的情况,楼主有没有考虑过?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-16 01:03:38 | 显示全部楼层

回帖奖励 +1 金钱

计算si还是有点复杂的。为什么不用fortran计算呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-19 10:54:20 | 显示全部楼层
往昔皇后 发表于 2016-4-15 16:08
在程序循环迭代部分,求饱和水汽压用的是迭代温度,求水汽压用的是500hpa温度,我在迭代过程中发现在还未满 ...

真的吗?那得好好研究一下,亲是否找出了问题所在
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-19 10:54:24 | 显示全部楼层
往昔皇后 发表于 2016-4-15 16:08
在程序循环迭代部分,求饱和水汽压用的是迭代温度,求水汽压用的是500hpa温度,我在迭代过程中发现在还未满 ...

真的吗?那得好好研究一下,亲是否找出了问题所在
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-19 10:55:07 | 显示全部楼层
海蓝蓝 发表于 2016-4-16 01:03
计算si还是有点复杂的。为什么不用fortran计算呢

当时fortran不熟,grads相对用的多谢,现在在研究fortran哒
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-19 10:55:15 | 显示全部楼层
海蓝蓝 发表于 2016-4-16 01:03
计算si还是有点复杂的。为什么不用fortran计算呢

当时fortran不熟,grads相对用的多谢,现在在研究fortran哒
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-21 09:16:50 | 显示全部楼层
sunflower 发表于 2016-4-19 10:55
当时fortran不熟,grads相对用的多谢,现在在研究fortran哒

Si的计算,关键在于计算气块温度。这里面很可能要用到假相当位温的计算。建议直接用Bolton公式。不要用定义式。假相当位温计算不准的话,会影响沙氏指数计算精度。具体见这篇文章:
http://bbs.06climate.com/forum.php?mod=viewthread&tid=42944
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-25 17:02:19 | 显示全部楼层

回帖奖励 +1 金钱

sunflower 发表于 2016-2-4 16:21
在家园各位亲的帮助下,终于算出来啦。目测家园里面还木有多少计算沙氏指数的gs哦,下面我把改正的问题说明 ...

美女好啊!
参考了你的SI程序,但是总卡在循环那出问题,菜鸟一个,想问下为什么sublin和subwrd分别选第二行和第一个数值呢,谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-31 21:44:48 | 显示全部楼层
海蓝蓝 发表于 2016-4-21 09:16
Si的计算,关键在于计算气块温度。这里面很可能要用到假相当位温的计算。建议直接用Bolton公式。不要用定 ...

非常感谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-6-29 10:39:16 | 显示全部楼层

回帖奖励 +1 金钱

楼主能共享下完整的程序么?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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