From ccbe99dfcead7d18f3c94ce80380f35ebbb6b252 2021-10-27 21:38:56
From: avaldez <avaldez@igp.gob.pe>
Date: 2021-10-27 21:38:56
Subject: [PATCH] LAST CHANGE WR

---

diff --git a/schainpy/model/data/jrodata.py b/schainpy/model/data/jrodata.py
index 7d7bec8..c1d2bb8 100644
--- a/schainpy/model/data/jrodata.py
+++ b/schainpy/model/data/jrodata.py
@@ -575,7 +575,6 @@ class Spectra(JROData):
         z = self.data_spc / factor
         z = numpy.where(numpy.isfinite(z), z, numpy.NAN)
         avg = numpy.average(z, axis=1)
-
         return 10 * numpy.log10(avg)
 
     def getCoherence(self, pairsList=None, phase=False):
diff --git a/schainpy/model/graphics/jroplot_parameters.py b/schainpy/model/graphics/jroplot_parameters.py
index 78e4c82..3f563a1 100644
--- a/schainpy/model/graphics/jroplot_parameters.py
+++ b/schainpy/model/graphics/jroplot_parameters.py
@@ -95,11 +95,9 @@ class PowerPlot(RTIPlot):
     colormap = 'jet'
 
     def update(self, dataOut):
-
         data = {
             'pow': 10*numpy.log10(dataOut.data_pow/dataOut.normFactor)
         }
-
         return data, {}
 
 class SpectralWidthPlot(RTIPlot):
@@ -204,9 +202,7 @@ class GenericRTIPlot(Plot):
         self.x = self.data.times
         self.y = self.data.yrange
         self.z = self.data['param']
-
         self.z = 10*numpy.log10(self.z)
-
         self.z = numpy.ma.masked_invalid(self.z)
 
         if self.decimation is None:
@@ -399,17 +395,11 @@ class WeatherPlot(Plot):
 
         data = {}
         meta = {}
-        print("aprox",dataOut.data_360[0])
         data['weather'] = 10*numpy.log10(dataOut.data_360[0]/(250.0))
-        #print(data['weather'])
         data['azi']     = dataOut.data_azi
-        print("UPDATE",data['azi'])
         return data, meta
 
     def const_ploteo(self,data_weather,data_azi,step,res):
-        #print("data_weather",data_weather)
-        print("data_azi",data_azi)
-        print("step",step)
         if self.ini==0:
             #------- AZIMUTH
             n     = (360/res)-len(data_azi)
diff --git a/schainpy/model/graphics/jroplot_spectra.py b/schainpy/model/graphics/jroplot_spectra.py
index 92542ce..b69bb3e 100644
--- a/schainpy/model/graphics/jroplot_spectra.py
+++ b/schainpy/model/graphics/jroplot_spectra.py
@@ -226,7 +226,6 @@ class RTIPlot(Plot):
     def setup(self):
         self.xaxis = 'time'
         self.ncols = 1
-        print("ch",self.data.channels)
         self.nrows = len(self.data.channels)
         self.nplots = len(self.data.channels)
         self.ylabel = 'Range [km]'
@@ -267,7 +266,6 @@ class RTIPlot(Plot):
                                        cmap=plt.get_cmap(self.colormap)
                                        )
                 if self.showprofile:
