grads如何读取txt文件

浮生岁月

ret=read('tj.txt')                                             *开始读取路径的txt
rc=sublin(ret,1)              *这句相当于读取头文件，或者是提取逻辑值rc用于判断文件是否存在，必须要的
if(rc>0)                                      *判断文件是否存在
say 'file error 1 'file
return
endif
rec=sublin(ret,2)                 *如果文件存在，开始往下读
y1=subwrd(rec,5)                  *读取上一步所读行的第5列（即，读取纬度）
x1=subwrd(rec,6)                   *读取该行的第6列（即，读取经度） ！！经度对应x1,纬度对应y1
dd=subwrd(rec,4)                  *读取该行

读也是可以读的。。就是麻烦

mikado

一般不会用GrADS直接处理txt文件的，太麻烦了。
如果你的txt是规律格点，用fortran转二进制文件很简单的，一般手册上都有程序例子
注意一下输出时候的顺序就行了

易小凯

兰溪之水 发表于 2012-7-31 07:56

                               
登录/注册后可看大图

会用也是可以的小凯哥，您最近跟爱人跑哪去玩了，都不来招呼兰溪了。。。。

你上次不是说8月份的么

billgates

浮生岁月 发表于 2012-7-31 13:15

                               
登录/注册后可看大图

ret=read('tj.txt')                                             *开始读取路径的txt
rc=sublin(ret,1)  ...

谢谢啊，请问这些步骤可以直接在grads中使用吗？

兰溪之水

易小凯 发表于 2012-7-31 21:03

                               
登录/注册后可看大图

你上次不是说8月份的么

7月底。。。

abd

billgates 发表于 2012-7-31 21:52

                               
登录/注册后可看大图

谢谢啊，请问这些步骤可以直接在grads中使用吗？

直接在grads中使用？大家都建议你不要用grads而用fortran，你还，。。。，看来你对grads有这么强烈的感情啊！呵呵！

Forcast

不会编程，用SURFER较好

wenrenrufeng

读取是可以的，比如用读取经纬度话台风路径。 读取少量数据还可以，读取数组就不划算了。
aa2=read(filename)     文件打开会有两行，第二行是打开的文件数据   第一行是文件打开的状况，若正确打开返回值是0
aa=sublin(aa2,2)

风飘摇

billgates 发表于 2012-7-31 10:24

                               
登录/注册后可看大图

我现在还不会用fortran转化为二进制，再写成ctl

会用fortran读取不？

billgates

wenrenrufeng 发表于 2012-8-1 15:31

                               
登录/注册后可看大图

读取是可以的，比如用读取经纬度话台风路径。 读取少量数据还可以，读取数组就不划算了。
aa2=read(filena ...

谢谢啊，我这里的数据量蛮大，估计要用fortran先转化一下了