Modeling the Subsurface Structure of Sunspots