-- &inputph tr2_ph = 1.0d-10 prefix = 'Pb' amass(1) = 207.2 outdir = './' fildyn = 'Pb.dyn' fildvscf = 'dvscf' / XQ1 XQ2 XQ3