Summary A totally conservative Eulerian 3D numerical code which conserves mass, momentum and energy both in the source and remap steps is developed. Mass, momentum and kinetic energy are conserved simultaneously during the remapping. The use of special form of linear viscosity makes the code more tolerant to the time step, leaving the second order of accuracy. Multimaterial flows including those contaminated by dust particles may be investigated using this program. The performance of the code is illustrated by modeling the Shoemaker-Levy 9 impact against Jupiter. Penetration of nonuniform fragments with complex structure into the Jovian atmosphere is investigated.