I am using the following code to try to sum the volume of a specific domain, but it returns me the sum of all domains. What is wrong with it? Thanks for you help =D
void volume (int ID)
Domain *domain = Get_Domain(1);
Thread *t = Lookup_Thread(domain,ID);
Message0("Volume integral: %f\n", VOLMOTC);
as it name implies thread_loop_c(t,domain) loops over all cell threads in your domain and thus neglects the previous specification with ID. Omitting the outer loop will do what you want.
|All times are GMT -4. The time now is 05:53.|