-                    print("test-------------------------------------1")
                     ax.plot_profile = self.pf_axes[n].plot(
                         data['rti'][n], self.y)[0]
                     ax.plot_noise = self.pf_axes[n].plot(numpy.repeat(data['noise'][n], len(self.y)), self.y,
diff --git a/schainpy/model/proc/jroproc_parameters.py b/schainpy/model/proc/jroproc_parameters.py
index 1b98f89..e6aaa10 100755
--- a/schainpy/model/proc/jroproc_parameters.py
+++ b/schainpy/model/proc/jroproc_parameters.py
@@ -1259,7 +1259,6 @@ class SpectralMoments(Operation):
         dataOut.data_pow = data_param[:,1]
         dataOut.data_dop = data_param[:,2]
         dataOut.data_width = data_param[:,3]
-
         return dataOut
 
     def __calculateMoments(self, oldspec, oldfreq, n0,
diff --git a/schainpy/model/proc/jroproc_voltage.py b/schainpy/model/proc/jroproc_voltage.py
index 5fb2314..87071df 100644
--- a/schainpy/model/proc/jroproc_voltage.py
+++ b/schainpy/model/proc/jroproc_voltage.py
@@ -1478,6 +1478,7 @@ class PulsePair(Operation):
         if self.__dataReady:
             dataOut.nCohInt        *= self.n
             dataOut.dataPP_POW      = data_intensity # S
+            print("help",data_power)
             dataOut.dataPP_POWER    = data_power     # P
             dataOut.dataPP_DOP      = data_velocity
             dataOut.dataPP_SNR      = data_snrPP
diff --git a/schainpy/scripts/USRP_ADQ_SPECM.py b/schainpy/scripts/USRP_ADQ_SPECM.py
index 01d2a23..5909e2f 100644
--- a/schainpy/scripts/USRP_ADQ_SPECM.py
+++ b/schainpy/scripts/USRP_ADQ_SPECM.py
@@ -32,12 +32,14 @@ controllerObj.setup(id = '191', name='Test_USRP', description=desc)
 ###path_pp = '/DATA_RM/TEST_HDF5_SPEC'
 
 #path = '/DATA_RM/USRP_22'
-path = '/DATA_RM/23/6v'
+###path = '/DATA_RM/23/6v'
+path = '/DATA_RM/TEST_19OCTUBRE/10MHZ'
 #path_pp = '/DATA_RM/TEST_HDF5'
+path_pp = '/DATA_RM/TEST_HDF5_19OCT' 
 # UTIMO TEST 22 DE SEPTIEMBRE
 #path_pp = '/DATA_RM/TEST_HDF5_SPEC_22'
 #path_pp = '/DATA_RM/TEST_HDF5_SPEC_3v'
-path_pp = '/DATA_RM/TEST_HDF5_SPEC_23/6v'
+###path_pp = '/DATA_RM/TEST_HDF5_SPEC_23/6v'
 
 
 #remotefolder = "/home/wmaster/graficos"
@@ -79,7 +81,7 @@ opObj11 = readUnitConfObj.addOperation(name='printInfo')
 #######################################################################
 
 
-V=6
+V=10
 IPP=400*1e-6
 n= int(1/(V*IPP))
 print("n numero de Perfiles a procesar con nFFTPoints ", n)
diff --git a/schainpy/scripts/basic_proc00001.py b/schainpy/scripts/basic_proc00001.py
index e0f478e..d8cad68 100644
--- a/schainpy/scripts/basic_proc00001.py
+++ b/schainpy/scripts/basic_proc00001.py
@@ -25,12 +25,15 @@ print("SETUP- RADAR METEOROLOGICO")
 V       = 10
 mode    = 1
 #path    = '/DATA_RM/23/6v'
-path    = '/DATA_RM/TEST_INTEGRACION_2M'
-path_ped='/DATA_RM/TEST_PEDESTAL/P20211012-082745'
+####path    = '/DATA_RM/TEST_INTEGRACION_2M'
+#path    = '/DATA_RM/TEST_19OCTUBRE/10MHZ'
+path    = '/DATA_RM/WR_20_OCT'
+#### path_ped='/DATA_RM/TEST_PEDESTAL/P20211012-082745'
+#### path_ped='/DATA_RM/TEST_PEDESTAL/P20211019-192244'
 figpath_pp  = "/home/soporte/Pictures/TEST_PP"
-figpath_mom = "/home/soporte/Pictures/TEST_MOM"
-plot        = 0
-integration = 1
+figpath_spec = "/home/soporte/Pictures/TEST_MOM"
+plot        = 1
+integration = 0
 save        = 0
 if save == 1:
     if mode==0:
@@ -70,11 +73,11 @@ time.sleep(15)
 #######################################################################
 ################# RANGO DE PLOTEO######################################
 dBmin = '1'
-dBmax = '85'
-xmin = '15'
-xmax = '15.25'
+dBmax = '65'
+xmin = '13.2'
+xmax = '13.5'
 ymin = '0'
-ymax = '600'
+ymax = '60'
 #######################################################################
 ########################FECHA##########################################
 str = datetime.date.today()
@@ -116,7 +119,7 @@ if mode ==0:
     procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
     if plot==1:
         opObj11 = procUnitConfObjB.addOperation(name='GenericRTIPlot',optype='external')
-        opObj11.addParameter(name='attr_data', value='dataPP_POW')
+        opObj11.addParameter(name='attr_data', value='dataPP_POWER')
         opObj11.addParameter(name='colormap', value='jet')
         opObj11.addParameter(name='xmin', value=xmin)
         opObj11.addParameter(name='xmax', value=xmax)
@@ -124,7 +127,7 @@ if mode ==0:
         opObj11.addParameter(name='zmax', value=dBmax)
         opObj11.addParameter(name='save', value=figpath_pp)
         opObj11.addParameter(name='showprofile', value=0)
-        opObj11.addParameter(name='save_period', value=50)
+        opObj11.addParameter(name='save_period', value=10)
 
     ####################### METODO ESCRITURA #######################################################################
     if save==1:
@@ -133,7 +136,7 @@ if mode ==0:
         #opObj10.addParameter(name='mode',value=0)
         opObj10.addParameter(name='blocksPerFile',value='100',format='int')
         opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list')
-        opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'
+        opObj10.addParameter(name='dataList',value='dataPP_POWER,dataPP_DOP,utctime',format='list')#,format='list'
     if integration==1:
         V=10
         blocksPerfile=360
@@ -175,9 +178,9 @@ else:
         opObj11.addParameter(name='xmax', value=xmax)
         opObj11.addParameter(name='zmin', value=dBmin)
         opObj11.addParameter(name='zmax', value=dBmax)
-        opObj11.addParameter(name='save', value=figpath_mom)
+        opObj11.addParameter(name='save', value=figpath_spec)
         opObj11.addParameter(name='showprofile', value=0)
-        opObj11.addParameter(name='save_period', value=100)
+        opObj11.addParameter(name='save_period', value=10)
 
     if save==1:
         opObj10 = procUnitConfObjC.addOperation(name='HDFWriter')
diff --git a/schainpy/scripts/basic_proc00002.py b/schainpy/scripts/basic_proc00002.py
new file mode 100644
index 0000000..4fa5ffb
--- /dev/null
+++ b/schainpy/scripts/basic_proc00002.py
@@ -0,0 +1,217 @@
+# Ing. AVP
+# 06/10/2021
+# ARCHIVO DE LECTURA
+import os, sys
+import datetime
+import time
+from schainpy.controller import Project
+#### NOTA###########################################
+# INPUT :
+# VELOCIDAD PARAMETRO :        V = 2°/seg
+# MODO PULSE PAIR O MOMENTOS:  0 : Pulse Pair ,1 : Momentos
+######################################################
+##### PROCESAMIENTO ##################################
+#####  OJO TENER EN CUENTA EL n= para el Pulse Pair ##
+#####                   O EL n= nFFTPoints         ###
+######################################################
+######## BUSCAMOS EL numero de IPP equivalente 1°#####
+######## Sea V la velocidad del Pedestal en °/seg#####
+######## 1° sera Recorrido en un tiempo de  1/V ######
+######## IPP del Radar 400 useg --> 60 Km ############
+######## n   = 1/(V(°/seg)*IPP(Km)) , NUMERO DE IPP ##
+######## n   = 1/(V*IPP) #############################
+######## VELOCIDAD DEL PEDESTAL ######################
+print("SETUP- RADAR METEOROLOGICO")
+V       = 10
+mode    = 1
+#path    = '/DATA_RM/23/6v'
+#path    = '/DATA_RM/TEST_INTEGRACION_2M'
+path    = '/DATA_RM/WR_20_OCT'
+
+#path_ped='/DATA_RM/TEST_PEDESTAL/P20211012-082745'
+path_ped='/DATA_RM/TEST_PEDESTAL/P20211020-131248'
+
+figpath_pp  = "/home/soporte/Pictures/TEST_PP"
+figpath_mom = "/home/soporte/Pictures/TEST_MOM"
+plot        = 0
+integration = 1
+save        = 0
+if save == 1:
+    if mode==0:
+        path_save = '/DATA_RM/TEST_HDF5_PP_23/6v'
+        path_save = '/DATA_RM/TEST_HDF5_PP'
+        path_save = '/DATA_RM/TEST_HDF5_PP_100'
+    else:
+        path_save = '/DATA_RM/TEST_HDF5_SPEC_23_V2/6v'
+
+print("* PATH data ADQ        :", path)
+print("* Velocidad Pedestal   :",V,"°/seg")
+############################ NRO Perfiles PROCESAMIENTO ###################
+V=V
+IPP=400*1e-6
+n= int(1/(V*IPP))
+print("* n - NRO Perfiles Proc:", n )
+################################## MODE ###################################
+print("* Modo de Operacion    :",mode)
+if mode ==0:
+    print("* Met. Seleccionado    : Pulse Pair")
+else:
+    print("* Met. Momentos        : Momentos")
+
+################################## MODE ###################################
+print("* Grabado de datos     :",save)
+if save ==1:
+    if mode==0:
+        ope= "Pulse Pair"
+    else:
+        ope= "Momentos"
+    print("* Path-Save Data -", ope , path_save)
+
+print("* Integracion de datos :",integration)
+
+time.sleep(15)
+#remotefolder = "/home/wmaster/graficos"
+#######################################################################
+################# RANGO DE PLOTEO######################################
+dBmin = '1'
+dBmax = '85'
+xmin = '15'
+xmax = '15.25'
+ymin = '0'
+ymax = '600'
+#######################################################################
+########################FECHA##########################################
+str = datetime.date.today()
+today = str.strftime("%Y/%m/%d")
+str2 = str - datetime.timedelta(days=1)
+yesterday = str2.strftime("%Y/%m/%d")
+#######################################################################
+########################SIGNAL CHAIN ##################################
+#######################################################################
+desc = "USRP_test"
+filename = "USRP_processing.xml"
+controllerObj = Project()
+controllerObj.setup(id = '191', name='Test_USRP', description=desc)
+#######################################################################
+######################## UNIDAD DE LECTURA#############################
+#######################################################################
+readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
+                                            path=path,
+                                            startDate="2021/01/01",#today,
+                                            endDate="2021/12/30",#today,
+                                            startTime='00:00:00',
+                                            endTime='23:59:59',
+                                            delay=0,
+                                            #set=0,
+                                            online=0,
+                                            walk=1,
+                                            ippKm = 60)
+
+opObj11 = readUnitConfObj.addOperation(name='printInfo')
+
+procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
+
+if mode ==0:
+    ####################### METODO PULSE PAIR ######################################################################
+    opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
+    opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO  QUE LA VELOCIDAD ES 10 GRADOS
+    #opObj11.addParameter(name='removeDC', value=1, format='int')
+    ####################### METODO Parametros ######################################################################
+    procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
+    if plot==1:
+        opObj11 = procUnitConfObjB.addOperation(name='GenericRTIPlot',optype='external')
+        opObj11.addParameter(name='attr_data', value='dataPP_POW')
+        opObj11.addParameter(name='colormap', value='jet')
+        opObj11.addParameter(name='xmin', value=xmin)
+        opObj11.addParameter(name='xmax', value=xmax)
+        opObj11.addParameter(name='zmin', value=dBmin)
+        opObj11.addParameter(name='zmax', value=dBmax)
+        opObj11.addParameter(name='save', value=figpath_pp)
+        opObj11.addParameter(name='showprofile', value=0)
+        opObj11.addParameter(name='save_period', value=50)
+
+    ####################### METODO ESCRITURA #######################################################################
+    if save==1:
+        opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')
+        opObj10.addParameter(name='path',value=path_save)
+        #opObj10.addParameter(name='mode',value=0)
+        opObj10.addParameter(name='blocksPerFile',value='100',format='int')
+        opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list')
+        opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'
+    if integration==1:
+        V=10
+        blocksPerfile=360
+        print("* Velocidad del Pedestal:",V)
+        tmp_blocksPerfile = 100
+        f_a_p= int(tmp_blocksPerfile/V)
+
+        opObj11 = procUnitConfObjB.addOperation(name='PedestalInformation')
+        opObj11.addParameter(name='path_ped', value=path_ped)
+        #opObj11.addParameter(name='path_adq', value=path_adq)
+        opObj11.addParameter(name='t_Interval_p', value='0.01', format='float')
+        opObj11.addParameter(name='blocksPerfile', value=blocksPerfile, format='int')
+        opObj11.addParameter(name='n_Muestras_p', value='100', format='float')
+        opObj11.addParameter(name='f_a_p', value=f_a_p, format='int')
+        opObj11.addParameter(name='online', value='0', format='int')
+
+        opObj11 = procUnitConfObjB.addOperation(name='Block360')
+        opObj11.addParameter(name='n', value='10', format='int')
+        opObj11.addParameter(name='mode', value=mode, format='int')
+
+        # este bloque funciona bien con divisores de 360 no olvidar 0 10 20 30 40 60 90 120 180
+
+        opObj11= procUnitConfObjB.addOperation(name='WeatherPlot',optype='other')
+
+
+else:
+    ####################### METODO SPECTROS ######################################################################
+    procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
+    procUnitConfObjB.addParameter(name='nFFTPoints', value=n, format='int')
+    procUnitConfObjB.addParameter(name='nProfiles' , value=n, format='int')
+
+    procUnitConfObjC = controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjB.getId())
+    procUnitConfObjC.addOperation(name='SpectralMoments')
+    if plot==1:
+        dBmin = '1'
+        dBmax = '65'
+        opObj11 = procUnitConfObjC.addOperation(name='PowerPlot',optype='external')
+        opObj11.addParameter(name='xmin', value=xmin)
+        opObj11.addParameter(name='xmax', value=xmax)
+        opObj11.addParameter(name='zmin', value=dBmin)
+        opObj11.addParameter(name='zmax', value=dBmax)
+        opObj11.addParameter(name='save', value=figpath_mom)
+        opObj11.addParameter(name='showprofile', value=0)
+        opObj11.addParameter(name='save_period', value=100)
+
+    if save==1:
+        opObj10 = procUnitConfObjC.addOperation(name='HDFWriter')
+        opObj10.addParameter(name='path',value=path_save)
+        #opObj10.addParameter(name='mode',value=0)
+        opObj10.addParameter(name='blocksPerFile',value='360',format='int')
+        #opObj10.addParameter(name='metadataList',value='utctimeInit,heightList,nIncohInt,nCohInt,nProfiles,channelList',format='list')#profileIndex
+        opObj10.addParameter(name='metadataList',value='utctimeInit,heightList,nIncohInt,nCohInt,nProfiles,channelList',format='list')#profileIndex
+        opObj10.addParameter(name='dataList',value='data_pow,data_dop,utctime',format='list')#,format='list'
+
+    if integration==1:
+       V=10
+       blocksPerfile=360
+       print("* Velocidad del Pedestal:",V)
+       tmp_blocksPerfile = 100
+       f_a_p= int(tmp_blocksPerfile/V)
+
+       opObj11 = procUnitConfObjC.addOperation(name='PedestalInformation')
+       opObj11.addParameter(name='path_ped', value=path_ped)
+       #opObj11.addParameter(name='path_adq', value=path_adq)
+       opObj11.addParameter(name='t_Interval_p', value='0.01', format='float')
+       opObj11.addParameter(name='blocksPerfile', value=blocksPerfile, format='int')
+       opObj11.addParameter(name='n_Muestras_p', value='100', format='float')
+       opObj11.addParameter(name='f_a_p', value=f_a_p, format='int')
+       opObj11.addParameter(name='online', value='0', format='int')
+
+       opObj11 = procUnitConfObjC.addOperation(name='Block360')
+       opObj11.addParameter(name='n', value='30', format='int')
+       opObj11.addParameter(name='mode', value=mode, format='int')
+
+       # este bloque funciona bien con divisores de 360 no olvidar 0 10 20 30 40 60 90 120 180
+       opObj11= procUnitConfObjC.addOperation(name='WeatherPlot',optype='other')
+controllerObj.start()
diff --git a/schainpy/scripts/check_pedestal.py b/schainpy/scripts/check_pedestal.py
index 17a6feb..25fc88a 100644
--- a/schainpy/scripts/check_pedestal.py
+++ b/schainpy/scripts/check_pedestal.py
@@ -6,7 +6,13 @@ import os,numpy,h5py
 import sys,time
 import matplotlib.pyplot as plt
 ####################################################
