如何将几百个站点资料输出成一个站点一个TXT

动力气象。

现在是800多个站点资料全在一个TXT力。我现在想让各个站点单独存放在以各自站点命名的TXT里。现在我按站点名输出TXT。但是没法将数据输出到各个TXT里。以下是我编的程序。想法比较简单。但是不知道该怎么修改。

program main
    character*5  sta(803),sta2(12155142)
    real lat(12155142),lon(12155142)
    integer v(12155142),year(12155142),mon(12155142),day(12155142)
    integer j,i
    OPEN(2,FILE='f:\803sta.txt',STATUS='old')
      
     do j=1,803,1
        read(2,*) sta(j)
     enddo
   
     open(3,file='f:\803newave.txt',status='old')
     do i=1,12155142
         read(3,*) sta2(i),lat(i),lon(i),year(i),mon(i),day(i),v(i)
     enddo
     
    do j=1,1,1
       do i=1,121
          open(4,file='f:\2\'//sta(j)//'.txt',status='unknown')
    if(sta(j)==sta2(i)) then
        write(4,'(a6,2x,f6.2,2x,f6.2,2x,i4,2x,i2,2x,i2,2x,i6)') sta2(i),lat(i),lon(i),year(i),mon(i),day(i),v(i)
        end if
     enddo
    enddo
   close(2)
   close(3)
close(4)
     end

月满西楼

写入文件应该用write吧。。。你用的是read。
do j=1,803
    open(4,file='f:\2\'//sta(j)//'.txt',status='unknown')
       do i=1,12155142   
         if(sta(j)==sta2(i)) then
           write(4,'(a6,2x,f6.2,2x,f6.2,2x,i4,2x,i2,2x,i2,2x,i6)') sta2(i),lat(i),lon(i),year(i),mon(i),day(i),v(i)
          end if
        enddo
    enddo

动力气象。

月满西楼 发表于 2016-11-7 16:10
写入文件应该用write吧。。。你用的是read。
do j=1,803
    open(4,file='f:\2\'//sta(j)//'. ...

编了一天。有点晕。。。。。

jolincai

循环啊 站点等于多少就输出成一个文件夹 批处理。

动力气象。

修改以后这程序对是对了。但是运算好慢。。。十几个小时都只是输出了一百多个站点。。这效率太低了有木有什么改进的方法的

言深深

在读取文件的时候就下功夫，增加数组的维度，不等于这个站点的，read到下一个维度里面去，这样程序快的多。

动力气象。

言深深 发表于 2016-11-8 17:12
在读取文件的时候就下功夫，增加数组的维度，不等于这个站点的，read到下一个维度里面去，这样程序快的多。

读文件倒不是重点。这个应该是输出的时候内循环太大了，导致时间一直很久，跑了两天了，还没跑完。不过您说的这种方法我不是很明白，可以的话还麻烦您稍微详细的解释下。不甚感激

言深深

动力气象。 发表于 2016-11-9 12:42
读文件倒不是重点。这个应该是输出的时候内循环太大了，导致时间一直很久，跑了两天了，还没跑完。不过您 ...

不对，你把一堆的东西呼噜噜read到一个数组里面，就不太好处理，循环增加了很多···
如果在read数据的时候，增加维度，就容易处理，比如下面这个数组，一维的存储成12个，二维的就是3*4，如果这时候需要输出等于23的，那么一维情况下你要比较6次得到结果，二维的比较5次就得到了，在数据量非常大的时候，高维数据更占优势
11
12
13
21
22
23
31
32
33
41
42
43

动力气象。

言深深 发表于 2016-11-9 15:55
不对，你把一堆的东西呼噜噜read到一个数组里面，就不太好处理，循环增加了很多···
如果在read数据的 ...

谢谢哈。我好好研究研究

Hilda_ZT

请问你是怎么下载到站点的数据的？
现在找数据找了很久还是到处找不到