| 查看: 1001 | 回复: 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
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
provels
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1534.9
- 散金: 81
- 红花: 1
- 帖子: 448
- 在线: 216.7小时
- 虫号: 639531
- 注册: 2008-10-28
- 专业: 凝聚态物性 II ：电子结构
2楼2021-09-11 14:40:02
20