schain Commit - r1405:49f004b59e80 · Jicamarca Repository

last 30-16:19

avaldezp -

r1405:49f004b59e80

parent child

Context file:

r1405:49f004b59e80

Collapse all files

schainpy/SOPHY_test0001.py created 644 +41 0

@@ -0,0 +1,41
	1	# Ing. AVP
	2	# 01/11/2021
	3	# ARCHIVO DE LECTURA
	4	import os, sys
	5	import datetime
	6	import time
	7	from schainpy.controller import Project
	8	print("----[Setup]-RadarMeteorologico--------")
	9	Vel = 6
	10	modo_proc = 1
	11	#-----------PATH DE DATOS-----------------------#
	12	path = "/DATA_RM/10"
	13	path_ped = "/DATA_RM/TEST_PEDESTAL/P20211111-173856"
	14	print("----[OPCIONES]------------------------")
	15	op_plot = 0
	16	op_integration = 0
	17	op_save = 0
	18	op_plot_spec = 0
	19
	20	########################SIGNAL CHAIN ##################################
	21	desc = "USRP_test"
	22	filename = "USRP_processing.xml"
	23	controllerObj = Project()
	24	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
	25	######################## UNIDAD DE LECTURA#############################
	26	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
	27	path=path,
	28	startDate="2021/11/11",#today,
	29	endDate="2021/12/30",#today,
	30	startTime='17:39:25',
	31	endTime='23:59:59',
	32	delay=0,
	33	#set=0,
	34	online=0,
	35	walk=1,
	36	ippKm = 60)
	37
	38	opObj11 = readUnitConfObj.addOperation(name='printInfo')
	39
	40	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc',inputId=readUnitConfObj.getId())
	41	controllerObj.start()

schainpy/scripts/IQ_digitalRF.py created 644 +72 0

@@ -0,0 +1,72
	1	import os, sys
	2	import datetime
	3	import time
	4	from schainpy.controller import Project
	5
	6	desc = "USRP_test"
	7	filename = "USRP_processing.xml"
	8	controllerObj = Project()
	9	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
	10
	11	############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############
	12	######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################
	13	path = '/home/alex/Downloads/test_rawdata'
	14	path = '/DATA_RM/WR_POT_09_2'
	15	figpath = '/home/alex/Downloads'
	16	figpath_pp = "/home/soporte/Pictures/TEST_POT"
	17	################# RANGO DE PLOTEO######################################
	18	dBmin = '30'
	19	dBmax = '60'
	20	xmin = '0'
	21	xmax ='24'
	22	ymin = '0'
	23	ymax = '600'
	24	########################FECHA##########################################
	25	str = datetime.date.today()
	26	today = str.strftime("%Y/%m/%d")
	27	str2 = str - datetime.timedelta(days=1)
	28	yesterday = str2.strftime("%Y/%m/%d")
	29	######################## UNIDAD DE LECTURA#############################
	30	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
	31	path=path,
	32	startDate="2021/01/01", #"2020/01/01",#today,
	33	endDate= "2021/12/01", #"2020/12/30",#today,
	34	startTime='00:00:00',
	35	endTime='23:59:59',
	36	delay=0,
	37	#set=0,
	38	online=0,
	39	walk =1,
	40	ippKm = 60 )
	41
	42	opObj11 = readUnitConfObj.addOperation(name='printInfo')
	43	#opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
	44	#######################################################################
	45	################ OPERACIONES DOMINIO DEL TIEMPO########################
	46	#######################################################################
	47
	48	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
	49	'''
	50	opObj11 = procUnitConfObjA.addOperation(name='PulsePairVoltage', optype='other')
	51	opObj11.addParameter(name='n', value='256', format='int')
	52	opObj11.addParameter(name='removeDC', value=1, format='int')
	53	'''
	54
	55	_type="iq"
	56	opObj10 = procUnitConfObjA.addOperation(name='ScopePlot', optype='external')
	57	#opObj10.addParameter(name='id', value='12')
	58	opObj10.addParameter(name='wintitle', value=_type )
	59	opObj10.addParameter(name='type', value=_type)
	60
	61
	62	'''
	63	type="WeatherPower"
	64	opObj10 = procUnitConfObjA.addOperation(name='PulsepairPowerPlot', optype='external')
	65	#opObj10.addParameter(name='id', value='12')
	66	opObj10.addParameter(name='wintitle', value=type )
	67
	68	opObj11 = procUnitConfObjA.addOperation(name='PulsepairVelocityPlot', optype='other')
	69	opObj11.addParameter(name='xmax', value=8)
	70	'''
	71
	72	controllerObj.start()

schainpy/scripts/SOPHY_test0001.py created 644 +118 0

@@ -0,0 +1,118
	1	# Ing. AVP
	2	# 01/11/2021
	3	# ARCHIVO DE LECTURA
	4	import os, sys
	5	import datetime
	6	import time
	7	from schainpy.controller import Project
	8	print("----[Setup]-RadarMeteorologico--------")
	9	Vel = 10
	10	modo_proc = 0 # 0 Pulse Pair 1 Spectros
	11	#-----------PATH DE DATOS-----------------------#
	12	#------------VERIFICAR SIEMPRE LA FECHA DE LA DATA
	13	path = "/DATA_RM/10"
	14	path_ped = "/DATA_RM/TEST_PEDESTAL/P20211110-171003"
	15	figpath_pp = "/home/soporte/Pictures/TEST_POT"
	16	figpath_ppi_pp = "/home/soporte/Pictures/ppi_PP_30DIC_4"
	17	#-------------------------------------------------------------------
	18	print("----[OPCIONES]------------------------")
	19	op_plot = 0
	20	op_integration = 1
	21	op_save = 1
	22	op_plot_spec = 0
	23	#-------------------------------------
	24	################# RANGO DE PLOTEO######################################
	25	dBmin = '1'
	26	dBmax = '85'
	27	xmin = '17.1'
	28	xmax = '17.25'
	29	ymin = '0'
	30	ymax = '600'
	31	#-------------------NRO Perfiles PROCESAMIENTO --------------------
	32	V=Vel
	33	IPP=400*1e-6
	34	n= int(1/(V*IPP))
	35	print("* n - NRO Perfiles Proc:", n )
	36	time.sleep(3)
	37
	38	#------------------------SIGNAL CHAIN ------------------------------
	39	desc = "USRP_test"
	40	filename = "USRP_processing.xml"
	41	controllerObj = Project()
	42	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
	43	######################## UNIDAD DE LECTURA#############################
	44	readUnitConfObj = controllerObj.addReadUnit(datatype = 'DigitalRFReader',
	45	path = path,
	46	startDate= "2021/11/10",#today,
	47	endDate = "2021/12/30",#today,
	48	startTime= '00:00:25',
	49	endTime = '23:59:59',
	50	delay = 0,
	51	online = 0,
	52	walk = 1,
	53	ippKm = 60)
	54
	55	opObj11 = readUnitConfObj.addOperation(name='printInfo')
	56	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc',inputId=readUnitConfObj.getId())
	57
	58	opObj11 = procUnitConfObjA.addOperation(name='selectHeights')
	59	opObj11.addParameter(name='minIndex', value='1', format='int')
	60	# opObj11.addParameter(name='maxIndex', value='10000', format='int')
	61	opObj11.addParameter(name='maxIndex', value='400', format='int')
	62
	63	if modo_proc ==0:
	64	#----------------------------------------PULSE PAIR --------------------------------------------------#
	65	opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
	66	opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO QUE LA VELOCIDAD ES 10 GRADOS
	67	#opObj11.addParameter(name='removeDC', value=1, format='int')
	68	#------------------------ METODO Parametros -----------------------------------------------------------
	69	procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
	70	if op_plot==1:
	71	opObj11 = procUnitConfObjB.addOperation(name='GenericRTIPlot',optype='external')
	72	opObj11.addParameter(name='attr_data', value='dataPP_POW')
	73	opObj11.addParameter(name='colormap', value='jet')
	74	#opObj11.addParameter(name='xmin', value=xmin)
	75	#opObj11.addParameter(name='xmax', value=xmax)
	76	opObj11.addParameter(name='zmin', value=dBmin)
	77	opObj11.addParameter(name='zmax', value=dBmax)
	78	opObj11.addParameter(name='save', value=figpath_pp)
	79	opObj11.addParameter(name='showprofile', value=0)
	80	opObj11.addParameter(name='save_period', value=50)
	81
	82	####################### METODO ESCRITURA #######################################################################
	83
	84	if op_integration==1:
	85	V=V
	86	blocksPerfile=100
	87	print("* Velocidad del Pedestal:",V)
	88	tmp_blocksPerfile = 100
	89	f_a_p= int(tmp_blocksPerfile/V)
	90
	91	opObj11 = procUnitConfObjB.addOperation(name='PedestalInformation')
	92	opObj11.addParameter(name='path_ped', value=path_ped)
	93	#opObj11.addParameter(name='path_adq', value=path_adq)
	94	opObj11.addParameter(name='t_Interval_p', value='0.01', format='float')
	95	opObj11.addParameter(name='blocksPerfile', value=blocksPerfile, format='int')
	96	opObj11.addParameter(name='n_Muestras_p', value='100', format='float')
	97	opObj11.addParameter(name='f_a_p', value=f_a_p, format='int')
	98	opObj11.addParameter(name='online', value='0', format='int')
	99
	100	opObj11 = procUnitConfObjB.addOperation(name='Block360')
	101	opObj11.addParameter(name='n', value='10', format='int')
	102	opObj11.addParameter(name='mode', value=modo_proc, format='int')
	103
	104	# este bloque funciona bien con divisores de 360 no olvidar 0 10 20 30 40 60 90 120 180
	105
	106	opObj11= procUnitConfObjB.addOperation(name='WeatherPlot',optype='other')
	107	opObj11.addParameter(name='save', value=figpath_ppi_pp)
	108	opObj11.addParameter(name='save_period', value=1)
	109
	110	if op_save==1:
	111	opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')
	112	opObj10.addParameter(name='path',value=figpath_ppi_pp)
	113	#opObj10.addParameter(name='mode',value=0)
	114	opObj10.addParameter(name='blocksPerFile',value='100',format='int')
	115	opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list')
	116	opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,azimuth,utctime',format='list')#,format='list'
	117
	118	controllerObj.start()

schainpy/scripts/avp_001.py created 644 +14 0

@@ -0,0 +1,14
	1	import numpy
	2	a= numpy.array([0,1,2,3,4,5,10,11,12,18,19,20,21,22,23,24,25,26,27,28])
	3	print(a)
	4	list=[]
	5	list2=[]
	6	for i in reversed(range(1,len(a))):
	7	dif=int(a[i])-int(a[i-1])
	8	print(i,a[i],dif )
	9	if dif>1:
	10	list.append(i-1)
	11	list2.append(dif-1)
	12	print("result")
	13	print(list)
	14	print(list2)

schainpy/scripts/pp__proc.py created 644 +147 0

