FORTRAN Module and Program

module ejemplo_modulo
    implicit none
    
    contains
    
    subroutine distancia (a, b)
        real, dimension(2) :: a, b
        real :: d
        real, dimension(2) :: aux
        
        aux = (b-a)**2
        d = SQRT(aux(1) + aux(2))
        
        write (*,*) d
        
    end subroutine distancia
end module ejemplo_modulo
program ejemplo
    use ejemplo_modulo
    
      implicit none
      real, dimension(2) :: a, b

      write(*,*) 'a, b= : '
      read(*,*) a,b

      call distancia(a, b)
        
end program ejemplo