The concept of 1D subsonic COIL model with a mixing length was generalized to include the influence of a variable magnetic field on the stimulated emission cross-section. Equations describing the chemical kinetics were solved taking into account together with the gas temperature also a simplified mixing model of oxygen and iodine molecules. With the external time variable magnetic field the model is no longer stationary. A transformation in the system moving with the mixture reduces partial differential equations to ordinary equations in time with initial conditions given either by the stationary flow at the moment when the magnetic field is switched on combined with the boundary conditions at the injector. Advantage of this procedure is a possibility to consider an arbitrary temporal dependence of the imposed magnetic field and to calculate directly the response of the laser output. The method was applied to model the experimental data measured with the subsonic version of the COIL device in the Institute of Physics, Prague, where the applied magnetic field had a saw-tooth dependence. We found that various values characterizing the laser performance, such as the power density distribution over the active zone cross-section, may have a fairly complicated structure given by combined effects of the delayed reaction to the magnetic switching and the flow velocity. This is necessarily translated in a time dependent spatial inhomogeneity of output beam intensity profile.