Some Fortran code