A new solvable two-matrix model and the BKP tau function

We present exactly solvable modifications of the two-matrix Zinn-Justin–Zuber model and write it as a tau function. The grand partition function of these matrix integrals is written as the fermion expectation value. The perturbation theory series is written explicitly in terms of a series in strict partitions. The related string equations are presented.