Analysis of the Fortran routines implementing the Arithmetic Mean method for the solution of block tridiagonal linear systems