| 查看: 2454 | 回复: 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
» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：用高斯计算，总是出现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
- 专业: 人工晶体
2楼2013-07-18 18:45:13
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 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
alvin510
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1784.3
- 散金: 423
- 红花: 4
- 帖子: 288
- 在线: 225.4小时
- 虫号: 1245360
- 注册: 2011-03-26
- 专业: 人工晶体
5楼2013-07-18 22:50:43
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
6楼2013-07-18 22:54:07
20