diff --git a/Dipolar-Gas-Simulator/+Scripts/editPhaseDiagram.m b/Dipolar-Gas-Simulator/+Scripts/editPhaseDiagram.m
index 8f26492..2a72f42 100644
--- a/Dipolar-Gas-Simulator/+Scripts/editPhaseDiagram.m
+++ b/Dipolar-Gas-Simulator/+Scripts/editPhaseDiagram.m
@@ -22,8 +22,8 @@ function editPhaseDiagram(M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST)
     yticks(1:ylen);
     xticklabels(string(NUM_ATOMS_LIST));
     yticklabels(string(SCATTERING_LENGTH_RANGE));
-    xlabel('Number of Atoms');
-    ylabel('Scattering Length (a₀)');
+    xlabel('Number of Atoms', 'Interpreter', 'tex', 'FontSize', 16);
+    ylabel('Scattering Length (\times a_o)', 'Interpreter', 'tex', 'FontSize', 16);
     grid on;
     title('Click a cell to modify value. Press "q" to quit.');
 
diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
index 91ed875..29414a9 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
@@ -613,8 +613,7 @@ for i = 1:length(SCATTERING_LENGTH_RANGE)
 end
 
 %% Visualize phase diagram
-% load('./Results/phase_diagram_matrix.mat')
-load('./Results/edited_phase_diagram.mat')
+load('./Results/PhaseDiagramUntilted.mat')
 PhaseDiagramMatrix   = M;
 ScatteringLengths    = SCATTERING_LENGTH_RANGE;
 NumberOfAtoms        = round(NUM_ATOMS_LIST * 1E-5,2);
@@ -623,32 +622,32 @@ ylen                 = length(ScatteringLengths);
 
 figure(1);
 set(gcf,'Position',[100 100 950 750])
+fig.WindowState = 'maximized';
+clf
+set(gcf,'Position', [100, 100, 1600, 900])
 set(gca,'FontSize',16,'Box','On','Linewidth',2);
 hImg = imagesc(M);
 set(gca, 'YDir', 'normal');
 colormap(parula);
-ax = gca;                % Get current axes
-ax.FontSize = 16;        % Set tick label font size (adjust as needed)
+colorbar;
 axis equal tight;
 xticks(1:xlen);
 yticks(1:ylen);
-xticklabels(string(NumberOfAtoms));
-yticklabels(string(ScatteringLengths));
-xlabel('Number of Atoms (x 1E5)', FontSize=16);
-ylabel('Scattering Length (a0)', FontSize=16);
-% title('Zero-temperature Phase Diagram for \theta = 0', 'Interpreter', 'tex', FontSize=16);
+xticklabels(string(NUM_ATOMS_LIST));
+yticklabels(string(SCATTERING_LENGTH_RANGE));
+xlabel('Number of Atoms', 'Interpreter', 'tex', 'FontSize', 16);
+ylabel('Scattering Length (\times a_o)', 'Interpreter', 'tex', 'FontSize', 16);
 grid on;
 
 %% Edit phase diagram
-% load('./Results/phase_diagram_matrix.mat')
-load('./Results/edited_phase_diagram.mat')
+load('./Results/PhaseDiagramUntilted.mat')
 PhaseDiagramMatrix   = M;
 ScatteringLengths    = SCATTERING_LENGTH_RANGE;
 NumberOfAtoms        = NUM_ATOMS_LIST;
 Scripts.editPhaseDiagram(M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST)
 
 %% Smoothen phase diagram
-load('./Results/edited_phase_diagram.mat');  % Load M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST
+load('./Results/PhaseDiagramUntilted.mat');  % Load M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST
 titleString = "[ \omega_x, \omega_y, \omega_z ] = 2 \pi \times [ 50, 20, 150 ] Hz; \theta = 0";
 Scripts.plotSmoothedPhaseDiagram(M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST, titleString);