- 积分
- 38418
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-20
- 最后登录
- 1970-1-1
成长值: 0
|
发表于 2012-4-21 10:06:49
|
显示全部楼层
孩子,你的文件名和读规则和你的数据不匹配,经过如下调整可以运行了,错误的地方在标记线之间
program main
implicit none
integer i,j,k
real x(10,12,33)
integer mydate !!! 读取的第一个文件!!!
character(len=80) filename,filedate !!! 用于存放文件名称!!!
mydate=1979
do k=1,33
write(filedate,'(i4.4)') mydate !!!!把mydate的值赋予filedate!!!
!--------------------------------------------------------------------
write(filename,*) 'cmap_mon_',trim(filedate),'_rain1.txt' !@~~~~~~漏了“rain1”~~~~~~@!
!--------------------------------------------------------------------
open(unit=k+50,file=filename)
do j=1,12
!--------------------------------------------------------------------
read(k+50,*) (x(i,j,k),i=1,10) !@~~~~~~此处读取方式需要和你的文件相匹配,隐循环表示从该行读取10个~~~~~~@!
!--------------------------------------------------------------------
end do
close(k+50)
mydate=mydate+1
end do
open(unit=100,file='station.txt',form="formatted",access="sequential",position="append")
do k=1,33
do j=1,12
write(100,'(10f8.2)') x(1,j,k),x(2,j,k),x(3,j,k),x(4,j,k),x(5,j,k),x(6,j,k),x(7,j,k),x(8,j,k),x(9,j,k),x(10,j,k)
end do
end do
stop
end |
|