爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5871|回复: 2

[求助] fortran 如何读入多个文件数据 将每个文件数据存一列

[复制链接]

新浪微博达人勋

发表于 2012-9-10 16:51:30 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
现在有大量的数据文件,批量读取都没有问题,原先是读取一个就对应输出一个文件,现在想将这些数据综合到一个文件当中,可是却遇到了问题,向高手请教!!

我想将每个文件里的数据存放于一列,如一个文件的数据放在第一列,第二个放在第二列,可是我输出了之后全部还是多行一列的数据,不知道哪里的逻辑关系搞错了,向大家请教,谢谢!

程序如下:希望m()的数据按列存放
open(10,file='../mref/0531.txt',form='formatted')
!input the filename
        
  read(10,*) hs
  
  do i=1,hs
        read(10,*) sourcename(i),stalon(i),stolon(i),stalat(i),stolat(i)
  enddo
!-----------------------------------------------------------------------
!output the file
        
  open(22,file='../mref/20.txt',form='formatted')

!-----------------------------------------------------------------------
!input the file
  do i1=1,hs
        close(11)
        open(11,file='../mref/'//sourcename(i1),form='binary')

。。。。。。。

    do k=1,RadarDataFileHeader%nZNumGrids
     m(i1)=0
     n(i1)=0
     l(i1)=0
  do i=ysta,ysto           
          do j=xsta,xsto         
         if (tt(k,i,j)>=5) then      
         m(i1)=m(i1)=+1                        
      endif                                 
      if (tt(k,i,j)>=30) then      
         n(i1)=n(i1)+1                        
      endif                                
      if (tt(k,i,j)>=40) then     
         l(i1)=l(i1)+1
      endif
   enddo
  enddo

  write(22,*) k,m(i1)
enddo
   
  enddo
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2012-9-10 16:57:55 | 显示全部楼层
用隐式循环就可以了
do i=1,m
write(1,"(<n>f10.2)")(x(i,j),j=1,n)
enddo
出来的就是m行n列
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-9-11 15:09:23 | 显示全部楼层
open里面文件号也可以用循环变量,循环打开 读取 就可以了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表