@@ -0,0 +1,147
	1	#!python
	2	'''
	3	'''
	4
	5	import os, sys
	6	import datetime
	7	import time
	8
	9	#path = os.path.dirname(os.getcwd())
	10	#path = os.path.dirname(path)
	11	#sys.path.insert(0, path)
	12
	13	from schainpy.controller import Project
	14
	15	desc = "USRP_test"
	16	filename = "USRP_processing.xml"
	17	controllerObj = Project()
	18	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
	19
	20	############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############
	21
	22	#######################################################################
	23	######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################
	24	#######################################################################
	25	#path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'
	26	#path = '/DATA_RM/TEST_INTEGRACION'
	27	#path = '/DATA_RM/TEST_ONLINE'
	28	#path = '/DATA_RM/TEST_INTEGRACION/ADQ_OFFLINE/'
	29	# ULTIMO TEST 22 DE SEPTIEMBRE
	30	path = '/DATA_RM/USRP_22'
	31	#path_pp = '/DATA_RM/TEST_HDF5'
	32	# UTIMO TEST 22 DE SEPTIEMBRE
	33	path_pp = '/DATA_RM/TEST_HDF5_PP_22'
	34	######################################################
	35	##### OJO TENER EN CUENTA EL n= para el Pulse Pair ###
	36	######################################################
	37	######## BUSCAMOS EL numero de IPP equivalente 1°#####
	38	######## Sea V la velocidad del Pedestal en °/seg#####
	39	######## 1° sera Recorrido en un tiempo de 1/V ######
	40	######## IPP del Radar 400 useg --> 60 Km ############
	41	######## n = 1/(V*IPP) , NUMERO DE IPP #############
	42	######## n = 1/(V*IPP) #############################
	43	V=2
	44	IPP=400*1e-6
	45	n= 1/(V*IPP)
	46	print("n numero de Perfiles a procesar con Pulse Pair: ", n)
	47
	48
	49
	50
	51	figpath = '/home/soporte/Pictures/TEST_INTEGRACION_IMG'
	52	#remotefolder = "/home/wmaster/graficos"
	53	#######################################################################
	54	################# RANGO DE PLOTEO######################################
	55	#######################################################################
	56	dBmin = '-5'
	57	dBmax = '20'
	58	xmin = '0'
	59	xmax ='24'
	60	ymin = '0'
	61	ymax = '600'
	62	#######################################################################
	63	########################FECHA##########################################
	64	#######################################################################
	65	str = datetime.date.today()
	66	today = str.strftime("%Y/%m/%d")
	67	str2 = str - datetime.timedelta(days=1)
	68	yesterday = str2.strftime("%Y/%m/%d")
	69	#######################################################################
	70	######################## UNIDAD DE LECTURA#############################
	71	#######################################################################
	72	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
	73	path=path,
	74	startDate="2021/01/01",#today,
	75	endDate="2021/12/30",#today,
	76	startTime='00:00:00',
	77	endTime='23:59:59',
	78	delay=0,
	79	#set=0,
	80	online=0,
	81	walk=1,
	82	ippKm = 60)
	83
	84	opObj11 = readUnitConfObj.addOperation(name='printInfo')
	85	#opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
	86	#######################################################################
	87	################ OPERACIONES DOMINIO DEL TIEMPO########################
	88	#######################################################################
	89
	90	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
	91
	92	#
	93	# codigo64='1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,0,0,1,0,1,1,1,0,0,0,1,0,'+\
	94	# '1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,0,0,0,1,0,0,1,0,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1'
	95
	96	#opObj11 = procUnitConfObjA.addOperation(name='setRadarFrequency')
	97	#opObj11.addParameter(name='frequency', value='70312500')
	98	opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
	99	opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO QUE LA VELOCIDAD ES 10 GRADOS
	100	opObj11.addParameter(name='removeDC', value=1, format='int')
	101	# Ploteo TEST
	102	'''
	103	opObj11 = procUnitConfObjA.addOperation(name='PulsepairPowerPlot', optype='other')
	104	opObj11 = procUnitConfObjA.addOperation(name='PulsepairSignalPlot', optype='other')
	105	opObj11 = procUnitConfObjA.addOperation(name='PulsepairVelocityPlot', optype='other')
	106	#opObj11.addParameter(name='xmax', value=8)
	107	opObj11 = procUnitConfObjA.addOperation(name='PulsepairSpecwidthPlot', optype='other')
	108	'''
	109	# OJO SCOPE
	110	#opObj10 = procUnitConfObjA.addOperation(name='ScopePlot', optype='external')
	111	#opObj10.addParameter(name='id', value='10', format='int')
	112	##opObj10.addParameter(name='xmin', value='0', format='int')
	113	##opObj10.addParameter(name='xmax', value='50', format='int')
	114	#opObj10.addParameter(name='type', value='iq')
	115	##opObj10.addParameter(name='ymin', value='-5000', format='int')
	116	##opObj10.addParameter(name='ymax', value='8500', format='int')
	117	#opObj11.addParameter(name='save', value=figpath, format='str')
	118	#opObj11.addParameter(name='save_period', value=10, format='int')
	119
	120	#opObj10 = procUnitConfObjA.addOperation(name='setH0')
	121	#opObj10.addParameter(name='h0', value='-5000', format='float')
	122
	123	#opObj11 = procUnitConfObjA.addOperation(name='filterByHeights')
	124	#opObj11.addParameter(name='window', value='1', format='int')
	125
	126	#codigo='1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,1,1,-1,-1,-1'
	127	#opObj11 = procUnitConfObjSousy.addOperation(name='Decoder', optype='other')
	128	#opObj11.addParameter(name='code', value=codigo, formatyesterday='floatlist')
	129	#opObj11.addParameter(name='nCode', value='1', format='int')
	130	#opObj11.addParameter(name='nBaud', value='28', format='int')
	131
	132	#opObj11 = procUnitConfObjA.addOperation(name='CohInt', optype='other')
	133	#opObj11.addParameter(name='n', value='100', format='int')
	134
	135	#######################################################################
	136	########## OPERACIONES ParametersProc########################
	137	#######################################################################
	138
	139	procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
	140	opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')
	141	opObj10.addParameter(name='path',value=path_pp)
	142	#opObj10.addParameter(name='mode',value=0)
	143	opObj10.addParameter(name='blocksPerFile',value='100',format='int')
	144	opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list')
	145	opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'
	146
	147	controllerObj.start()

schainpy/scripts/pp_proc.py created 644 +115 0

@@ -0,0 +1,115
	1	#!python
	2	'''
	3	'''
	4
	5	import os, sys
	6	import datetime
	7	import time
	8
	9	#path = os.path.dirname(os.getcwd())
	10	#path = os.path.dirname(path)
	11	#sys.path.insert(0, path)
	12
	13	from schainpy.controller import Project
	14
	15	desc = "USRP_test"
	16	filename = "USRP_processing.xml"
	17	controllerObj = Project()
	18	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
	19
	20	############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############
	21
	22	#######################################################################
	23	######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################
	24	#######################################################################
	25	#path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'
	26	#path = '/DATA_RM/TEST_INTEGRACION'
	27	#path = '/DATA_RM/TEST_ONLINE'
	28	#path = '/DATA_RM/TEST_INTEGRACION/ADQ_OFFLINE/'
	29	# ULTIMO TEST 22 DE SEPTIEMBRE
	30	path = '/DATA_RM/USRP_22'
	31	#path_pp = '/DATA_RM/TEST_HDF5'
	32	# UTIMO TEST 22 DE SEPTIEMBRE
	33	path_pp = '/DATA_RM/TEST_HDF5_PP_22'
	34	######################################################
	35	##### OJO TENER EN CUENTA EL n= para el Pulse Pair ###
	36	######################################################
	37	######## BUSCAMOS EL numero de IPP equivalente 1°#####
	38	######## Sea V la velocidad del Pedestal en °/seg#####
	39	######## 1° sera Recorrido en un tiempo de 1/V ######
	40	######## IPP del Radar 400 useg --> 60 Km ############
	41	######## n = 1/(V*IPP) , NUMERO DE IPP #############
	42	######## n = 1/(V*IPP) #############################
	43	V=2
	44	IPP=400*1e-6
	45	n= 1/(V*IPP)
	46	print("n numero de Perfiles a procesar con Pulse Pair: ", n)
	47
	48
	49
	50
	51	figpath = '/home/soporte/Pictures/TEST_INTEGRACION_IMG'
	52	#remotefolder = "/home/wmaster/graficos"
	53	#######################################################################
	54	################# RANGO DE PLOTEO######################################
	55	#######################################################################
	56	dBmin = '-5'
	57	dBmax = '20'
	58	xmin = '0'
	59	xmax ='24'
	60	ymin = '0'
	61	ymax = '600'
	62	#######################################################################
	63	########################FECHA##########################################
	64	#######################################################################
	65	str = datetime.date.today()
	66	today = str.strftime("%Y/%m/%d")
	67	str2 = str - datetime.timedelta(days=1)
	68	yesterday = str2.strftime("%Y/%m/%d")
	69	#######################################################################
	70	######################## UNIDAD DE LECTURA#############################
	71	#######################################################################
	72	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
	73	path=path,
	74	startDate="2021/01/01",#today,
	75	endDate="2021/12/30",#today,
	76	startTime='00:00:00',
	77	endTime='23:59:59',
	78	delay=0,
	79	#set=0,
	80	online=0,
	81	walk=1,
	82	ippKm = 60)
	83
	84	opObj11 = readUnitConfObj.addOperation(name='printInfo')
	85	#opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
	86	#######################################################################
	87	################ OPERACIONES DOMINIO DEL TIEMPO########################
	88	#######################################################################
	89
	90	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
	91
	92	#
	93	# codigo64='1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,0,0,1,0,1,1,1,0,0,0,1,0,'+\
	94	# '1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,0,0,0,1,0,0,1,0,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1'
	95
	96	#opObj11 = procUnitConfObjA.addOperation(name='setRadarFrequency')
	97	#opObj11.addParameter(name='frequency', value='70312500')
	98	opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
	99	opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO QUE LA VELOCIDAD ES 10 GRADOS
	100	opObj11.addParameter(name='removeDC', value=1, format='int')
	101
	102
	103	#######################################################################
	104	########## OPERACIONES ParametersProc########################
	105	#######################################################################
	106
	107	procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
	108	opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')
	109	opObj10.addParameter(name='path',value=path_pp)
	110	#opObj10.addParameter(name='mode',value=0)
	111	opObj10.addParameter(name='blocksPerFile',value='100',format='int')
	112	opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list')
	113	opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'
	114
	115	controllerObj.start()

schainpy/scripts/spec_momentos_proc.py created 644 +88 0

