| 查看: 2738 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
alvin510银虫 (小有名气)
|
[求助]
Fortran编程过程中遇到错误，求大侠帮忙看下
|
||
|
我在用FORTRAN计算NLO过程中，编的一段程序遇到错误，请高手帮忙看下，万分感谢。 其中的一段程序为： module def implicit none integer :: nk,nkpts,ndum,ns,nspins,nb,locat_str character (len=80) :: dummy real num_e real fermi_energy integer num_bands real,allocatable :: weight(  real,allocatable :: kpoints(:, end module program main !读入 kpoints, weights 和 eigenvalues use def implicit none integer,parameter :: fileid=10 character (len=80) :: filename="F:\FORTRAN\NLOC\CGA\unit_bands.txt" logical alive inquire (file=filename,exist=alive) if(.not. alive) then write(*,*) trim(filename), "doesn't exist!" stop end if open(unit=10,file=filename) read(unit=10,"(19XI6,/,26XI2,/,20XF7.1,/,22XI7,/,31XF13.8)"  nkpts,spins,num_e,num_bands,fermi_energyallocate(kpoints(1,nkpts)) allocate(kpoints(2,nkpts)) allocate(kpoints(3,nkpts)) allocate(weight(nkpts)) do nk = 1, nkpts read(fileid, '(a)') dummy locat_str = index(dummy,'K-point') read(dummy(locat_str+7 ,*) ndum, kpoints(1,nk), kpoints(2,nk), kpoints(3,nk), weight(nk)do ns = 1, nspins read(fileid,*) dummy do nb = 1, num_bands read(fileid,*) eigenvalues(nb, nk, ns) enddo enddo enddo close(fileid) stop end program 其中调用了一个数据文件unit_bands.txt。这个文件我已经上传了上来。 再次感谢你的帮助！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : unit_bands.txt
- 附件 2 : CGA-01.txt
2013-07-18 18:42:01, 39.01 K
2013-07-18 18:43:19, 1.17 K
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有87人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：用高斯计算，总是出现2070错误代码，请大侠们帮帮忙
已经有5人回复
- MATLAB, 代码错误，请大侠帮忙看看。
已经有5人回复
- 大侠们~本人编译vasp.5.lib过程中遇到如下错误 敢问何解？不胜感激~
已经有7人回复
- win8安装Endnote失败，显示错误1723，求高手大侠解决啊
已经有4人回复
- VS20？？+intel visual fortran2011XE做并行计算的，能介绍一下经验吗？
已经有19人回复
- 急求fortran运行错误原因，在线等
已经有7人回复
- 求大侠帮忙解决一个B类错误！！谢谢
已经有5人回复
- Perl 生物信息学编程求助，运行出现问题，烦请帮忙查找一下错误，谢谢！
已经有3人回复
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 有会看fortran程序的，能帮忙看看这个程序吗？
已经有10人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 求助Fortran编程
已经有8人回复
- MAC ox 下怎么用fortran 编程？
已经有9人回复
- 求助fortran编程时怎么使用数据库
已经有4人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】fortran编程问题
已经有3人回复
- 【求助】运行fortran程序出现错误 : error LNK2019 怎么回事啊
已经有8人回复
- 【求助】同样一个问题用matlab和fortran编程计算，结果相差比较大，不知道为什么
已经有4人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】ANSYS 中命令流编程（fortran）
已经有11人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【分享】C/C++常用算法设计方法,fortran其他语言编程也可以参考下
已经有8人回复
alvin510
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1784.3
- 散金: 423
- 红花: 4
- 帖子: 288
- 在线: 225.4小时
- 虫号: 1245360
- 注册: 2011-03-26
- 专业: 人工晶体
5楼2013-07-18 22:50:43
alvin510
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1784.3
- 散金: 423
- 红花: 4
- 帖子: 288
- 在线: 225.4小时
- 虫号: 1245360
- 注册: 2011-03-26
- 专业: 人工晶体
2楼2013-07-18 18:45:13
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
3楼2013-07-18 21:10:42
alvin510
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1784.3
- 散金: 423
- 红花: 4
- 帖子: 288
- 在线: 225.4小时
- 虫号: 1245360
- 注册: 2011-03-26
- 专业: 人工晶体
|
谢谢你发现这个问题。后来我自己又发现了几个问题一并做了更改，发现还是出现问题了。我的代码如下： module def implicit none integer :: nk,nkpts,ndum,ns,nspins,nb,locat_str character (len=80) :: dummy real num_e real fermi_energy integer num_bands real,allocatable :: weight( real,allocatable :: kpoints(:, real,allocatable :: eigenvalues(:,:, end module program main !读入 kpoints, weights 和 eigenvalues use def implicit none integer,parameter :: fileid=10 character (len=80) :: filename="F:\FORTRAN\NLOC\CGA\unit_bands.txt" logical alive inquire (file=filename,exist=alive) if(.not. alive) then write(*,*) trim(filename), "doesn't exist!" stop end if open(unit=fileid,file=filename) read(unit=fileid,"(19XI6,/,26XI2,/,20XF7.1,/,22XI7,/,31XF13.8)" nkpts,nspins,num_e,num_bands,fermi_energy!read(unit=fileid,"(19XI6)" nkpts!read(unit=fileid,"(26XI2)" nspins!read(unit=fileid,"(20XF7.1)" num_e!read(unit=fileid,"(22XI7)" num_bands!read(unit=fileid,"(31XF13.8)" fermi_energyallocate(kpoints(3,nkpts)) allocate(weight(nkpts)) allocate(eigenvalues(num_bands,nkpts,nspins)) do nk = 1, nkpts read(fileid, '(a)') dummy locat_str = index(dummy,'K-point') read(dummy(locat_str+7 ,*) ndum, kpoints(1,nk), kpoints(2,nk), kpoints(3,nk), weight(nk)do ns = 1, nspins read(fileid,*) dummy do nb = 1, num_bands read(fileid,*) eigenvalues(nb, nk, ns) enddo enddo enddo close(fileid) stop end program |
4楼2013-07-18 22:47:56
