Hi everyone
In calc_force_dem.f file, the interaction between particle LL and its neighbour particle I is modeled. However, ‘force_coh’ is only assigned to particle LL through postcohesive(LL). Why ‘force_coh’ is not assigned to particle I at the same time, i.e.
PostCohesive(LL) = PostCohesive(LL) + FORCE_COH / PMASS(LL)
PostCohesive(I) = PostCohesive(I) + FORCE_COH / PMASS(I)
Particle LL and I may not be combined repeatedly during neighbour searching.