Add DEMderiv.R containing methods to compute surface attributes