计算沙氏指数的程序

astiny · 发表于 2019-6-3 10:53:42

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 astiny 于 2021-1-26 15:40 编辑

基本上是照抄的

李社宏.用C语言开发的气象常用参数和物理量计算函数库（一）[J].陕西气象,1994(03):42-45.

复制代码

然鹅几十年前的期刊了。目测代码完全是靠人工排版的（辛苦编辑了），所以一眼望上去，满眼都是bug。然后凭脑补修正了一下。

用1975年5月30日07时南京探空数据验证了一下。差不了太多。那我就当自己是对的了吧。。然后也抄了ncl里面的一些函数。
最后希望懂行的大佬批评指正。
然后github
https://github.com/IncubatorShokuhou/pyMeteo

有空的同学可以帮一下忙。

eeeee · 发表于 2019-6-3 16:16:02

感谢楼主提供的资料，其实可以直接把C代码抄过来然后调用Python的C-API来在Python里面使用的，速度也会快很多，在计算格点的时候特别有用，我这两天尝试一下。

astiny · 发表于 2019-6-3 21:09:38

eeeee 发表于 2019-6-3 16:16
感谢楼主提供的资料，其实可以直接把C代码抄过来然后调用Python的C-API来在Python里面使用的，速度也会快很 ...

其实不是特别用得着。因为对网格数据做这种只有加减乘除的运算的话，用numba编成机械码速度就可以很快。另外人眼ocr辛苦了（

eeeee · 发表于 2019-6-3 21:14:13

astiny 发表于 2019-6-3 21:09
其实不是特别用得着。因为对网格数据做这种只有加减乘除的运算的话，用numba编成机械码速度就可以很快。 ...

我以前试过各种加速方法，包括numba，cython，swig啥的，最后觉得还是cython好用(
numba在有些情况下效果不明显

sujuantao123 · 发表于 2019-7-17 15:55:07

学习一下

		自动登录	找回密码
密码			立即注册

[源代码] 计算沙氏指数的程序

登录后查看更多精彩内容~