starting work on discrete energy distributions

completely untested: DTD looks risky
fortran code and all relevant src updated

compiles fine but untested
Merge branch 'bugfixes' into discrete-distributions
Work looks finished

Everything seems to work fine
Don't go wild with the scale parameters though!