diamond phonon &inputph tr2_ph=1.0d-14, prefix='diamond', outdir='./work/' fildyn='diamond.dyn', amass(1)=12 ldisp=.true. nq1 = 2 nq2 = 2 nq3 = 2 /