fortran版diamond4转grads的源程序

霹雳 · 发表于 2013-10-29 12:43:47

下载来看看

Figo · 发表于 2013-11-4 12:43:37

x=int(head(13)+0.0001)
y=int(head(14)+0.0001)
请问这两句是什么意思？
谢谢！

mofangbao · 发表于 2013-11-4 14:20:13

Figo 发表于 2013-11-4 12:43

登录/注册后可看大图

x=int(head(13)+0.0001)
y=int(head(14)+0.0001)
请问这两句是什么意思？

强制类型转换至于为何要加0.0001，是因为我当时用的编译器在不加的情况下会有些精度问题你可以试着去掉

Figo · 发表于 2013-11-4 17:34:35

mofangbao 发表于 2013-11-4 14:20

登录/注册后可看大图

强制类型转换至于为何要加0.0001，是因为我当时用的编译器在不加的情况下会有些精度问题你可以试着去掉 ...

谢谢回复！

风云翩跹 · 发表于 2013-11-5 10:10:37

感谢楼主！

风云翩跹 · 发表于 2013-11-5 11:05:47

清风大大，我运行了你的程序出现了和这个人一样的问题，麻烦问下怎么解决？
http://bbs.06climate.com/forum.php?mod=viewthread&tid=12408

mofangbao · 发表于 2013-11-5 12:24:19

风云翩跹发表于 2013-11-5 11:05

登录/注册后可看大图

清风大大，我运行了你的程序出现了和这个人一样的问题，麻烦问下怎么解决？
http://bbs.06climate.com/for ...

这个程序写的通用性并不是很强，主要是文件头的读取上有的4类文件头不是这样行列分布的，就要自己重写读取的文件头的部分，通过参数个数来判断会更好，有错误调试就行

风云翩跹 · 发表于 2013-11-5 12:49:46

mofangbao 发表于 2013-11-5 12:24

登录/注册后可看大图

这个程序写的通用性并不是很强，主要是文件头的读取上有的4类文件头不是这样行列分布的，就要自己重写读 ...

大大，我这个500hpa预报资料是这样的，您看怎么设置相关参数才能让我的grd文件不出y方向的错误啊

13110420.024 (48.47 KB, 下载次数: 8)

mofangbao · 发表于 2013-11-5 14:17:24

风云翩跹发表于 2013-11-5 12:49

登录/注册后可看大图

大大，我这个500hpa预报资料是这样的，您看怎么设置相关参数才能让我的grd文件不出y方向的错误啊

通过文件头可以看出来，你要自己分析，然后试着自己写吧

249557221 · 发表于 2013-11-14 16:37:06

正遇到这个问题，感谢大神

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

[源代码] fortran版diamond4转grads的源程序