diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster_in_plane_trap.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster_in_plane_trap.m
index 6fd87a8..e6584c2 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster_in_plane_trap.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster_in_plane_trap.m
@@ -17,7 +17,7 @@ OptionsStruct.TrapFrequencies         = [vz/AspectRatio, vz/AspectRatio, 500];
 OptionsStruct.TrapPotentialType       = 'Harmonic';
 
 OptionsStruct.NumberOfGridPoints      = [128, 128];
-OptionsStruct.Dimensions              = [10, 10];      
+OptionsStruct.Dimensions              = [20, 20];      
 OptionsStruct.TimeStepSize            = 0.005;           % in s
 OptionsStruct.MinimumTimeStepSize     = 1E-5;            % in s
 OptionsStruct.TimeCutOff              = 2E6;             % in s
@@ -26,14 +26,14 @@ OptionsStruct.ResidualTolerance       = 1E-05;
 OptionsStruct.NoiseScaleFactor        = 0.05;
 
 OptionsStruct.MaxIterations           = 10;      
-OptionsStruct.VariationalWidth        = 1.2;       
+OptionsStruct.VariationalWidth        = 1.3;       
 OptionsStruct.WidthLowerBound         = 0.01;
 OptionsStruct.WidthUpperBound         = 12;
 OptionsStruct.WidthCutoff             = 5e-3;
 
-OptionsStruct.PlotLive                = true;
+OptionsStruct.PlotLive                = false;
 OptionsStruct.JobNumber               = 1;
-OptionsStruct.RunOnGPU                = false;
+OptionsStruct.RunOnGPU                = true;
 OptionsStruct.SaveData                = true;
 OptionsStruct.SaveDirectory           = './Results/Data_TiltingOfDipoles/HarmonicTrap/Hz500';
 options                               = Helper.convertstruct2cell(OptionsStruct);
@@ -70,14 +70,14 @@ OptionsStruct.ResidualTolerance       = 1E-05;
 OptionsStruct.NoiseScaleFactor        = 0.05;
 
 OptionsStruct.MaxIterations           = 10;      
-OptionsStruct.VariationalWidth        = 1.2;       
+OptionsStruct.VariationalWidth        = 1.3;       
 OptionsStruct.WidthLowerBound         = 0.01;
 OptionsStruct.WidthUpperBound         = 12;
 OptionsStruct.WidthCutoff             = 5e-3;
 
-OptionsStruct.PlotLive                = true;
+OptionsStruct.PlotLive                = false;
 OptionsStruct.JobNumber               = 2;
-OptionsStruct.RunOnGPU                = false;
+OptionsStruct.RunOnGPU                = true;
 OptionsStruct.SaveData                = true;
 OptionsStruct.SaveDirectory           = './Results/Data_TiltingOfDipoles/HarmonicTrap/Hz500';
 options                               = Helper.convertstruct2cell(OptionsStruct);
@@ -88,97 +88,4 @@ pot                                   = VariationalSolver2D.Potentials(options{:
 solver.Potential                      = pot.trap();
 
 %-% Run Solver %-%
-[Params, Transf, psi, V, VDk]         = solver.run();
-
-%% v_z = 500, theta = 0: a_s = 76.41
-
-OptionsStruct = struct;
-
-OptionsStruct.NumberOfAtoms           = 101250;     
-OptionsStruct.DipolarPolarAngle       = 0;
-OptionsStruct.DipolarAzimuthAngle     = 0;
-OptionsStruct.ScatteringLength        = 76.41;        
-
-vz                                    = 500;
-AspectRatio                           = 5;
-OptionsStruct.TrapFrequencies         = [vz/AspectRatio, vz/AspectRatio, 500];
-OptionsStruct.TrapPotentialType       = 'Harmonic';
-
-OptionsStruct.NumberOfGridPoints      = [128, 128];
-OptionsStruct.Dimensions              = [10, 10];      
-OptionsStruct.TimeStepSize            = 0.005;           % in s
-OptionsStruct.MinimumTimeStepSize     = 1E-5;            % in s
-OptionsStruct.TimeCutOff              = 2E6;             % in s
-OptionsStruct.EnergyTolerance         = 5E-10;
-OptionsStruct.ResidualTolerance       = 1E-05;
-OptionsStruct.NoiseScaleFactor        = 0.05;
-
-OptionsStruct.MaxIterations           = 10;      
-OptionsStruct.VariationalWidth        = 1.2;       
-OptionsStruct.WidthLowerBound         = 0.01;
-OptionsStruct.WidthUpperBound         = 12;
-OptionsStruct.WidthCutoff             = 5e-3;
-
-OptionsStruct.PlotLive                = true;
-OptionsStruct.JobNumber               = 3;
-OptionsStruct.RunOnGPU                = false;
-OptionsStruct.SaveData                = true;
-OptionsStruct.SaveDirectory           = './Results/Data_TiltingOfDipoles/HarmonicTrap/Hz500';
-options                               = Helper.convertstruct2cell(OptionsStruct);
-clear OptionsStruct
-
-solver                                = VariationalSolver2D.DipolarGas(options{:});
-pot                                   = VariationalSolver2D.Potentials(options{:});
-solver.Potential                      = pot.trap();
-
-%-% Run Solver %-%
-[Params, Transf, psi, V, VDk]         = solver.run();
-
-%% v_z = 500, theta = 15: a_s = 77.45
-
-OptionsStruct = struct;
-
-OptionsStruct.NumberOfAtoms           = 101250;     
-OptionsStruct.DipolarPolarAngle       = deg2rad(15);
-OptionsStruct.DipolarAzimuthAngle     = 0;
-OptionsStruct.ScatteringLength        = 77.45;        
-
-vz                                    = 500;
-AspectRatio                           = 5;
-OptionsStruct.TrapFrequencies         = [vz/AspectRatio, vz/AspectRatio, 500];
-OptionsStruct.TrapPotentialType       = 'Harmonic';
-
-OptionsStruct.NumberOfGridPoints      = [128, 128];
-OptionsStruct.Dimensions              = [10, 10];      
-OptionsStruct.TimeStepSize            = 0.005;           % in s
-OptionsStruct.MinimumTimeStepSize     = 1E-5;            % in s
-OptionsStruct.TimeCutOff              = 2E6;             % in s
-OptionsStruct.EnergyTolerance         = 5E-10;
-OptionsStruct.ResidualTolerance       = 1E-05;
-OptionsStruct.NoiseScaleFactor        = 0.05;
-
-OptionsStruct.MaxIterations           = 10;      
-OptionsStruct.VariationalWidth        = 1.2;       
-OptionsStruct.WidthLowerBound         = 0.01;
-OptionsStruct.WidthUpperBound         = 12;
-OptionsStruct.WidthCutoff             = 5e-3;
-
-OptionsStruct.PlotLive                = true;
-OptionsStruct.JobNumber               = 4;
-OptionsStruct.RunOnGPU                = false;
-OptionsStruct.SaveData                = true;
-OptionsStruct.SaveDirectory           = './Results/Data_TiltingOfDipoles/HarmonicTrap/Hz500';
-options                               = Helper.convertstruct2cell(OptionsStruct);
-clear OptionsStruct
-
-solver                                = VariationalSolver2D.DipolarGas(options{:});
-pot                                   = VariationalSolver2D.Potentials(options{:});
-solver.Potential                      = pot.trap();
-
-%-% Run Solver %-%
-[Params, Transf, psi, V, VDk]         = solver.run();
-%% - Analysis
-SaveDirectory = './Results/Data_TiltingOfDipoles/HarmonicTrap/Hz500';
-JobNumber     = 4;
-% Plotter.visualizeGSWavefunction2D(SaveDirectory, JobNumber)
-[contrast, period_X, period_Y] = Scripts.analyzeGSWavefunction(SaveDirectory, JobNumber);
\ No newline at end of file
+[Params, Transf, psi, V, VDk]         = solver.run();
\ No newline at end of file