From 4b7d1ed94ccad1ad5a167ad21df1030f3d2f8cca Mon Sep 17 00:00:00 2001
From: Cristian C Lalescu <Cristian.Lalescu@ds.mpg.de>
Date: Sun, 10 Jan 2016 11:26:18 +0100
Subject: [PATCH] bugfix --- code now works for multiple interpolators
---
bfps/NavierStokes.py | 1 +
1 file changed, 1 insertion(+)
diff --git a/bfps/NavierStokes.py b/bfps/NavierStokes.py
index 70a4426d..901cbb05 100644
--- a/bfps/NavierStokes.py
+++ b/bfps/NavierStokes.py
@@ -551,6 +551,7 @@ class NavierStokes(bfps.fluid_base.fluid_particle_base):
'if (ps0->iteration % niter_part == 0)\n' +
'{\n')
for s in range(nspecies):
+ neighbours = self.parameters[interpolator[s] + '_neighbours']
self.particle_start += 'sprintf(fname, "tracers{0}");\n'.format(s0 + s)
self.particle_end += ('ps{0}->write(stat_file);\n' +
'delete ps{0};\n').format(s0 + s)
--
GitLab