Issue with dust.py area

[ledm]

There's an issue with area calculation and masking in dust.py.
Needs a bit of work.

This issue is an extension of #144