diff --git a/script/jitter_scan.py b/script/jitter_scan.py
index 2468bb1..8c189b8 100644
--- a/script/jitter_scan.py
+++ b/script/jitter_scan.py
@@ -141,12 +141,18 @@ class all_statistics_data (object):
  m5    = statistics_data()
  m10   = statistics_data()
 
- def return_string(self):
+ def return_string(self, fmt_type):
      # give results back as formatted string
+     # fmt_type 1 = phase / 0 = amplt
      ret_str = "Statistics:\tmin\tmean\tmax\tstdev"
-     ret_str = ret_str + "\nSaturation    :\t" + "%.3e" %  (self.sat.min) + "\t" + "%.3e" %  (self.sat.mean) + "\t" + "%.3e" %  (self.sat.max) + "\t" + "%.3e" %  (self.sat.stdev)
-     ret_str = ret_str + "\n5% out of sat :\t" + "%.3e" %  (self.m5.min)  + "\t" + "%.3e" %  (self.m5.mean)  + "\t" + "%.3e" %  (self.m5.max)  + "\t" + "%.3e" %  (self.m5.stdev)
-     ret_str = ret_str + "\n10% out of sat:\t" + "%.3e" %  (self.m10.min) + "\t" + "%.3e" %  (self.m10.mean) + "\t" + "%.3e" %  (self.m10.max) + "\t" + "%.3e" %  (self.m10.stdev)
+     if fmt_type = 1:
+         ret_str = ret_str + "\nSaturation    :\t" + "%.3e" %  (self.sat.min) + "\t" + "%.3e" %  (self.sat.mean) + "\t" + "%.3e" %  (self.sat.max) + "\t" + "%.3e" %  (self.sat.stdev)
+         ret_str = ret_str + "\n5% out of sat :\t" + "%.3e" %  (self.m5.min)  + "\t" + "%.3e" %  (self.m5.mean)  + "\t" + "%.3e" %  (self.m5.max)  + "\t" + "%.3e" %  (self.m5.stdev)
+         ret_str = ret_str + "\n10% out of sat:\t" + "%.3e" %  (self.m10.min) + "\t" + "%.3e" %  (self.m10.mean) + "\t" + "%.3e" %  (self.m10.max) + "\t" + "%.3e" %  (self.m10.stdev)
+     else
+         ret_str = ret_str + "\nSaturation    :\t" + "%.3f" %  (self.sat.min) + "\t" + "%.3f" %  (self.sat.mean) + "\t" + "%.3f" %  (self.sat.max) + "\t" + "%.3f" %  (self.sat.stdev)
+         ret_str = ret_str + "\n5% out of sat :\t" + "%.3f" %  (self.m5.min)  + "\t" + "%.3f" %  (self.m5.mean)  + "\t" + "%.3f" %  (self.m5.max)  + "\t" + "%.3f" %  (self.m5.stdev)
+         ret_str = ret_str + "\n10% out of sat:\t" + "%.3f" %  (self.m10.min) + "\t" + "%.3f" %  (self.m10.mean) + "\t" + "%.3f" %  (self.m10.max) + "\t" + "%.3f" %  (self.m10.stdev)
      return ret_str
 
      
@@ -204,7 +210,20 @@ vsum_jit_amplt = all_statistics_calc(scan_result.getReadable(8))
 vsum_jit_phase = all_statistics_calc(scan_result.getReadable(9))
 
 
-
+# adjust plot ranges for the color bar
+(p0,p1,p2,p3,p4,p5,p6,p7,p8,p9) = get_plots(”Jitter Scan”)             # get_plots() return the plots of the default panel
+a_lim = 2.0e-4
+p_lim = 0.004
+p0.setColorScale(0.0, a_lim)
+p1.setColorScale(0.0, p_lim)
+p2.setColorScale(0.0, a_lim)
+p3.setColorScale(0.0, p_lim)
+p4.setColorScale(0.0, a_lim)
+p5.setColorScale(0.0, p_lim)
+p6.setColorScale(0.0, a_lim)
+p7.setColorScale(0.0, p_lim)
+p8.setColorScale(0.0, a_lim)
+p9.setColorScale(0.0, p_lim)
 
 
 #Setting attributes to the scan group
@@ -233,20 +252,20 @@ set_attribute(path, "HVPS", hvps)
 
 msg = "HVPS command = " + str(hvps)[:6] + "V)"
 msg = msg + "\n---------------------------------------------------------"
-msg = msg + "\nREF amplt\n" + ref_jit_amplt.return_string()
-msg = msg + "\n\nREF phase\n" + ref_jit_phase.return_string()
+msg = msg + "\nREF amplt\n" + ref_jit_amplt.return_string(1)
+msg = msg + "\n\nREF phase\n" + ref_jit_phase.return_string(0)
 msg = msg + "\n---------------------------------------------------------"
-msg = msg + "\nIQM amplt\n" + iqm_jit_amplt.return_string()
-msg = msg + "\n\nIQM phase\n" + iqm_jit_phase.return_string()
+msg = msg + "\nIQM amplt\n" + iqm_jit_amplt.return_string(1)
+msg = msg + "\n\nIQM phase\n" + iqm_jit_phase.return_string(0)
 msg = msg + "\n---------------------------------------------------------"
-msg = msg + "\nPRE amplt\n" + pre_jit_amplt.return_string()
-msg = msg + "\n\nPRE phase\n" + pre_jit_phase.return_string()
+msg = msg + "\nPRE amplt\n" + pre_jit_amplt.return_string(1)
+msg = msg + "\n\nPRE phase\n" + pre_jit_phase.return_string(0)
 msg = msg + "\n---------------------------------------------------------"
-msg = msg + "\nKLY amplt\n" + kly_jit_amplt.return_string()
-msg = msg + "\n\nKLY phase\n" + kly_jit_phase.return_string()
+msg = msg + "\nKLY amplt\n" + kly_jit_amplt.return_string(1)
+msg = msg + "\n\nKLY phase\n" + kly_jit_phase.return_string(0)
 msg = msg + "\n---------------------------------------------------------"
-msg = msg + "\nVSUM amplt\n" + vsum_jit_amplt.return_string()
-msg = msg + "\n\nVSUM phase\n" + vsum_jit_phase.return_string()
+msg = msg + "\nVSUM amplt\n" + vsum_jit_amplt.return_string(1)
+msg = msg + "\n\nVSUM phase\n" + vsum_jit_phase.return_string(0)
 msg = msg + "\n---------------------------------------------------------"
 msg = msg + "\n\nPlots axis: X-axis = Three amplitude setpoints from left to right: 10% below sat / 5% below sat / sat"
 msg = msg + "\nY-axis = Readback phase (with feedback closed = setpoint phase)"