diff --git a/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m b/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m
index 404d771..e68a59d 100644
--- a/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m
+++ b/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m
@@ -1,6 +1,6 @@
 function [evals, modes] = solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem)
 
-Size         = length(psi(:));
+Size         = length(psi(:)); % 2-D matrix will be unwrapped and the corresponding matrix of (N^2)^2 elements solved for.
 Neigs        = length(psi(:));
 opts.tol     = 1e-16;
 opts.disp    = 1;
diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
index 4f08fb3..683d5d0 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
@@ -128,6 +128,9 @@ muchem = solver.Calculator.calculateChemicalPotential(psi,Params,VParams,Transf,
 
 [evals, modes] = BdGSolver2D.solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem);
 
+% Save the eigenvalues and eigenvectors to a .mat file
+save(sprintf(strcat(solver.SaveDirectory, '/Run_%03i/bdg_eigen_data.mat'),solver.JobNumber), 'evals', 'modes');
+
 %% - Create Variational2D and Calculator object with specified options
 
 OptionsStruct = struct;
diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m
index ff0cf5e..afc632c 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m
@@ -204,3 +204,6 @@ V = 0.0*(Params.gx.*X.^2+Params.gy.*Y.^2);
 muchem = solver.Calculator.calculateChemicalPotential(psi,Params,VParams,Transf,VDk,V);
 
 [evals, modes] = BdGSolver2D.solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem);
+
+% Save the eigenvalues and eigenvectors to a .mat file
+save(sprintf(strcat(solver.SaveDirectory, '/Run_%03i/bdg_eigen_data.mat'),Params.njob), 'evals', 'modes');