Added reference fortran code