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\)