| 查看: 5720 | 回复: 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» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
| 顶一下，感谢分享！ |
20楼2013-08-05 16:46:21
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
