Refactor IceSheet to use array parallelism

[mcgibbon]

Currently the IceSheet object operates on each column within a Python for loop. Many of the operations could be written using array parallelism, which would run a lot faster. This would be purely a performance enhancement.