Generalize the direct inversion algorithm to handle data of arbitrary lon/lat resolution

[csyhuang]

The additional direct inversion algorithm assumes 1°-resolution data (since we were in a hurry submitting Neal et al with the code). Generalize the code to handle data of arbitrary latitudinal resolution (uniform grid).

See pull request from Chris for details #45