【兰溪金盆洗手帖】将GrADS读取的数据写到ASCII文件（txt等）

backham060

谢谢分享，收藏了

MD_ELC

跨专业做ArcGIS需要用到grd文件，从下数据装&学软件打开图已经历时两天了_(:з」∠)_现在就差转ASCII了~谢谢楼主分享！！！

suixinxiu

大神，收下我的膝盖吧

damien1987

谢谢兰溪，我也遇到个问题，我是用第二版的grads2ascii读WRF后处理之后产生的ctl dat文件，ctl描述是附件这样的

输出变量之后得到的txt文件里是504*436的数据，实际有效数据只有201*201，但是txt里包含了大量的缺省值，请问怎样才能只把201*201的有效数据导出到txt里？

兰溪之水

damien1987 发表于 2015-1-16 19:52
谢谢兰溪，我也遇到个问题，我是用第二版的grads2ascii读WRF后处理之后产生的ctl dat文件，ctl描述是附件这 ...

WRF不是等经纬度的，所以GrADS是经过插值，你可以把504和436都改成201，并把pdef那行注释掉，这样就可以输出201*201的等距离数据。

damien1987

兰溪之水 发表于 2015-1-16 20:02
WRF不是等经纬度的，所以GrADS是经过插值，你可以把504和436都改成201，并把pdef那行注释掉，这样就可以 ...

好的，我去试一下，谢谢兰溪！

alexclimate

必须给兰溪多多加分~

alexclimate

必须给兰溪多多加分~

zhuangzi

兰溪大大，有个问题想请您指导下，谢谢ing~~~
我下载了一个日尺度的最低气温的全球数据，文件名为
tasmin_day_GFDL-CM3_rcp45_r1i1p1_20060101-20101231.nc

现在想将此nc文件中的每日的数据都提取出来。已下载您开发的grads2ascii.gs文件。
编写此nc文件的ctl如下：

dset E:\tasmin_day_GFDL-CM3_rcp45_r1i1p1_20060101-20101231.nc

title GFDL-CM3 model output prepared for CMIP5 RCP4.5
undef 1e+20
dtype netcdf
xdef 144 linear 0 2.5
ydef 90 linear -90 2.0
zdef 1 linear 0 1
tdef 1826 linear 12Z1JAN2006 1440mn
vars 1
tasmin=>tasmin 0 t,y,x Daily Minimum Near-Surface Air Temperature
endvars

用这个ctl文件就可以打开nc文件。下一步编写了调用grads2ascii.gs函数的一个名为test.gs的文件，具体内容为：
'reinit '
'open E:\tasmin.ctl'
'E:\grads2ascii tasmin  E:\tasmin.txt [%g 90 1]'
'reinit'
;

用这个test.gs文件运行产生了text文件，可是不能将“tasmin_day_GFDL-CM3_rcp45_r1i1p1_20060101-20101231.nc”文件的每日温度，都生成text 文件。只产生了一个text,我理解的应该产生1826天的1826个text 文件。
已纠结很久，求指导~~~~~谢谢谢谢

yuenai

楼楼真是太棒了！！！~~经常用grads， 但很多功能还不太懂