Output format issue with PARCHG

Problems running VASP, runtime errors, technical questions.
Please check whether the answer to your question is given in
the VASP online manual or has been discussed in this forum
previously!

Moderators: Moderator, Global Moderator

Output format issue with PARCHG

Postby rubel » Mon Mar 26, 2018 12:10 am

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
rubel
Newbie
Newbie
 
Posts: 5
Joined: Fri Jan 27, 2017 11:25 pm

Re: Output format issue with PARCHG

Postby admin » Thu Apr 05, 2018 8:01 am

Thank you. The code will be modified according to your suggestion.
admin
Administrator
Administrator
 
Posts: 2907
Joined: Tue Aug 03, 2004 8:18 am


Return to Bugreports

Who is online

Users browsing this forum: No registered users and 1 guest