24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dwma

铁杆木虫 (正式写手)

应助: 1 (幼儿园)
金币: 8740.1
散金: 94
红花: 7
帖子: 892
在线: 327.8小时
虫号: 429825
注册: 2007-08-05
性别: GG
专业: 中国近代史、现代史

[交流] 【求助】Fortran程序

我的程序如下：（用来对两个文件的数据求和，问题是主程序不能调用子程序但编译没问题）.哪位高手给看看？多谢了！

Program Main
Implicit none
Real,Dimension(5,3):: Value1,Value2,Value3
Integer Dimen1, Dimen2
Integer I,J,K
Dimen1=5
Dimen2=3
Call ReadValue(Dimen1,Dimen2,1000,'file1.txt')
Call ReadValue(Dimen1,Dimen2,2000,'file2.txt')
Do I=1,5
Value3(I,1)=(Value1(I,1)+Value2(I,1))/2
end do
Do J=1,5
Do K=2,3
Value3(J,K)=Value1(J,K)+Value2(J,K)
end do
end do
Open(unit=3000,file='file3')
write (3000,*) ((Value3(J,K),K=1,3),J=1,5)
end

Subroutine ReadValue(Dimen1,Dimen2,UnitNum,Filename)
Integer Dimen1,Dimen2,UnitNum,I1,I2
Character*(*) Filename
Dimension Value(Dimen1,Dimen2)
Open (unit=UnitNum,file=Filename,status='old')
do I1=1,Dimen1
do I2=1,Dimen2
Read(UnitNum,*) Value(I1,I2)
end do
end do
end

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

1楼 2008-04-02 21:26:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xujc1983

木虫 (著名写手)

苦逼青椒一枚

应助: 71 (初中生)
金币: 2065.3
散金: 2373
红花: 18
帖子: 1096
在线: 575小时
虫号: 285496
注册: 2006-10-14
性别: GG
专业: 凝聚态物性I:结构、力学和

★
zzgyb(金币+1,VIP+0):谢谢你的帮助，欢迎继续关注计算模拟版！

引用回帖:

Originally posted by dwma at 2008-4-3 08:39:
多谢了！楼上哥们写的和我原来的一样，就是不行。测试结果程序执行到第一个Call前就终止了。哪位再看看？

是不是你的Fortran程序的问题，还是两个txt文件内数据的格式问题？
这段程序我刚刚调试过，一切正常啊
我的file1是：
1
2
.
.
.
15
file2是：
1
1
.
.
.
1
最后file3是：
1.000000    3.000000    4.000000    2.500000    6.000000
7.000000    4.000000    9.000000    10.00000    5.500000
12.00000    13.00000    7.000000    15.00000    16.00000

赞一下(10人)

回复此楼

5楼2008-04-03 13:15:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 5 个回答

xujc1983

木虫 (著名写手)

苦逼青椒一枚

应助: 71 (初中生)
金币: 2065.3
散金: 2373
红花: 18
帖子: 1096
在线: 575小时
虫号: 285496
注册: 2006-10-14
性别: GG
专业: 凝聚态物性I:结构、力学和

★ ★ ★ ★
zzgyb(金币+4,VIP+0):谢谢你的参与，欢迎继续关注计算模拟版！

你的子程序写的有问题，子程序将文件中的数据保存在Value数组中，但是主程序调用的时候并没有返回该数组，应该调用的时候加上这个参数：
Program Main
Implicit none
Real,Dimension(5,3):: Value1,Value2,Value3
Integer Dimen1, Dimen2
Integer I,J,K
Dimen1=5
Dimen2=3
Call ReadValue(Dimen1,Dimen2,1000,'file1.txt',Value1)
Call ReadValue(Dimen1,Dimen2,2000,'file2.txt',Value2)
Do I=1,5
Value3(I,1)=(Value1(I,1)+Value2(I,1))/2
end do
Do J=1,5
Do K=2,3
Value3(J,K)=Value1(J,K)+Value2(J,K)
end do
end do
Open(unit=3000,file='file3')
write (3000,*) ((Value3(J,K),K=1,3),J=1,5)
end

Subroutine ReadValue(Dimen1,Dimen2,UnitNum,Filename,Value)
Integer Dimen1,Dimen2,UnitNum,I1,I2
Character*(*) Filename
Dimension Value(Dimen1,Dimen2)
Open (unit=UnitNum,file=Filename,status='old')
do I1=1,Dimen1
do I2=1,Dimen2
Read(UnitNum,*) Value(I1,I2)
end do
end do
end

赞一下(10人)

回复此楼

2楼2008-04-03 00:37:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