Create custom output file
Hello all,
I have created a function in Gadget2 that calculates a velocity increment (dv), and I add this increment on the file timestep.c every 100 timesteps.
However, I'd like to know if my function is actually doing what I want, so I want to output a file that writes "dv1, dv2, dv3".
What I did was to add this immediately after I calculate dv on timestep.c, on the advance_and_find_timesteps function:
if (ThisTask ==0)
{
fprintf( t_out, "%e, %e, %e, %e, %e, %e\n", P[i].Pos[0], P[i].Pos[1], P[i].Pos[2], dv1, dv2, dv3 ) ;
}
The file is created and the values are written, but it only writes 1400 values (out of 20,000 particles). I am not sure how to fix this so that all the 20,000 values are written. If I remove the if statement, the code crashes, saying "Segmentation fault: 11. Signal code: Address not mapped"
If anyone has any suggestions, please let me know.
Thank you very much.
-Jesus
Received on 2017-02-25 03:48:32
This archive was generated by hypermail 2.3.0
: 2023-01-10 10:01:32 CET