Flowchart of our MOKA^3D kinematic model. The method consists in assigning the model parameters and assuming a coordinate system and a velocity field. Each cloud is assigned the coordinates in the source reference frame (green), which is later transformed to the observer reference frame (blue), through Euler rotation matrices. If necessary, the method requires the convolution with the PSF and/or the spectral LSF. Then, the model is binned in the space sampled by the observed data cube (red). The model emission line profile percentiles at 1% and 99% are compared to the observed ones and, if they are comparable, each particle is weighted according to the observed emission in the corresponding bin. The model is evaluated by means of the κ estimator (Eq. (6)). The free parameters are varied and the loop repeats, until finding the suitable parameter configuration that minimises κ and reproduces the observed spectrum in each spaxel.