diff --git a/wk/fast_analysis.m b/wk/fast_analysis.m index 9dcf610fc9eef151688a87fd9964da13fbaeffc5..e7c14fd1e1bd85c049d1d305e97fc81568f5b002 100644 --- a/wk/fast_analysis.m +++ b/wk/fast_analysis.m @@ -5,13 +5,15 @@ addpath(genpath([gyacomodir,'matlab/compute'])) % ... add addpath(genpath([gyacomodir,'matlab/load'])) % ... add default_plots_options % Partition of the computer where the data have to be searched -PARTITION=''; - +% PARTITION='/Users/ahoffmann/gyacomo/results/paper_3/'; +PARTITION='/misc/gyacomo23_outputs/paper_3/'; %% Paper 3 -resdir = '/Users/ahoffmann/gyacomxlio/results/paper_3/DTT_rho85/3x2x192x48x32'; -% resdir = '/Users/ahoffmann/gyacomo/results/paper_3/DTT_rho85/3x2x192x48x32_NT'; -% resdir = '/Users/ahoffmann/gyacomo/results/paper_3/DTT_rho98/3x2x192x48x32'; -resdir = '/Users/ahoffmann/gyacomo/results/paper_3/DTT_rho98/3x2x192x48x32_0.25grad'; +% resdir = 'DTT_rho85/3x2x192x48x32'; +% resdir = 'DTT_rho85/3x2x192x48x32_NT'; +% resdir = 'DTT_rho98/3x2x192x48x32'; +% resdir = 'DTT_rho98/3x2x192x48x32_0.25grad'; +% resdir = 'LM_DIIID_rho95/5x3x512x92x32'; +resdir = 'LM_DIIID_rho95/3x2x512x92x32'; %% J0 = 00; J1 = 02; diff --git a/wk/lin_run_script.m b/wk/lin_run_script.m index acb62f80a1fd2ab976f78e2ab808ac7253f95876..97644a1790f9d1067e49b2f388db569bb269b40b 100644 --- a/wk/lin_run_script.m +++ b/wk/lin_run_script.m @@ -31,7 +31,8 @@ run lin_JET_rho97 % run lin_KBM %% Change parameters NY = 2; -NX = 4; +NX = 6; +NZ = 32; PMAX = 2; JMAX = PMAX/2; ky = 0.5; diff --git a/wk/lin_scan_script.m b/wk/lin_scan_script.m index ffdd64ced3af101493570f2d21519387fc696c9e..8831740be18598ee3c075df149aef21272dbafd6 100644 --- a/wk/lin_scan_script.m +++ b/wk/lin_scan_script.m @@ -6,8 +6,8 @@ %% Set up the paths for the necessary Matlab modules gyacomodir = pwd; gyacomodir = gyacomodir(1:end-2); -% mpirun = 'mpirun'; -mpirun = '/opt/homebrew/bin/mpirun'; % for macos +mpirun = 'mpirun'; +% mpirun = '/opt/homebrew/bin/mpirun'; % for macos addpath(genpath([gyacomodir,'matlab'])) % Add matlab folder addpath(genpath([gyacomodir,'matlab/plot'])) % Add plot folder addpath(genpath([gyacomodir,'matlab/compute'])) % Add compute folder @@ -18,8 +18,8 @@ addpath(genpath([gyacomodir,'wk/parameters'])) % Add parameters folder RUN = 1; % To run or just to load RERUN = 0; % rerun if the data does not exist default_plots_options -% EXECNAME = 'gyacomo23_sp'; % single precision -EXECNAME = 'gyacomo23_dp'; % double precision +EXECNAME = 'gyacomo23_sp'; % single precision +%EXECNAME = 'gyacomo23_dp'; % double precision %% Setup parameters % run lin_DTT_AB_rho85 @@ -35,6 +35,7 @@ SIMID = [SIMID,'_scan']; P_a = [2 4 6 8]; % P_a = 2; ky_a = logspace(-1.5,1.5,30); +CO = 'SG'; %% Scan loop % arrays for the result g_ky = zeros(numel(ky_a),numel(P_a)); @@ -47,7 +48,7 @@ for PMAX = P_a i = 1; for ky = ky_a LY = 2*pi/ky; - DT = 1e-3*log(ky/0.05);%min(1e-2,1e-3/ky); + DT = 1e-5;%/(1+log(ky/0.05));%min(1e-2,1e-3/ky); TMAX = 20;%min(10,1.5/ky); DTSAVE0D = 0.1; DTSAVE3D = 0.1;