Comment by basementcat

Note this is free form FORTRAN 90 code; not your grandparent's f77 optimized for card punch entry.