Here is the comparison for Zeldovich, all cases use PPMP HLLC and VL, the Dual Energy parameters are:

Also, I found that Enzo can use the 4th order method to get the gradient of the potential, here I am using that method.

$\eta_1 = 0.001$ $\eta_2 = 0.010$

Now also using the truncated error condition like Ramses using $\beta=0.05$ and same $\eta_2$