added fortran code