GRIDobj.toposhielding computes topographic shielding factors for cosmogenic nuclide production rates.
This composes several functions including gradient8 (#65), aspect (#102) and castshadow (#117) and applies a kind of map-reduce pattern, which might be more efficiently computed in libtopotoolbox.
