| 查看: 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» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有20人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有48人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
21楼2014-02-25 17:43:50
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
