diff --git a/pysqm/plot.py b/pysqm/plot.py
index f411e28c9b31f3d3eca4fe86831db40f233e3289..379ee6dfbb97a5e2a287587326efae7ad4f16bb7 100644
--- a/pysqm/plot.py
+++ b/pysqm/plot.py
@@ -449,10 +449,10 @@ class Plot(object):
 
 
         self.thegraph_sunalt.set_title(\
-         'Sky Brightness ('+config._device_shorttype+'-'+\
+         'Sky Brightness ('+config._device_shorttype+'_'+\
          config._observatory_name+')\n',fontsize='x-large')
         self.thegraph_sunalt.set_xlabel('Solar altitude (deg)',fontsize='large')
-        self.thegraph_sunalt.set_ylabel('Sky Brightness (mag/arcsec2)',fontsize='medium')
+        self.thegraph_sunalt.set_ylabel('Sky Brightness (mag/$\mathregular{arcsec^2}$)',fontsize='medium')
 
         # Auxiliary plot (Temperature)
         '''
@@ -495,7 +495,7 @@ class Plot(object):
         #self.thegraph_time.set_title('Sky Brightness (SQM-'+config._observatory_name+')',\
         # fontsize='x-large')
         self.thegraph_time.set_xlabel('Time (UTC'+UTC_offset_label+')',fontsize='large')
-        self.thegraph_time.set_ylabel('Sky Brightness (mag/arcsec2)',fontsize='medium')
+        self.thegraph_time.set_ylabel('Sky Brightness (mag/$\mathregular{arcsec^2}$)',fontsize='medium')
 
         # Auxiliary plot (Temperature)
         '''
@@ -596,7 +596,7 @@ class Plot(object):
         aftermidnight_label = str(Data.aftermidnight.label_dates).replace('[','').replace(']','')
 
         self.thegraph_sunalt.text(0.00,1.015,\
-         config._device_shorttype+'-'+config._observatory_name+' '*5+'Serial #'+str(Data.serial_number),\
+         config._device_shorttype+'_'+config._observatory_name+' '*5+'Serial #'+str(Data.serial_number),\
          color='0.25',fontsize='small',fontname='monospace',\
          transform = self.thegraph_sunalt.transAxes)
 
@@ -677,7 +677,7 @@ class Plot(object):
         aftermidnight_label = str(Data.aftermidnight.label_dates).replace('[','').replace(']','')
 
         self.thegraph_time.text(0.00,1.015,\
-         config._device_shorttype+'-'+config._observatory_name+' '*5+'Serial #'+str(Data.serial_number),\
+         config._device_shorttype+'_'+config._observatory_name+' '*5+'Serial #'+str(Data.serial_number),\
          color='0.25',fontsize='small',fontname='monospace',\
          transform = self.thegraph_time.transAxes)
 
@@ -711,10 +711,10 @@ def save_stats_to_file(Night,NSBData,Ephem):
      '# Col 1: Date\n'+\
      '# Col 2: Total measures\n'+\
      '# Col 3: Number of Best NSB measures\n'+\
-     '# Col 4: Median of best N NSBs (mag/arcsec2)\n'+\
-     '# Col 5: Err in the median of best N NSBs (mag/arcsec2)\n'+\
+     '# Col 4: Median of best N NSBs (mag/arcsec^2)\n'+\
+     '# Col 5: Err in the median of best N NSBs (mag/arcsec^2)\n'+\
      '# Col 6: Window size for the smoothing function\n'+\
-     '# Col 7: Mean of Abs diff of NSBs data - fourier model (mag/arcsec2)\n'+\
+     '# Col 7: Mean of Abs diff of NSBs data - fourier model (mag/arcsec^2)\n'+\
      '# Col 8: Min Temp (C) between astronomical twilights\n'+\
      '# Col 9: Max Temp (C) between astronomical twilights\n\n'
      #'# Col 6: Number of terms of the low-freq fourier model\n'+\