Add new script run_phonons for Phonon calculations.