SiC
&inputph
outdir='./',
prefix='sic',
tr2_ph=1.0d-18,
lmultipole = .true.
fildrho = 'drho.dat.02'
fildvscf = 'drhodv.dat.02'
/
0.02000000 -0.00000000 0.00000000