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

rdf = RadialDistribution(md_trajectory,
                         cutoff_radius=8.0*Angstrom,
                         pair_selection=[Aluminum, Oxygen])

# Get the bin_centers and the histogram of the radial distribution.
distances = rdf.distances().inUnitsOf(Angstrom)
histogram = rdf.data()

# Plot the data using pylab.
import pylab

pylab.plot(distances, histogram, label='Al-O RDF')
pylab.xlabel('r (Ang)')
pylab.ylabel('g(r)')
pylab.legend()

pylab.show()