q = np.tile(q, (expparams.shape[0], 1, 1, 1)).transpose((1,2,3,0))