09-10-26 Updated Iterative energy balance solution for soil surface

After discover of numerical problems with some parameter combinations causing errors in the energy balance close variable a new default setting have been to estimate fluxes when 40 iterations are not enough to find an energy balance close. Both old and new are based on estimating fluxes from previous periods with a successful energy balance closure. However, new values are also forcing soil temperature and soil vapour at soil surface during such conditions to be equal to air temperature and saturation vapour pressure at air temperature. To trace such conditions when the energy balance is not find by the iterative solution a new auxiliary output : EBalance Iterations has been introduces as part of the soil evaporation module. The change of the model is not expected to produce very small deviations from previous version except for rare combinations of soil surface energy balances.