用GrADS读取日尺度的nc数据，转换成每日的text文件

zhuangzi

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

现在想将此nc文件中的每日的数据都提取出来。已下载一个能将将GrADS读取的数据写到ASCII文件grads2ascii.gs文件。但是只能转换出一个text，和1826天不符。我觉得转换出的是一个年均值，我需要的是每日的值。首先，我编写此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文件。下一步编写了一个名为test.gs 的文件来调用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 文件。
已纠结很久，求指导~~~~~谢谢谢谢

传说中的谁

你让grads写到一个文件他当然就只写一个啦

zhuangzi

传说中的谁 发表于 2015-3-12 15:35
你让grads写到一个文件他当然就只写一个啦

啊，有道理，我试试去，谢谢您~~
其实我还在考虑是不是不用在ctl文件中设置下t 的时间了是吗？

muguangyuze

不懂啊，也同问

冷水舟

麻烦你帮我看看我提问的grads转化为nc数据的问题~~~~谢谢{:eb303:}{:eb303:}

猫的他

用ncl的asciiwrite岂不是更简单