code works in fortran 95 with GNU