用grads把nc转为grd遇到的问题

bigbingr

用文件grads把e降水资料的nc文件写成grd出现的错误

Efwrite.gs文件如下：
'reinit'
'sdfopen d:\paper\data\e008543m1980_2012_pre_chi_remapcon_24h_muldpm_seas(1).nc'
'set gxout fwrite'
'set fwrite d:\paper\redata\epre.grd'

'set lat 18 53.5'
'set lon 72.5 136'
'set t 1 133'
'd var142'
'disable fwrite'
;
得到了一个epre.grd文件，用grads打开时，出现如下问题：
file:///C:/Users/hp/AppData/Local/Temp/ksohtml/wpsC86D.tmp.jpg
不知道怎么办了。。。。。。。。,
求大神解救，谢谢啦

赤炼修

给epre.grd写个ctl文件，再打开这个ctl文件

bigbingr

赤炼修 发表于 2016-4-16 10:00
给epre.grd写个ctl文件，再打开这个ctl文件

用过了，还是不行
*the ctl of efwrite's epre
dset d:\paper\redata\epre.grd
title the average data of season from 1980 to 2012
dtype netcdf
undef -9e+33
xdef 128 linear 72.5 0.5
ydef 72 linear 18 0.5
zdef 1 linear 0 1
tdef 133 linear 00:33Z01JAN2001 3mo
vars 1
var142 0 0 rain
endvars

洪hong

你先直接在GRADS黑色框框里面输入语句：sdfopen （你的nc文件）（注意路径)
然后用语句：q ctlinfo 查询nc的ctl描述文件，从中知道变量名，而不是你写的'd var142',根据查询到的变量名和维数设定，再写一个新的gs文件提取NC文件变成grd文件。

赤炼修

参照楼上。如果打不开nc文件，可用ncl 的filedump功能查看nc文件信息。

bigbingr

洪hong 发表于 2016-4-16 22:52
你先直接在GRADS黑色框框里面输入语句：sdfopen （你的nc文件）（注意路径)
然后用语句：q ctlinfo 查询nc ...

变量名就是我用了q ctlinfo查出来的，确实是var142,然后ctl文件也是我照着内置ctl写的，可是就是不行，我就不知道了

bigbingr

赤炼修 发表于 2016-4-17 00:17
参照楼上。如果打不开nc文件，可用ncl 的filedump功能查看nc文件信息。

ncl不会啊

洪hong

洪hong 发表于 2016-4-16 22:52
你先直接在GRADS黑色框框里面输入语句：sdfopen （你的nc文件）（注意路径)
然后用语句：q ctlinfo 查询nc ...

你nc的提取，我建议这样提取：
'sdfopen D:\paper\data\*****(你的文件名).nc'
'set gxout fwrite'
'set fwrite d:\paper\redata\epre.grd'
'set lon 72.5 136'
'set lat 18 53.5'
'set lev (你需要的层次，有850,500,200hPa等等，如果只有一层就写：‘set z 1’)'
'set t 1 133'
'd var142'
'disable fwrite'
;

lehuoxiaoyu

学习学习{:eb502:}{:eb502:}

淡定冷静

一模一样的问题，怎么解决