Also change integer(4) to integer(i4). This makes the code more portable
and adhering to the Fortran standard.