@@ -0,0 +1,88
	1	#!python
	2
	3	import os, sys
	4	import datetime
	5	import time
	6
	7	from schainpy.controller import Project
	8
	9	desc = "USRP_test"
	10	filename = "USRP_processing.xml"
	11	controllerObj = Project()
	12	controllerObj.setup(id = '191', name='Test_USRP', description=desc)
	13
	14	############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############
	15
	16	#######################################################################
	17	######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################
	18	#######################################################################
	19
	20	path = '/DATA_RM/TEST_19OCTUBRE/10MHZ'
	21	#path_pp = '/DATA_RM/TEST_HDF5'
	22	path_pp = '/DATA_RM/TEST_HDF5_19OCT'
	23
	24	#######################################################################
	25	################# RANGO DE PLOTEO######################################
	26	#######################################################################
	27	dBmin = '-5'
	28	dBmax = '20'
	29	xmin = '0'
	30	xmax ='24'
	31	ymin = '0'
	32	ymax = '600'
	33	#######################################################################
	34	########################FECHA##########################################
	35	#######################################################################
	36	str = datetime.date.today()
	37	today = str.strftime("%Y/%m/%d")
	38	str2 = str - datetime.timedelta(days=1)
	39	yesterday = str2.strftime("%Y/%m/%d")
	40	#######################################################################
	41	######################## UNIDAD DE LECTURA#############################
	42	#######################################################################
	43	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
	44	path=path,
	45	startDate="2021/01/01",#today,
	46	endDate="2021/12/30",#today,
	47	startTime='00:00:00',
	48	endTime='23:59:59',
	49	delay=0,
	50	#set=0,
	51	online=0,
	52	walk=1,
	53	ippKm = 60)
	54
	55	opObj11 = readUnitConfObj.addOperation(name='printInfo')
	56	#opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
	57	#######################################################################
	58	################ OPERACIONES DOMINIO DEL TIEMPO########################
	59	#######################################################################
	60	V=10 # aca se coloca la velocidad
	61	IPP=400*1e-6
	62	n= int(1/(V*IPP))
	63	print("n numero de Perfiles a procesar con nFFTPoints ", n)
	64
	65	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
	66
	67	procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
	68	procUnitConfObjB.addParameter(name='nFFTPoints', value=n, format='int')
	69	procUnitConfObjB.addParameter(name='nProfiles' , value=n, format='int')
	70
	71	#######################################################################
	72	########## OPERACIONES ParametersProc########################
	73	#######################################################################
	74
	75	procUnitConfObjC= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjB.getId())
	76
	77	procUnitConfObjC.addOperation(name='SpectralMoments')
	78
	79	opObj10 = procUnitConfObjC.addOperation(name='HDFWriter')
	80	opObj10.addParameter(name='path',value=path_pp)
	81	#opObj10.addParameter(name='mode',value=0)
	82	opObj10.addParameter(name='blocksPerFile',value='100',format='int')
	83	#opObj10.addParameter(name='metadataList',value='utctimeInit,heightList,nIncohInt,nCohInt,nProfiles,channelList',format='list')#profileIndex
	84	opObj10.addParameter(name='metadataList',value='utctimeInit,heightList,nIncohInt,nCohInt,nProfiles,channelList',format='list')#profileIndex
	85
	86	opObj10.addParameter(name='dataList',value='data_pow,data_dop,utctime',format='list')#,format='list'
	87
	88	controllerObj.start()

schainpy/scripts/test_wradlib.py created 644 +36 0

@@ -0,0 +1,36
	1	import numpy as np
	2	import matplotlib.pyplot as pl
	3	import wradlib
	4	import warnings
	5	#export WRADLIB_DATA=/path/to/wradlib-data
	6	warnings.filterwarnings('ignore')
	7	'''
	8	try:
	9	get_ipython().magic('matplotlib inline')
	10	except:
	11	pl.ion()
	12	'''
	13	filename = wradlib.util.get_wradlib_data_file("/home/soporte/Downloads/raa00-dx_10908-0806021735-fbg---bin.gz")
	14	img, meta = wradlib.io.read_dx(filename)
	15	print("Shape of polar array: %r\n" % (img.shape,))
	16	print("Some meta data of the DX file:")
	17	print("\tdatetime: %r" % (meta["datetime"],))
	18	print("\tRadar ID: %s" % (meta["radarid"],))
	19
	20	img[200:250,:]= np.ones([50,img.shape[1]])*np.nan
	21
	22	img[300:360,:]= np.ones([60,img.shape[1]])*np.nan
	23
	24	cgax, pm= wradlib.vis.plot_ppi(img)
	25	txt = pl.title('Simple PPI')
	26	print("coordenada angular",img[:,0],len(img[:,0]))
	27	print("COORDENADA 0",img[0],len(img[0]))
	28	cbar = pl.gcf().colorbar(pm, pad=0.075)
	29
	30	#r = np.arange(40, 80)
	31	#az = np.arange(200, 250)
	32	#ax, pm = wradlib.vis.plot_ppi(img[200:250, 40:80], r, az, autoext=False)
	33	#ax, pm = wradlib.vis.plot_ppi(img[200:250, 40:80], r, az)
	34
	35	#txt = pl.title('Sector PPI')
	36	pl.show()

schainpy/sophy_test0001.py created 644 +17 0

@@ -0,0 +1,17
	1	# Ing. AVP
	2	# 01/11/2021
	3	# ARCHIVO DE LECTURA
	4	import os, sys
	5	import datetime
	6	import time
	7	from schainpy.controller import Project
	8	print("----[Setup]-RadarMeteorologico--------")
	9	Vel = 6
	10	modo_proc = 1
	11	path = "/DATA_RM/10"
	12	path_ped = "/DATA_RM/TEST_PEDESTAL/P20211111-173856"
	13	print("----[OPCIONES]------------------------")
	14	op_plot = 0
	15	op_integration = 0
	16	op_save = 0
	17	op_plot_spec = 0

schainpy/model/graphics/jroplot_parameters.py +116 -28

@@ -390,6 +390,8 class WeatherPlot(Plot):
390	self.plots_adjust.update({'wspace': 0.4, 'hspace':0.4, 'left': 0.1, 'right': 0.9, 'bottom': 0.08})	390	self.plots_adjust.update({'wspace': 0.4, 'hspace':0.4, 'left': 0.1, 'right': 0.9, 'bottom': 0.08})
391	self.flag =0	391	self.flag =0
392	self.indicador= 0	392	self.indicador= 0
		393	self.last_data_azi = None
		394	self.val_mean = None
393		395
394	def update(self, dataOut):	396	def update(self, dataOut):
395		397
@@ -402,41 +404,115 class WeatherPlot(Plot):
402	factor = dataOut.normFactor	404	factor = dataOut.normFactor
403		405
404	####print("factor",factor)	406	####print("factor",factor)
405	data['weather'] = 10*numpy.log10(dataOut.data_360[0]/(factor))	407	data['weather'] = 10*numpy.log10(dataOut.data_360[1]/(factor))
406	print("weather",data['weather'])	408	####print("weather",data['weather'])
407	data['azi'] = dataOut.data_azi	409	data['azi'] = dataOut.data_azi
408	return data, meta	410	return data, meta
409		411
		412	def analizeDATA(self,data_azi):
		413	list1 = []
		414	list2 = []
		415	dat = data_azi
		416	for i in reversed(range(1,len(dat))):
		417	if dat[i]>dat[i-1]:
		418	diff = int(dat[i])-int(dat[i-1])
		419	else:
		420	diff = 360+int(dat[i])-int(dat[i-1])
		421	if diff > 1:
		422	list1.append(i-1)
		423	list2.append(diff-1)
		424	return list1,list2
		425
		426	def fixDATANEW(self,data_azi,data_weather):
		427	list1,list2 = self.analizeDATA(data_azi)
		428	if len(list1)== 0:
		429	return data_azi,data_weather
		430	else:
		431	resize = 0
		432	for i in range(len(list2)):
		433	resize= resize + list2[i]
		434	new_data_azi = numpy.resize(data_azi,resize)
		435	new_data_weather= numpy.resize(date_weather,resize)
		436
		437	for i in range(len(list2)):
		438	j=0
		439	position=list1[i]+1
		440	for j in range(list2[i]):
		441	new_data_azi[position+j]=new_data_azi[position+j-1]+1
		442
		443	return new_data_azi
		444
		445	def fixDATA(self,data_azi):
		446	data=data_azi
		447	for i in range(len(data)):
		448	if numpy.isnan(data[i]):
		449	data[i]=data[i-1]+1
		450	return data
		451
		452	def replaceNAN(self,data_weather,data_azi,val):
		453	print("----------------activeNEWFUNCTION")
		454	data= data_azi
		455	data_T= data_weather
		456	print("data_azi",data_azi)
		457	print("VAL:",val)
		458	print("SHAPE",data_T.shape)
		459	for i in range(len(data)):
		460	if numpy.isnan(data[i]):
		461	print("NAN")
		462	data_T[i,:]=numpy.ones(data_T.shape[1])*val
		463	#data_T[i,:]=numpy.ones(data_T.shape[1])*numpy.nan
		464	return data_T
		465
410	def const_ploteo(self,data_weather,data_azi,step,res):	466	def const_ploteo(self,data_weather,data_azi,step,res):
411	if self.ini==0:	467	if self.ini==0:
412	#------- AZIMUTH	468	#------- AZIMUTH
413	n = (360/res)-len(data_azi)	469	n = (360/res)-len(data_azi)
414	start = data_azi[-1] + res	470	##### new
415	~~end~~ = data_azi[0] - ~~res~~	471	data_azi_old = data_azi
		472	data_azi_new = self.fixDATA(data_azi)
		473	#ata_azi_new = self.fixDATANEW(data_azi)
		474
		475	start = data_azi_new[-1] + res
		476	end = data_azi_new[0] - res
		477	##### new
		478	self.last_data_azi = end
416	if start>end:	479	if start>end:
417	end = end + 360	480	end = end + 360
418	azi_vacia = numpy.linspace(start,end,int(n))	481	azi_vacia = numpy.linspace(start,end,int(n))
419	azi_vacia = numpy.where(azi_vacia>360,azi_vacia-360,azi_vacia)	482	azi_vacia = numpy.where(azi_vacia>360,azi_vacia-360,azi_vacia)
420	data_azi = numpy.hstack((data_azi,azi_vacia))	483	data_azi = numpy.hstack((data_azi_new,azi_vacia))
421	# RADAR	484	# RADAR
422	val_mean = numpy.mean(data_weather[:,0])	485	val_mean = numpy.mean(data_weather[:,-1])
		486	self.val_mean = val_mean
423	data_weather_cmp = numpy.ones([(360-data_weather.shape[0]),data_weather.shape[1]])*val_mean	487	data_weather_cmp = numpy.ones([(360-data_weather.shape[0]),data_weather.shape[1]])*val_mean
		488	data_weather = self.replaceNAN(data_weather=data_weather,data_azi=data_azi_old,val=self.val_mean)
424	data_weather = numpy.vstack((data_weather,data_weather_cmp))	489	data_weather = numpy.vstack((data_weather,data_weather_cmp))
425	else:	490	else:
426	# azimuth	491	# azimuth
427	flag=0	492	flag=0
428	start_azi = self.res_azi[0]	493	start_azi = self.res_azi[0]
		494	#### new
		495	data_azi_old = data_azi
		496	### weather ###
		497	data_weather = self.replaceNAN(data_weather=data_weather,data_azi=data_azi_old,val=self.val_mean)
		498
		499	if numpy.isnan(data_azi[0]):
		500	data_azi[0]=self.last_data_azi+1
		501	data_azi = self.fixDATA(data_azi)
		502	####
		503
429	start = data_azi[0]	504	start = data_azi[0]
430	end = data_azi[-1]	505	end = data_azi[-1]
431	print("start",start)	506	self.last_data_azi= end
432	print("~~end"~~,~~end~~)	507	####print("start",start)
		508	####print("end",end)
433	if start< start_azi:	509	if start< start_azi:
434	start = start +360	510	start = start +360
435	if end <start_azi:	511	if end <start_azi:
436	end = end +360	512	end = end +360
437		513
438	print("start",start)	514	####print("start",start)
439	print("end",end)	515	####print("end",end)
440	#### AQUI SERA LA MAGIA	516	#### AQUI SERA LA MAGIA
441	pos_ini = int((start-start_azi)/res)	517	pos_ini = int((start-start_azi)/res)
442	len_azi = len(data_azi)	518	len_azi = len(data_azi)
@@ -448,9 +524,10 class WeatherPlot(Plot):
448	dif= 360-pos_ini	524	dif= 360-pos_ini
449	comp= len_azi-dif	525	comp= len_azi-dif
450		526
451	print(pos_ini)	527	#-----------------
452	print(~~len_az~~i)	528	####print(pos_ini)
453	print(~~"shape"~~,~~self~~.~~res_azi~~.~~shape~~)	529	####print(len_azi)
		530	####print("shape",self.res_azi.shape)
