爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 4051|回复: 3

[脚本编辑] fortran批量读取数据

[复制链接]
发表于 2018-5-3 09:20:52 | 显示全部楼层 |阅读模式

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

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

x
每年一个文件,每个文件是153个数的一维时间序列,我想用FORTRAN把数据都读到一个A数组里,可不知为什么读出来的数都是重复的,哪位大神能帮我看看程序出了什么问题?不胜感激!!

PARAMETER (yrn=34,nt=153)
        real A(yrn,nt),B(yrn,nt),sum,SUMA,SUMB
        integer yr,t
      character*5::year(yrn)=(/'1981','1982','1983','1984','1985',
     &'1986','1987','1988','1989','1990','1991','1992','1993','1994',
     &'1995','1996','1997','1998','1999','2000','2001','2002','2003',
     &'2004','2005','2006','2007','2008','2009','2010','2011','2012',
     &'2013','2014'/)
       
      do yr=1,yrn
         do t=1,nt
           A(yr,t)=0.0
           B(yr,t)=0.0
         enddo
        enddo

      do yr=1,34
        OPEN (1,FILE='O:\qj\850uaave'//trim(year(yr))//'.txt')
      OPEN (2,FILE='O:\qj\dp850uaave'//trim(year(yr))//'.txt')

      read(1,*) (A(yr,t),t=1,nt)
      read(2,*) (B(yr,t),t=1,nt)

        close(1)
        close(2)
      enddo

        write(*,*) A
        write(*,*)B

end[img=0,1]file:///C:\Documents and Settings\Administrator\Application Data\Tencent\Users\760928732\QQ\WinTemp\RichOle\G]YKW)[1EPGLWUII[PNLO35.png[/img] G]YKW)[1EPGLWUII[PNLO35.png

密码修改失败请联系微信:mofangbao
发表于 2018-5-3 12:36:49 | 显示全部楼层
看着很高大上啊
密码修改失败请联系微信:mofangbao
发表于 2018-5-3 12:37:31 | 显示全部楼层
下次好好的操作一下
密码修改失败请联系微信:mofangbao
发表于 2018-5-3 15:36:01 | 显示全部楼层
嗯,默认write数组是按列优先的,所以 这个结果里表示的是 第1年的第1个数,第2年的第1个数……第34年的第1个数……
可以试试改成按行优先输出
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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