<div dir="ltr">Dear YT users,<div><br></div><div>I'm trying to make a density projection, in the following lines I show you the lines that I have.</div><div><br></div><div>#Ramses output</div><div><div>ds = yt.load(path+"output_%05d"%snap+"/info_%05d.txt" %snap)</div><div>disk = ds.disk('max',[0,0,1],(20,'pc'),(5,'pc'))</div><div>proj=ds.proj('density',axis='z',weight_field='density',data_source=disk,method='integrate')</div></div><div><br></div><div>However when I print proj["density"] I get a list of nan, however when I use the whole domain as my data_source I get reasonable numbers.</div><div><br></div><div>Thanks,</div><div>Matías Suazo</div><div>MSc student at Universidad de Chile</div></div>