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