fortran90 code, compiles with gfortran