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