A numerical model is constructed using the PLASIMO toolkit to simulate a microwave configuration which is similar to that which is in use for optical glass fibre production. The simulations offer the flow patterns and the electromagnetic (EM) energy incoupling of a two-temperature argon plasma. The Yee algorithm was used for the EM module, whereas the Semi implicit method for pressure linked equations algorithm was used to calculate the pressure and velocity field. It is found that at 400 W coupling of 2.46 GHz EM radiation in 1000 Pa argon results in a plasma that is not in local thermodynamic equilibrium, in the sense that Te/Th ≥ 4 and the ionization degree is smaller than that predicted by the Saha equations. The model results are subjected to various sanity checks.