[Argen]

I am using parallel computing for OpenFOAM and some FORTRAN subroutines are called by OF during computation. But, it's found that the errors occur when any Fortran subroutine (say my_test(a)) is called where a is defined as "double a[nCells]". The interesting is that it's no problem for calling it if a single processor is used. Any idea for this problem?

[Argen]

any suggestion to sort it out?

[wyldckat]

Greetings Argen,

OK, the reason why it works in serial mode, but not in parallel, is very likely because "nCells" refers to the whole mesh and not the local processor's mesh only.

Uhm... I know there are applications in OpenFOAM that have similar needs, for both serial and parallel execution...
Sorry, I can't find any in a short period of time. You'll have to check the existing "applications" folder and see directly what code already exists that can do this.

The "applications/test" folder has many times true pearls of wisdom!
Here's an example that might be a good hint for what you're looking for - "applications/test/router/Test-processorRouter.C":

// Collect centres of individual meshes (for visualization only)
Gather<point> meshCentres(meshCentre(mesh));

The "Gather" template provides a way for gathering information from all processors. Nonetheless, this would imply that whenever you need to use the FORTRAN code, only the master processor could do the calculation...

Best regards and good luck!
Bruno