In this work, we provide a viscous solution of the wave field generated by librating a disk (harmonic oscillation of the rotation rate) in a stably stratified rotating fluid. The zonal flow (mean flow correction) generated by the nonlinear interaction of the wave field is also calculated in the weakly nonlinear framework. We focus on the low dissipative limit relevant for geophysical applications and for which the wave field and the zonal flow exhibit generic features (Ekman scaling, universal structures, etc.). General expressions are obtained which depend on the disk radius a * , the libration frequency ω * , the rotation rate Ω * of the frame, the buoyancy frequency N * of the fluid, its kinematic diffusion ν * and its thermal diffusivity κ * . When the libration frequency is in the inertia-gravity frequency interval (min(Ω * , N * ) < ω * < max(Ω * , N * )), the presence of conical internal shear layers is observed in which the spatial structures of the harmonic response and of the mean flow correction are provided. At the point of focus of these internal shear layers on the rotation axis, the largest amplitudes are obtained: the angular velocity of the harmonic response and the mean flow correction are found to be O(εE −1/3 ) and (ε 2 E −2/3 ) respectively, where ε is the libration amplitude and E = ν * /(Ω * a * 2 ) is the Ekman number. We show that the solution in the internal shear layers and in the focus region is at leading order the same as that generated by an oscillating source of axial flow localized at the edge of the disk (oscillating Dirac ring source).