md_trajectory = nlread('alumina_trajectory.nc')[-1]

t_profile = TemperatureProfile(md_trajectory)

# Get the bin_centers and the temperature values.
z            = t_profile.zValues().inUnitsOf(Angstrom)
temperatures = t_profile.data().inUnitsOf(Kelvin)

# Plot the data using pylab.
import pylab

pylab.plot(z, temperatures, label='Temperature profile')
pylab.xlabel('z (Ang)')
pylab.ylabel('T(z) (K)')
pylab.legend()

pylab.show()