| 查看: 843 | 回复: 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
» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
provels
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1534.9
- 散金: 81
- 红花: 1
- 帖子: 448
- 在线: 216.7小时
- 虫号: 639531
- 注册: 2008-10-28
- 专业: 凝聚态物性 II ：电子结构
2楼2021-09-11 14:40:02