【已解决】fwrite提取wrfout后用grads画图问题求助

陌小默的 · 发表于 2016-10-25 17:43:18

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

x

本帖最后由陌小默的于 2016-10-26 14:03 编辑

用fwrite函数提取wrfout数据，直接输代码画出来的图是正常的，想提取多个层次的数据，写成dat后再用grads画图，结果出来的和原数据差异很大，脚本和结果贴在下面，希望各位大神抽空帮忙看一下，万分感谢！！！

陌小默的 · 发表于 2016-10-25 17:45:44

求各位大神帮忙看看

river · 发表于 2016-10-26 11:43:41

把缺测值改改吧，应该是缺测值的问题。改成-9.99e+33之类的试试吧

陌小默的 · 发表于 2016-10-26 13:07:22

river 发表于 2016-10-26 11:43
把缺测值改改吧，应该是缺测值的问题。改成-9.99e+33之类的试试吧

您好，是提取出来数据的ctl中缺测值麽？我修改了还是一样错误(>﹏<)缺测修改前后没有变化，不知道该怎么解决？十分感谢

陌小默的 · 发表于 2016-10-26 13:07:27

river 发表于 2016-10-26 11:43
把缺测值改改吧，应该是缺测值的问题。改成-9.99e+33之类的试试吧

您好，是提取出来数据的ctl中缺测值麽？我修改了还是一样错误(>﹏<)缺测修改前后没有变化，不知道该怎么解决？十分感谢

陌小默的 · 发表于 2016-10-26 14:03:01

river 发表于 2016-10-26 11:43
把缺测值改改吧，应该是缺测值的问题。改成-9.99e+33之类的试试吧

您好，十分感谢您的提醒，我的问题解决了，就是缺测的问题，提取数据的时候用q fwrite查看output的缺省值，然后对应修改就可以了，十分感谢。

river · 发表于 2016-10-26 18:45:15

陌小默的发表于 2016-10-26 14:03
您好，十分感谢您的提醒，我的问题解决了，就是缺测的问题，提取数据的时候用q fwrite查看output的缺省值 ...

很聪明的办法啊，我都不知道，每次都是自己瞎猜，我也学了一招！谢谢

zys1022 · 发表于 2016-10-26 20:41:15

我用fortran读取grads数据文件dat，怎么相同的经纬度下气温等参数与grads环境下gs文件抽取的数据有差别啊，请大侠们指点迷津。fortran读取的量级、参数都一一对应。不会是gs自动插值引起的吧。不知毛病出在哪

BLUEANGEL2012 · 发表于 2017-7-15 15:43:42

想请教一下楼主提取数据的时候为什么是set x 1 xdef， set y 1 ydef而不是pdef后面的格点数？

lgz · 发表于 2017-7-15 16:03:22

提取数据的时候用q fwrite查看output的缺省值

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

[图形美化] 【已解决】fwrite提取wrfout后用grads画图问题求助