* ThermalEquilibrium.f90: check forces only for
   potential SLAB-POTENTIAL
 * IndependentOscillatorsModel.f90: in case the
   interpolating function for the SD has small negative
   values, an absolute value has been added