In this paper we investigate several classes of harmonic functions with varying argument of coefficients which are defined by means of the principle of subordination between harmonic functions. Properties such as the coefficient estimates, distortion theorems, convolution properties, radii of convexity, starlikeness and the closure properties of these classes under the generalized Bernardi-Libera-Livingston integral operators are investigated