quad2fe_3d is now working properly for any mesh
XX call it properly
Use quad2fe_3d()
XX remove debug stuff