Hi, you can use for instance the approach posted by the OP here - Is this a correct way to calculate one-site and two-site entanglement spectrum and entropy?. Then from the eigenvals. of the reduced density matrix you can easily get the entanglement entropy and other such quantities.
This works very well as long as the subsystem isn’t very large (a few sites or so, and also depends on the number of indices involved).