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

shenyuan09

赞！确实是环境变量的问题，我这次添加了“GASCRP，变量值为：c:\OpenGrADS\Contents\Resources\Scripts”
就可以了！
不过，还有一个问题：昨天后来我就将函数拆开来输入，发现程序并不识别result这个参数，但是在兰大的这个程序中却能用，这是为啥呢？

兰溪之水

shenyuan09 发表于 2013-1-31 09:46

                               
登录/注册后可看大图

赞！确实是环境变量的问题，我这次添加了“GASCRP，变量值为：c:\OpenGrADS\Contents\Resources\Scripts”
...

result要在gs脚本里面用~

shenyuan09

非常感谢！我是入门新手，兰大真是帮大忙了！
等把数据读入c++后，我再琢磨琢磨批量的方法，再请教兰大！

shenyuan09

想问一下，为啥脚本读出的TRMM 3B42数据是576800个，不是1440*400个数？然后，数据的排序时从左到右是经度增加，从上到下是纬度减小么，感觉和经纬度对不上的样子，和matlab也对不上啊

shenyuan09

我现在觉得按照读入数据的顺序应该是纬度增加的

shenyuan09

呵呵，写了一个批处理的文件：
'reinit'
i=2
*将需要处理的所有文件名写在filename.txt文本中，空格相连，不回车
readline=read('d:\filename.txt')
*？？为总文件数+2
while (i<？？)
bb = subwrd (readline,i)
rc = write ('d:\filenames.txt',bb)
i = i + 1
'sdfopen d:\misc1\data2\'%bb%'.nc'
*？？为需要输出的变量
'grads2ascii ？？ d:\'%bb%'.txt %5.3f 144 1'
endwhile
say "good bye!"

shenyuan09

发现一个错误，得加上两个reinit，不然，只能打开第一个文件：
位置如下：
'grads2ascii r d:\'%bb%'.txt %5.3f 144 1'
'reinit'
endwhile
'reinit'
say "good bye!"

红河岸

问下版主，这只能处理nc格式么？那像模式输出的数据如grd和dat格式的，我试过了都不行啊

兰溪之水

红河岸 发表于 2013-2-22 10:10

                               
登录/注册后可看大图

问下版主，这只能处理nc格式么？那像模式输出的数据如grd和dat格式的，我试过了都不行啊

只要是grads能正确读取的都行，这就是为啥叫grads2ascii，不叫nc2ascii的原因~

兰溪之水

红河岸 发表于 2013-2-22 10:10

                               
登录/注册后可看大图

问下版主，这只能处理nc格式么？那像模式输出的数据如grd和dat格式的，我试过了都不行啊

只要是grads能正确读取的都行，这就是为啥叫grads2ascii，不叫nc2ascii的原因~