From cf1d40f4ae22d8b2cf66e83eb3ed7c6ec75f2787 Mon Sep 17 00:00:00 2001
From: Cristian C Lalescu <Cristian.Lalescu@ds.mpg.de>
Date: Mon, 30 Apr 2018 17:40:09 +0200
Subject: [PATCH] force divfree on test vorticity
---
bfps/TEST.py | 24 +++++++++++------------
bfps/cpp/full_code/test_interpolation.cpp | 1 +
2 files changed, 13 insertions(+), 12 deletions(-)
diff --git a/bfps/TEST.py b/bfps/TEST.py
index 96090076..d9b3dedd 100644
--- a/bfps/TEST.py
+++ b/bfps/TEST.py
@@ -349,18 +349,18 @@ class TEST(_code):
write_to_file = False,
spectra_slope = 1.0,
amplitude = 0.05)
- data[:] = 0.0
- # ABC
- data[0, 0, 1, 1] = -0.5*(1j)
- data[0, 0, 1, 2] = 0.5*(1j)
- data[0, 1, 0, 0] = -0.5*(1j)
- data[0, self.parameters['nz'] - 1, 0, 0] = 0.5*(1j)
- data[0, 1, 0, 1] = 0.5
- data[0, self.parameters['nz'] - 1, 0, 1] = 0.5
- data[1, 0, 0, 0] = 0.5
- data[self.parameters['ny'] - 1, 0, 0, 0] = 0.5
- data[1, 0, 0, 2] = -0.5*(1j)
- data[self.parameters['ny'] - 1, 0, 0, 2] = 0.5*(1j)
+ #data[:] = 0.0
+ ## ABC
+ #data[0, 0, 1, 1] = -0.5*(1j)
+ #data[0, 0, 1, 2] = 0.5*(1j)
+ #data[0, 1, 0, 0] = -0.5*(1j)
+ #data[0, self.parameters['nz'] - 1, 0, 0] = 0.5*(1j)
+ #data[0, 1, 0, 1] = 0.5
+ #data[0, self.parameters['nz'] - 1, 0, 1] = 0.5
+ #data[1, 0, 0, 0] = 0.5
+ #data[self.parameters['ny'] - 1, 0, 0, 0] = 0.5
+ #data[1, 0, 0, 2] = -0.5*(1j)
+ #data[self.parameters['ny'] - 1, 0, 0, 2] = 0.5*(1j)
ofile['vorticity/complex/{0}'.format(0)] = data
with h5py.File(os.path.join(self.work_dir, self.simname + '_output.h5'), 'a') as ofile:
ofile.require_group('tracers0')
diff --git a/bfps/cpp/full_code/test_interpolation.cpp b/bfps/cpp/full_code/test_interpolation.cpp
index 4cfdca66..e702e708 100644
--- a/bfps/cpp/full_code/test_interpolation.cpp
+++ b/bfps/cpp/full_code/test_interpolation.cpp
@@ -108,6 +108,7 @@ int test_interpolation<rnumber>::do_work()
this->simname + std::string("_input.h5"),
"vorticity",
0, true);
+ this->kk->template force_divfree<rnumber>(this->vorticity->get_cdata());
// compute velocity
invert_curl(this->kk, this->vorticity, this->velocity);
--
GitLab