Add calc_DThetaDsig
Add calcuration method for buoyancy freqency