We describe a new hybrid N -body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equati...