Made all exponents in pow functions in surface budget and diffusion
integer to increase performance
Add copyright statement to JPM for heptadiagonal matrix solver
Fixed out-of-bounds read in Brunt vaisala calculation