关于grads软件如何读取文本文件的问题

pengyulong

最近，在做台风数值预报的时候，遇到一个问题。我已经知道了一些点（大概有8000多个）的经纬度信息和时间信息，如何根据这些信息，求取fnl资料的的信息。如果用循环的话，这些经纬度信息保存在2文本文件中，如何利用gs编程将其读取存储到一个变量中？求各位大侠赐教。

四叶草

求取fnl资料的的信息？这没看懂是什么意思？gs脚本读取文本文档里的信息，大致就是用到read(name),sublin,subwrd等命令。

river

read (filename)
This functions reads individual records from file filename. Repeated calls must be made to read consecutive records. The result is a string containing two lines: the first line is the return code, the 2nd line is the record read from the file. The record may be a maximum of 80 characters. Use the sublin function to separate the result. Return codes are:
0 - ok
1 - open error
2 - end of file
8 - file open for write
9 - I/O error
Files are opened when the first call to read is made for a particular file name. Files are closed when the execution of the script file terminates (note that files remain open between function calls, etc).
write (filename, record <, append>)
This functions writes records to output file filename. On the first call to write for a particular file, the file is opened in write mode. This will destroy an existing file! If you use the optional append flag, the file will be opened in append mode, and all writes will be appended to the end of the file. Return codes are:
0 - ok
1 - open error
8 - file open for read

pengyulong

我要读的文件里包含的时间，经度，纬度的信息有8000多条，怎么分别读取到三个字符串变量中啊？麻烦高手贴个代码？

pengyulong

四叶草 发表于 2015-12-9 20:18
求取fnl资料的的信息？这没看懂是什么意思？gs脚本读取文本文档里的信息，大致就是用到read(name),sublin,s ...

是读取，打错了。我想用grads按给定的时间，经纬度信息读取grib1数据,但是这些经纬度信息都是杂乱无章的（8000多条）.所以想利用read()函数读取这些经纬度信息到字符串，但是不知道怎么操作？

pengyulong

现在，我想到一种思路，但是代码老是报错。我把代码贴在下面：希望大家给我指点一下：

1. 'reinit'
   
2. 'c'
   
3. 'open fnl_2000.ctl'
   
4. lonstr='\t'
   
5. while(1)
   
6. r=read('longstr.txt')
   
7. temp=subwrd(r,2)
   
8. if(temp!=0);break;endif;
   
9. line=sublin(r,1);
   
10. lonstr=lonstr line;
    
11. endwhile
    
12. ;
    
13. 
    

复制代码

四叶草

pengyulong 发表于 2015-12-10 15:34
是读取，打错了。我想用grads按给定的时间，经纬度信息读取grib1数据,但是这些经纬度信息都是杂乱无章的 ...

不太明白楼主的意思！你的杂乱无章的数据是台风路径吗？你最终目的是要干什么？grib1数据做什么用？

pengyulong

四叶草 发表于 2015-12-10 17:11
不太明白楼主的意思！你的杂乱无章的数据是台风路径吗？你最终目的是要干什么？grib1数据做什么用？

是的。我要提取台风经过的点的一些类似U风，V风的实时数据。做它台风强度预测用的。

四叶草

pengyulong 发表于 2015-12-10 17:31
是的。我要提取台风经过的点的一些类似U风，V风的实时数据。做它台风强度预测用的。

这样说就明白楼主的意思了！不过处理起来有点儿复杂。个人认为，根据路径点提取点上的U、V实测数据，可以根据绘制台风路径的gs脚本来。楼主，你要通过时间去匹配FNL数据，并打开相应的文件？

pengyulong

四叶草 发表于 2015-12-10 19:29
这样说就明白楼主的意思了！不过处理起来有点儿复杂。个人认为，根据路径点提取点上的U、V实测数据，可以 ...

是啊，我手里有台风经过的经纬度以及时间的信息，现在想提取fnl实测代码数据，现在那个gs代码，还是有点问题。可以帮忙看一下吗？

1. "reinit"
   
2. "set display color white";*设置背景色为白色,默认为黑色
   
3. "c"
   
4. "open fnl_20.ctl"
   
5. FileName='lonStr.txt'
   
6. rf=readfile(FileName)
   
7. str=subwrd(rf,2)*这样，能提取成功！但是我想用while循环依次提取str中的数字，可我一加上while循环就报错。很无语。
   
8. say str
   
9. ;*定义函数
   
10. function readfile(FileName)
    
11. *读取整个ASCII文本文件,最多读10行
    
12. Lines=''
    
13. while(1);
    
14.     r=read(FileName)
    
15.     temp=subwrd(r,1);*第一行仅有一个单词,为状态码,读取无错误应状态码为0
    
16.     if (temp != 0);break;endif;*如果读取有误或读到EOF就退出
    
17.     line=sublin(r,2);*提取该行字符串
    
18.     Lines=Lines' ' line ' '
    
19. endwhile
    
20. *返回文本字符串
    
21. return(Lines)
    
22. ;
    

复制代码