aim
function get_bath(this) result(bath_sites) class(aim) :: this integer :: bath_sites(this%n_bath) bath_sites = this%bath_sites end function get_bath