Some bugs have been corrected. Erroneous calculation of heat flow between layers. Previous version was enhancing heat flow between layer 1 and 2 making the variation to large in the subsoil. A new analytical solution is now an option in the soil heat flow module PEJ