SPECFEM1D comparison script