diff --git a/doc/examples/time_stepping/3d/submit_1simul.py b/doc/examples/time_stepping/3d/submit_1simul.py
index a74f6a2eeda8b09c5f46c087eac7e5bcc68804ef_ZG9jL2V4YW1wbGVzL3RpbWVfc3RlcHBpbmcvM2Qvc3VibWl0XzFzaW11bC5weQ==..49384adc3740924c8c882e6a9aeb409299edd7be_ZG9jL2V4YW1wbGVzL3RpbWVfc3RlcHBpbmcvM2Qvc3VibWl0XzFzaW11bC5weQ== 100644
--- a/doc/examples/time_stepping/3d/submit_1simul.py
+++ b/doc/examples/time_stepping/3d/submit_1simul.py
@@ -1,3 +1,12 @@
 from util_submit import submit_simul
 
-submit_simul(1.0, 800, "RK4", cfl_coef=None, nb_proc=None)
+nb_proc = 10
+nx = 400
+
+submit_simul(2 / 3, nx, "RK4", nb_proc=nb_proc, truncation_shape="cubic")
+submit_simul(1.0, nx, "RK4", nb_proc=nb_proc, truncation_shape="cubic")
+submit_simul(2 / 3, nx, "RK4", nb_proc=nb_proc)
+submit_simul(1.0, nx, "RK4", nb_proc=nb_proc)
+
+submit_simul(1.0, nx, "RK2_phaseshift_random", nb_proc=nb_proc)
+submit_simul(1.0, nx, "RK2_phaseshift_exact", nb_proc=nb_proc)
diff --git a/doc/examples/time_stepping/3d/util_submit.py b/doc/examples/time_stepping/3d/util_submit.py
index a74f6a2eeda8b09c5f46c087eac7e5bcc68804ef_ZG9jL2V4YW1wbGVzL3RpbWVfc3RlcHBpbmcvM2QvdXRpbF9zdWJtaXQucHk=..49384adc3740924c8c882e6a9aeb409299edd7be_ZG9jL2V4YW1wbGVzL3RpbWVfc3RlcHBpbmcvM2QvdXRpbF9zdWJtaXQucHk= 100644
--- a/doc/examples/time_stepping/3d/util_submit.py
+++ b/doc/examples/time_stepping/3d/util_submit.py
@@ -13,7 +13,12 @@
 
 
 def submit_simul(
-    coef_dealiasing, nx, type_time_scheme, cfl_coef=None, nb_proc=None
+    coef_dealiasing,
+    nx,
+    type_time_scheme,
+    cfl_coef=None,
+    nb_proc=None,
+    truncation_shape=None,
 ):
     nb_nodes = 1
     if nb_proc is None:
@@ -35,9 +40,6 @@
         f'--max-elapsed "{max_elapsed}"'
     )
 
-    if cfl_coef:
-        command += f" -cfl {cfl_coef}"
-
     name_run = f"{type_time_scheme}_trunc{coef_dealiasing:.3f}"
 
     if cfl_coef:
@@ -41,5 +43,6 @@
     name_run = f"{type_time_scheme}_trunc{coef_dealiasing:.3f}"
 
     if cfl_coef:
+        command += f" -cfl {cfl_coef}"
         name_run += f"_cfl{cfl_coef}"
 
@@ -44,5 +47,9 @@
         name_run += f"_cfl{cfl_coef}"
 
+    if truncation_shape is not None:
+        command += f" --truncation_shape {truncation_shape}"
+        name_run += f"_{truncation_shape}"
+
     print(f"submitting:\npython {command}")
 
     if not cluster: