diff --git a/E_corrs.m b/E_corrs.m
index 473f734..34222d3 100644
--- a/E_corrs.m
+++ b/E_corrs.m
@@ -7,7 +7,7 @@
 
 
 %%
-split='psilo' 
+split='main' 
 load(fullfile(['data/',split,'.mat']))
 numClusters = 4;
 
diff --git a/T_sweep_sps.m b/T_sweep_sps.m
index d7e8dfa..4504c59 100644
--- a/T_sweep_sps.m
+++ b/T_sweep_sps.m
@@ -12,7 +12,7 @@
 
 %% set inputs
 numClusters = 4;
-split='psilo' 
+split='main' 
 load(fullfile(['data/',split,'.mat']))
 
 
diff --git a/ami_calc.m b/ami_calc.m
index bf2052f..96b254b 100644
--- a/ami_calc.m
+++ b/ami_calc.m
@@ -5,7 +5,7 @@
 basedir = '/Users/sps253/Documents/energy_landscape';
 cd(basedir);
 
-split = 'psilo'
+split = 'main'
 
 load(fullfile(['data/',split,'.mat']))
 
diff --git a/complexity_measures.m b/complexity_measures.m
index e6755dd..d582ab4 100644
--- a/complexity_measures.m
+++ b/complexity_measures.m
@@ -8,7 +8,7 @@
 savedir = fullfile(basedir,'results','example');mkdir(savedir);	
 %% set inputs
 
-split='psilo' 
+split='main' 
 load(fullfile(['data/',split,'.mat']))
 numClusters = 4;
 
diff --git a/countclusters.m b/countclusters.m
index ee905f0..6c9c81b 100644
--- a/countclusters.m
+++ b/countclusters.m
@@ -7,7 +7,7 @@
 
 %% set inputs
 
-split='psilo'
+split='main'
 load(fullfile(['data/',split,'.mat']))
 savedir = fullfile(basedir,'results','example');mkdir(savedir);		% set save directory
 
diff --git a/elbow_sps.m b/elbow_sps.m
index 0af3e8c..59974e6 100644
--- a/elbow_sps.m
+++ b/elbow_sps.m
@@ -9,7 +9,7 @@
 basedir = '/Users/sps253/Documents/energy_landscape';
 cd(basedir);
 
-split = 'psilo'
+split = 'main'
 
 savedir = fullfile(basedir,'results','example');mkdir(savedir);
 load(fullfile(basedir,['data/',split,'.mat'])) % make sure this file matches the split you want to run below
diff --git a/m5HTcomparisons.m b/m5HTcomparisons.m
index c43648a..b2cc5b3 100644
--- a/m5HTcomparisons.m
+++ b/m5HTcomparisons.m
@@ -5,7 +5,7 @@
 
 %% set inputs
 numClusters = 4;
-split='psilo' 
+split='main' 
 load(fullfile(['data/',split,'.mat']))
 
 T = 0.001  % set time scale parameters based on values from paper
diff --git a/plotcentroidsSPS.m b/plotcentroidsSPS.m
index bbeeb0a..bf82479 100644
--- a/plotcentroidsSPS.m
+++ b/plotcentroidsSPS.m
@@ -9,7 +9,7 @@
 %% plot 
 
 numClusters = 4;
-split='psilo' 
+split='main' 
 
 load(fullfile(savedir,['Partition_bp',num2str(split),'_k',num2str(numClusters),'.mat']))
 
diff --git a/repeatkmeans_sps.m b/repeatkmeans_sps.m
index 210c735..66c9159 100644
--- a/repeatkmeans_sps.m
+++ b/repeatkmeans_sps.m
@@ -6,7 +6,7 @@
 
 %% load data
 
-split = 'psilo' %%<- options: main, gsr, sch, music, psilo
+split = 'main' %%<- options: main, gsr, sch, music, psilo
 
 
 load(fullfile(['data/',split,'.mat']))
diff --git a/si_scripts/energy_as_increasing_k.m b/si_scripts/energy_as_increasing_k.m
index f23bba6..5bd7d58 100644
--- a/si_scripts/energy_as_increasing_k.m
+++ b/si_scripts/energy_as_increasing_k.m
@@ -28,7 +28,7 @@
 %% load partitions from elbow plot
 
 maxk = round(sqrt(TR)) -1; %k^2 must be less than TR to capture all transitions
-k_rng=2:14;
+k_rng=2:maxk;
 
 C_sub = [repmat([1 0 0],nsubjs,1); repmat([0 0 0],nsubjs,1)];
 C_avg = [[1 0 0];[0 0 0]];
diff --git a/spin_receptor_map.m b/spin_receptor_map.m
index 601bb09..e71da0d 100644
--- a/spin_receptor_map.m
+++ b/spin_receptor_map.m
@@ -8,7 +8,7 @@
 savedir = fullfile(basedir,'results','example');mkdir(savedir);		% set save directory
 
 %%
-split='psilo' 
+split='main' 
 load(fullfile(['data/',split,'.mat']))
 
 numClusters=4; T=0.001; 
diff --git a/subcentroids.m b/subcentroids.m
index 6f86160..04fb217 100644
--- a/subcentroids.m
+++ b/subcentroids.m
@@ -8,7 +8,7 @@
 %% load BOLD data
 
 
-split='psilo'
+split='main'
 load(fullfile(['data/',split,'.mat']))
 
 numClusters=4;
diff --git a/subj_energy.m b/subj_energy.m
index b314b7c..24a9104 100644
--- a/subj_energy.m
+++ b/subj_energy.m
@@ -10,7 +10,7 @@
 T = 0.001 % set time scale parameters based on values from T_sweep_sps.m
 
 numClusters=4;
-split='psilo'
+split='main'
 load(fullfile(['data/',split,'.mat']))
 
 load(['Partition_bp',num2str(split),'_k',num2str(numClusters),'.mat'],'partition','centroids','clusterNames'); %make sure file matches centroids you want to initialize from
diff --git a/systems_plot_sps.m b/systems_plot_sps.m
index 4f62130..6cc8888 100644
--- a/systems_plot_sps.m
+++ b/systems_plot_sps.m
@@ -5,7 +5,7 @@
 addpath(genpath('code'))
 %% set inputs
 numClusters = 4;
-split='psilo'
+split='main'
 
 savedir = fullfile(basedir,'results','example');mkdir(savedir);		% set save directory
 load(fullfile(savedir,['Partition_bp',num2str(split),'_k',num2str(numClusters),'.mat']))
diff --git a/transProbs.m b/transProbs.m
index 51d9fcf..c317ae2 100644
--- a/transProbs.m
+++ b/transProbs.m
@@ -6,7 +6,7 @@
 
 %% set inputs
 numClusters = 4;
-split='psilo'
+split='main'
 load(fullfile(['data/',split,'.mat']))
 savedir = fullfile(basedir,'results','example');mkdir(savedir);		% set save directory
 load(fullfile(savedir,['Partition_bp',num2str(split),'_k',num2str(numClusters),'.mat']))