[yt-dev] inline profiles

[Nathan Goldbaum]

Hey Stephen,

I think you have two options and in principle both routes should give you the same answer.  The first option is to do as you suggest and use yt to do the cloud in cell interpolation of the dark matter particles onto the grid.  The second option would be to to make sure GravitatingMassFieldParticles in enzo is valid by calling the appropriate enzo grid member functions.  Enzo has a ton of machinery to make sure the CIC operation happens as fast as possible using the non-blocking communication machinery, however, that means you will have a lot of programming overhead to get enzo to produce the field in the middle of a timestep.

The 'best' answer probably depends on your needs, the size of the simulation, etc.

Cheers,

Nathan

On Aug 6, 2012, at 2:56 PM, Stephen Skory wrote:

> Hi all,
> 
> this is about inline Enzo+yt, but since I think the solution might be
> on the yt end, I'm sending it here.
> 
> I am attempting to profile halos located inline, and I am running into
> trouble when I attempt to profile over the 'TotalMassMsun' field. It
> is dying when it tries to access the 'Dark_Matter_Density' field,
> which I think it expects to be available, but isn't because Enzo only
> produces that when it writes to disk (AFAIK?). If this sounds correct
> to those in the know, I'd like to ask if anyone has a good idea how to
> get around this? Should I make a derived field (perhaps over-writing
> the built-in 'Dark_Matter_Density') using CICDeposit? I suppose I
> could try this solution first without emailing y'all, but I was
> thinking that if I/we find a working solution, it could be added to yt
> so the next person doesn't have to figure it out, so I'd like to
> minimize the false starts of poor/incorrect implementations.
> 
> Thanks for your comments.
> 
> -- 
> Stephen Skory
> s at skory.us
> http://stephenskory.com/
> 510.621.3687 (google voice)
> _______________________________________________
> yt-dev mailing list
> yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org