Working scripts for pmode modelling