Circulation of cerebrospinal fluid (CSF) plays an important role during development. In zebrafish embryo, the flow of CSF has been found to be bidirectional in the central canal of the spinal cord. In order to compare conditions and genetic mutants across each other, we recently automated the quantification of the velocity profile of exogenous fluorescent particles in the CSF. We demonstrated that the beating of motile and tilted cilia localized on the ventral side of the central canal was sufficient to generate locally such bidirectionality. Our approach can easily be extended to characterize CSF flow in various genetic mutants. We provide here a detailed protocol and a user interface program to quantify CSF dynamics. In order to interpret potential changes in CSF flow profiles, we provide additional tools to measure the central canal diameter, characterize cilia dynamics and compare experimental data with our theoretical model in order to estimate the impact of cilia in generating a volume force in the central canal. Our approach can also be of use for measuring particle velocity in vivo and modeling flow in diverse biological solutions.