C calling Fortran routines works for full finite element solution