I did some simple volume rendering with the following script:

    volume2 = AMRKDTree(pf, fields=["Dark_Matter_Density"],
                no_ghost=False, tree_type="domain",
                le=c-0.5*WW, re=c+0.5*WW)
    cam = pf.h.camera(c, L, W, N, tf, volume=volume2, no_ghost=False,
                north_vector=L, steady_north=True)
    cam.snapshot(fn="%s_iso-DMdensity-%3.3d.png" % (filenameTHIS, j))

I got rather strange results in that the pictures look symmetric, which I am pretty
sure can not be true.
I attach the obtained plot.
Note that I am using KD tree and using 32 cores. 

Your help at your earliest convenience is appreciated.