This calculation takes into account corrections for real gases. It uses the same polynomial fit as Subsurface. Note that this calculation assumes room temperature conditions. Fill slowly and wait for it to cool down.
This script is by Robert C. Helling (helling@atdotde.de). Find its source on Github and more information at TheTheoreticalDiver.