Add latitude dependence to the Coriolis parameter

[timspainNERSC]

Currently, the Coriolis parameter is the fixed value of 2ω_rot with no latitude dependence. Due to the cos θ dependence of the Coriolis force, this is pretty accurate in the polar regions. But latitude dependence should be included in the updateMomentum() functions of the dynamics classes.