Clean up of the harmonics.f90
The output can now handle second derivatives