New solvable two-matrix model and 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 out explicitly in terms of series in strict partitions. The related string equations are presented.