Skip to content
Snippets Groups Projects
Commit 0a0786eb authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

store old scripts

parent c2ec0ff1
Branches
Tags
No related merge requests found
...@@ -9,12 +9,12 @@ CLUSTER.TIME = '99:00:00'; % allocation time hh:mm:ss ...@@ -9,12 +9,12 @@ CLUSTER.TIME = '99:00:00'; % allocation time hh:mm:ss
%% %%
SIMID = 'p2_linear_new'; % Name of the simulation SIMID = 'p2_linear_new'; % Name of the simulation
RERUN = 0; % rerun if the data does not exist RERUN = 0; % rerun if the data does not exist
RUN = 1; RUN = 0;
% KT_a = [3:0.5:6.5 6.96]; KT_a = [3:0.5:6.5 6.96];
KT_a = [6.96]; % KT_a = [6.96];
P_a = [2]; % P_a = [2];
% P_a = [3:1:29]; % P_a = [3:1:29];
% P_a = 2:2:10; P_a = 2:2:10;
J_a = floor(P_a/2); J_a = floor(P_a/2);
% collision setting % collision setting
CO = 'DG'; CO = 'DG';
......
...@@ -10,6 +10,8 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0' ...@@ -10,6 +10,8 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0'
% datafname = 'p2_linear/8x24_ky_0.3_P_8_J_4_kT_3_6.96_nu_0_1_DGDK.mat'; % datafname = 'p2_linear/8x24_ky_0.3_P_8_J_4_kT_3_6.96_nu_0_1_DGDK.mat';
% datafname = 'p2_linear/8x24_ky_0.3_P_16_J_8_kT_3_6.96_nu_0_1_DGDK.mat'; % datafname = 'p2_linear/8x24_ky_0.3_P_16_J_8_kT_3_6.96_nu_0_1_DGDK.mat';
%% Scans over KT and PJ, keeping ky, CO constant (CBC_kT_PJ_scan.m) %% Scans over KT and PJ, keeping ky, CO constant (CBC_kT_PJ_scan.m)
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_30_DGDK_0.001.mat';
% datafname = 'p2_linear_check/8x24_ky_0.3_kT_3_9_P_2_30_DGDK_0.001.mat';
% datafname = 'p2_linear/8x24_ky_0.3_kT_3_6.96_P_2_40_DGDK_0.05.mat'; % datafname = 'p2_linear/8x24_ky_0.3_kT_3_6.96_P_2_40_DGDK_0.05.mat';
% datafname = 'p2_linear/8x24_ky_0.3_kT_3_6.96_P_2_30_DGDK_0.05.mat'; % datafname = 'p2_linear/8x24_ky_0.3_kT_3_6.96_P_2_30_DGDK_0.05.mat';
% datafname = 'p2_linear/8x24_ky_0.3_kT_3_6.96_P_2_40_DGDK_0.025.mat'; % datafname = 'p2_linear/8x24_ky_0.3_kT_3_6.96_P_2_40_DGDK_0.025.mat';
...@@ -46,7 +48,7 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0' ...@@ -46,7 +48,7 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0'
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_10_DGDK_0.05.mat'; % datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_10_DGDK_0.05.mat';
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_30_DGGK_0.05.mat'; % datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_30_DGGK_0.05.mat';
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_20_SGGK_0.05.mat'; % datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_20_SGGK_0.05.mat';
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_10_LDGK_0.05.mat'; datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_10_LDGK_0.05.mat';
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_20_DGDK_0.1.mat'; % datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_20_DGDK_0.1.mat';
% datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_20_DGGK_0.1.mat'; % datafname = 'p2_linear_new/8x24_ky_0.3_kT_3_6.96_P_2_20_DGGK_0.1.mat';
...@@ -74,7 +76,7 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0' ...@@ -74,7 +76,7 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0'
% datafname = 'p2_linear_new/8x24_P_4_ky_0.05_1_SGGK_nu_0.005_0.05_kT_5.3.mat'; % datafname = 'p2_linear_new/8x24_P_4_ky_0.05_1_SGGK_nu_0.005_0.05_kT_5.3.mat';
% datafname = 'p2_linear_new/8x24_P_8_ky_0.05_1_SGGK_nu_0.005_0.05_kT_5.3.mat'; % datafname = 'p2_linear_new/8x24_P_8_ky_0.05_1_SGGK_nu_0.005_0.05_kT_5.3.mat';
datafname = 'p2_linear_new/8x24_P_16_ky_0.05_1_SGGK_nu_0.005_0.05_kT_5.3.mat'; % datafname = 'p2_linear_new/8x24_P_16_ky_0.05_1_SGGK_nu_0.005_0.05_kT_5.3.mat';
% datafname = 'p2_linear_new/8x24_P_4_ky_0.05_1_LDGK_nu_0.005_0.05_kT_5.3.mat'; % datafname = 'p2_linear_new/8x24_P_4_ky_0.05_1_LDGK_nu_0.005_0.05_kT_5.3.mat';
% datafname = 'p2_linear_new/8x24_P_8_ky_0.05_1_LDGK_nu_0.005_0.05_kT_5.3.mat'; % datafname = 'p2_linear_new/8x24_P_8_ky_0.05_1_LDGK_nu_0.005_0.05_kT_5.3.mat';
...@@ -91,13 +93,13 @@ if FILTERGAMMA ...@@ -91,13 +93,13 @@ if FILTERGAMMA
d.data = d.data.*(d.data>0.025); d.data = d.data.*(d.data>0.025);
d.err = d.err.*(d.data>0.025); d.err = d.err.*(d.data>0.025);
end end
if 0 if 1
%% Pcolor of the peak %% Pcolor of the peak
figure; figure;
% [XX_,YY_] = meshgrid(d.s1,d.s2); % [XX_,YY_] = meshgrid(d.s1,d.s2);
[XX_,YY_] = meshgrid(1:numel(d.s1),1:numel(d.s2)); [XX_,YY_] = meshgrid(1:numel(d.s1),1:numel(d.s2));
pclr=imagesc_custom(XX_,YY_,d.data'.*(d.data>0)'); % pclr=imagesc_custom(XX_,YY_,d.data'.*(d.data>0)');
% pclr=contourf(1:numel(d.s1),1:numel(d.s2),d.data'.*(d.data>0)'); pclr=contourf(1:numel(d.s1),1:numel(d.s2),d.data'.*(d.data>0)');
% pclr=surf(1:numel(d.s1),1:numel(d.s2),d.data'.*(d.data>0)'); % pclr=surf(1:numel(d.s1),1:numel(d.s2),d.data'.*(d.data>0)');
title(d.title); title(d.title);
xlabel(d.s1name); ylabel(d.s2name); xlabel(d.s1name); ylabel(d.s2name);
...@@ -110,7 +112,7 @@ clb.Label.String = '$\gamma c_s/R$'; ...@@ -110,7 +112,7 @@ clb.Label.String = '$\gamma c_s/R$';
clb.Label.Interpreter = 'latex'; clb.Label.Interpreter = 'latex';
clb.Label.FontSize= 18; clb.Label.FontSize= 18;
end end
if 1 if 0
%% Scan along first dimension %% Scan along first dimension
figure figure
colors_ = jet(numel(d.s2)); colors_ = jet(numel(d.s2));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment