  The new code for numerical simulation of magnetic hydrodynamical
astrophysical flows with consideration of chemical reactions is given in the
paper. At the heart of the code - the new original low-dissipation numerical
method based on a combination of operator splitting approach and
piecewise-parabolic method on the local stencil. The details of the numerical
method are described; the main tests and the scheme of parallel implementation
are shown. The chemodynamics of the hydrogen while the turbulent formation of
molecular clouds is modeled.
