An effective method to compute closure ordering for nilpotent orbits of θ-representations