Examine performance of loops over soil layers

As @swensosc pointed out today, a fair amount of time is spent in loops over soil layers. It would probably be worth looking at some of these more closely to see if we can identify some general ways the performance of these multi-level loops can be improved. For example, can we improve the cache efficiency or vectorizability? Can we move some code (conditionals, etc.) to outside these loops?