-path_ped = '/DATA_RM/TEST_PEDESTAL/P20211012-082745'
+#################################################################
+# LA FECHA 21-10-20 CORRESPONDE A LAS PRUEBAS DEL DIA MIERCOLES
+# 1:15:51 pm hasta 3:49:32 pm
+#################################################################
+
+#path_ped = '/DATA_RM/TEST_PEDESTAL/P20211012-082745'
+path_ped = '/DATA_RM/TEST_PEDESTAL/P20211020-131248'
 # Metodo para verificar numero
 def isNumber(str):
     try:
@@ -34,6 +40,7 @@ t=  0
 time_        = numpy.zeros([m])
 # creacion de
 for i in range(m):
+    print("order:",i)
     tmp_azi_pos = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="azi_pos")
     tmp_ele_pos = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="ele_pos")
     tmp_azi_vel = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="azi_vel")
diff --git a/schainpy/scripts/pp_spec_comparation_pow.txt b/schainpy/scripts/pp_spec_comparation_pow.txt
new file mode 100644
index 0000000..3c0e1f0
--- /dev/null
+++ b/schainpy/scripts/pp_spec_comparation_pow.txt
@@ -0,0 +1,13 @@
+
+PULSE PAIR
+help [[1.364 1.376 1.28  ... 1.352 1.372 1.332]
+ [2.912 3.012 3.06  ... 3.056 2.6   2.936]]
+ help [[1.256 1.304 1.212 ... 1.228 1.328 1.528]
+ [2.744 2.604 2.492 ... 2.544 2.916 2.644]]
+ help [[1.176 1.248 1.48  ... 1.388 1.396 1.216]
+  [2.564 2.524 2.756 ... 2.772 2.7   2.684]]
+
+  help [[1.312 1.328 1.28  ... 1.28  1.444 1.532]
+ [3.    2.964 3.092 ... 3.104 3.016 2.984]]
+
+MOMENTOS
diff --git a/schainpy/scripts/pulsepair_test0001.py b/schainpy/scripts/pulsepair_test0001.py
new file mode 100644
index 0000000..101c68a
--- /dev/null
+++ b/schainpy/scripts/pulsepair_test0001.py
@@ -0,0 +1,69 @@
+# Ing. AVP
+# 06/10/2021
+# ARCHIVO DE LECTURA
+import os, sys
+import datetime
+import time
+from schainpy.controller import Project
+
+print("SETUP- RADAR METEOROLOGICO")
+V       = 10
+#######################################################################
+################# RANGO DE PLOTEO######################################
+dBmin = '1'
+dBmax = '65'
+xmin = '13.2'
+xmax = '13.5'
+ymin = '0'
+ymax = '60'
+
+path    = '/DATA_RM/WR_20_OCT'
+figpath_pp  = "/home/soporte/Pictures/TEST_PP"
+
+
+IPP=400*1e-6
+n= int(1/(V*IPP))
+#n=250
+print("* n - NRO Perfiles Proc:", n )
+
+desc = "USRP_test"
+filename = "USRP_processing.xml"
+controllerObj = Project()
+controllerObj.setup(id = '191', name='Test_USRP', description=desc)
+#######################################################################
+######################## UNIDAD DE LECTURA#############################
+#######################################################################
+readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
+                                            path=path,
+                                            startDate="2021/01/01",#today,
+                                            endDate="2021/12/30",#today,
+                                            startTime='00:00:00',
+                                            endTime='23:59:59',
+                                            delay=0,
+                                            #set=0,
+                                            online=0,
+                                            walk=1,
+                                            ippKm = 60)
+
+opObj11 = readUnitConfObj.addOperation(name='printInfo')
+
+procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
+
+opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
+opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO  QUE LA VELOCIDAD ES 10 GRADOS
+#opObj11.addParameter(name='removeDC', value=1, format='int')
+
+procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
+
+opObj11 = procUnitConfObjB.addOperation(name='GenericRTIPlot',optype='external')
+opObj11.addParameter(name='attr_data', value='dataPP_POWER')
+opObj11.addParameter(name='colormap', value='jet')
+opObj11.addParameter(name='xmin', value=xmin)
+opObj11.addParameter(name='xmax', value=xmax)
+opObj11.addParameter(name='zmin', value=dBmin)
+opObj11.addParameter(name='zmax', value=dBmax)
+opObj11.addParameter(name='save', value=figpath_pp)
+opObj11.addParameter(name='showprofile', value=0)
+#opObj11.addParameter(name='save_period', value=10)
+
+controllerObj.start()
diff --git a/schainpy/scripts/specmom_test0001.py b/schainpy/scripts/specmom_test0001.py
new file mode 100644
index 0000000..e8c73ff
--- /dev/null
+++ b/schainpy/scripts/specmom_test0001.py
@@ -0,0 +1,80 @@
+# Ing. AVP
+# 06/10/2021
+# ARCHIVO DE LECTURA
+import os, sys
+import datetime
+import time
+from schainpy.controller import Project
+
+print("SETUP- RADAR METEOROLOGICO")
+V       = 10
+#######################################################################
+################# RANGO DE PLOTEO######################################
+dBmin = '1'
+dBmax = '65'
+xmin = '13.2'
+xmax = '13.5'
+ymin = '0'
+ymax = '60'
+
+path    = '/DATA_RM/WR_20_OCT'
+figpath_spec = "/home/soporte/Pictures/TEST_MOM"
+
+
+IPP=400*1e-6
+n= int(1/(V*IPP))
+print("* n - NRO Perfiles Proc:", n )
+time.sleep(5)
+desc = "USRP_test"
+filename = "USRP_processing.xml"
+controllerObj = Project()
+controllerObj.setup(id = '191', name='Test_USRP', description=desc)
+#######################################################################
+######################## UNIDAD DE LECTURA#############################
+#######################################################################
+readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
+                                            path=path,
+                                            startDate="2021/01/01",#today,
+                                            endDate="2021/12/30",#today,
+                                            startTime='00:00:00',
+                                            endTime='23:59:59',
+                                            delay=0,
+                                            #set=0,
+                                            online=0,
+                                            walk=1,
+                                            ippKm = 60)
+
+opObj11 = readUnitConfObj.addOperation(name='printInfo')
+
+procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
+
+procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
+procUnitConfObjB.addParameter(name='nFFTPoints', value=n, format='int')
+procUnitConfObjB.addParameter(name='nProfiles' , value=n, format='int')
+'''
+opObj11 = procUnitConfObjB.addOperation(name='RTIPlot', optype='external')
+#.addParameter(name='id', value='2', format='int')
+opObj11.addParameter(name='wintitle', value='RTIPlot', format='str')
+opObj11.addParameter(name='xmin', value=xmin)
+opObj11.addParameter(name='xmax', value=xmax)
+opObj11.addParameter(name='zmin', value=dBmin, format='int')
+opObj11.addParameter(name='zmax', value=dBmax, format='int')
+'''
+#opObj13 = procUnitConfObjB.addOperation(name='removeDC')
+#opObj13.addParameter(name='mode', value='2', format='int')
+
+procUnitConfObjC = controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjB.getId())
+procUnitConfObjC.addOperation(name='SpectralMoments')
+
+dBmin = '1'
+dBmax = '65'
+opObj11 = procUnitConfObjC.addOperation(name='PowerPlot',optype='external')
+opObj11.addParameter(name='xmin', value=xmin)
+opObj11.addParameter(name='xmax', value=xmax)
+opObj11.addParameter(name='zmin', value=dBmin)
+opObj11.addParameter(name='zmax', value=dBmax)
+opObj11.addParameter(name='save', value=figpath_spec)
+opObj11.addParameter(name='showprofile', value=0)
+#opObj11.addParameter(name='save_period', value=10)
+
+controllerObj.start()
diff --git a/schainpy/scripts/test_sim00011b.py b/schainpy/scripts/test_sim00011b.py
index 2794af6..927f1a9 100644
--- a/schainpy/scripts/test_sim00011b.py
+++ b/schainpy/scripts/test_sim00011b.py
@@ -28,7 +28,8 @@ ymax = '600'
 #######################################################################
 #######################################################################
 #path       = '/DATA_RM/TEST_HDF5_SPEC'
-path       = '/DATA_RM/TEST_HDF5_SPEC_23/6v/'
+#path       = '/DATA_RM/TEST_HDF5_SPEC_23/6v/'
+path       = '/DATA_RM/TEST_HDF5_19OCT'
 figpath    = '/home/soporte/Downloads/23/6v'
 desc       = "Simulator Test"
 desc_data  =  {