use wind3d_horiz_coordinates global parameter to determine dpar('datacoords')
heights method was not negating depths
Note that due to the bug in heights, where depths were returned as positive
numbers, the results returned by dat.Ssoil have been erroneous.

In dat.Ssoil, if sum==TRUE, only use Tsoils that are above the depth of Gsoil.

Cosmetic cleanup, converting = to <-.
Removing unneeded extra parameters to apply functions.
Merge branch 'master' of github.com:ncareol/R-eol