initial_old fortran code