In this paper we consider a simultaneous system of spatially interrelated cross sectional equations. Our specification incorporates spatial lags in the endogenous and exogenous variables. In modeling the disturbance process we allow for both spatial correlation as well as correlation across equations. The data set is taken to be a single cross section of observations. The model may be viewed as an extension of the widely used single equation Cliff-Ord model. We suggest computationally simple limited and full information instrumental variable estimators for the parameters of the system and give formal large sample results.