def _calc_taub(w, aod700, p):
p0 = <NUM_LIT><EOL>tb1 = <NUM_LIT> + <NUM_LIT>*np.log(w) + <NUM_LIT>*np.log(w)**<NUM_LIT:2><EOL>tb0 = <NUM_LIT> + <NUM_LIT>*np.log(w) + <NUM_LIT>*np.log(w)**<NUM_LIT:2><EOL>tbp = <NUM_LIT>*w + <NUM_LIT><EOL>taub = tb1*aod700 + tb0 + tbp*np.log(p/p0)<EOL>return taub<EOL>
Calculate the taub coefficient