Theoretical modeling of coexisting tetragonal and rhombohedral heterophase polydomain structures in lead zirconate titanate ferroelectric films near the morphotropic phase boundary