Triaxial mass models producing ellipticity variations and isophote twists are investigated. The spherical mass distributions are made triaxial by adding two more terms, each consisting of a radial function muliplied by a ...