-- &inputph tr2_ph = 1.0d-12 prefix = 'diam' amass(1) = 12.01078 outdir = './' fildyn = 'diam.dyn' fildvscf = 'dvscf' / XQ1 XQ2 XQ3