diff --git a/schainpy/model/jrodata.py b/schainpy/model/jrodata.py
index f2334f3..c8c4530 100644
--- a/schainpy/model/jrodata.py
+++ b/schainpy/model/jrodata.py
@@ -479,7 +479,7 @@ class Spectra(JROData):
         pwcode = 1
         if self.flagDecodeData:
             pwcode = numpy.sum(self.code[0]**2)
-        normFactor = min(self.nFFTPoints,self.nProfiles)*self.nIncohInt*self.nCohInt*self.windowOfFilter*pwcode
+        normFactor = min(self.nFFTPoints,self.nProfiles)*self.nIncohInt*self.nCohInt*pwcode
         
         return normFactor
     
diff --git a/schainpy/model/jrodataIO.py b/schainpy/model/jrodataIO.py
index 3ac33c5..ff7b3d6 100644
--- a/schainpy/model/jrodataIO.py
+++ b/schainpy/model/jrodataIO.py
@@ -1032,7 +1032,7 @@ class JRODataWriter(JRODataIO, Operation):
         if self.fp != None:
             self.fp.close()
         
-        timeTuple = time.localtime( self.dataOut.dataUtcTime)
+        timeTuple = time.localtime( self.dataOut.utctime)
         subfolder = 'D%4.4d%3.3d' % (timeTuple.tm_year,timeTuple.tm_yday)
 
         fullpath = os.path.join( path, subfolder )
@@ -1382,7 +1382,7 @@ class VoltageReader(JRODataReader):
             if not( self.readNextBlock() ):
                 return 0
         
-            self.dataOut.dtype = self.dtype
+            self.dataOut.dtype = numpy.dtype([('real','<f8'),('imag','<f8')]) #self.dtype
             
             self.dataOut.nProfiles = self.processingHeaderObj.profilesPerBlock
             
@@ -1966,7 +1966,7 @@ class SpectraReader(JRODataReader):
                 #desplaza a la derecha en el eje 2 determinadas posiciones
                 cspc = numpy.roll( cspc, shift, axis=2 )
             
-            self.processingHeaderObj.shif_fft = True
+#            self.processingHeaderObj.shif_fft = True
 
         spc = numpy.transpose( spc, (0,2,1) )
         self.data_spc = spc
@@ -2040,7 +2040,7 @@ class SpectraReader(JRODataReader):
     
         self.dataOut.flagNoData = False
 
-        self.dataOut.dtype = self.dtype
+        self.dataOut.dtype = numpy.dtype([('real','<f8'),('imag','<f8')])#self.dtype
 
 #        self.dataOut.nChannels = self.nRdChannels
         
@@ -2280,9 +2280,9 @@ class SpectraWriter(JRODataWriter):
         if self.flagIsNewFile == 0:
             self.getBasicHeader()
         
-        self.data_spc = self.dataOut.data_spc
-        self.data_cspc = self.dataOut.data_cspc
-        self.data_dc = self.dataOut.data_dc
+        self.data_spc = self.dataOut.data_spc.copy()
+        self.data_cspc = self.dataOut.data_cspc.copy()
+        self.data_dc = self.dataOut.data_dc.copy()
         
         # #self.processingHeaderObj.dataBlocksPerFile)
         if self.hasAllDataInBuffer():
diff --git a/schainpy/model/jroheaderIO.py b/schainpy/model/jroheaderIO.py
index 948ef42..2a1b51f 100644
--- a/schainpy/model/jroheaderIO.py
+++ b/schainpy/model/jroheaderIO.py
@@ -89,6 +89,7 @@ class BasicHeader(Header):
         return 1
     
     def write(self, fp):
+        self.utc -= self.__LOCALTIME
         headerTuple = (self.size,self.version,self.dataBlock,self.utc,self.miliSecond,self.timeZone,self.dstFlag,self.errorCount)
         header = numpy.array(headerTuple,self.struct)        
         header.tofile(fp)
@@ -390,7 +391,7 @@ class ProcessingHeader(Header):
             if ((self.processFlags & PROCFLAG.DEFINE_PROCESS_CODE) == PROCFLAG.DEFINE_PROCESS_CODE):
                 self.nCode = int(numpy.fromfile(fp,'<u4',1))
                 self.nBaud = int(numpy.fromfile(fp,'<u4',1))
-                self.code = numpy.fromfile(fp,'<f4',self.nCode*self.nBaud).reshape(self.nBaud,self.nCode)
+                self.code = numpy.fromfile(fp,'<f4',self.nCode*self.nBaud).reshape(self.nCode,self.nBaud)
             
             if ((self.processFlags & PROCFLAG.SHIFT_FFT_DATA) == PROCFLAG.SHIFT_FFT_DATA):
                 self.shif_fft = True
@@ -449,13 +450,14 @@ class ProcessingHeader(Header):
 
             
         if self.processFlags & PROCFLAG.DEFINE_PROCESS_CODE == PROCFLAG.DEFINE_PROCESS_CODE:
