fn =data[i]
title =titles[i]
pro =fn.h.proj(axis,field,center=fn.domain_center)
proj =fn.h.proj(0,field,weight_field=field,center=fn.domain_center)
pro_frb =pro.to_frb((25./h, "Mpc"), 2000)
axs1[i] =plt.subplot(2,2,i+1)
Image=axs1[i].imshow(pro_frb[field],origin='lower',norm=LogNorm(),extent=(-12.5/h,12.5/h,-12.5/h,12.5/h))
vmin_0,vmax_0 = Image.get_clim()
***********************************************************************
I tried to find out something similar with the new yt.ProjectionPlot( ) without success, Could you please help me?
Kind Regards,