Added script to do multiple runs over mu
Scaled output flux to 11.2 km because plotting routines expect 11.2
Added Lscale parameter to lcplot to simulate a small area near the polar cap
Looking at crust heating at different depths
Turned off perpendicular K components