请教fortran编程问题

11楼2009-07-05 12:02:16

maomao1210

金虫 (正式写手)

程序强帖: 5
应助: 2 (幼儿园)
金币: 1443.3
散金: 242
红花: 16
沙发: 1
帖子: 990
在线: 441.4小时
虫号: 253215
注册: 2006-05-20
性别: MM
专业: 考古理论

看了一下，总结的不错，支持。

赞一下

12楼2009-07-06 14:59:34

xiejieshuo

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 63.8
帖子: 253
在线: 24小时
虫号: 793394
注册: 2009-06-13
专业: 传热传质学

很好,学习下

13楼2009-07-23 11:16:37

duanyan

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3011.2
散金: 58
红花: 1
帖子: 198
在线: 174.5小时
虫号: 702930
注册: 2009-02-17
性别: MM
专业: 凝聚态物性 II ：电子结构

请教fortran编程问题

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

楼主你好，我现在需要用fortran编个小程序来处理一个数据文件，但我刚接触fortran,看了一周也没搞定，能否请楼主帮忙！这个数据文件名是yl.trace，文件中的数据如下，我的目的是挑出和第二列300对应的第一列和第5列中的数。先谢谢啦！
#    Ef    T [K]          N       DOS(Ef)          S          s/t
0.19535 50.0000    7.18654871  0.10367155E+03  0.35288872E-06  0.44791871E+21
0.19535  100.0000    7.18654516  0.10386861E+03  0.95161005E-06  0.44650443E+21
0.19535  150.0000    7.18644396  0.10418505E+03  0.12966957E-05  0.44624494E+21
0.19535  200.0000    7.18620694  0.10449336E+03  0.16732701E-05  0.44613906E+21
0.19535  250.0000    7.18583102  0.10478008E+03  0.21028510E-05  0.44625119E+21
0.19535  300.0000    7.18533997  0.10499602E+03  0.25607245E-05  0.44658290E+21
0.19535  350.0000    7.18477461  0.10511645E+03  0.30247535E-05  0.44712185E+21
0.19535  400.0000    7.18417905  0.10514146E+03  0.34766052E-05  0.44784913E+21
0.19535  450.0000    7.18359153  0.10508490E+03  0.39010528E-05  0.44874165E+21
0.19535  500.0000    7.18304064  0.10496560E+03  0.42859995E-05  0.44977496E+21
0.19535  550.0000    7.18254516  0.10480245E+03  0.46225413E-05  0.45092558E+21
0.19535  600.0000    7.18211541  0.10461204E+03  0.49048114E-05  0.45217217E+21
0.19535  650.0000    7.18175517  0.10440788E+03  0.51296242E-05  0.45349605E+21
0.19535  700.0000    7.18146358  0.10420040E+03  0.52960040E-05  0.45488108E+21
0.19535  750.0000    7.18123663  0.10399723E+03  0.54046868E-05  0.45631346E+21
0.19535  800.0000    7.18106839  0.10380373E+03  0.54576528E-05  0.45778138E+21
0.19585 50.0000    7.16067850  0.10336246E+03  0.82244420E-06  0.44470893E+21
0.19585  100.0000    7.16056763  0.10398434E+03  0.10391732E-05  0.44505273E+21
0.19585  150.0000    7.16037815  0.10435867E+03  0.13564653E-05  0.44493263E+21
0.19585  200.0000    7.16005609  0.10472610E+03  0.16882630E-05  0.44491950E+21
0.19585  250.0000    7.15960340  0.10505058E+03  0.20613994E-05  0.44506583E+21
0.19585  300.0000    7.15905607  0.10528193E+03  0.24602573E-05  0.44540437E+21
0.19585  350.0000    7.15846067  0.10539912E+03  0.28675326E-05  0.44594088E+21
0.19585  400.0000    7.15786068  0.10540740E+03  0.32678437E-05  0.44666449E+21
0.19585  450.0000    7.15729023  0.10532568E+03  0.36474104E-05  0.44755581E+21
0.19585  500.0000    7.15677269  0.10517706E+03  0.39944573E-05  0.44859217E+21
0.19585  550.0000    7.15632171  0.10498352E+03  0.42997278E-05  0.44975089E+21
0.19585  600.0000    7.15594318  0.10476370E+03  0.45566892E-05  0.45101088E+21
0.19585  650.0000    7.15563736  0.10453236E+03  0.47614020E-05  0.45235331E+21
0.19585  700.0000    7.15540068  0.10430053E+03  0.49121665E-05  0.45376169E+21
0.19585  750.0000    7.15522718  0.10407609E+03  0.50090737E-05  0.45522173E+21
0.19585  800.0000    7.15510961  0.10386438E+03  0.50535506E-05  0.45672105E+21
0.19635 50.0000    7.13474799  0.10426602E+03  0.34658479E-06  0.44397901E+21
0.19635  100.0000    7.13453727  0.10428378E+03  0.10713758E-05  0.44351661E+21
0.19635  150.0000    7.13425552  0.10464092E+03  0.13736112E-05  0.44358722E+21
0.19635  200.0000    7.13383702  0.10504105E+03  0.16814473E-05  0.44369293E+21
0.19635  250.0000    7.13330055  0.10538304E+03  0.20106252E-05  0.44389742E+21
0.19635  300.0000    7.13269548  0.10561000E+03  0.23557447E-05  0.44426331E+21
0.19635  350.0000    7.13207300  0.10570628E+03  0.27080397E-05  0.44481345E+21
0.19635  400.0000    7.13147448  0.10568369E+03  0.30569238E-05  0.44554530E+21
0.19635  450.0000    7.13092868  0.10556650E+03  0.33908719E-05  0.44644373E+21
0.19635  500.0000    7.13045269  0.10538156E+03  0.36988342E-05  0.44748855E+21
0.19635  550.0000    7.13005434  0.10515327E+03  0.39714431E-05  0.44865849E+21
0.19635  600.0000    7.12973470  0.10490165E+03  0.42016593E-05  0.44993325E+21
0.19635  650.0000    7.12949022  0.10464210E+03  0.43848946E-05  0.45129430E+21
0.19635  700.0000    7.12931455  0.10438586E+03  0.45187920E-05  0.45272516E+21
0.19635  750.0000    7.12919977  0.10414070E+03  0.46028358E-05  0.45421131E+21
0.19635  800.0000    7.12913731  0.10391167E+03  0.46379162E-05  0.45574005E+21

