This is a follow-up to https://ask.sagemath.org/question/73771/ where the OP requested help at MSE forum and got a solution in the form of a Magma code.
I wonder if Sage has enough functionality to get the same things done, without calling Magma. Could someone experienced in algebraic geometry please take a look at that code and convert it to Sage if possible?
