| 2 | 1/1 | 返回列表 |
| 查看: 878 | 回復(fù): 1 | |||
zll201101金蟲 (正式寫手)
|
[求助]
VASP輸出結(jié)果JCARBX、JCARBY、JCARBZ怎樣作圖?
|
|
VASP計算時輸入文件INCAR加入如下參數(shù)時: 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 輸出結(jié)果會產(chǎn)生JCARBX、JCARBY、JCARBZ三個文件(見附件)。文件內(nèi)容可能是電流密度矢量圖類似的內(nèi)容。請問這結(jié)果用什么軟件作圖? 產(chǎn)生這三個文件的源代碼應(yīng)該是下面這段: ! 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 |
金蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|