added prepare_obs (again), change_parametrisation and plot_model_diff. Models can now be calculated into different parametrisations, so that an inversion based on another parametrisation can be carried out. (for now, only rhovsvp <--> rhomulambda). plot_model_diff makes it easier to visualise the difference between two inversion results, e.g. with and without hard constraints