You can solve the problem like this:
using Compat, Compat.Random
const twisters = [MersenneTwister() for i ∈ 1:Threads.nthreads()];
#Multithreaded.
function tf!(x::Vector{Matrix{Float64}},N::Int64)
Threads.@threads for ii=1:N
id = Threads.threadid()
twister = twisters[id]
@inbounds x_thd = x[id]
for nn=1:100
for mm=1:100
@inbounds x_thd[mm,nn] += randn(twister)
end
end
end
return nothing
end