| 查看: 1035 | 回复: 1 | ||
zll201101金虫 (正式写手)
|
[求助]
VASP输出结果JCARBX、JCARBY、JCARBZ怎样作图？
|
|
VASP计算时输入文件INCAR加入如下参数时： LCHIMAG=.TRUE. !use linear response to calculate shielding tensor,perturbation theory with respect to B field LWRTCUR=.TRUE. !write current response to B field to file LBFCONST=.TRUE. !constant B field appiled BCONST= 10.0000 10.0000 10.0000 !applied magnetic field 输出结果会产生JCARBX、JCARBY、JCARBZ三个文件（见附件）。文件内容可能是电流密度矢量图类似的内容。请问这结果用什么软件作图？ 产生这三个文件的源代码应该是下面这段： ! we might want to look at the current IF (LWRTCUR) THEN IF (IO%IU6>=0) THEN OPEN(UNIT=99,FILE=DIR_APP(1  IR_LEN)//'JCAR'//DIR_TEXT(BDIR),STATUS='UNKNOWN')! write header CALL OUTPOS(99,.FALSE.,INFO%SZNAM1,T_INFO,LATT_CUR%SCALE,LATT_CUR%A,.FALSE.,T_INFO%POSION) WRITE(99,'(3I5)') GRID_SOFT%NGX,GRID_SOFT%NGY,GRID_SOFT%NGZ WRITE(99,'(5(1X,E17.11))') (0._q, I=1,GRID_SOFT%NGX*GRID_SOFT%NGY*GRID_SOFT%NGZ) DO J=1,3 WRITE(99,'(3I5)') GRID_SOFT%NGX,GRID_SOFT%NGY,GRID_SOFT%NGZ WRITE(99,'(5(1X,E17.11))') & & (REAL(S_BARE(I,J),q),I=1,GRID_SOFT%NGX*GRID_SOFT%NGY*GRID_SOFT%NGZ) ENDDO CLOSE(99) ENDIF ENDIF |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : JCARBX
- 附件 2 : JCARBY
- 附件 3 : JCARBZ
2019-08-24 15:27:39, 3.24 M
2019-08-24 15:28:06, 3.24 M
2019-08-24 15:28:33, 3.24 M
» 猜你喜欢
极限拉扯，如何抉择
已经有10人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
provels
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1534.9
- 散金: 81
- 红花: 1
- 帖子: 448
- 在线: 216.7小时
- 虫号: 639531
- 注册: 2008-10-28
- 专业: 凝聚态物性 II ：电子结构
2楼2021-09-11 14:40:02
