MOD: density gradient calculation

Computing the density lapse rate (rate of change with altitude) can be
computed on a non-uniform grid using "basic", low-order centered,
forward, and backward finite difference methods. For the time being I
am leaving the high-order methods as a to-do.