| 查看: 1350 | 回复: 2 | |||
damao4361556铁虫 (正式写手)
|
[交流]
PWSCF 算完能带 和 dos 后的数据 如何转换成origin 能打开的 格式？ 已有1人参与
|
|
PWSCF 算完能带 和 dos 后的数据 如何转换成origin 能打开的 格式？ 求分享 有关的小代码 自带的 plotband.x 等 感觉 不好用。 |
回复此楼» 猜你喜欢
急需调剂
已经有3人回复
-
22408 312求调剂
已经有21人回复
-
恳请有学校收留
已经有4人回复
-
297，工科调剂?
已经有9人回复
-
300求调剂
已经有4人回复
-
291求调剂
已经有5人回复
-
085404 22408 309分求调剂
已经有10人回复
-
296求调剂
已经有13人回复
-
通信工程求调剂！！！
已经有6人回复
-
294求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何选择能带计算软件和入门能带计算
已经有7人回复
- 学习PWSCF(二)：态密度的计算
已经有7人回复
- pwscf-5.0得到的声子谱很乱是怎么回事？
已经有8人回复
- QE-PWscf计算声子谱流程。
已经有10人回复
- Problem - nscf - c_bands: eigenvalues not converged
已经有20人回复
- [软件交流]您在使用哪些免费视图、分析软件呢？
已经有23人回复
- pwscf4.3.2画能带图
已经有14人回复
- 电子态密度怎么画图啊？
已经有4人回复
- 【求助】ms能带图
已经有11人回复
- 【求助】band.dat用origin画能带图
已经有5人回复
- 【求助】pwscf的声子态密度图和色散曲线图能用origin画出来吗？
已经有5人回复
- 【重点讨论】计算能带时如何确定高对称点和k-points path?
已经有40人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
dos算完的结果是可以直接导入origin的，band算完的程序要用一个小程序转换下。下面这个程序是我从小木虫找到的。我改了下，有两个输入文件一个band.dat,还有一个ev写费米能级。 program prog real, allocatable :: e(:,  real, allocatable :: k(:, ,dk( ,x(real, dimension(3) ::k0,a real ev character(len=32):: xx, yy, zz,evv integer :: nbands, i, j, n, nk open(10, file='band.dat', status='old') open(22, file='ev', status='old') read(10,*) xx, yy, nbands, zz, nk read(22,*) evv read(evv,*) ev write(*,*)"nk=",nk allocate(k(nk,3)) allocate(e(nk,nbands)) allocate(dk(0:nk)) allocate(x(1:nk)) open(11,file='phoband.dat') do i=1,nk read(10,*) (k(i,j),j=1,3) read(10,*) (e(i,n),n=1,nbands) enddo do n=1,nk dk(0)=0 if (n.eq.1) then k0=k(n, endif a=k(n, -k0 dk(n)=dk(n-1)+sqrt(dot_product(a,a)) x(n)=dk(n) k0=k(n, enddo do i=1,nk write(11, 8888) x(i), (e(i,j)-ev,j=1,nbands) enddo 8888 format (1x, f9.4, 2x ,1000(f12.6,2x)) end program prog |
2楼2015-05-01 20:58:08
damao4361556
铁虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 38.8
- 红花: 3
- 帖子: 479
- 在线: 435.8小时
- 虫号: 917440
- 注册: 2009-12-01
- 专业: 理论和计算化学
|
你好 非常感谢您的回复 请问 为什么我这个 编译总出错呢？ 是直接 编译 gfortran ***.f 吗？ [wl_903@hpc-login3 band]$ gfortran Getband.f Getband.f:1.2: program prog 1 Error: Non-numeric character in statement label at (1) Getband.f:1.2: program prog 1 Error: Unclassifiable statement at (1) Getband.f:2.31: real, allocatable :: e(:, 1 Error: Missing ')' in statement at or before (1) Getband.f:3.42: real, allocatable :: k(:, ,dk( ,x( 1 Error: Missing ')' in statement at or before (1) Getband.f:14.18: allocate(k(nk,3)) 1 Error: Syntax error in ALLOCATE statement at (1) Getband.f:15.18: allocate(e(nk,nbands)) 1 Error: Syntax error in ALLOCATE statement at (1) Getband.f:16.19: allocate(dk(0:nk)) 1 Error: Syntax error in ALLOCATE statement at (1) Getband.f:17.18: allocate(x(1:nk)) 1 Error: Syntax error in ALLOCATE statement at (1) Getband.f:20.16: read(10,*) (k(i,j),j=1,3) 1 Error: Expected variable in READ statement at (1) Getband.f:21.16: read(10,*) (e(i,n),n=1,nbands) 1 Error: Expected variable in READ statement at (1) Getband.f:26.11: dk(0)=0 1 Error: Unclassifiable statement at (1) Getband.f:28.72: k0=k(n, 1 Error: Syntax error in argument list at (1) Getband.f:30.72: a=k(n,-k0 1 Error: Syntax error in argument list at (1) Getband.f:31.12: dk(n)=dk(n-1)+sqrt(dot_product(a,a)) 1 Error: Statement function at (1) is recursive Getband.f:32.72: x(n)=dk(n) 1 Error: Unexpected STATEMENT FUNCTION statement at (1) Getband.f:33.72: k0=k(n, 1 Error: Syntax error in argument list at (1) Getband.f:38.6: 8888 format (1x, f9.4, 2x ,1000(f12.6,2x)) 1 Error: Bad continuation line at (1) Getband.f:38.6: 8888 format (1x, f9.4, 2x ,1000(f12.6,2x)) 1 Error: Unclassifiable statement at (1) Getband.f:39.22: end program prog 1 Error: Syntax error in END PROGRAM statement at (1) Error: Unexpected end of file in 'Getband.f' 需要当前目录存在两个文件是吗？ 一个 是ev 里面的数据是费米能级 另外一个是 band.dat 里面的内容是 QE 算能带后得到的？ |
3楼2015-05-04 06:13:03
