set nrowmax 72 load gs.xml set xc PBE set wf_dyn PSDA set ecutprec 4 set atoms_dyn MD set dt 10 set thermostat BDP set th_temp 400 set th_time 10000 run 10 5 save md.xml