454	if flag==0:	531	if flag==0:
455	# AZIMUTH	532	# AZIMUTH
456	self.res_azi[pos_ini:pos_ini+len_azi] = data_azi	533	self.res_azi[pos_ini:pos_ini+len_azi] = data_azi
@@ -470,17 +547,28 class WeatherPlot(Plot):
470	return data_weather,data_azi	547	return data_weather,data_azi
471		548
472	def plot(self):	549	def plot(self):
473	print("--------------------------------------",self.ini,"-----------------------------------")	550	#print("--------------------------------------",self.ini,"-----------------------------------")
474	#numpy.set_printoptions(suppress=True)	551	#numpy.set_printoptions(suppress=True)
475	#print(self.data.times)	552	####print("times: ",self.data.times)
476	thisDatetime = datetime.datetime.utcfromtimestamp(self.data.times[-1])	553	thisDatetime = datetime.datetime.utcfromtimestamp(self.data.times[-1]).strftime('%Y-%m-%d %H:%M:%S')
		554	#print("times: ",thisDatetime)
477	data = self.data[-1]	555	data = self.data[-1]
478	# ALTURA altura_tmp_h	556	####ALTURA altura_tmp_h
479	altura_h = (data['weather'].shape[1])/10.0	557	###print("Y RANGES",self.data.yrange,len(self.data.yrange))
480	stoprange = float(altura_h1.5)#stoprange = float(331.5) por ahora 400	558	###altura_h = (data['weather'].shape[1])/10.0
481	rangestep = float(0.15)	559	###stoprange = float(altura_h0.3)#stoprange = float(331.5) por ahora 400
482	r = numpy.arange(0, stoprange, rangestep)	560	###rangestep = float(0.03)
		561	###r = numpy.arange(0, stoprange, rangestep)
		562	###print("r",r,len(r))
		563	#-----------------------------update----------------------
		564	r= self.data.yrange
		565	delta_height = r[1]-r[0]
		566	#print("1",r)
		567	r_mask= numpy.where(r>=0)[0]
		568	r = numpy.arange(len(r_mask))*delta_height
		569	#print("2",r)
483	self.y = 2*r	570	self.y = 2*r
		571	######self.y = self.data.yrange
484	# RADAR	572	# RADAR
485	#data_weather = data['weather']	573	#data_weather = data['weather']
486	# PEDESTAL	574	# PEDESTAL
@@ -491,27 +579,27 class WeatherPlot(Plot):
491	#print("shape wr_data", wr_data.shape)	579	#print("shape wr_data", wr_data.shape)
492	#print("shape wr_azi",wr_azi.shape)	580	#print("shape wr_azi",wr_azi.shape)
493	#print("step",step)	581	#print("step",step)
494	print("Time---->",self.data.times[-1],thisDatetime)	582	####print("Time---->",self.data.times[-1],thisDatetime)
495	#print("alturas", len(self.y))	583	#print("alturas", len(self.y))numpy.where(r>=0)
496	self.res_weather, self.res_azi = self.const_ploteo(data_weather=data['weather'],data_azi=data['azi'],step=step,res=res)	584	self.res_weather, self.res_azi = self.const_ploteo(data_weather=data['weather'][:,r_mask],data_azi=data['azi'],step=step,res=res)
497	#numpy.set_printoptions(suppress=True)	585	#numpy.set_printoptions(suppress=True)
498	#print("resultado",self.res_azi)	586	#print("resultado",self.res_azi)
499	##########################################################	587	###########################/DATA_RM/10_tmp/ch0###############################
500	################# PLOTEO ###################	588	################# PLOTEO ###################
501	##########################################################	589	##########################################################
502		590
503	for i,ax in enumerate(self.axes):	591	for i,ax in enumerate(self.axes):
504	if ax.firsttime:	592	if ax.firsttime:
505	plt.clf()	593	plt.clf()
506	cgax, pm = wrl.vis.plot_ppi(self.res_weather,r=r,az=self.res_azi,fig=self.figures[0], proj='cg', vmin=1, vmax=60)	594	cgax, pm = wrl.vis.plot_ppi(self.res_weather,r=r,az=self.res_azi,fig=self.figures[0], proj='cg', vmin=8, vmax=35)
507	else:	595	else:
508	plt.clf()	596	plt.clf()
509	cgax, pm = wrl.vis.plot_ppi(self.res_weather,r=r,az=self.res_azi,fig=self.figures[0], proj='cg', vmin=1, vmax=60)	597	cgax, pm = wrl.vis.plot_ppi(self.res_weather,r=r,az=self.res_azi,fig=self.figures[0], proj='cg', vmin=8, vmax=35)
510	caax = cgax.parasites[0]	598	caax = cgax.parasites[0]
511	paax = cgax.parasites[1]	599	paax = cgax.parasites[1]
512	cbar = plt.gcf().colorbar(pm, pad=0.075)	600	cbar = plt.gcf().colorbar(pm, pad=0.075)
513	caax.set_xlabel('x_range [km]')	601	caax.set_xlabel('x_range [km]')
514	caax.set_ylabel('y_range [km]')	602	caax.set_ylabel('y_range [km]')
515	plt.text(1.0, 1.05, 'azimuth '+str(thisDatetime)+"step"+str(self.ini), transform=caax.transAxes, va='bottom',ha='right')	603	plt.text(1.0, 1.05, 'azimuth '+str(thisDatetime)+" step "+str(self.ini), transform=caax.transAxes, va='bottom',ha='right')
516		604
517	self.ini= self.ini+1	605	self.ini= self.ini+1

schainpy/model/io/jroIO_digitalRF.py +1 -1

                     start_index, end_index = self.digitalReadObj.get_bounds(
                         channelNameList[channelList[0]])
-                    print("*****",start_index,end_index)
+                    ##print("*****",start_index,end_index)
                     if not startUTCSecond:
                         startUTCSecond = start_index / self.__sample_rate

