diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
index 3f5fb3e..eb1778f 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
@@ -553,50 +553,6 @@ SaveDirectory = './Results/Data_3D/GradientDescent/aS_079_theta_000_phi_000_N_50
 JobNumber     = 0;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
 %%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_080_theta_000_phi_000_N_50000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_081_theta_000_phi_000_N_50000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_082_theta_000_phi_000_N_60000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_083_theta_000_phi_000_N_60000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_084_theta_000_phi_000_N_70000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_085_theta_000_phi_000_N_70000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_086_theta_000_phi_000_N_80000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_087_theta_000_phi_000_N_100000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_089_theta_000_phi_000_N_105000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_084_theta_000_phi_000_N_65000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
-SaveDirectory = './Results/Data_3D/GradientDescent/aS_090_theta_000_phi_000_N_70000';
-JobNumber     = 0;
-Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
-%%
 SaveDirectory = './Results/Data_3D/GradientDescent/';
 JobNumber     = 1;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
diff --git a/Dipolar-Gas-Simulator/submit_jobs.sh b/Dipolar-Gas-Simulator/submit_jobs.sh
index e086fa2..7cbe973 100644
--- a/Dipolar-Gas-Simulator/submit_jobs.sh
+++ b/Dipolar-Gas-Simulator/submit_jobs.sh
@@ -2,9 +2,9 @@
 
 # Use space-separated floating-point/integer values
 SCATTERING_LENGTH_RANGE="[79.0 80.0 81.0 82.0 83.0 84.0 85.0 86.0 87.0 88.0 89.0 90.0]"
-POLAR_ANGLE_RANGE="[0.0]"
+POLAR_ANGLE_RANGE="[20.0 30.0 40.0 50.0]"
 AZIMUTHAL_ANGLE_RANGE="[0.0]"
-NUM_ATOMS_LIST="[50000 60000 70000 80000]"
+NUM_ATOMS_LIST="[50000 55000 60000 65000 70000 75000 80000 85000 90000 95000 100000 105000]"
 CHUNK_SIZE=4
 
 # ----------- Count total combinations for SLURM array -----------
@@ -32,7 +32,7 @@ sbatch --export=SCATTERING_LENGTH_RANGE="$SCATTERING_LENGTH_RANGE",POLAR_ANGLE_R
 #SBATCH --cpus-per-task=8
 #SBATCH --gres=gpu:1
 #SBATCH --mem=16G
-#SBATCH --time=03:30:00
+#SBATCH --time=04:00:00
 #SBATCH --job-name=simulation
 #SBATCH --error=simulation_%A_%a.err
 #SBATCH --output=simulation_%A_%a.out