Now parallelized in both theta and zeta using 2D DMDA