多要素，多时次站点资料转格点资料

chnlhttom

参考mofangbao的单时次站点资料处理方法http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903，修改得到多时次，多要素站点资料处理方法。

1. Program sta2gr
   
2. Implicit none
   
3. !这里是程序的变量声明
   
4. Character*8 stid
   
5. Real lon,lat, ave,max,min,tim
   
6. integer nlev,flag,i,j,sta,time
   
7. !变量声明结束
   
8. !程序开始,输入时次，站点数
   
9. time=6
   
10. sta=37
    
11. 
    
12. tim=0.0
    
13. flag=1
    
14. 
    
15. Open(1,file='test.txt',status='old')
    
16. open(2,file='test.grd',status='replace',form='binary')
    
17. ！一行一行读取数据，写入数据
    
18. j=1
    
19. 100 if(j<=time)then
    
20. i=1
    
21. 10 if(i<=sta)then
    
22. nlev=1
    
23. Read(1,*)stid,lat,lon,ave,max,min
    
24. Write(2)stid,lat,lon,tim,nlev,flag,ave,max,min
    
25. 
    
26. i=i+1
    
27. Goto 10
    
28. 
    
29. else
    
30. nlev=0
    
31. Write(2)stid,lat,lon,tim,nlev,flag
    
32. j=j+1
    
33. end if
    
34. goto 100
    
35. else
    
36. close(1)
    
37. close(2)
    
38. end if
    
39. !程序结束
    
40. End

复制代码

多时次数据并未使用数组的读取方式，而是采用判断语句编写了一个计数器，每读取完一个时次的所有站点数据后写结束语句然后继续读取下一个站点的数据。
测试数据排列方式是站号，纬度，经度，变量1（Ave），变量2（max），变量3（min），一个时次所有数据然后接下一个时次所有数据。

lulu77

参考下。。。。

lgz

参考下。。。。

Elric

楼主，如果一个txt文件为单时次，但是需要读取多时次（也就是多个txt文件，怎么做啊）program test
    implicit none
    integer i
    character(40)::filename='F:\bishe\1998-2016\1998\98010108.000'
    character*8 stid
    real lon,lat,var,rain,time
    integer nlev,flag
    do i=101,131
    write(filename(28:30),'(i3)')i
    time=0.0
    nlev=1
    flag=1
open(1,file=filename,status='old')
open(2,file='f:\bishe\sta.grd',form='binary')
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
10  read(1,*,end=100)stid,lon,lat,var,rain
    write(2)stid,lon,lat,var,rain
    goto 10
100 continue
    close(1)
    nlev=0
    print*,stid,lon,lat,var,rain
    write(2)stid,lat,lon,time,nlev,flag
    close(2)
    end do
    end
    我是这么写的的，感觉循环不对啊，可否指导一下

Elric

=。=求指导

爱喝柠檬汁

50136        52.58        122.31        2016        1        1        4        0        4                               
50136        52.58        122.31        2016        1        2        0        0        0                               
50136        52.58        122.31        2016        1        3        3        0        3                               
50136        52.58        122.31        2016        1        4        0        0        0                               
50136        52.58        122.31        2016        1        5        8        5        13                               
50136        52.58        122.31        2016        1        6        0        0        0                               
50136        52.58        122.31        2016        1        7        0        0        0                               
50136        52.58        122.31        2016        1        8        0        0        0                               
50136        52.58        122.31        2016        1        9        2        0        2                               
50136        52.58        122.31        2016        1        10        0        0        0
50136        52.58        122.31        2016        1        11        0        0        0                               
50136        52.58        122.31        2016        1        12        0        0        0                               
50136        52.58        122.31        2016        1        13        0        0        0                               
50136        52.58        122.31        2016        1        14        0        32700        32700                               
50136        52.58        122.31        2016        1        15        1        32700        1                               
50136        52.58        122.31        2016        1        16        0        0        0                               
50136        52.58        122.31        2016        1        17        3        2        5                               
50136        52.58        122.31        2016        1        18        32700        1        1                               
50136        52.58        122.31        2016        1        19        0        0        0                               
50136        52.58        122.31        2016        1        20        0        0        0                               
50136        52.58        122.31        2016        1        21        0        0        0                               
50136        52.58        122.31        2016        1        22        0        0        0                               
50136        52.58        122.31        2016        1        23        0        0        0                               
50136        52.58        122.31        2016        1        24        0        0        0                               
50136        52.58        122.31        2016        1        25        32700        5        5                               
50136        52.58        122.31        2016        1        26        0        5        5                               
50136        52.58        122.31        2016        1        27        2        0        2                               
50136        52.58        122.31        2016        1        28        0        0        0                               
50136        52.58        122.31        2016        1        29        0        0        0                               
50136        52.58        122.31        2016        1        30        0        0        0                               
50136        52.58        122.31        2016        1        31        0        0        0

爱喝柠檬汁

楼主请问这样的数据应该怎样编写fortran程序，我的资料是全国的，有839个站，而且每个站都是一天一行

gaoyuan215

谢谢分享。正在试验。