Fortran code