From 237d3fff169cac1c493ee845044d049b7f81e632 Mon Sep 17 00:00:00 2001
From: Antoine <antoine.hoffmann@epfl.ch>
Date: Fri, 4 Aug 2023 09:54:56 +0200
Subject: [PATCH] cleaning
---
wk/load_metadata_scan.m | 14 +++++---------
1 file changed, 5 insertions(+), 9 deletions(-)
diff --git a/wk/load_metadata_scan.m b/wk/load_metadata_scan.m
index 6615ee53..f107d4af 100644
--- a/wk/load_metadata_scan.m
+++ b/wk/load_metadata_scan.m
@@ -15,13 +15,9 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0'
% datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_1_P_2_8_DGGK_0.05_be_0.0039.mat';
% datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_0.9_P_2_8_kN_1.33_DGGK_0.56901_be_0.0039.mat';
% datafname = 'lin_DTT_AB_rho85_PT_scan/8x24_ky_0.05_1.5_P_2_8_kN_1.33_DGGK_0.1_be_0.0039.mat';
-% datafname = 'lin_JET_rho97_scan/4x32_ky_0.01_10_P_2_8_kN_10_DGGK_0.1_be_0.0031.mat';
-% datafname = 'lin_JET_rho97_scan/4x32_ky_0.01_10_P_2_8_kN_10_DGGK_0.05_be_0.0031.mat';
-% datafname = 'lin_JET_rho97_scan/8x32_ky_0.01_10_P_2_8_kN_10_DGGK_0.05_be_0.0031.mat';
-% datafname = 'lin_JET_rho97_scan/8x32_ky_0.01_10_P_2_2_kN_10_DGGK_0.05_be_0.0031.mat';
-datafname = 'lin_JET_rho97_scan/8x32_ky_0.031623_31.6228_P_2_2_kN_10_DGGK_0.1_be_0.0031.mat';
+datafname = 'lin_JET_rho97_scan/8x32_ky_0.031623_31.6228_P_1_1_kN_31.4286_DGGK_0.1_be_0.0031.mat';
%% Chose if we filter gamma>0.05
-FILTERGAMMA = 0;
+FILTERGAMMA = 1;
%% Load data
fname = ['../results/',datafname];
@@ -29,7 +25,7 @@ d = load(fname);
gamma = real(d.data); g_err = real(d.err);
omega = imag(d.data); w_err = imag(d.err);
if FILTERGAMMA
- d.data = real(gamma).*(real(gamma)>0.025) + imag(gamma);
+ gamma = gamma.*(gamma>0.025);
end
if 0
%% Pcolor of the peak
@@ -46,7 +42,7 @@ set(gca,'YTick',1:numel(d.s2),'YTicklabel',d.s2)
colormap(jet)
colormap(bluewhitered)
clb=colorbar;
-clb.Label.String = '$\gamma c_s/R$';
+clb.Label.String = '$\gamma R/c_s$';
clb.Label.Interpreter = 'latex';
clb.Label.FontSize= 18;
end
@@ -83,7 +79,7 @@ for i = 1:numel(d.s2)
'color',colors_(i,:));
hold on;
end
-xlabel(d.s1name); ylabel('$\omega c_s/R$');title(d.title);
+xlabel(d.s1name); ylabel('$\omega R/c_s$');title(d.title);
xlim([d.s1(1) d.s1(end)]);
colormap(colors_);
--
GitLab