Networks with as many mechanical constraints as degrees of freedom and no redundant constraints are minimally rigid or isostatic. Isostatic networks are relevant in the study of network glasses, soft matter, and sphere packings. Because of being at the verge of mechanical collapse, they have anomalous elastic and dynamical properties not found in the more commonly occurring hyperstatic networks. In particular, while hyperstatic networks are only slightly affected by geometric disorder, the elastic properties of isostatic networks are dramatically altered by it. In this paper, we show how disorder and system size strongly affect the ability of isostatic networks to sustain a compressive load. We develop an analytic method to calculate the bulk compressive modulus B for various boundary conditions as a function of disorder strength and system size. For simplicity, we consider square and cubic lattices with L^{d} sites, each having d mechanical degrees of freedom, and dL^{d} rotatable springs in the presence of hot-solid disorder of magnitude ε. Additionally, ∼L^{θ} sites may be fixed, thus introducing a nonextensive number of redundancies, either in the bulk or on the boundaries of the system. In all cases, B is analytically and numerically shown to decay as L^{-μ} with μ_{large}=d-θ for large disorder and μ_{small}=max{(d-θ-1),0} for small disorder. Furthermore B(L,ε)L^{μ_{small}}=g(λ) with λ=L^{(μ_{large}-μ_{small})}ε^{2} a scaling variable such that λ<<1 is small disorder and λ>1 is large disorder. The faster decay to zero of B in the large disorder regime results from a broad distribution of spring tensions, including tensions of both signs in equal proportions, which is remarkable since the system is under a purely compressive load. Notably, the bulk modulus is discontinuous at ε=0, a consequence of the fact that the regular network sits at an unstable degenerate configuration.