| 查看: 2963 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助成功】Band－procar 和Dos－procar 程序执行出错 已有8人参与
|
|||
|
我从版内下载的分析PROCAR的程序Band－procar.f 和Dos－procar.f 通过g77 编译成可执行文件。但是执行后出现segmentation fault (core dumped)错误提示。这是为什么啊？我上网查了一下这表示指针使用出错。我的编程水平太差实在查不出问题在哪。麻烦哪位高人指点一下，十分感谢！以下是源程序 implicit real*8(a-h,o-z) parameter (nbd = 200) parameter (nkd = 500) parameter (nxd = 300) parameter (natmd = 20) dimension a(3,3),b(3,3),c(3),e(nkd,nbd),sk(nkd,3) dimension xx(nxd) ,wei(nkd) dimension dump(20),oc(nkd,nbd,natmd,4) open(7,file='PROCAR',form='FORMATTED',status='OLD') pi = 3.141592654 read(7,103) dump write(*,*) 'Spin polarized calculation? (no=1,yes=2):' read (*,*) ispin if ((ispin.ne.1).and.(ispin.ne.2)) then write(*,*) ' INPUT ERROR, ispin must equal to 1 or 2 ' stop endif c write(*,*) 'Enter # of interval (npoints) and division (ndiv):' c read (*,*) npoints,ndiv open(9,file='KPOINTS',form='FORMATTED',status='OLD') read(9,100) temp read(9,*) ndiv write(*,*) 'Enter the range of energy to plot:' read (*,*) er1,er2 emin=min(er1,er2) emax=max(er1,er2) write(*,*) 'Enter the value of fermi energy:' read(*,*) ef if (ispin.eq.1) then open(11,file='band.dat') elseif (ispin.eq.2) then open(11,file='band-up.dat') open(12,file='band-dn.dat') endif open(8,file='POSCAR',form='FORMATTED',status='OLD') read(8,100) temp c write(6,100) temp 100 format(20a4) read (8,*) aa c WRITE(6,*) aa c c *** read lattice constant from POSCAR** c do i=1,3 read (8,*) (a(i,j),j=1,3) c WRITE(6,500) (a(i,j),j=1,3) 500 format (3f12.8) enddo do i=1,3 do j=1,3 a(i,j)=aa*a(i,j) enddo c WRITE(6,500) (a(i,j),j=1,3) enddo c c *** read lattice vector from POSCAR*** c volume=a(1,1)*a(2,2)*a(3,3)+a(1,2)*a(2,3)*a(3,1) &+a(1,3)*a(2,1)*a(3,2)-a(1,1)*a(2,3)*a(3,2) &-a(1,2)*a(2,1)*a(3,3)-a(1,3)*a(2,2)*a(3,1) do i=1,3 if (i .eq. 1) then j=2 k=3 else if (i .eq. 2) then j=3 k=1 else j=1 k=2 endif c(1)=a(j,2)*a(k,3)-a(j,3)*a(k,2) c(2)=a(j,3)*a(k,1)-a(j,1)*a(k,3) c(3)=a(j,1)*a(k,2)-a(j,2)*a(k,1) do j=1,3 b(i,j)=2*pi*c(j)/volume c WRITE (6,*) b(i,j) enddo enddo do 9000 isp=1,ispin read(7,104) nk,nband,nion do 1000 k = 1,nk read(7,103) dump read(7,105) kp,(sk(k,j),j=1,3),wei(k) c write(6,105) kp,(sk(k,j),j=1,3),wei(k) read(7,103) dump do nb = 1,nband read(7,106) nb1,e(k,nb),occ c write(6,106) nb1,e(k,nb),occ read(7,103) dump read(7,103) dump c write(6,*) 'nion=',nion niont = nion +1 if (nion .eq. 1) niont = 1 do ion = 1,niont read(7,107) (oc(k,nb,ion,j),j=1,4) c write(6,107) (oc(k,nb,ion,j),j=1,4) enddo read(7,103) dump c write(6,103) dump enddo 1000 continue weight = 0.0 do k = 1, nk weight = weight + wei(k) enddo do k = 1,nk wei(k) = wei(k) / weight enddo 101 format(10x,f9.5) 102 format(f10.5) 103 format(20a4) 104 format(16x,i3,20x,i5,19x,i4) 105 format(10x,i3,5x,3f11.8,13x,f11.8) 106 format(4x,i4,9x,f14.8,7x,f12.8) 107 format(3x,4f7.3) c c *** find reciprocal lattice vector *** xx(1) = 0.0 nn = 1 do k = 1,nk-1 dkx=(sk(k+1,1)-sk(k,1))*b(1,1) + (sk(k+1,2)-sk(k,2))*b(2,1) & + (sk(k+1,3)- sk(k,3))*b(3,1) dky=(sk(k+1,1)-sk(k,1))*b(1,2) + (sk(k+1,2)-sk(k,2))*b(2,2) & + (sk(k+1,3)- sk(k,3))*b(3,2) dkz=(sk(k+1,1)-sk(k,1))*b(1,3) + (sk(k+1,2)-sk(k,2))*b(2,3) & + (sk(k+1,3)- sk(k,3))*b(3,3) del = sqrt ( dkx**2 + dky**2 + dkz**2 ) nn = nn +1 xx(nn) = xx(nn-1) + del enddo do n=1,nband if (mod(n,2).ne.0) then do k=1,nk ee = e(k,n) - ef if ( ee .gt. emax ) ee = emax if ( ee .lt. emin ) ee = emin write (10+isp,300) xx(k),ee enddo elseif (mod(n,2).eq.0) then do i=nk,1,-1 ee = e(i,n) - ef if ( ee .gt. emax ) ee = emax if ( ee .lt. emin ) ee = emin write (10+isp,300) xx(i),ee enddo endif enddo 300 format (f12.8,2x,f12.8) if (mod(nband,2) .ne. 0) then write (10+isp,300) xx(nk),emin write (10+isp,300) xx(1),emin else write (10+isp,300) xx(1),emin endif c c *** write xx-ee *** c npoints=nk/ndiv do n=2,npoints kk=(n-1)*ndiv write (10+isp,300) xx(kk),emin write (10+isp,300) xx(kk),emax write (10+isp,300) xx(kk),emin enddo write (10+isp,300) xx(nk),emin write (10+isp,300) xx(nk),emax write (10+isp,300) xx(1),emax write (10+isp,300) xx(1),emin zero=0.0 write (10+isp,300) xx(1),zero write (10+isp,300) xx(nk),zero 9000 continue stop end [ Last edited by chk0521 on 2010-6-25 at 16:36 ] |
回复此楼» 猜你喜欢
271材料工程求调剂
已经有5人回复
-
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
-
材料工程专硕调剂
已经有6人回复
-
一志愿天大材料与化工（085600）总分338
已经有4人回复
-
085700资源与环境308求调剂
已经有3人回复
-
求材料调剂
已经有8人回复
-
294求调剂材料与化工专硕
已经有5人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有4人回复
-
一志愿吉林大学材料学硕321求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 二硫苏糖醇加入还原性胶束中的具体操作方法和加入量。
已经有10人回复
- 请问高分子和有机药物合成对身体影响大吗
已经有8人回复
- 求助：岩石薄片岩矿分析和XRD得到的各矿物成分及含量有较大差异，怎么解释？
已经有6人回复
- 循环伏安 循环极化和动电位再活化曲线有什么分别？
已经有9人回复
- 水文局和海事局，知道的来说一下
已经有13人回复
- 得不到纯产物，哪出错了，亲们帮忙看看啊
已经有3人回复
- 一个二倍体植物的全基因组测序程序及花费？
已经有18人回复
- 用c语言编程怎来表达时滞！谢谢！
已经有6人回复
- 请问学位论文查重会不会和自己的文章重复呀？
已经有11人回复
- Fortran的subroutine和function有什么区别，是在于输入输出的特点吗？
已经有13人回复
- 关于CASTEP 中计算时出错
已经有6人回复
- 出错 warning#2 #180
已经有7人回复
- MS一打开就出错，任何任务都不能运行
已经有12人回复
- DOS图中怎么判断非金属性，半金属性，金属性
已经有7人回复
- 请教：orca 中用CCSD(T)计算单点能出错--提示内存不足 该怎么办？
已经有10人回复
- MS软件老出错，求指导
已经有3人回复
- 制备丙烯酸乳液，十二烷基磺酸钠和十二烷基苯磺酸钠能否通用，大侠们帮帮忙啊~~~
已经有8人回复
repent03
金虫 (小有名气)
- 1ST强帖: 1
- 应助: 3 (幼儿园)
- 金币: 1339.3
- 散金: 6
- 帖子: 87
- 在线: 40小时
- 虫号: 1068065
- 注册: 2010-08-02
- 专业: 半导体材料
9楼2012-03-22 00:45:22
lvjian8596
木虫 (小有名气)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 0.014
- 金币: 3531.7
- 红花: 2
- 帖子: 296
- 在线: 209.1小时
- 虫号: 545009
- 注册: 2008-04-13
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2010-06-24 08:44:53
3楼2010-06-24 15:51:51
calos818
木虫 (著名写手)
新人菜鸟
- 应助: 11 (小学生)
- 金币: 3125.7
- 散金: 77
- 红花: 3
- 帖子: 1201
- 在线: 730.8小时
- 虫号: 698507
- 注册: 2009-02-08
- 专业: 凝聚态物性 II ：电子结构
4楼2010-10-24 18:23:05
