-- &inputph prefix = 'diam', epsil = .false., fildyn = 'diam.dyn', ldisp = .true. fildvscf = 'dvscf' nq1=6, nq2=6, nq3=6, tr2_ph = 1.0d-12, /