24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2149  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Gina88

木虫 (正式写手)

[求助] FORTRAN中如何实现重复读取一行?

如题:fortran中如何实现重复读取一行啊?
第一次读取的时候是将其作为字符串读取的,进行统计,看看数据的个数,然后好定义数组的大小,
第二次读取的时候是将其作为数来读取的,并存入所定义的数组中。
先谢谢大侠了!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by snoopyzhao at 2012-12-16 22:40:57
查一下 Fortran 中有关“内部文件”的用法:

这里举一个小例子:


program test
character(len = 27) :: a = '1.0 2.0 3.0 4.0 5.0 6.0 7.0'
real, dimension(7) :: b

read(a, *) b(1:7)
write(*,*) b ...

您好,我又来找您帮忙了:我写了下面一个fortran的程序。也用您原来教我的write(*,*)查看了,查看了里面的ispin的确是2,因为在91行的循环中我用了i=1,ispin。可是怎么每次i只能等于1,而不能i=2啊,也就是91行的循环怎么不起作用啊?
多次麻烦您,谢谢您了!


dos-procar1.f90(5.78KB)
http://kuai.xunlei.com/d/LnU4DYxJXv3aUAQA6ad?p=130497


DOSCAR(810.18KB)
http://kuai.xunlei.com/d/LnU4DZBJsv3aUAQA876?p=130497


EIGENVAL(27.2KB)
http://kuai.xunlei.com/d/LnU4DZZJ-P3aUAQA9ec?p=130497


PROCAR(1.69MB)
http://kuai.xunlei.com/d/LnU4DZhJHv7aUAQA1d3?p=130497
21楼2012-12-26 21:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Gina88: 金币+5, ★★★★★最佳答案 2012-12-16 23:31:48
jjdg: 金币+1, 感谢参与 2012-12-17 02:25:06
查一下 Fortran 中有关“内部文件”的用法:

这里举一个小例子:
CODE:
program test
character(len = 27) :: a = '1.0 2.0 3.0 4.0 5.0 6.0 7.0'
real, dimension(7) :: b

read(a, *) b(1:7)
write(*,*) b

end program test

» 本帖已获得的红花(最新10朵)

2楼2012-12-16 22:40:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

送鲜花一朵
引用回帖:
2楼: Originally posted by snoopyzhao at 2012-12-16 22:40:57
查一下 Fortran 中有关“内部文件”的用法:

这里举一个小例子:


program test
character(len = 27) :: a = '1.0 2.0 3.0 4.0 5.0 6.0 7.0'
real, dimension(7) :: b

read(a, *) b(1:7)
write(*,*) b ...

嗯,您说的这个方法确实可行,我已经用上了。
呵呵,谢谢啦!

不过那个call system的还没有解决,不知道您会不?
还有就是在用read(28,'(a)')line判断一个很长很长的行的时候,好像判断不出来。不知道您有没有好的办法?比如处理下面的这个文件得到的结果就不正确。

DOSCAR(2.11MB)
http://kuai.xunlei.com/d/LnU4DcS3cujNUAQA891?p=130497
3楼2012-12-16 23:31:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by snoopyzhao at 2012-12-16 22:40:57
查一下 Fortran 中有关“内部文件”的用法:

这里举一个小例子:


program test
character(len = 27) :: a = '1.0 2.0 3.0 4.0 5.0 6.0 7.0'
real, dimension(7) :: b

read(a, *) b(1:7)
write(*,*) b ...

感觉读取上面的那个文件之所以出错,就是因为没有将那个文件中的一行完整的读出来一样。
就不知道要怎么样才能完整的读取一行?
谢谢了!
4楼2012-12-16 23:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见