Subroutine example fortran 90