fortran code