Observations of Globular Cluster ellipticity distributions related to some fundamental parameters give strong evidence for a decay of rotational energy in these systems with time. In order to study the effectiveness of angular momentum transport (or loss, resp.) a code has been written which solves the Fokker-Planck equation in (E, Jz)-space and follows the evolution from some initial conditions through core collapse (and possibly gravothermal oscillations) up to the post-collapse phase. For the purpose of comparability with N-body simulations rotating initial model configurations according to the prescriptions of Lupton & Gunn (1987) have been constructed. These models are intended to continue previous work by Goodman (1983, Fokker-Planck) and Akiyama & Sugimoto (1989, N-Body). In this contribution the derivation of the flux coefficients is given.