The shear viscosity of the quark-gluon plasma (QGP) plays a crucial role in interpreting current measurements from heavy-ion collisions and is a key input to hydro-dynamical models. The interest in shear viscosity also lies in the fact that QGP is the most ideal fluid ever observed and has the shear viscosity to entropy ratio ($\eta / s$) close to the theoretical bound $\eta / s \geq 1/ 4 \pi$ in the strong coupling region within AdS/CFT formalism. The lattice determination of $\eta / s$ has been explored for the pure gauge case, but its determination in full QCD remains unexplored, despite its significant importance. In this proceeding, we present updates on extending our quenched findings to full QCD. Specifically, we focus on the renormalization of the energy-momentum tensor with the gradient flow method and provide a progress update on determining the relevant renormalization coefficients for shear viscosity. For this purpose, we have used an imaginary isospin chemical potential.