We study the statistical mechanics of the KDP, F and modified KDP models defined on a 4 x4 lattice. Direct fields that are parallel and perpendicular to the crystal axis as well as staggered field are considered. Numerical results on the energy, specific heat, polarization and polarizability are given.