added final section n the main code, for printing the new initial conditions for a given photon energy
slight adjustments have been made in the rdinput, fourier and spline routines