if len(modelparams.shape) == 1:
            modelparams = modelparams[:, np.newaxis]