We present an algorithm to compute the primary decomposition of a submodule N of the free module ℤ[x1,...,xn]m. For this purpose we use algorithms for primary decomposition of ideals in the polynomial ring over the integers. The idea is to compute first the minimal associated primes of N, i.e. the minimal associated primes of the ideal Ann (ℤ[x1,...,xn]m/N) in ℤ[x1,...,xn] and then compute the primary components using pseudo-primary decomposition and extraction, following the ideas of Shimoyama-Yokoyama. The algorithms are implemented in Singular.