# Read self-consistent calculations
molecule_configuration_list = nlread(
    'nh3_scf.nc',
    MoleculeConfiguration
)

# Choose the first configuration in the list
molecule_configuration = molecule_configuration_list[0]