| 查看: 5722 | 回复: 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人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
青基如何破局
已经有4人回复
-
基金不中，共勉
已经有8人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
刘仕晨
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3276.8
- 散金: 6
- 红花: 11
- 帖子: 319
- 在线: 179.2小时
- 虫号: 2201333
- 注册: 2012-12-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
13楼2013-04-01 22:51:16
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
