The fast azimuthal advection upon which FARGO is based is presented in the following publications:

 A&AS paper describing the algorithm "Fast Advection in Rotating Gaseous Objects"
 Conference proceedings which include a discussion of the timestep limitation
(it specifically shows that it is correct for the embedded protoplanet case).