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

abd

兰溪之水 发表于 2015-5-5 23:56
你看我帖子的具体例子在输出的文本里面也是没有输出这一行文字的，只是在grads的终端显示而已，你把你具 ...

我那个gs文件里输出对应的内容是如下：

*确定需要取平均的区域的经纬度
latb1=38.7
latb2=39.13
lonb1=116.48
lonb2=116.69

'set t 1 90'
'set z 'zi
'/home/grads2ascii.gs  aave(div.1*1e+3,lon='lonb1',lon='lonb2',lat='latb1',lat='latb2')  /home/div/a-div-'z.zi'.txt %f '

最后输出的文件里是这样的：

-0.059037
Printing Grid -- 1 Values -- Undef = 1e+30
-0.054918
Printing Grid -- 1 Values -- Undef = 1e+30
-0.049154
Printing Grid -- 1 Values -- Undef = 1e+30
-0.049711
Printing Grid -- 1 Values -- Undef = 1e+30
-0.049791
Printing Grid -- 1 Values -- Undef = 1e+30
-0.049915
Printing Grid -- 1 Values -- Undef = 1e+30
。。。。。。。

凡不凡

给力，收藏在学习

兰溪之水

abd 发表于 2015-5-6 10:51
我那个gs文件里输出对应的内容是如下：

*确定需要取平均的区域的经纬度

哦，我知道了，你输出时间序列就带了这字符串，你可以用grads2ascii第二版循环输出，
http://bbs.06climate.com/forum.php?mod=viewthread&tid=11515

'set z 'zi
it=1
while(it<=90)
    'set t 'it
    '/home/grads2ascii.gs  aave(div.1*1e+3,lon='lonb1',lon='lonb2',lat='latb1',lat='latb2')  /home/div/a-div-'z.zi'.txt a %f '
    it = it + 1
endwhile

abd

兰溪之水 发表于 2015-5-6 12:54
哦，我知道了，你输出时间序列就带了这字符串，你可以用grads2ascii第二版循环输出，
http://bbs.06clim ...

嗯嗯，好的， 我改了下，现在好了，非常感谢您了！！！！！

abd

兰溪之水 发表于 2015-5-6 12:54
哦，我知道了，你输出时间序列就带了这字符串，你可以用grads2ascii第二版循环输出，
http://bbs.06clim ...

问题好像出在循环那儿，时间和高度层都分别用while循环来输出就可以了，好像不能用简单的“‘set t 1 90’”， 呵呵，多谢了！ 好人一生平安！ O(∩_∩)O~  

edwardouc

太厉害了  赞一个

daizi

好                           

grace639

                                

272861619

攒钱下载一下

grace639

顶