Skip to content
Snippets Groups Projects
Commit 5162d638 authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

fixes typo, adds explanation

parent dd135920
No related branches found
No related tags found
1 merge request!35Feature/convergence test
Pipeline #120903 passed
......@@ -72,7 +72,7 @@ The full function is given below:
.. code:: python
def generate_initial_conditions():
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......@@ -121,7 +121,7 @@ NSVE solvers):
.. code:: python
def run_simulations_field(err_vs_t = False):
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......
......@@ -63,7 +63,7 @@ run:
.. code:: python
def generate_initial_conditions():
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......@@ -110,11 +110,26 @@ To put it simply, for the simulation "nsvep_base" the dataset
However, the "tracers0/rhs/8" dataset for simulation "nsvep_1x" contains
meaningful values (that **are** used because we request that `iter0` is
different from 0).
In order to generate consistent values of the "rhs" datasets, we first
run a simulation at the highest resolution (smallest time-step), and
then we extract the appropriate snapshots --- see also the second part
of the trajectory plot function.
Finally, the `launch` method is called in an unusual way for TurTLE.
Firstly, we demand a DNS that starts at an iteration different from 0.
Secondly, we manually copy the forcing parameters from the source
simulation; since the source cannot be specified explicitly, the `DNS`
class cannot search for the source parameters itself.
Thirdly, I will note that one may in fact break this call by modifying
the `iter0` and `checkpoints_per_file` values, since we have explicitly
placed `iter0` in the `checkpoint` 0 file.
However, a full discussion of the `iter0` and `checkpoint` subtleties is
beyond the scope of this tutorial.
.. code:: python
def run_simulations_particles():
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......@@ -147,7 +162,8 @@ different from 0).
'--tracers0_neighbours', '{0}'.format(neighbours),
'--tracers0_smoothness', '{0}'.format(smoothness),
'--tracers0_integration_steps', '4',
'--cpp_random_particles', '{0}'.format(particle_random_seed)])
'--cpp_random_particles', '{0}'.format(particle_random_seed),
'--checkpoints_per_file', '{0}'.format(16)])
# 2. Prepare initial conditions
for neighbours, smoothness in neighbours_smoothness_list:
......@@ -216,7 +232,8 @@ different from 0).
'--tracers0_neighbours', '{0}'.format(neighbours),
'--tracers0_smoothness', '{0}'.format(smoothness),
'--tracers0_integration_steps', '4',
'--cpp_random_particles', '0'], # turn off C++ particle initialization
'--cpp_random_particles', '0', # turn off C++ particle initialization
'--checkpoints_per_file', '{0}'.format(16)],
iter0 = 8*factor)
source_file.close()
return None
......
......@@ -21,7 +21,7 @@ def main():
return None
def generate_initial_conditions():
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......@@ -65,7 +65,7 @@ def generate_initial_conditions():
return None
def run_simulations_field(err_vs_t = False):
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......
......@@ -26,7 +26,7 @@ def main():
return None
def generate_initial_conditions():
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......@@ -55,7 +55,7 @@ def generate_initial_conditions():
return None
def run_simulations_particles():
# change these two values as neded.
# change these two values as needed.
# number of MPI processes to use
nprocesses = 8
# number of OpenMP threads per MPI process to use
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment