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