hello everyone.
Based on the Discrete element model(DEM), I am trying to output the total force of each particle (FC(:, in MFIX code) to see the force direction on particles, just like the particle velocity.
And I want to output it by using the array ORIENTATION, so what i do is replacing the array ORIENTATION by FC in file vtp_mod.f, as is shown as below.
while the problem is the output data of FC is zero.
So how can i modify the code and realize this operation?