Classic single-grid Fortran routines