lattice
subroutine set_ndim(this, n_dim) class(lattice) :: this integer, intent(in) :: n_dim this%n_dim = n_dim end subroutine set_ndim