From 3e99d1187336e5b05000f9aab3b6cd4d92ac2898 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Thu, 26 Jan 2023 13:41:07 +0100
Subject: [PATCH] script to store meta results

---
 wk/p2_heatflux_salpha_convergence.m | 19 +++++++++++++++++++
 1 file changed, 19 insertions(+)
 create mode 100644 wk/p2_heatflux_salpha_convergence.m

diff --git a/wk/p2_heatflux_salpha_convergence.m b/wk/p2_heatflux_salpha_convergence.m
new file mode 100644
index 00000000..2d850670
--- /dev/null
+++ b/wk/p2_heatflux_salpha_convergence.m
@@ -0,0 +1,19 @@
+figure
+%% KT 6.96, nuDGDK = 0.05, 128x64x24, Nexc 5
+P   = [4     6     8     10   ];
+Qx  = [67.62 67.50 59.21 64.17];
+std = [15.42 20.32 17.25 16.05];
+    errorbar(P,Qx,std/2,'s-',...
+    'LineWidth',2.0,...
+    'DisplayName','KT 6.96, nuDGDK 0.05'); hold on
+xlabel('$P$, $J=P/2$'); ylabel('$Q_x$');
+%% KT 5.3, nuDGDK = 0.05, 128x64x24, Nexc 5
+P   = [4     6     8     10   ];
+Qx  = [44.10 21.61 16.04 0.489];
+std = [10.61 6.952 4.166 0.061];
+    errorbar(P,Qx,std/2,'s-',...
+    'LineWidth',2.0,...
+    'DisplayName','KT 5.3, nuDGDK 0.05');
+xlabel('$P$, $J=P/2$'); ylabel('$Q_x$');
+
+title('GYAC, turb. heat flux conv.');
\ No newline at end of file
-- 
GitLab