We present a new, versatile numerical code PLUTO, suitable for the solution of high-mach number flow in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multiphysics (newtonian or relativistic HD or MHD), multialgorithm modular environment particularly oriented toward the treatment of astrophysical flows in th presence of discontinuities. The modular structure exploits a general framework for integrating a system of conservation laws, built on modern Godunov-type shock-capturing schemes. Recently implemented features include Adaptive Mesh Refinement, radiative cooling, resistivity and heat conduction.