schainpy/model/io/jroIO_param.py +63 -41

                 This unit reads HDF5 files created with `HDFWriter` operation contains
                 by default two groups Data and Metadata all variables would be saved as `dataOut`
                 attributes.
                 It is possible to read any HDF5 file by given the structure in the `description`
                 parameter, also you can add extra values to metadata with the parameter `extras`.
                     Dictionary with the description of the HDF5 file
                 extras : dict, optional
                     Dictionary with extra metadata to be be added to `dataOut`
                 Examples
                 --------
                 desc = {
                     'Data': {
                         'data_output': ['u', 'v', 'w'],
                 extras = {
                     'timeZone': 300
                 }
                 reader = project.addReadUnit(
                     name='HDFReader',
                     path='/path/to/files',
                     self.set_kwargs(**kwargs)
                     if not self.ext.startswith('.'):
                         self.ext = '.{}'.format(self.ext)
                     if self.online:
                         log.log("Searching files in online mode...", self.name)
                         for nTries in range(self.nTries):
                             fullpath = self.searchFilesOnLine(self.path, self.startDate,
                                 self.endDate, self.expLabel, self.ext, self.walk,
                                 self.filefmt, self.folderfmt)
                             try:
                                 fullpath = next(fullpath)
                             except:
                                 fullpath = None
                             if fullpath:
                                 break
                             log.warning(
                                 'Waiting {} sec for a valid file in {}: try {} ...'.format(
                                     self.delay, self.path, nTries + 1),
                                 self.name)
                             time.sleep(self.delay)
                         if not(fullpath):
                             raise schainpy.admin.SchainError(
                                 'There isn\'t any valid file in {}'.format(self.path))
                         pathname, filename = os.path.split(fullpath)
                         self.year = int(filename[1:5])
                         self.doy = int(filename[5:8])
                         self.set = int(filename[8:11]) - 1
                     else:
                         log.log("Searching files in {}".format(self.path), self.name)
                         self.filenameList = self.searchFilesOffLine(self.path, self.startDate,
                             self.endDate, self.expLabel, self.ext, self.walk, self.filefmt, self.folderfmt)
                     self.setNextFile()
                     return
                 def readFirstHeader(self):
                     '''Read metadata and data'''
                     self.__readMetadata()
                     self.__readData()
                     self.__setBlockList()
                     if 'type' in self.meta:
                         self.dataOut = eval(self.meta['type'])()
                     for attr in self.meta:
                         setattr(self.dataOut, attr, self.meta[attr])
                     self.blockIndex = 0
                     return
                 def __setBlockList(self):
                 def __readData(self):
                     data = {}
                     if self.description:
                         for key, value in self.description['Data'].items():
                             if isinstance(value, str):
                                 array = numpy.array(array)
                             else:
                                 log.warning('Unknown type: {}'.format(name))
                             if name in self.description:
                                 key = self.description[name]
                             else:
                     self.data = data
                     return
                 def getData(self):
                     for attr in self.data:
                 The HDF5 file contains by default two groups Data and Metadata where
                 you can save any `dataOut` attribute specified by `dataList` and `metadataList`
                 parameters, data attributes are normaly time dependent where the metadata
                 are not.
                 It is possible to customize the structure of the HDF5 file with the
                 optional description parameter see the examples.
                 Parameters:
                     If True the name of the files corresponds to the timestamp of the data
                 description : dict, optional
                     Dictionary with the desired description of the HDF5 file
                 Examples
                 --------
                 desc = {
                     'data_output': {'winds': ['z', 'w', 'v']},
                     'utctime': 'timestamps',
                         'heightList': 'heights'
                     }
                 }
                 writer = proc_unit.addOperation(name='HDFWriter')
                 writer.addParameter(name='path', value='/path/to/file')
                 writer.addParameter(name='blocksPerFile', value='32')
                 metadataList = None
                 currentDay = None
                 lastTime = None
+                last_Azipos = None
+                last_Elepos = None
+                mode       = None
                 def __init__(self):
                     Operation.__init__(self)
                     return
+                def generalFlag(self):
+                    print("GENERALFLAG")
+                    if self.mode== "weather":
+                        if self.last_Azipos == None:
+                            tmp = self.dataOut.azimuth
+                            print("ang azimuth writer",tmp)
+                            self.last_Azipos = tmp
+                            flag = False
+                            return flag
+                        print("ang_azimuth writer",self.dataOut.azimuth)
+                        result = self.dataOut.azimuth - self.last_Azipos
+                        self.last_Azipos = self.dataOut.azimuth
+                        if result<0:
+                            flag = True
+                            return flag
                 def setup(self, path=None, blocksPerFile=10, metadataList=None, dataList=None, setType=None, description=None):
                     self.path = path
                     self.blocksPerFile = blocksPerFile
                             dsDict['shape'] = dataAux.shape
                             dsDict['dsNumber'] = dataAux.shape[0]
                             dsDict['dtype'] = dataAux.dtype
                         dsList.append(dsDict)
                     self.dsList = dsList
                         self.lastTime = currentTime
                         self.currentDay = dataDay
                         return False
                     timeDiff = currentTime - self.lastTime
                     #Si el dia es diferente o si la diferencia entre un dato y otro supera la hora
                         return False
                 def run(self, dataOut, path, blocksPerFile=10, metadataList=None,
-                        dataList=[], setType=None, description={}):
+                        dataList=[], setType=None, description={},mode= None):
                     self.dataOut = dataOut
+                    self.mode    = mode
                     if not(self.isConfig):
                         self.setup(path=path, blocksPerFile=blocksPerFile,
                                    metadataList=metadataList, dataList=dataList,
                                    setType=setType, description=description)
                     self.putData()
                     return
                 def setNextFile(self):
                     ext = self.ext
                     path = self.path
                     setFile = self.setFile
                             return 'pair{:02d}'.format(x)
                         else:
                             return 'channel{:02d}'.format(x)
                 def writeMetadata(self, fp):
                     if self.description:
                     return
                 def writeData(self, fp):
                     if self.description:
                         if 'Data' in self.description:
                             grp = fp.create_group('Data')
                     dtsets = []
                     data = []
                     for dsInfo in self.dsList:
                         if dsInfo['nDim'] == 0:
                             ds = grp.create_dataset(
                                 self.getLabel(dsInfo['variable']),
                                 (self.blocksPerFile, ),
                                 chunks=True,
                                 dtype=numpy.float64)
                             dtsets.append(ds)
                             data.append((dsInfo['variable'], -1))
                                 sgrp = grp
                             for i in range(dsInfo['dsNumber']):
                                 ds = sgrp.create_dataset(
                                     self.getLabel(dsInfo['variable'], i),
                                     (self.blocksPerFile, ) + dsInfo['shape'][1:],
                                     chunks=True,
                                     dtype=dsInfo['dtype'])
                     fp.flush()
                     log.log('Creating file: {}'.format(fp.filename), self.name)
                     self.ds = dtsets
                     self.data = data
                     self.firsttime = True
                     return
                 def putData(self):
+                    print("**************************PUT DATA***************************************************")
-                    if (self.blockIndex == self.blocksPerFile) or self.timeFlag():
+                    if (self.blockIndex == self.blocksPerFile) or self.timeFlag() or self.generalFlag():
                         self.closeFile()
                         self.setNextFile()

schainpy/model/proc/jroproc_parameters.py +201 -66

@@ -23,6 +23,8 from numpy import NaN
23	from scipy.optimize.optimize import OptimizeWarning	23	from scipy.optimize.optimize import OptimizeWarning
24	warnings.filterwarnings('ignore')	24	warnings.filterwarnings('ignore')
25		25
		26	from time import sleep
		27
26	import matplotlib.pyplot as plt	28	import matplotlib.pyplot as plt
27		29
28	SPEED_OF_LIGHT = 299792458	30	SPEED_OF_LIGHT = 299792458
@@ -1323,8 +1325,8 class SpectralMoments(Operation):
1323	if (ss1 > m):	1325	if (ss1 > m):
1324	ss1 = m	1326	ss1 = m
1325		1327
1326	valid = numpy.arange(int(m + bb0 - ss1 + 1)) + ss1	1328	#valid = numpy.arange(int(m + bb0 - ss1 + 1)) + ss1
1327	#valid = numpy.arange(1,oldspec.shape[0])# valid perfil completo igual pulsepair	1329	valid = numpy.arange(1,oldspec.shape[0])# valid perfil completo igual pulsepair
1328	signal_power = ((spec2[valid] - n0) * fwindow[valid]).mean() # D. Scipión added with correct definition	1330	signal_power = ((spec2[valid] - n0) * fwindow[valid]).mean() # D. Scipión added with correct definition
1329	total_power = (spec2[valid] * fwindow[valid]).mean() # D. Scipión added with correct definition	1331	total_power = (spec2[valid] * fwindow[valid]).mean() # D. Scipión added with correct definition
1330	power = ((spec2[valid] - n0) * fwindow[valid]).sum()	1332	power = ((spec2[valid] - n0) * fwindow[valid]).sum()
@@ -4003,6 +4005,39 class PedestalInformation(Operation):
4003	def __init__(self):	4005	def __init__(self):
4004	Operation.__init__(self)	4006	Operation.__init__(self)
4005		4007
		4008
		4009	def getAnguloProfile(self,utc_adq,list_pedestal):
		4010	##print("NEW-METHOD")
		4011	utc_adq = utc_adq
		4012	list_pedestal = list_pedestal
		4013	utc_ped_list = []
		4014
		4015	for i in range(len(list_pedestal)):
		4016	#print(i)# OJO IDENTIFICADOR DE SINCRONISMO
		4017	utc_ped_list.append(self.gettimeutcfromDirFilename(path=self.path_ped,file=list_pedestal[i]))
		4018
		4019	nro_file,utc_ped,utc_ped_1 =self.getNROFile(utc_adq,utc_ped_list)
		4020	#print("utc_adq",utc_adq)
		4021	print("utc_ped",utc_ped)
		4022	print("DIFF",utc_adq-utc_ped)
		4023	print("nro_file",nro_file)
		4024	if nro_file < 0:
		4025	#print("-----------------------------------------------------------------")
		4026	#print("INSERTANDO ANGULO NAN")
		4027	return numpy.NaN
		4028	else:
		4029	nro_key_p = int((utc_adq-utc_ped)/self.t_Interval_p)-1 # ojito al -1 estimado alex
		4030	ff_pedestal = list_pedestal[nro_file]
		4031	#angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azimuth")
		4032	angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azi_pos")
		4033	if 99>=nro_key_p>0:
		4034	##print("angulo_array :",angulo[nro_key_p])
		4035	return angulo[nro_key_p]
		4036	else:
		4037	#print("-----------------------------------------------------------------")
		4038	return numpy.NaN
		4039
		4040
4006	def getfirstFilefromPath(self,path,meta,ext):	4041	def getfirstFilefromPath(self,path,meta,ext):
4007	validFilelist = []	4042	validFilelist = []
4008	#print("SEARH",path)	4043	#print("SEARH",path)
@@ -4068,9 +4103,9 class PedestalInformation(Operation):
4068		4103
4069	def getNROFile(self,utc_adq,utc_ped_list):	4104	def getNROFile(self,utc_adq,utc_ped_list):
4070	c=0	4105	c=0
4071	print("insidegetNROFile")	4106	#print("insidegetNROFile")
4072	print(utc_adq)	4107	#print(utc_adq)
4073	print(len(utc_ped_list))	4108	#print(len(utc_ped_list))
4074	for i in range(len(utc_ped_list)):	4109	for i in range(len(utc_ped_list)):
4075	if utc_adq>utc_ped_list[i]:	4110	if utc_adq>utc_ped_list[i]:
4076	#print("mayor")	4111	#print("mayor")
@@ -4105,11 +4140,11 class PedestalInformation(Operation):
4105	###utc_adq = self.gettimeutcadqfromDirFilename(path=self.path_adq,file=self.list_adq[0])	4140	###utc_adq = self.gettimeutcadqfromDirFilename(path=self.path_adq,file=self.list_adq[0])
4106	print("dios existe donde esta")	4141	print("dios existe donde esta")
4107		4142
4108	#print("utc_ped_list",utc_ped_list)	4143	print("utc_ped_list",utc_ped_list)
4109	###print("utc_adq",utc_adq)	4144	#print("utc_adq",utc_adq)
4110	# utc_adq_dataOut	4145	# utc_adq_dataOut
4111	utc_adq_dataOut =dataOut.utctime	4146	utc_adq_dataOut =dataOut.utctime
4112	print("O~~ffline-utc_adq~~_dataout",utc_adq_dataOut)	4147	print("OFFLINE-UTC_ADQ_dataout",utc_adq_dataOut)
4113		4148
4114	nro_file,utc_ped,utc_ped_1 = self.getNROFile(utc_adq=utc_adq_dataOut, utc_ped_list= utc_ped_list)	4149	nro_file,utc_ped,utc_ped_1 = self.getNROFile(utc_adq=utc_adq_dataOut, utc_ped_list= utc_ped_list)
4115		4150
@@ -4118,12 +4153,14 class PedestalInformation(Operation):
4118	nro_key_p = int((utc_adq_dataOut-utc_ped)/self.t_Interval_p)-1 # ojito al -1 estimado alex	4153	nro_key_p = int((utc_adq_dataOut-utc_ped)/self.t_Interval_p)-1 # ojito al -1 estimado alex
4119	print("nro_key_p",nro_key_p)	4154	print("nro_key_p",nro_key_p)
4120		4155
4121	ff_pedestal = self.list_~~pedestal~~[~~nro_file~~]	4156	ff_pedestal = self.list_
		4157	pedestal[nro_file]
4122	#angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azimuth")	4158	#angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azimuth")
4123	angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azi_pos")	4159	angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azi_pos")
4124		4160
4125	print("utc_pedestal_init :",utc_ped+nro_key_p*self.t_Interval_p)	4161	print("utc_pedestal_init :",utc_ped+nro_key_p*self.t_Interval_p)
4126	print("angulo_array :",angulo[nro_key_p])	4162	if nro_key_p>0:
		4163	print("angulo_array :",angulo[nro_key_p])
4127	self.nro_file = nro_file	4164	self.nro_file = nro_file
4128	self.nro_key_p = nro_key_p	4165	self.nro_key_p = nro_key_p
4129		4166
@@ -4160,6 +4197,7 class PedestalInformation(Operation):
4160	self.t_Interval_p = t_Interval_p	4197	self.t_Interval_p = t_Interval_p
4161	self.n_Muestras_p = n_Muestras_p	4198	self.n_Muestras_p = n_Muestras_p
4162	self.blocksPerfile= blocksPerfile	4199	self.blocksPerfile= blocksPerfile
		4200	#print("dasdadasdasda",self.blocksPerfile)
4163	self.f_a_p = f_a_p	4201	self.f_a_p = f_a_p
4164	self.online = online	4202	self.online = online
4165	self.angulo_adq = numpy.zeros(self.blocksPerfile)	4203	self.angulo_adq = numpy.zeros(self.blocksPerfile)
@@ -4172,14 +4210,14 class PedestalInformation(Operation):
4172	print("LIST NEW", self.list_pedestal[:20])	4210	print("LIST NEW", self.list_pedestal[:20])
4173	#self.list_adq = self.getfirstFilefromPath(path=self.path_adq,meta="D",ext=".hdf5")	4211	#self.list_adq = self.getfirstFilefromPath(path=self.path_adq,meta="D",ext=".hdf5")
4174	print("***********Longitud list pedestal**************",len(self.list_pedestal))	4212	print("***********Longitud list pedestal**************",len(self.list_pedestal))
4175		4213	'''
4176	if self.online:	4214	if self.online:
4177	print("Enable Online")	4215	print("Enable Online")
4178	self.setup_online(dataOut)	4216	self.setup_online(dataOut)
4179	else:	4217	else:
4180	#self.setup_offline(dataOut,list_pedestal=self.list_pedestal,list_adq=self.list_adq)	4218	#self.setup_offline(dataOut,list_pedestal=self.list_pedestal,list_adq=self.list_adq)
4181	self.setup_offline(dataOut,list_pedestal=self.list_pedestal)	4219	self.setup_offline(dataOut,list_pedestal=self.list_pedestal)
4182		4220	'''
4183		4221
4184	def setNextFileP(self,dataOut):	4222	def setNextFileP(self,dataOut):
4185	if self.online:	4223	if self.online:
@@ -4190,52 +4228,109 class PedestalInformation(Operation):
4190	return data_pedestal	4228	return data_pedestal
4191		4229
4192		4230
		4231
		4232	def checkPedFile(self,path,nro_file):
		4233	try:
		4234	utc_pedestal_now = self.gettimeutcfromDirFilename(path=path,file=self.list_pedestal[nro_file])
		4235	utc_pedestal_past = self.gettimeutcfromDirFilename(path=path,file=self.list_pedestal[nro_file-1])
		4236	diff_utc = utc_pedestal_now - utc_pedestal_past
		4237	except:
		4238	diff_utc = -1
		4239	return diff_utc
		4240
