Taylor couette flow

I m trying to simulate Taylor -couette flow in fortran.I consider the flow in the annular space between the cylinder as an axisymmetric flow. I read some articles but I can not understand how the azimuthal speed is they calculated.
Help me please if you have a simple code example on this subject.