Dear VASP community:
I used LPARD=.TRUE. to write band and k-point decomposed charge density (vasp 5.4.4). When plotting generated PARCHG files with VESTA, I noticed that some plots have a bug that appears as a constant 3D charge density after a certain coordinate. The reason is a mismatch between the value and the format that results in "***********" output as shown below:
oleg$ grep "*" PARCHG*
PARCHG: 0.61416E-01 -2.3446 -4.8390 -7.1151 -8.8709 -9.8715 *********** -9.2820 -7.8766 -6.0362
To solve the problem, I edited "pardens.F" file and replaced in line 1157
CALL OUTCHG(GRIDC,iunit,.FALSE.,CHTOT)
with
CALL OUTCHG(GRIDC,iunit,.TRUE.,CHTOT)
This enables "LLONG" output format for the charge densities in the file "fileio.F"
IF (LLONG) THEN
FORM='(1(1X,E17.11))'
NWRITE=5
ELSE
FORM='(1(1X,G11.5))'
NWRITE=10
ENDIF
It seems that G11.5 may not be sufficient in some cases.
I hope this will be helpful
Oleg
Output format issue with PARCHG
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 12
- Joined: Fri Jan 27, 2017 11:25 pm
- License Nr.: 5-2663
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
Re: Output format issue with PARCHG
Thank you. The code will be modified according to your suggestion.