4193	def setNextFileoffline(self,dataOut):	4241	def setNextFileoffline(self,dataOut):
4194	##tmp=0	4242	print("error")
		4243	print("no entiendo")
		4244	flag_NOPedfile = False
		4245	cont_NOPedFile = 0
		4246	if self.nro_file<0:
		4247	print("adq empieza antes")
		4248	return numpy.ones(self.blocksPerfile)*numpy.nan
		4249	print("INICIO-----------------------------------------")
		4250	sleep(3)
4195	for j in range(self.blocksPerfile):	4251	for j in range(self.blocksPerfile):
4196	~~###~~print(~~"NUMERO DEL BLOQUE---->",~~j)	4252	print(j)
4197	###print("nro_key_p",self.nro_key_p)	4253	iterador = self.nro_key_p + self.f_a_p*self.c_ped
		4254	self.c_ped = self.c_ped + 1
		4255	print("Iterador-->", iterador)
		4256	if iterador < self.n_Muestras_p:
		4257	self.nro_file = self.nro_file
		4258	else:
		4259	if flag_NOPedfile==False:
		4260	###########################
		4261	self.nro_file = self.nro_file +1
		4262	###########################
		4263	print("nro_file",self.nro_file)
		4264	diff_utc = self.checkPedFile(path=self.path_ped,nro_file=self.nro_file)
		4265	print("diff_utc",diff_utc)
		4266	if diff_utc==1:
		4267	utc_ped_setnext=self.gettimeutcfromDirFilename(path=self.path_ped,file=self.list_pedestal[self.nro_file])
		4268	utc_adq_setnext=dataOut.utctime
		4269	print("utc_pedestal",utc_ped_setnext)
		4270	print("utc_adq",utc_adq_setnext)
		4271	print("self.c_ped",self.c_ped)
		4272	#dif = self.blocksPerfile-(self.nro_key_p+self.f_a_p*(self.c_ped-2))
		4273	dif = self.n_Muestras_p-(self.nro_key_p+self.f_a_p*(self.c_ped-2))
		4274	self.c_ped = 1
		4275	##tmp = j
		4276	##print("tmp else",tmp)
		4277	self.nro_key_p= self.f_a_p-dif
		4278	iterador = self.nro_key_p
		4279	print("iterador else",iterador)
		4280	if diff_utc >1:
		4281	print("FALTAN DATOS DEL PEDESTAL AMIGO WAIT")
		4282	sleep(1)
		4283	#aqui no hay DATA pero creo el nro_key_p "como si existiera" y reinicio el c_ped
		4284	print("continua bro")
		4285	print("self.c_ped",self.c_ped)
		4286	dif = self.n_Muestras_p-(self.nro_key_p+self.f_a_p*(self.c_ped-2))
		4287	self.c_ped = 1
		4288	self.nro_key_p = self.f_a_p-dif
		4289	#
		4290	iterador = self.nro_key_p
		4291	print("iterador else",iterador)
		4292	#
		4293	utc_ped_nofile = self.gettimeutcfromDirFilename(path=self.path_ped,file=self.list_pedestal[self.nro_file-1])+1
		4294	#iterador = self.nro_key_p #no hay iterador pero importa el nro_key_p
		4295	flag_NOPedfile= True
		4296	print("q paso")
		4297
		4298	if diff_utc <0:
		4299	print("No hay mas archivos")
		4300	self.angulo_adq = numpy.ones(self.blocksPerfile)*numpy.nan
		4301	else:
		4302	print("NO EXISTE DATA - IMAGINA")
		4303	dif = self.n_Muestras_p-(self.nro_key_p+self.f_a_p*(self.c_ped-2))
		4304	self.c_ped = 1
		4305	##tmp = j
		4306	##print("tmp else",tmp)
		4307	self.nro_key_p= self.f_a_p-dif
4198		4308
4199	#iterador = self.nro_key_p +self.f_a_p*(j-tmp)
4200	iterador = self.nro_key_p +self.f_a_p*self.c_ped
4201	self.c_ped = self.c_ped +1
4202		4309
4203	print("iterador------------->",iterador)
4204	if iterador < self.n_Muestras_p:
4205	self.nro_file = self.nro_file
4206	else:
4207	self.nro_file = self.nro_file+1
4208	print("PRUEBA-------------")
4209	utc_ped_setnext=self.gettimeutcfromDirFilename(path=self.path_ped,file=self.list_pedestal[self.nro_file])
4210	utc_adq_setnext=dataOut.utctime
4211	print("utc_pedestal",utc_ped_setnext)
4212	print("utc_adq",utc_adq_setnext)
4213
4214	print("self.c_ped",self.c_ped)
4215	#dif = self.blocksPerfile-(self.nro_key_p+self.f_a_p*(self.c_ped-2))
4216	dif = self.n_Muestras_p-(self.nro_key_p+self.f_a_p*(self.c_ped-2))
4217
4218	self.c_ped = 1
4219	##tmp = j
4220	##print("tmp else",tmp)
4221	self.nro_key_p= self.f_a_p-dif
4222	iterador = self.nro_key_p
4223	print("iterador else",iterador)
4224	#self.c_ped = self.c_ped +1
4225		4310
4226	print("nro_file",self.nro_file)	4311
4227	#print("tmp",tmp)	4312	if flag_NOPedfile == False:
4228	try:	4313	# AQUI NECESITO DOS COSAS "nro_file" y el "iterador"
4229	ff_pedestal = self.list_pedestal[self.nro_file]	4314	ff_pedestal = self.list_pedestal[self.nro_file]
4230	print("ff_pedestal",ff_pedestal)	4315	print("ff_pedestal",ff_pedestal)
4231	except:	4316	angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azi_pos")
4232	print("############# EXCEPCION ######################")	4317	save_angulo = angulo[iterador]
4233	return numpy.ones(self.blocksPerfile)*numpy.nan	4318	self.angulo_adq[j]= angulo[iterador]
		4319	else:
		4320	print(" DATA SIMULADA------------------------------------",j)
		4321	self.angulo_adq[j]= numpy.nan
		4322	cont_NOPedFile = cont_NOPedFile + 1
		4323	print("contador:",cont_NOPedFile)
		4324	sleep(2)
		4325	if cont_NOPedFile==self.f_a_p:
		4326	cont_NOPedFile = 0
		4327	utc_ped_nofile = utc_ped_nofile+1
		4328	utc_ped_setnext= self.gettimeutcfromDirFilename(path=self.path_ped,file=self.list_pedestal[self.nro_file])
		4329	if utc_ped_setnext-utc_ped_nofile==1:
		4330	fladata_pg_NOPedfile=False
		4331
4234		4332
4235	#angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azimuth")
4236	angulo = self.getDatavaluefromDirFilename(path=self.path_ped,file=ff_pedestal,value="azi_pos")
4237		4333
4238	self.angulo_adq[j]= angulo[iterador]
4239		4334
4240	return self.angulo_adq	4335	return self.angulo_adq
4241		4336
@@ -4298,32 +4393,68 class PedestalInformation(Operation):
4298		4393
4299	#def run(self, dataOut,path_ped,path_adq,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online):	4394	#def run(self, dataOut,path_ped,path_adq,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online):
4300	def run(self, dataOut,path_ped,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online):	4395	def run(self, dataOut,path_ped,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online):
		4396	if not self.isConfig:
		4397	self.setup(dataOut, path_ped,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online)
		4398	self.__dataReady = True
		4399	self.isConfig = True
		4400	#print("config TRUE")
		4401
		4402	utc_adq = dataOut.utctime
		4403	print("utc_adq---------------",utc_adq)
		4404	list_pedestal = self.list_pedestal
		4405	angulo = self.getAnguloProfile(utc_adq=utc_adq,list_pedestal=list_pedestal)
		4406	print("angulo**********",angulo)
		4407	dataOut.flagNoData = False
		4408	if numpy.isnan(angulo):
		4409	dataOut.flagNoData = True
		4410	return dataOut
		4411	dataOut.azimuth = angulo
		4412	return dataOut
4301		4413
		4414	'''
4302	if not self.isConfig:	4415	if not self.isConfig:
4303	print("######################SETUP#########################################")	4416	print("######################SETUP#########################################")
		4417	print("CONFIG dataOut.utctime---",dataOut.utctime)
4304	#self.setup( dataOut, path_ped,path_adq,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online)	4418	#self.setup( dataOut, path_ped,path_adq,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online)
4305	self.setup( dataOut, path_ped,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online)	4419	self.setup( dataOut, path_ped,t_Interval_p,n_Muestras_p,blocksPerfile,f_a_p,online)
4306	self.isConfig = True	4420	print("aqui bro")
		4421	if self.nro_file>=0:
		4422	self.isConfig = True
		4423	else:
		4424	self.isConfig = False
		4425	print("RUN---dataOut.utctime---",dataOut.utctime)
4307		4426
4308	dataOut.flagNoData = True	4427	dataOut.flagNoData = True
4309	###print("profIndex",self.__profIndex)	4428	#sleep(0.5)
4310		4429	#print("--------------------------------------------------------------------------------")
		4430	#print("profIndex",self.__profIndex)
		4431	'''
		4432	'''
4311	if self.__profIndex==0:	4433	if self.__profIndex==0:
4312	angulo_adq = self.setNextFileP(dataOut)	4434	angulo_adq = self.setNextFileP(dataOut)
4313	dataOut.azimuth = angulo_adq	4435	dataOut.azimuth = angulo_adq
4314	~~######~~print("TIEMPO:",dataOut.utctime)	4436	print("TIEMPO--------:",dataOut.utctime)
4315	##print("####################################################################")	4437	print("####################################################################")
4316	######print("angulos",dataOut.azimuth,len(dataOut.azimuth))	4438	####print("angulos",dataOut.azimuth,len(dataOut.azimuth))
4317	self.__dataReady = True	4439	self.__dataReady = True
		4440	'''
		4441	'''
4318	self.__profIndex += 1	4442	self.__profIndex += 1
4319	####print("TIEMPO_bucle:",dataOut.utctime)	4443	####print("TIEMPO_bucle:",dataOut.utctime)
4320	~~###~~#print("profIndex",self.__profIndex)	4444	#print("profIndex",self.__profIndex)
4321	if self.__profIndex== blocksPerfile:	4445	if self.__profIndex== blocksPerfile:
		4446	self.__dataReady = True
4322	self.__profIndex = 0	4447	self.__profIndex = 0
4323	if self.__dataReady:	4448	if self.__dataReady:
4324	#print(self.__profIndex,dataOut.azimuth[:10])	4449	#print(self.__profIndex,dataOut.azimuth[:10])
		4450	angulo_adq = self.setNextFileP(dataOut)
		4451	dataOut.azimuth = angulo_adq
		4452	print("TIEMPO--------:",dataOut.utctime)
		4453	print("####################################################################")
		4454	####print("angulos",dataOut.azimuth,len(dataOut.azimuth))
4325	dataOut.flagNoData = False	4455	dataOut.flagNoData = False
4326	return dataOut	4456	'''
		4457	#return dataOut
4327		4458
4328		4459
4329	class Block360(Operation):	4460	class Block360(Operation):
@@ -4392,9 +4523,10 class Block360(Operation):
4392	#print(data.azimuth[self.index])	4523	#print(data.azimuth[self.index])
4393	#print("index",self.index)	4524	#print("index",self.index)
4394		4525
4395	self.__buffer2[self.__profIndex] = data.azimuth[self.index]	4526	#####self.__buffer2[self.__profIndex] = data.azimuth[self.index]
		4527	self.__buffer2[self.__profIndex] = data.azimuth
