Dear All,

I want to use two functions which are basically distribution functions in momentum space. I want to compute the normalization constant for those functions and after normalization, I will get the intersection point of those distribution functions, say p_inter_point.

But now the issue is that the distribution functions will depend on the temperature, density, and some other derived fields of their respective grid. And for each grid, I have to compute the normalization and intersection point (p_inter_point) for the two distributions.

I used the yt-data enzo_tiny_cosmology and this is the sample code I run for example.

Please have a look and suggest me the right way to achieve the goal.

Many thanks and regards

-Prateek Gupta