-- &inputph tr2_ph = 1.0d-12 prefix = 'SiC' amass(1) = 28.0855 amass(2) = 1.0d12 ! amass(2) = 12.0107 outdir = './' fildyn = 'SiC.dyn' fildvscf = 'dvscf' ! nat_todo = 1 / XQ1 XQ2 XQ3 1