24小时热门版块排行榜    

查看: 2246  |  回复: 6

bjwang

铁杆木虫 (正式写手)

[求助] 三维数组换行输入到文件(intel fortran 编译器)

我有一个三维数组r2i(imax,jmx.kmax),想把它写入到一个文件,要求如下:
第一行排r2i(1-imax,1,1)对应的数组元素,然后换行排r2i(1-imax,2,1)等等。
我是按下面来执行的,但达不到目的,输入到文件的列数好像大于imax的值,不知错在何处?哪位高手帮忙解决一下,谢谢!
open(10,file='/home/omega-i.dat',status='new')
do k=1,kmax
do j=1,jmax
write(17,"(1X,E15.6,\)" (r2i(i,j,k),i=1,imax)
enddo
enddo
close(10)

[ Last edited by jjdg on 2011-7-13 at 21:59 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

It is a shame for the soul to be first to give way in this life,when thy body does not give way.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
jjdg(金币+2): 感谢参与 2011-07-22 01:40:29
bjwang(金币+10): 谢谢您的热情回复,其实我已解决问题 2011-07-22 08:25:26
主要是你在 write 那行的格式不正确,你那个一行只能写一个数据,呵呵……
试一下下面这个看看没有什么问题
CODE:
character (len=20) :: fm

open(10,file='/home/omega-i.dat',status='new')

write(fm,'(a,i0,a)') '(1X,',imax,'(E15.6))'

do k=1,kmax
   do j=1,jmax
      write(17,fm) (r2i(i,j,k),i=1,imax)
   enddo
enddo
close(10)

2楼2011-07-21 15:40:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)

★
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 13:13:35
我是利用
  write(17,300) (r2is(i,j,k),i=1,imax)
300 FORMAT (107E15.6)    !107为imax值
来达到目的的!!!!!
  但最后发现,对imax值较小时(目前验证为107)是可行的,但为165时就不行了,照样可以运行和输出,但另一个读取这个输出的软件(软件应该没问题)出错。这和Fortran 90/95标准中规定,当程序代码中的一行超过132个字符时,需要用&续行有关吗?
我修改如下(但编译说无法终止):
if(imax.le.132)then

do k=1,kmax
  do j=1,jmax
   write(17,300) (r2is(i,j,k),i=1,imax)
   300 FORMAT (165E15.6)   !HERE, 107=imax for alter row per imax
   enddo
enddo
close(17)

else

do k=1,kmax
  do j=1,jmax

  write(17,400) (r2is(i,j,k),i=1,130)
  400 FORMAT (130E15.6)

  write(17,500)'&'
   500 FORMAT (A1)

  write(17,600) (r2is(i,j,k),i=131,imax)
  600 FORMAT (39E15.6)  ! imax-130
   enddo
enddo
close(17)

这里两个问题:
1)是Fortran 90/95标准中规定132列造成的问题吗?
2)如果是,下面的程序哪点有问题?
谢谢各位兄弟了!!!
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
3楼2011-07-31 12:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★
jjdg(金币+1): 感谢参与 2011-07-31 14:51:57
问题是你能不能开一个 r2is(165,165,165) 的矩阵?

你的代码中没有哪一行超过 132 字符啊……
4楼2011-07-31 13:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)

ben_ladeng: 楼下有提示哦 2011-07-31 16:29:48
谢谢snoopyzhao,我的代码没有超过132字符的。但我写入的imax列数为165时,这个程序也可正常编译和执行,但输出的文件被另一个程序(要求imax列排在一行的格式)读取时出错,我怀疑是那个输出文件在列数大时格式不对。难道fortran矩阵的大小有限制?
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
5楼2011-07-31 14:58:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

给个完整的程序和测试文件吧,否则搞不清楚你想干咐,以及问题出在哪里……
6楼2011-07-31 16:18:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)

jjdg: 欢迎常来程序语言版 2011-08-04 13:01:26
谢谢兄弟,是我的另一个读取程序的一个bug所致,和换行没有关系。
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
7楼2011-08-03 10:13:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bjwang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 11:46 by YHiw1gou1Eda
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 11:35 by YHiw1gou1Eda
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0TEMgxdT1eIw 2026-11-09 4/200 2026-11-10 11:02 by YHiw1gou1Eda
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 0TEMgxdT1eIw 2026-11-09 4/200 2026-11-10 10:18 by YHiw1gou1Eda
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-10 09:23 by 1yv6QACPoY1E
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-10 01:40 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见