Bulk calculation

[Muhsin]

Hi all,

I'm trying to calculate the bulk temperature and bulk velocity to the axial direction. I ran code saturne in parallel with 4 processor. I found that my calculation is stopped after 'call sum of parameters...' is written and suppose the problem occured when it's calling parsom (as you could see in the subroutine and listing attached).

Code: Select all

if (irangp.ge.0) then
write (nfecra,*) 'call sum of parameters...'
call parsom(velave)
write (nfecra,*) 'velave called...'
call parsom(tempave)
write (nfecra,*) 'tempave called...'
call parsom(totvol)
write (nfecra,*) 'totvol called...'
endif

I wonder what exactly is wrong in the script. Would anyone care to help and show me the problem caused. Thank you

[Yvan Fournier]

Hello,

There is one error in your code (which I did not see immediately):

parsom needs to be called on all ranks, but you call it inside a test on (nlelt.gt.0), which depends on the mesh and partition.

Try simply removing the test on "ilelt.gt.0" (always run the code inside): when ilelt is equal to 0; the loop on ilelt will not execute, so the code should work fine (i.e. the test on "ilelt.gt.0" is redundant and breaks parallelism).

Regards,

Yvan