ADD: gradient caclulation

Adding Fortran routine for calculating the gradient of a scalar
function. This routine currently contains one subroutine capable of
calculating the density lapse rate from sounding data on a non-uniform
grid. The reason why np.gradient is not used is due to the non-uniform
nature of this problem.