The orientation of a fluid of anisotropic particles at both solid (impenetrable) and vapour-like (penetrable) interfaces was studied using classical molecular dynamics (MD) simulations. A range of uniaxial oblate Gay-Berne (GB) ellipsoids, representative of typical small OSC molecules, were studied to elucidate the general effects of shape and interaction anisotropy. These simple models greatly increase computational efficiency compared with all-atom simulations, allowing a wide variety of parameters and conditions to be examined. Simulations were carried out using LAMMPS (version 3 March 2020), and visualisation and analysis of simulation trajectories were conducted using OVITO. The GB potential is an anisotropic form of the widely used Lennard-Jones (LJ) pair potential, and captures the short-ranged excluded-volume repulsion and longer ranged van der Waals attraction between uncharged anisotropic molecules. The potential is characterised by parameters σ and ε that define the length and energy scales of the potential, respectively, parameters ν and µ that tune the shape of the potential, and, for uniaxial particles, the short and long principal diameters, σ F and σ S , and dimensionless relative well depths, ε F and ε S , along the corresponding axes (the subscripts "F" and "S" denote "face" and "side", respectively). The anisotropy of the uniaxial particles is described by the shape anisotropy parameter κ = σ F /σ S and the interaction anisotropy parameter κ = ε S /ε F . A full description of the GB potential and parameters is given in the ESI, section S1.