How to handle the list in an openmp parallelism

gonski

I met across a big problem in parallelizing my code. In my code, it needs to create a list to consider interactions between particles, see the following pseudo code. To parallelize the part of code, How can I avoid race conditions for k, f1,f2 in multi-thread calculations in the parallel loop?

Cheers, Gonski

Pseudo code to create neighbor list for particles

k=1

do i=1, n

Note: here many lines have been snipped

f1(k)=c1

if(something )then

Note: here many lines have been snipped

f2(k)=c2

endif

Note: here many lines have been snipped

if(something)then

k=k+1

endif

enddo