| 查看: 880 | 回复: 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
» 猜你喜欢
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
-
268求调剂
已经有3人回复
-
311求调剂
已经有10人回复
-
被我言中：新模板不强调格式了，假专家开始管格式了
已经有4人回复
provels
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1534.9
- 散金: 81
- 红花: 1
- 帖子: 448
- 在线: 216.7小时
- 虫号: 639531
- 注册: 2008-10-28
- 专业: 凝聚态物性 II ：电子结构
2楼2021-09-11 14:40:02
