We present a comprehensive near‐infrared study of two molecular bow shocks in two protostellar outflows, HH 99 in R Coronae Australis and VLA 1623A (HH 313) in Rho Ophiuchi. New, high‐resolution, narrow‐band images reveal the well‐defined bow shock morphologies of both sources. These are compared with two‐dimensional MHD modelling of molecular bows from which we infer flow inclination angles, shock speeds and the magnetic field in the pre‐shock gas in each system. With combined echelle spectroscopy and low‐resolution K‐band spectra we further examine the kinematics and excitation of each source. Bow shock models are used to interpret excitation (CDR) diagrams and estimate the extinction and, in the case of VLA 1623, the ortho–para ratio associated with the observed H2 population. For the first time, morphology, excitation and kinematics are fitted with a single bow shock model.
Specifically, we find that HH 99 is best fitted by a C‐type bow shock model (although a J‐type cap is probably responsible for the [Fe ii] emission). The bow is flowing away from the observer (at an angle to the line of sight of ∼45°) at a speed of roughly 100 km s−1. VLA 1623A is interpreted in terms of a C‐type bow moving towards the observer (at an angle to the line of sight of ∼75°) at a speed of ∼80 km s−1. The magnetic field associated with HH 99 is thought to be orientated parallel to the flow axis; in VLA 1623A the field is probably oblique to the flow axis, since this source is clearly asymmetric in our H2 images.