grads提取nc文件失败

曦懿-叶子 · 发表于 2019-9-14 16:24:14

最近在学习grads软件，使用是ECMWF上下载的nc格式的数据，在论坛里看了许多的帖子，自己想简单的试试看，但是一直没有成功，又不知道出现的是什么问题，请各位大佬帮一下忙。谢谢！想要提取500hpa的数据，编写了一个gs文件，但是执行时一直说错误。
gs文件如下：

'reinit'
'sdfopen D:\DLW-shuju\shum_1998-2018-2030.nc'
'set gxout fwrite'
'set fwrite D:\DLW-shuju\lianxi\shum-z-500.grd'
'set lon 95 110'
'set lat 20 30'
'set lev 500'
'set t 1 252'
'd q'
'disable fwrite'
'reinit'

QQ截图20190914162239.png

river · 发表于 2019-9-14 21:01:08

你的那个脚本确定在那个路径下？

曦懿-叶子 · 发表于 2019-9-15 10:40:12

river 发表于 2019-9-14 21:01
你的那个脚本确定在那个路径下？

恩恩在的我已经试过很多次所以这个不会错

river · 发表于 2019-9-15 10:51:50

曦懿-叶子发表于 2019-9-15 10:40
恩恩在的我已经试过很多次所以这个不会错

你的这个gs文件的图标是不是现实的是个txt文档的图标？

曦懿-叶子 · 发表于 2019-9-15 15:27:59

对哦我刚刚看了一下是文本文档

帅帅f · 发表于 2019-9-15 20:33:40

直接把D盘的grads文件夹移动到c盘再试试

曦懿-叶子 · 发表于 2019-9-16 13:53:51

帅帅f 发表于 2019-9-15 20:33
直接把D盘的grads文件夹移动到c盘再试试

好的，我试试看行不行

D_Fora · 发表于 2019-9-19 09:44:31

曦懿-叶子发表于 2019-9-16 13:53
好的，我试试看行不行

还是建议换成NCL或者Python、Matlab吧，Grads比较适合初步看图。个人经验而言，花这么多时间折腾Grads有些浪费。

曦懿-叶子 · 发表于 2019-9-19 15:05:35

D_Fora 发表于 2019-9-19 09:44
还是建议换成NCL或者Python、Matlab吧，Grads比较适合初步看图。个人经验而言，花这么多时间折腾Grads有 ...

嗯嗯我也觉得花了很多时间但是结果不会很理想呀

lightmoon · 发表于 2020-1-9 12:21:51

很简单的，你作为数据输出，并不能一次输出一个三维变量，你只需要把时间变成循环就好了。
eg：
'reinit'
'set gxout fwrite'
'set fwrite F:\Data\2019\uwnd.2019.grd'
'sdfopen F:\ERA-Interim\2019.nc'
it=1
'set x 1 144'
while(it<=360)
      'set t 'it
      iz=1
      while(iz<=17)
         'set z 'iz
         'd ave(u,t='it',t='it'+3)'
         iz=iz+1
      endwhile
      it=it+4
endwhile
'disable fwrite'
say 'end'

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

[图形美化] grads提取nc文件失败