[ Last edited by duanyan on 2009-8-5 at 09:47 ]

努力.努力.再努力！！

14楼2009-08-05 09:22:27

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18678.2
红花: 29
帖子: 3803
在线: 1420.7小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
woshilsh(金币+3,VIP+0):辛苦了，常来哦 8-5 19:21

其实大概就是这个样子吧……
read(5,*)a, b, c, d,e,f
if (b .eq. 300) write(6,*) a, e

自己写个循环就可以了……

赞一下(22人)

15楼2009-08-05 17:21:26

duanyan

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3011.2
散金: 58
红花: 1
帖子: 198
在线: 174.5小时
虫号: 702930
注册: 2009-02-17
性别: MM
专业: 凝聚态物性 II ：电子结构

谢谢

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by snoopyzhao at 2009-8-5 17:21:
其实大概就是这个样子吧……
read(5,*)a, b, c, d,e,f
if (b .eq. 300) write(6,*) a, e

自己写个循环就可以了……

不好意思，我是用着fortran才开始学习，看了几天了，现在还在郁闷中，能否麻烦你说的更详细一些，比如如何读取yl.trace文件中的数据，如何把提取的数据存为另外一个文件。。。，

努力.努力.再努力！！

16楼2009-08-05 18:02:37

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18678.2
红花: 29
帖子: 3803
在线: 1420.7小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by duanyan at 2009-8-5 18:02:

不好意思，我是用着fortran才开始学习，看了几天了，现在还在郁闷中，能否麻烦你说的更详细一些，比如如何读取yl.trace文件中的数据，如何把提取的数据存为另外一个文件。。。，

我只会一点 fortran77，下面是个简单的程序，你凑合着看吧。

CODE:

C234567

      PROGRAM YL

      CHARACTER COMMENT*128

      REAL*8 EF, T, N, DOS, S, ST

      OPEN(4,FILE="YL.TRACE",STATUS="OLD")

      OPEN(7,FILE="YL.RESULT",STATUS="NEW")

      READ(4,'(A)') COMMENT

10    READ(4,20,END=40) EF, T, N, DOS, S, ST

      IF (T .EQ. 300.0D0) WRITE(7,30) EF, ST

      GOTO 10

20    FORMAT (F10.5,F10.4,F16.8,3E16.8E2)

30    FORMAT (F10.5,E16.8E2)

40    CONTINUE

      END

17楼2009-08-05 21:10:26

duanyan

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3011.2
散金: 58
红花: 1
帖子: 198
在线: 174.5小时
虫号: 702930
注册: 2009-02-17
性别: MM
专业: 凝聚态物性 II ：电子结构

谢谢！

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by snoopyzhao at 2009-8-5 21:10:

我只会一点 fortran77，下面是个简单的程序，你凑合着看吧。

[code]
C234567
   PROGRAM YL
   CHARACTER COMMENT*128
   REAL*8 EF, T, N, DOS, S, ST
   OPEN(4,FILE="YL.TRACE ...

对我很有帮助，真的非常感谢！

努力.努力.再努力！！

18楼2009-08-06 08:37:39

calos818

木虫 (著名写手)

新人菜鸟

应助: 11 (小学生)
金币: 3293.2
散金: 77
红花: 3
帖子: 1196
在线: 728.6小时
虫号: 698507
注册: 2009-02-08
专业: 凝聚态物性 II ：电子结构

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

FFT里的码位倒置该怎么理解？

19楼2009-08-09 13:19:01

xiejieshuo

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 63.8
帖子: 253
在线: 24小时
虫号: 793394
注册: 2009-06-13
专业: 传热传质学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

我第一次用fortran，花了两三个月编了一个七八百行的程序，现在结果还行

赞一下(2人)

20楼2009-08-10 15:18:26