Thank you for this helpful package! it helped me a lot.
Is there is a possibility to add to your package, or alternatively use your results and write a different code to calculate the Kretschmann scalar which is given by:
K = R_{abcd}R^{abcd}
Where R_{abcd} is the Riemann tensor which obtained by:
R_{abcd}=g_{ce}R_{abd}^{e}