finalised scripts for release

runs/idealised_delaunay.py

@@ -120,8 +120,8 @@ def autoreload():
 
 params = var.params()
 params.plot = True
-params.lmbda_fa = 0.0
 params.lmbda_sa = 0.1
+params.lmbda_sa = 1e-6
 dplot = diagnostics.diag_plotter(params, nhi, nhj)
 dplot.show((0, 1), sols, kls=kls_ref, v_extent=v_extent, dfft_plot=True)
 
@@ -391,7 +391,7 @@ def autoreload():
         v_extent=v_extent,
         ir_args=ir_args,
         phys_lbls=["", ""],
-        fn=fn,
+        fn="",
     )
     print(uw_01_rc.sum())
 
@@ -413,7 +413,7 @@ def autoreload():
         v_extent=v_extent,
         ir_args=ir_args,
         phys_lbls=["", ""],
-        fn=fn,
+        fn="",
     )
     print(uw_02_rc.sum())
 

runs/idealised_isosceles.py

@@ -274,6 +274,7 @@ def csam_run(cell, n_modes, lmbda_fg, lmbda_sg):
 
 # %%
 #### plot the errors
+print("amplitudes:")
 print(sums)
 plotter.error_bar_abs_plot(
     selected_errs,
@@ -302,7 +303,10 @@ def csam_run(cell, n_modes, lmbda_fg, lmbda_sg):
 # %%
 #### print the errors
 np.set_printoptions(suppress=True)
+print("percentage error in amplitude:")
 print(np.around(sum_errs, 5) * 100)
+print("")
+print("L2-errors:")
 print(errs)
 
 # %%