In this paper we provide a framework for constructing general complex geometrical optics solutions for several systems of two variables that can be reduced to a system with the Laplacian as the leading order term. We apply these special solutions to the problem of reconstructing inclusions inside a domain filled with known conductivity from local boundary measurements. Computational results demonstrate the versatility of these solutions to determine electrical inclusions.