if (ts < tns + 1):
                                        # exploiting the symmetry
                                        pdf[i, 1:Nx1, ts] += fxind[IJ:J] * dt
                                        err[i, 1:Nx1, ts] += err0[IJ:J] * dt2
                                        terr[i, 1:Nx1, ts] += terr0[IJ:J] * dt