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.

Now also using the truncated error condition like Ramses using and same