VLW code 
VLW code includes two parts, the first called VLW A with which the detonation performance of all kinds of explosives can be computed accurately, the second called VLW B with which the combustion properties of rocket and gun propellants can be computed accurately too.
For the computations required the input data which
are elemental composition, heat of formation, and density.
VLW code includes an equilibrium subroutine that can solve the
equilibrium composition of a system of 10 elements, 20 gaseous species, and five
other phases subject to the limitation that only one phase may disappear.
The examples of calculations have been put in the 
form of calculation.html
It may be said that the VLW code is versatile.
Now it is on sale. If you are interested,
please mail me here.