爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 22591|回复: 11

[分享资料] 求助:nc数据多变量转grd文件出错

[复制链接]
回帖奖励 1 金钱 回复本帖可获得 1 金钱奖励! 每人限 1 次(中奖概率 10%)

新浪微博达人勋

发表于 2021-7-6 17:35:15 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 lily2012_01 于 2021-7-7 11:19 编辑

最近在用nc风场数据,发现多变量写入GRD文件报错。不知道是GRD生成有问题,还是给GRD写的ctl文件写的有问题,请各位大拿看看哪里有问题,谢谢!
gs文件如下:
'sdfopen g:/ERA5/input-nc/20200517-21.nc'
'set gxout fwrite'
'set fwrite g:/ERA5/output-wind/out20200517-21.grd'
'set x 1 1440'
'set y 1 721'

tt=1
while (tt <= 144)
'set t 'tt''
'set z 1'
'd u10'
'd v10'
'd fg10'
'd i10fg'

tt=tt+1
endwhile
'disable fwrite'
'c'

给GRD文件写的ctl信息如下:
dset g:/ERA5/output-wind/out20200517-21.grd
title wind 20200517-21
undef -9.99e+33
dtype netcdf
xdef 1440 linear 0 0.25
ydef 721 linear -90 0.25
zdef 1 linear 0 1
tdef 144 linear 00Z17MAY2020 60mn
vars 4
u10=>u10  0  t,y,x  10 metre U wind component
v10=>v10  0  t,y,x  10 metre V wind component
fg10=>fg10  0  t,y,x  10 metre wind gust since previous post-processing
i10fg=>i10fg  0  t,y,x  Instantaneous 10 metre wind gust
endvars

用GRADS打开生成的GRD文件报错如下:


报错PNG.PNG
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-7-6 23:54:41 | 显示全部楼层
为啥要写grd呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-7-7 09:38:02 | 显示全部楼层
fortran读取方便
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-7 10:05:28 | 显示全部楼层
问题是什么呢?报错内容是什么?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-7-7 11:18:38 | 显示全部楼层
lqouc 发表于 2021-7-7 10:05
问题是什么呢?报错内容是什么?

ga-> open g:/ERA5/output-wind/out20200517-21.ctl
Scanning description file:  g:/ERA5/output-wind/out20200517-21.ctl
Error: nc_open failed to open file g:/ERA5/output-wind/out20200517-21.grd
NetCDF: Unknown file format
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-7-7 14:24:16 | 显示全部楼层
lqouc 发表于 2021-7-7 10:05
问题是什么呢?报错内容是什么?

打开GRD文件的ctl报错如下:
报错PNG.PNG
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-7-7 15:11:59 | 显示全部楼层
找到原因了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-17 16:39:27 | 显示全部楼层
dtype netcdf  
很明显,你写了个普通二进制文件,非要在描述文件里告诉GrADS它是NC格式的,所以就打不开
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-22 09:24:23 | 显示全部楼层
我在数据转换时,出现如下报错信息,怎么解决大神。
Traceback (most recent call last):
  File "grd2nc.py", line 3, in <module>
    open_CtlDataset('E:/Study/Data_assimilation/paper/model2data/DATA_back45/nature/t30.ctl').to_netcdf('output.nc')
  File "/mnt/e/study/Data_assimilation/paper/model2data/DATA_back45/result/xgrads/io.py", line 103, in open_CtlDataset
    ctl = CtlDescriptor(encoding=encoding, file=desfile)
  File "/mnt/e/study/Data_assimilation/paper/model2data/DATA_back45/result/xgrads/core.py", line 107, in __init__
    if os.path.getsize(abspath) / (1024.0*1024.0) > 2:
  File "/home/changmh/anaconda3/lib/python3.8/genericpath.py", line 50, in getsize
    return os.stat(filename).st_size
FileNotFoundError: [Errno 2] No such file or directory: 'E:/Study/Data_assimilation/paper/model2data/DATA_back45/nature/t30.ctl'
(base) changmh@changmh:/mnt/e/study/Data_assimilation/paper/model2data/DATA_back45/result$
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-25 22:12:14 | 显示全部楼层
学习了
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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