pr0 = 0.5*(1 + np.sum(modelparams*expparams['axis'],1))