Elliptic pfaffians and solvable lattice models

We introduce and study twelve multivariable theta functions defined by pfaffians with elliptic function entries. We show that, when the crossing parameter is a cubic root of unity, the domain wall partition function for the eight-vertex-solid-on-solid model can be written as a sum of two of these pfaffians. As a limit case, we express the domain wall partition function for the three-colour model as a sum of two Hankel determinants. We also show that certain solutions of the TQ-equation for the supersymmetric eight-vertex model can be expressed in terms of elliptic pfaffians.