Major update

- separation of mass and energy balance calculations
- new main subroutine which needs to be called (calls in other routines need to change that!)
- loop over energy balance with included energy resulting from melting/refreezing
Update to account for changes in 'surface_physics.f90'