We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from /fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses , a common mass for the first-and second-generation squarks and a distinct third-generation squark mass , a common mass for the first-and second-generation sleptons and a distinct third-generation slepton mass , a common trilinear mixing parameter A, the Higgs mixing parameter , the pseudoscalar Higgs mass and . In the fit including , a Bino-like is preferred, whereas a Higgsino-like is mildly favoured when the constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, , into the range indicated by cosmological data. In the fit including , coannihilations with and the Wino-like or with nearly-degenerate first- and second-generation sleptons are active, whereas coannihilations with the and the Higgsino-like or with first- and second-generation squarks may be important when the constraint is dropped. In the two cases, we present functions in two-dimensional mass planes as well as their one-dimensional profile projections and best-fit spectra. Prospects remain for discovering strongly-interacting sparticles at the LHC, in both the scenarios with and without the constraint, as well as for discovering electroweakly-interacting sparticles at a future linear collider such as the ILC or CLIC.