Modeling short-range substitution order and disorder in crystals: Application to the Ga/Si distribution in a natrolite zeolite.