Correct the sign of Shear Production
Move helper functions into subdirectory
calculate diff between model and observations
move palete and ticbtic setting into ..limits.plt