The MHD simulations of stellar jets recently included complex models of radiative emission computation, allowing for better predictions in terms of emission line ratios. Employing also Adaptive Mesh Refinement, the large-scale propagation of jets could be followed. The simulation of multiple shockwaves originating in perturbations close to the jet origin and travelling along the jet beam allows for the construction of synthetic emission maps at various wavelengths, to be directly compared to observations. We apply this procedure for the jets originating from RW Aurigae.