Switch to plyr + doMC and sparse matrices

1. If memory is an issue (see #2), then sparse matrices should be better. I
tested with a 100x100 matrix with 100 ones.

2. I replaced snow/snowfall by doMC and plyr. This will make it easier
to run the analysis on clusters. There is an example at the end.