  We present a new modeling tool for planetary nebulae (PNe), based on 3D
photoionization calculations. From models for two theoretical PNe, we show that
the enhancement in the equatorial zone observed in several PNe is not
necessarily due to a density gradient, as usually interpreted. We study the
morpho-kinematic properties of the PN NGC 3132 and show that a bipolar Diabolo
shape successfully reproduces the observed images, as well as the high
resolution observations.
