lattice
integer elemental function get_ndim(this) class(lattice), intent(in) :: this get_ndim = this%n_dim end function get_ndim