  We present estimates of Single Spin Asymmetry (SSA) in the electroproduction
of $J/\psi$ taking into account the transverse momentum dependent (TMD)
evolution of the gluon Sivers function and using Color Evaporation Model of
charmonium production. We estimate SSA for JLab, HERMES, COMPASS and eRHIC
energies using recent parameters for the quark Sivers functions which are
fitted using an evolution kernel in which the perturbative part is resummed up
to next-to-leading logarithms (NLL) accuracy. We find that these SSAs are much
smaller as compared to our first estimates obtained using DGLAP evolution but
are comparable to our estimates obtained using TMD evolution where we had used
approximate analytical solution of the TMD evolution equation for the purpose.
