changed exponent in ra calculation from 0.945 to 0.495 following Greek et al. 1989. Changed denominator order of operations in Te.calc and added epsilon, following Greek et al. 1989