def interface_temperature(Ts, Tatm, **kwargs):
<EOL>lev = Tatm.domain.axes['<STR_LIT>'].points<EOL>lev_bounds = Tatm.domain.axes['<STR_LIT>'].bounds<EOL>f = interp1d(lev, Tatm, axis=-<NUM_LIT:1>)  <EOL>Tinterp = f(lev_bounds[<NUM_LIT:1>:-<NUM_LIT:1>])<EOL>Ttoa = Tatm[...,<NUM_LIT:0>]<EOL>Tinterp = np.concatenate((Ttoa[..., np.newaxis], Tinterp, Ts), axis=-<NUM_LIT:1>)<EOL>return Tinterp<EOL>
Compute temperature at model layer interfaces.