There are also some FDM and spectral-based solvers in Python:
Also several powerful Python FEM frameworks, such as FEniCs and FireDrake:
The UFM language these are based on isn’t too far from Fortran.
One more recent one, again using a Python code generator that outputs C++ (or even SIMD-heavy assembly):
The authors went all the way to 10^12 degrees of freedom.