-            nCode = self.nCode #Probar con un dato que almacene codigo, hasta el momento no se hizo la prueba
+            nCode = numpy.array([self.nCode], numpy.dtype('u4')) #Probar con un dato que almacene codigo, hasta el momento no se hizo la prueba
             nCode.tofile(fp)
 
-            nBaud = self.nBaud
+            nBaud = numpy.array([self.nBaud], numpy.dtype('u4'))
             nBaud.tofile(fp)
 
-            code = self.code.reshape(nCode*nBaud)
+            code = self.code.reshape(self.nCode*self.nBaud)
+            code = code.astype(numpy.dtype('<f4'))
             code.tofile(fp)
             
         return 1
diff --git a/schainpy/model/jroplot.py b/schainpy/model/jroplot.py
index 5ec9fac..b3a82ec 100644
--- a/schainpy/model/jroplot.py
+++ b/schainpy/model/jroplot.py
@@ -53,7 +53,7 @@ class CrossSpectraPlot(Figure):
                 counter += 1
     
     def run(self, dataOut, idfigure, wintitle="", pairsList=None, showprofile='True',
-            xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None,
+            xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, normalize=True,
             save=False, figpath='./', figfile=None,
             power_cmap='jet', coherence_cmap='jet', phase_cmap='RdBu_r'):
         
@@ -87,11 +87,14 @@ class CrossSpectraPlot(Figure):
         
         if len(pairsIndexList) > 4:
             pairsIndexList = pairsIndexList[0:4]
-        factor = dataOut.normFactor
+        
+        factor = 1
+        if normalize:
+            factor = dataOut.normFactor
         x = dataOut.getVelRange(1)
         y = dataOut.getHeiRange()
         z = dataOut.data_spc[:,:,:]/factor
-#        z = numpy.where(numpy.isfinite(z), z, numpy.NAN)
+
         avg = numpy.average(z, axis=1)
         noise = dataOut.getNoise()/factor
         
@@ -238,7 +241,7 @@ class RTIPlot(Figure):
                 counter += 1
     
     def run(self, dataOut, idfigure, wintitle="", channelList=None, showprofile='True',
-            xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None,
+            xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, normalize=True,
             timerange=None,
             save=False, figpath='./', figfile=None):
         
@@ -272,7 +275,9 @@ class RTIPlot(Figure):
         
         tmin = None
         tmax = None
-        factor = dataOut.normFactor
+        factor = 1
+        if normalize:
+            factor = dataOut.normFactor
         x = dataOut.getTimeRange()
         y = dataOut.getHeiRange()
         
@@ -418,7 +423,7 @@ class SpectraPlot(Figure):
                 counter += 1
     
     def run(self, dataOut, idfigure, wintitle="", channelList=None, showprofile='True',
-            xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None,
+            xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, normalize=True, 
             save=False, figpath='./', figfile=None):
         
         """
@@ -445,7 +450,9 @@ class SpectraPlot(Figure):
                 if channel not in dataOut.channelList:
                     raise ValueError, "Channel %d is not in dataOut.channelList"
                 channelIndexList.append(dataOut.channelList.index(channel))
-        factor = dataOut.normFactor
+        factor = 1
+        if normalize:
+            factor = dataOut.normFactor
         x = dataOut.getVelRange(1)
         y = dataOut.getHeiRange()
         
@@ -651,7 +658,7 @@ class ProfilePlot(Figure):
                 self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan, colspan, 1)
     
     def run(self, dataOut, idfigure, wintitle="", channelList=None,
-            xmin=None, xmax=None, ymin=None, ymax=None,
+            xmin=None, xmax=None, ymin=None, ymax=None, normalize=True,
             save=False, figpath='./', figfile=None):
         
         if channelList == None:
@@ -664,7 +671,9 @@ class ProfilePlot(Figure):
                     raise ValueError, "Channel %d is not in dataOut.channelList"
                 channelIndexList.append(dataOut.channelList.index(channel))
                 
-        factor = dataOut.normFactor
+        factor = 1
+        if normalize:
+            factor = dataOut.normFactor
         y = dataOut.getHeiRange()        
         x = dataOut.data_spc[channelIndexList,:,:]/factor
         x = numpy.where(numpy.isfinite(x), x, numpy.NAN) 
@@ -954,7 +963,7 @@ class RTIfromNoise(Figure):
         
                         
     def run(self, dataOut, idfigure, wintitle="", channelList=None, showprofile='True',
-            xmin=None, xmax=None, ymin=None, ymax=None,
+            xmin=None, xmax=None, ymin=None, ymax=None, normalize=True,
             timerange=None,
             save=False, figpath='./', figfile=None):
         
@@ -975,7 +984,9 @@ class RTIfromNoise(Figure):
         tmax = None
         x = dataOut.getTimeRange()
         y = dataOut.getHeiRange()
-        factor = dataOut.normFactor
+        factor = 1
+        if normalize:
+            factor = dataOut.normFactor
         noise = dataOut.getNoise()/factor
         noisedB = 10*numpy.log10(noise)