4396	#print("q pasa")	4528	#print("q pasa")
4397	self.index+=1	4529	#####self.index+=1
4398	#print("index",self.index,data.azimuth[:10])	4530	#print("index",self.index,data.azimuth[:10])
4399	self.__profIndex += 1	4531	self.__profIndex += 1
4400	return #················· Remove DC···································	4532	return #················· Remove DC···································
@@ -4424,7 +4556,7 class Block360(Operation):
4424	data_p = None	4556	data_p = None
4425	#print("dataOu",dataOut.dataPP_POW)	4557	#print("dataOu",dataOut.dataPP_POW)
4426	self.putData(data=dataOut,mode = self.mode)	4558	self.putData(data=dataOut,mode = self.mode)
4427	#print("profIndex",self.__profIndex)	4559	##### print("profIndex",self.__profIndex)
4428	if self.__profIndex == self.n:	4560	if self.__profIndex == self.n:
4429	data_360,n,data_p = self.pushData(data=dataOut)	4561	data_360,n,data_p = self.pushData(data=dataOut)
4430	self.__dataReady = True	4562	self.__dataReady = True
@@ -4441,32 +4573,35 class Block360(Operation):
4441	if data_360 is None:	4573	if data_360 is None:
4442	return None, None,None	4574	return None, None,None
4443		4575
		4576
4444	avgdatatime = self.__initime	4577	avgdatatime = self.__initime
		4578	if self.n==1:
		4579	avgdatatime = datatime
4445	deltatime = datatime - self.__lastdatatime	4580	deltatime = datatime - self.__lastdatatime
4446	self.__initime = datatime	4581	self.__initime = datatime
4447	#print(data_360.shape,avgdatatime,data_p.shape)	4582	#print(data_360.shape,avgdatatime,data_p.shape)
4448	return data_360,avgdatatime,data_p	4583	return data_360,avgdatatime,data_p
4449		4584
4450	def run(self, dataOut,n = None,mode=None,**kwargs):	4585	def run(self, dataOut,n = None,mode=None,**kwargs):
4451	~~###~~#print("BLOCK 360 HERE WE GO MOMENTOS")	4586	#print("BLOCK 360 HERE WE GO MOMENTOS")
4452	if not self.isConfig:	4587	if not self.isConfig:
4453	self.setup(dataOut = dataOut, n = n ,mode= mode ,**kwargs)	4588	self.setup(dataOut = dataOut, n = n ,mode= mode ,**kwargs)
4454	self.index = 0	4589	####self.index = 0
4455	#print("comova",self.isConfig)	4590	#print("comova",self.isConfig)
4456	self.isConfig = True	4591	self.isConfig = True
4457	if self.index==dataOut.azimuth.shape[0]:	4592	####if self.index==dataOut.azimuth.shape[0]:
4458	self.index=0	4593	#### self.index=0
4459	data_360, avgdatatime,data_p = self.blockOp(dataOut, dataOut.utctime)	4594	data_360, avgdatatime,data_p = self.blockOp(dataOut, dataOut.utctime)
4460	dataOut.flagNoData = True	4595	dataOut.flagNoData = True
4461		4596
4462	if self.__dataReady:	4597	if self.__dataReady:
4463	dataOut.data_360 = data_360 # S	4598	dataOut.data_360 = data_360 # S
4464	##print("---------------------------------------------------------------------------------")	4599	##print("---------------------------------------------------------------------------------")
4465	##print("---------------------------DATAREADY---------------------------------------------")	4600	#####print("---------------------------DATAREADY---------------------------------------------")
4466	##print("---------------------------------------------------------------------------------")	4601	##print("---------------------------------------------------------------------------------")
4467	##print("data_360",dataOut.data_360.shape)	4602	##print("data_360",dataOut.data_360.shape)
4468	dataOut.data_azi = data_p	4603	dataOut.data_azi = data_p
4469	##print("azi: ",dataOut.data_azi)	4604	#####print("azi: ",dataOut.data_azi)
4470	#print("jroproc_parameters",data_p[0],data_p[-1])#,data_360.shape,avgdatatime)	4605	#print("jroproc_parameters",data_p[0],data_p[-1])#,data_360.shape,avgdatatime)
4471	dataOut.utctime = avgdatatime	4606	dataOut.utctime = avgdatatime
4472	dataOut.flagNoData = False	4607	dataOut.flagNoData = False

schainpy/model/utils/scp.py +1 -1

             class SCPException(Exception):
                 """SCP exception class"""
-                pass
  No newline at end of file
+                pass

schainpy/scripts/basic_proc00001.py +89 -23

@@ -22,19 +22,38 from schainpy.controller import Project
22	######## n = 1/(V*IPP) #############################	22	######## n = 1/(V*IPP) #############################
23	######## VELOCIDAD DEL PEDESTAL ######################	23	######## VELOCIDAD DEL PEDESTAL ######################
24	print("SETUP- RADAR METEOROLOGICO")	24	print("SETUP- RADAR METEOROLOGICO")
25	V = 10	25	V = 6
26	mode = 1	26	mode = 1
27	#path = '/DATA_RM/23/6v'	27	#--------------------------PATH -----------------------------
28	#~~###~~path = '/DATA_RM/~~TEST_INTEGRACION_2M~~'	28	# path = '/DATA_RM/23/6v'
29	#path = '/DATA_RM/TEST_~~19OCTUBRE/10MHZ~~'	29	# path = '/DATA_RM/TEST_INTEGRACION_2M'
30	path = '/DATA_RM/~~WR_20_OCT~~'	30	# path = '/DATA_RM/TEST_19OCTUBRE/10MHZ'
		31	# path = '/DATA_RM/WR_20_OCT'
31	#### path_ped='/DATA_RM/TEST_PEDESTAL/P20211012-082745'	32	#### path_ped='/DATA_RM/TEST_PEDESTAL/P20211012-082745'
32	####path_ped='/DATA_RM/TEST_PEDESTAL/P20211019-192244'	33	####path_ped='/DATA_RM/TEST_PEDESTAL/P20211019-192244'
33	figpath_pp = "/home/soporte/Pictures/TEST_PP"	34	#path = '/DATA_RM/WR_POT_09_1'
34	figpath_spec = "/home/soporte/Pictures/TEST_MOM"	35	#path = '/DATA_RM/WR_POT_09_2'
		36	#path ="/DATA_RM/10"
		37	#path ="/DATA_RM/11"
		38	#path ="/DATA_RM/DRONE/2MHZ_6V"
		39	path ="/DATA_RM/DRONE/2MHZ_10V_ELEVACION"
		40	path ="/DATA_RM/DRONE/2MHZ_6V_AZIMUTH"
		41	#-------------------------PATH-PLOTEO------------------------------------
		42	#figpath_pp = "/home/soporte/Pictures/TEST_PP"
		43	#figpath_spec = "/home/soporte/Pictures/TEST_MOM_15"
		44	#figpath_spec = "/home/soporte/Pictures/11_"
		45	figpath_spec = "/home/soporte/Pictures/23NOV"
		46	figpath_spec = "/home/soporte/Pictures/23NOV_TEST_AZI"
		47	figpath_pp = "/home/soporte/Pictures/TEST_POT"
		48	figpath_pp = "/home/soporte/Pictures/TEST_POT_10_"
		49
		50
		51	#--------------------------OPCIONES-----------------------------------
35	plot = 0	52	plot = 0
36	integration = 1	53	integration = 0
37	save = 0	54	save = 0
		55	plot_spec = 1
		56	#-------------------------------------------------------------------------
38	if save == 1:	57	if save == 1:
39	if mode==0:	58	if mode==0:
40	path_save = '/DATA_RM/TEST_HDF5_PP_23/6v'	59	path_save = '/DATA_RM/TEST_HDF5_PP_23/6v'
@@ -68,16 +87,16 if save ==1:
68		87
69	print("* Integracion de datos :",integration)	88	print("* Integracion de datos :",integration)
70		89
71	time.sleep(15)	90	time.sleep(3)
72	#remotefolder = "/home/wmaster/graficos"	91	#remotefolder = "/home/wmaster/graficos"
73	#######################################################################	92	#######################################################################
74	################# RANGO DE PLOTEO######################################	93	################# RANGO DE PLOTEO######################################
75	dBmin = '1'	94	dBmin = '8'
76	dBmax = '65'	95	dBmax = '35'
77	xmin = '13.2'	96	xmin = '12.1'#17.1,17.5
78	xmax = '1~~3.5~~'	97	xmax = '12.2'#17.2,17.8
79	ymin = '0'	98	ymin = '0'#### PONER A 0
80	ymax = '60'	99	ymax = '8'#### PONER A 8
81	#######################################################################	100	#######################################################################
82	########################FECHA##########################################	101	########################FECHA##########################################
83	str = datetime.date.today()	102	str = datetime.date.today()
@@ -96,13 +115,13 controllerObj.setup(id = '191', name='Test_USRP', description=desc)
96	#######################################################################	115	#######################################################################
97	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',	116	readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
98	path=path,	117	path=path,
99	startDate="2021/01/01",#today,	118	startDate="2021/11/23",#today,
100	endDate="2021/12/30",#today,	119	endDate="2021/12/23",#today,
101	startTime='00:00:00',	120	startTime='15:00:00',#'17:39:25',
102	endTime='~~23:59~~:59',	121	endTime='16:00:59',#23:59:59',
103	delay=0,	122	delay=0,
104	#set=0,	123	#set=0,
105	online=0,	124	online=1,
106	walk=1,	125	walk=1,
107	ippKm = 60)	126	ippKm = 60)
108		127
@@ -110,6 +129,22 opObj11 = readUnitConfObj.addOperation(name='printInfo')
110		129
111	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())	130	procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
112		131
		132	opObj11 = procUnitConfObjA.addOperation(name='setH0')
		133	opObj11.addParameter(name='h0', value='-3.3', format='float')
		134	#-2.8 5MHZ
		135	#-3.3 2MHZ
		136
		137	###opObj11 = procUnitConfObjA.addOperation(name='selectChannels')
		138	###opObj11.addParameter(name='channelList', value='[1]',format='intList')
		139
		140
		141	opObj11 = procUnitConfObjA.addOperation(name='selectHeights')
		142	opObj11.addParameter(name='minIndex', value='1', format='int')
		143	###opObj11.addParameter(name='maxIndex', value='10000', format='int')
		144	opObj11.addParameter(name='maxIndex', value='160', format='int')
		145	## 400 PARA 5MHZ 12KM 2000
		146	## 160 ṔARA 2MHZ 12KM 800
		147
113	if mode ==0:	148	if mode ==0:
114	####################### METODO PULSE PAIR ######################################################################	149	####################### METODO PULSE PAIR ######################################################################
115	opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')	150	opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
@@ -123,6 +158,8 if mode ==0:
123	opObj11.addParameter(name='colormap', value='jet')	158	opObj11.addParameter(name='colormap', value='jet')
124	opObj11.addParameter(name='xmin', value=xmin)	159	opObj11.addParameter(name='xmin', value=xmin)
125	opObj11.addParameter(name='xmax', value=xmax)	160	opObj11.addParameter(name='xmax', value=xmax)
		161	#opObj11.addParameter(name='ymin', value=ymin)
		162	#opObj11.addParameter(name='ymax', value=ymax)
