To make my life easier, I wrote python wrappers to improve the usability of old FORTRAN finite element code.
You can see FEM1D here and FEM2D here.