this is before I tested vs the Fortran code