126	opObj11.addParameter(name='zmin', value=dBmin)	163	opObj11.addParameter(name='zmin', value=dBmin)
127	opObj11.addParameter(name='zmax', value=dBmax)	164	opObj11.addParameter(name='zmax', value=dBmax)
128	opObj11.addParameter(name='save', value=figpath_pp)	165	opObj11.addParameter(name='save', value=figpath_pp)
@@ -167,20 +204,49 else:
167	procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())	204	procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
168	procUnitConfObjB.addParameter(name='nFFTPoints', value=n, format='int')	205	procUnitConfObjB.addParameter(name='nFFTPoints', value=n, format='int')
169	procUnitConfObjB.addParameter(name='nProfiles' , value=n, format='int')	206	procUnitConfObjB.addParameter(name='nProfiles' , value=n, format='int')
		207	if plot_spec==1:
		208	'''
		209	opObj11 = procUnitConfObjB.addOperation(name='SpectraPlot', optype='external')
		210	#.addParameter(name='id', value='2', format='int')
		211	opObj11.addParameter(name='wintitle', value='SpectraPlot', format='str')
		212	opObj11.addParameter(name='channelList', value='1', format='intlist')
		213	#opObj11.addParameter(name='xmin', value=xmin)
		214	#opObj11.addParameter(name='xmax', value=xmax)
		215	opObj11.addParameter(name='ymin', value=ymin)
		216	opObj11.addParameter(name='ymax', value=ymax)
		217	opObj11.addParameter(name='zmin', value=dBmin, format='int')
		218	opObj11.addParameter(name='zmax', value=dBmax, format='int')
		219	opObj11.addParameter(name='save', value=figpath_spec)
		220	opObj11.addParameter(name='showprofile', value=0)
		221	opObj11.addParameter(name='save_period', value=1)
		222	'''
		223	opObj11 = procUnitConfObjB.addOperation(name='RTIPlot', optype='external')
		224	#.addParameter(name='id', value='2', format='int')
		225	opObj11.addParameter(name='wintitle', value='RTIPlot', format='str')
		226	opObj11.addParameter(name='channelList', value='1', format='intlist')
		227	opObj11.addParameter(name='xmin', value=xmin)
		228	opObj11.addParameter(name='xmax', value=xmax)
		229	opObj11.addParameter(name='ymin', value=ymin)
		230	opObj11.addParameter(name='ymax', value=ymax)
		231	opObj11.addParameter(name='zmin', value=dBmin, format='int')
		232	opObj11.addParameter(name='zmax', value=dBmax, format='int')
		233	opObj11.addParameter(name='save', value=figpath_spec)
		234	opObj11.addParameter(name='showprofile', value=0)
		235	opObj11.addParameter(name='save_period', value=1)
170		236
171	procUnitConfObjC = controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjB.getId())	237	procUnitConfObjC = controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjB.getId())
172	procUnitConfObjC.addOperation(name='SpectralMoments')	238	procUnitConfObjC.addOperation(name='SpectralMoments')
173	if plot==1:	239	if plot==1:
174	dBmin = '1'
175	dBmax = '65'
176	opObj11 = procUnitConfObjC.addOperation(name='PowerPlot',optype='external')	240	opObj11 = procUnitConfObjC.addOperation(name='PowerPlot',optype='external')
177	opObj11.addParameter(name='xmin', value=xmin)	241	opObj11.addParameter(name='xmin', value=xmin)
178	opObj11.addParameter(name='xmax', value=xmax)	242	opObj11.addParameter(name='xmax', value=xmax)
		243	opObj11.addParameter(name='ymin', value=ymin)
		244	opObj11.addParameter(name='ymax', value=ymax)
179	opObj11.addParameter(name='zmin', value=dBmin)	245	opObj11.addParameter(name='zmin', value=dBmin)
180	opObj11.addParameter(name='zmax', value=dBmax)	246	opObj11.addParameter(name='zmax', value=dBmax)
181	opObj11.addParameter(name='save', value=figpath_spec)	247	opObj11.addParameter(name='save', value=figpath_spec)
182	opObj11.addParameter(name='showprofile', value=0)	248	opObj11.addParameter(name='showprofile', value=0)
183	opObj11.addParameter(name='save_period', value=10)	249	opObj11.addParameter(name='save_period', value=1)
184		250
185	if save==1:	251	if save==1:
186	opObj10 = procUnitConfObjC.addOperation(name='HDFWriter')	252	opObj10 = procUnitConfObjC.addOperation(name='HDFWriter')

schainpy/scripts/basic_proc00002.py +89 -37

             ######## n   = 1/(V*IPP) #############################
             ######## VELOCIDAD DEL PEDESTAL ######################
             print("SETUP- RADAR METEOROLOGICO")
+            V       = 6
             V       = 10
-            mode    = 0
+            mode    = 1
-            #path    = '/DATA_RM/23/6v'
+            #--------------------------PATH -----------------------------
-            #path    = '/DATA_RM/TEST_INTEGRACION_2M'
+            # path    = '/DATA_RM/23/6v'
-            path    = '/DATA_RM/WR_20_OCT'
+            # 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'
+            path ="/DATA_RM/10"
+            # path  = '/DATA_RM/WR_POT_09_1'
+            #path ="/DATA_RM/11"
+            #path ="/DATA_RM/11"
+            #-------------------------PATH-PLOTEO------------------------------------
+            #figpath_pp  = "/home/soporte/Pictures/TEST_PP"
+            #figpath_spec = "/home/soporte/Pictures/TEST_MOM"
+            figpath_spec = "/home/soporte/Pictures/ppi"
+            figpath_ppi  = "/home/soporte/Pictures/ppi_SPEC_10DIC"
+            figpath_ppi_pp  = "/home/soporte/Pictures/ppi_2_10"
+            figpath_pp  = "/home/soporte/Pictures/TEST_POT"
+            figpath_pp  = "/home/soporte/Pictures/TEST_POT2"
             #path_ped='/DATA_RM/TEST_PEDESTAL/P20211012-082745'
-            path_ped='/DATA_RM/TEST_PEDESTAL/P20211020-131248'
+            #path_ped='/DATA_RM/TEST_PEDESTAL/P20211020-131248'
+            #path_ped='/DATA_RM/TEST_PEDESTAL/P20211110-171003'
+            #path_ped='/DATA_RM/TEST_PEDESTAL/P20211111-173856'
+            path_ped  =  "/DATA_RM/TEST_PEDESTAL/P20211110-171003"
-            figpath_pp  = "/home/soporte/Pictures/TEST_PP"
+            figpath_pp  = "/home/soporte/Pictures/25TEST_PP"
             figpath_mom = "/home/soporte/Pictures/TEST_MOM"
+            #--------------------------OPCIONES-----------------------------------
             plot        = 0
             integration = 1
-            save        = 0
+            save        = 1
             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'
+                    path_save = '/DATA_RM/TEST_EMPTHY'
                 else:
                     path_save = '/DATA_RM/TEST_HDF5_SPEC_23_V2/6v'
+                    path_save = '/DATA_RM/TEST_EMPTHY_SPEC'
+                    path_save = '/DATA_RM/LAST_TEST_16_VACIO3'
+                    path_save = '/DATA_RM/LAST_TEST_30_360'
             print("* PATH data ADQ        :", path)
             print("* Velocidad Pedestal   :",V,"°/seg")
             print("* Integracion de datos :",integration)
-            time.sleep(5)
+            time.sleep(3)
             #remotefolder = "/home/wmaster/graficos"
             #######################################################################
             ################# RANGO DE PLOTEO######################################
             dBmin = '1'
             dBmax = '85'
-            xmin = '15'
+            xmin = '17'
-            xmax = '15.25'
+            xmax = '17.25'
             ymin = '0'
             ymax = '600'
             #######################################################################
             #######################################################################
             readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
                                                         path=path,
-                                                        startDate="2021/01/01",#today,
+                                                        startDate="2021/11/10",#today,
                                                         endDate="2021/12/30",#today,
-                                                        startTime='00:00:00',
+                                                        startTime='17:10:25',
                                                         endTime='23:59:59',
                                                         delay=0,
                                                         #set=0,
             opObj11 = readUnitConfObj.addOperation(name='printInfo')
-            procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
+            procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc',inputId=readUnitConfObj.getId())
+            #opObj11 = procUnitConfObjA.addOperation(name='setH0')
+            #opObj11.addParameter(name='h0', value='-2.8', format='float')
+            opObj11 = procUnitConfObjA.addOperation(name='selectHeights')
+            opObj11.addParameter(name='minIndex', value='1', format='int')
+            #    opObj11.addParameter(name='maxIndex', value='10000', format='int')
+            opObj11.addParameter(name='maxIndex', value='400', format='int')
             if mode ==0:
                 ####################### METODO PULSE PAIR ######################################################################
                     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
+                    V=V
-                    blocksPerfile=360
+                    blocksPerfile=100
                     print("* Velocidad del Pedestal:",V)
                     tmp_blocksPerfile = 100
                     f_a_p= int(tmp_blocksPerfile/V)
                     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')
+                    opObj11.addParameter(name='save', value=figpath_ppi_pp)
+                    opObj11.addParameter(name='save_period', value=1)
+                    '''
+                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,azimuth,utctime',format='list')#,format='list'
             else:
                 ####################### METODO SPECTROS ######################################################################
                     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
+                   V=V
-                   blocksPerfile=360
+                   blocksPerfile=100
                    print("* Velocidad del Pedestal:",V)
                    tmp_blocksPerfile = 100
                    f_a_p= int(tmp_blocksPerfile/V)
                    opObj11.addParameter(name='online', value='0', format='int')
                    opObj11 = procUnitConfObjC.addOperation(name='Block360')
-                   opObj11.addParameter(name='n', value='10', format='int')
+                   opObj11.addParameter(name='n', value='1', 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')
+                   opObj11.addParameter(name='save', value=figpath_ppi)
+                   opObj11.addParameter(name='save_period', value=1)
+                   '''
+                if save==1:
+                   opObj10 = procUnitConfObjC.addOperation(name='HDFWriter')
+                   opObj10.addParameter(name='path',value=path_save)
+                   opObj10.addParameter(name='mode',value='weather')
+                   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,azimuth,utctime',format='list')#,format='list'
+                   opObj10.addParameter(name='dataList',value='data_360,data_azi,utctime',format='list')#,format='list'
+                   '''
                    # 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')
+                   opObj11.addParameter(name='save', value=figpath_ppi)
+                   opObj11.addParameter(name='save_period', value=1)
+                   '''
             controllerObj.start()

schainpy/scripts/check_pedestal.py +4 -1

             #path_ped = '/DATA_RM/TEST_PEDESTAL/P20211012-082745'
             path_ped = '/DATA_RM/TEST_PEDESTAL/P20211020-131248'
+            path_ped = '/DATA_RM/TEST_PEDESTAL/P20211110-171003'
+            path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173856'
+            #path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173409'
             # Metodo para verificar numero
             def isNumber(str):
                 try:
             plt.plot(t,ele_pos,label='ELEVATION_POS',color='red')#*10
             #ax.set_xlim(0, 9)
-            ax.set_ylim(-5, 400)
+            ax.set_ylim(-5, 20)
             plt.ylabel("Azimuth Position")
             plt.xlabel("Muestra")
             plt.title('Azimuth Position vs Muestra ', fontsize=20)

General Comments 0

Write
Preview

You need to be logged in to leave comments. Login now

No TODOs yet

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages