24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3458  |  回复: 21

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
6楼: Originally posted by physliebe at 2014-02-09 19:33:13
WRITE(12,"(3f10.5)"; ) A,这种输出矩阵的方法试了一下,报错...

那是因为这个论坛的什么表情的缘故,我是直接抄你的代码的。我的意思就是那句话给你的Fortran是按列读取的错觉。
11楼2014-02-10 02:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
8楼: Originally posted by duhb655 at 2014-02-09 23:52:07
不好意思,BB没了,就送红花一朵。这个我好像解决了,但有个问题,原TXT中有的列是整型有的是实型,我这样一下全给别人弄成实型了,麻烦您看一下7楼,不知道您能不能给点意见,谢谢啦!...

你要想整型的读成整型,实型的读成实型,那就得用类似 2I5,2F7.5 这种方式来读取了……
12楼2014-02-10 02:27:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
8楼: Originally posted by duhb655 at 2014-02-09 23:52:07
不好意思,BB没了,就送红花一朵。这个我好像解决了,但有个问题,原TXT中有的列是整型有的是实型,我这样一下全给别人弄成实型了,麻烦您看一下7楼,不知道您能不能给点意见,谢谢啦!...

前面的回复不对

如果你用一个矩阵来存贮前面的那个文本,你只能指定一种类型,在你这种情况下自然只能是实型……
13楼2014-02-10 02:35:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


jjdg: 金币+1, 辛苦了 2014-02-13 11:52:05
引用回帖:
7楼: Originally posted by duhb655 at 2014-02-09 23:47:30
zhao老师,要不您帮帮忙,我给出一个小例子,然后我说出想要做什么(很简单),然后您帮我写一下主要的句子我学习学习,行不。TXT数据如下:
50     40     28.30186          54.27303
51     41     28.28920  ...

一个 A 是没有办法既存整弄,又存实型的。分成两个可以:
CODE:
PROGRAM test

   IMPLICIT NONE
   INTEGER :: i, ios
   INTEGER, DIMENSION(11,2) :: c1_2
   REAL, DIMENSION(11,2) :: c3_4
   
   OPEN(UNIT = 12, FILE = 'aa.txt')
   
   DO i=1,11
      READ(12,'(I2,TR5,I2,TR5,F8.5,TR10,F8.5)') c1_2(i,:), c3_4(i,:)
   END DO
   
   DO i=1,11
      WRITE(*, '(2(TR2,I2),2(TR2,F8.5))') c1_2(i,:), c3_4(i,:)
   END DO

END PROGRAM test

14楼2014-02-10 02:45:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duhb655

金虫 (著名写手)

引用回帖:
14楼: Originally posted by snoopyzhao at 2014-02-10 02:45:07
一个 A 是没有办法既存整弄,又存实型的。分成两个可以:


PROGRAM test

   IMPLICIT NONE
   INTEGER :: i, ios
   INTEGER, DIMENSION(11,2) :: c1_2
   REAL, DIMENSION(11,2) :: c3_4
   
   OPEN ...

谢谢大神,那么晚还没休息啊!按您的指导我分别存在了A和B,那么能把A,B分别按他们的类型一起写到一个文件里吗?我试了一下,错的啊,有一半对的,一半会有*******号;
15楼2014-02-11 00:09:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
15楼: Originally posted by duhb655 at 2014-02-11 00:09:27
谢谢大神,那么晚还没休息啊!按您的指导我分别存在了A和B,那么能把A,B分别按他们的类型一起写到一个文件里吗?我试了一下,错的啊,有一半对的,一半会有*******号;...

把写文件的代码贴出来……

******的这个,多半是位置不够用
16楼2014-02-11 00:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duhb655

金虫 (著名写手)

引用回帖:
16楼: Originally posted by snoopyzhao at 2014-02-11 00:16:24
把写文件的代码贴出来……

******的这个,多半是位置不够用...

program test
implicit none
integer  :: a(11,2),i
real     :: b(11,2)

open(1,file='Result.txt')
do i=1,11
read(1,*) a(i,,b(i,
enddo
open(2,file='zre.txt')
print*,a(:,1)
write(2,"(2i6,2f10.5)" (a(i,,i=1,11),(b(i,,i=1,11)
end
结果:
    50    40   0.00000   0.00000
    52    42   0.00000   0.00000
    54    44   0.00000   0.00000
    56    46   0.00000   0.00000
    58    48   0.00000   0.00000
    60    50  28.30186  54.27303
************  28.27396  54.26752
************  28.23691  54.28119
************  28.22947  54.30813
************  28.20773  54.31368
************  28.15063  54.28722
17楼2014-02-11 00:48:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
17楼: Originally posted by duhb655 at 2014-02-11 00:48:43
program test
implicit none
integer  :: a(11,2),i
real     :: b(11,2)

open(1,file='Result.txt')
do i=1,11
read(1,*) a(i,,b(i,
enddo
open(2,file='zre.txt')
print*,a(:,1)
write(2," ...

你这个程序运行的时候没有报错?好奇怪……虽然 gfortran 编译没有问题,但运行会报错

Fortran runtime error: Expected REAL for item 3 in formatted transfer, got INTEGER

你最后那个 write 语句是不对的。那句话的意思是先输出 a,再输出 b,所以……
18楼2014-02-11 01:23:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duhb655

金虫 (著名写手)

引用回帖:
18楼: Originally posted by snoopyzhao at 2014-02-11 01:23:56
你这个程序运行的时候没有报错?好奇怪……虽然 gfortran 编译没有问题,但运行会报错

Fortran runtime error: Expected REAL for item 3 in formatted transfer, got INTEGER

你最后那个 write 语句是不对的 ...

刚从家里回到学校,没有报错哈,嘿嘿。那我想得到的那种结果的WRITE应该怎么写啊?
19楼2014-02-11 21:51:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bennett_14

木虫 (正式写手)

引用回帖:
5楼: Originally posted by polypro at 2014-02-09 15:07:37
接3L,将write那一句替换就行了
DO I=1,11
   write(13,*) A(I,1:3)
ENDDO

高手在这儿,此楼完全正确。解答及理由充分清晰。顶
君子当自强不息,死而后已!
20楼2014-02-11 23:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 duhb655 的主题更新
信息提示
请填处理意见