| 查看: 5664 | 回复: 23 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
PWSCF 具体实例 详细解答各种问题 已有13人参与
|
|||
|
PWSCF 具体实例 详细解答各种问题，摘自网络一个blogblog 地址 http://zhaoyunawuli.blog.163.com ... 976200811147143279/ 网易的 网名： 随梦想而飞 1,附 处理PWSCF中计算得到本征频率的程序 自己修改了 用的时候需要把*.freq 的文件 手动改成input.freq program prog real, allocatable :: e(:,  real, allocatable :: k(:, ,dk( ,x(real, dimension(3) ::k0,a character(len=32):: xx, yy, zz integer :: nbands, i, j, n, nk open(10, file='input.freq', status='old') read(10,*) xx, yy, nbands, zz, nk 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),j=1,nbands) enddo 8888 format (1x, f9.4, 2x ,1000(f12.6,2x)) end program prog 2,处理能带的程序中在 do i=1, nhighk-1 前一行 加上 ntkp=0 其它的我暂时还没用到，有什么问题希望共同交流[ Last edited by lu_wei07 on 2013-3-5 at 16:37 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : PWSCF教程，从入门到精通.docx
- 附件 2 : phfreq.exe
2013-03-04 20:04:39, 45.86 K
2013-03-05 16:37:30, 349 K
» 本帖已获得的红花（最新10朵）
renren123123» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有34人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有78人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有12人回复
-
FileCode能看出啥？
已经有30人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
关于Filecode分析方法
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
23楼2016-03-07 17:47:34
0.8
| 顶一下，感谢分享！ |
2楼2013-03-04 21:23:36
hakuna
木虫 (知名作家)
- 1ST强帖: 3
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
3楼2013-03-05 06:49:10
0.8
| 谢谢分享！！！！！！！！！ |
4楼2013-03-05 06:56:35
