Ph
 &inputph 
    prefix='calc'
    outdir='si_nano',
    fildyn='dyn',
    tr2_ph=1.0d-14
    trans=.true.,
    epsil=.false.,
    ldisp=.true.,
    nq1=1,
    nq2=1,
    nq3=8,
    drho_star%open = .true.
    drho_star%ext  = 'drho',
    drho_star%dir  = 'FILDRHO'
/
