schain Commit - r1355:cfdf8000eba5 · Jicamarca Repository

ultimas actualizaciones y pruebas

avaldez -

r1355:cfdf8000eba5 test_WRM

parent child

Context file:

r1355:cfdf8000eba5

schainpy/scripts/readFileconfig.py created 644 +55 0

@@ -0,0 +1,55
	1	print("LECTURA DE ARCHIVOS DE CONFIGURACION")
	2	class ReadfileWR():
	3	def __init__(self,filename):
	4	f = open(filename, "r")
	5	i=0
	6	self.dict={'resolution':None,'vel_ped_azi':None,'pos_ped_ele':None,'ipp':None,'n':None,'len_ped':None,\
	7	't_s_ped':None,'t_f_ped':None,'b_f_adq':None,'t_f_adq':None,'mode':None}
	8	while(True):
	9	##print(i)
	10	linea = f.readline()
	11	if i==4:
	12	resolution=float(linea)
	13	self.dict['resolution']=resolution
	14	if i==6:
	15	vel_pedestal_a=float(linea)
	16	self.dict['vel_ped_azi']=vel_pedestal_a
	17	if i==8:
	18	vel_pedestal_e=float(linea)
	19	self.dict['pos_ped_ele']=vel_pedestal_e
	20	if i==10:
	21	ipp = float(linea)
	22	self.dict['ipp']= round(ipp,5)
	23	if i==12:
	24	n = float(linea)
	25	self.dict['n']= n
	26	if i==14:
	27	len_pedestal= float(linea)
	28	self.dict['len_ped']= len_pedestal
	29	if i==16:
	30	time_x_sample_ped=float(linea)
	31	self.dict['t_s_ped']= time_x_sample_ped
	32	if i==18:
	33	time_x_file_ped = float(linea)
	34	self.dict['t_f_ped']= time_x_file_ped
	35	if i==20:
	36	bloques_x_file_adq= float(linea)
	37	self.dict['b_f_adq']=bloques_x_file_adq
	38	if i==22:
	39	time_x_file_adq = float(linea)
	40	self.dict['t_f_adq'] = time_x_file_adq
	41	if i==24:
	42	mode= int(linea)
	43	self.dict['mode'] = mode
	44	#print(linea)
	45	if not linea:
	46	break
	47	i+=1
	48	f.close()
	49	def getDict(self):
	50	return self.dict
	51
	52
	53	filename= "/home/developer/Downloads/config_WR.txt"
	54	dict= ReadfileWR(filename).getDict()
	55	print(dict)

schainpy/model/graphics/jroplot_parameters.py +12 -10

                     return data, meta
                 def plot(self):
+                    thisDatetime = datetime.datetime.utcfromtimestamp(self.data.times[-1])
+                    print("thisDatetime",thisDatetime)
                     stoprange = float(33*1.5)
                     rangestep = float(0.15)
                     r      = numpy.arange(0, stoprange, rangestep)
                     step   = (360/(res*tmp_v.shape[0]))
                     mode   = 1
                     if mode==0:
-                        print("self.ini",self.ini)
+                        #print("self.ini",self.ini)
                         val             = numpy.mean(tmp_v[:,0])
                         self.len_azi    = len(tmp_z)
                         ones            = numpy.ones([(360-tmp_v.shape[0]),tmp_v.shape[1]])*val
                         self.ini            = self.ini+1
                     if mode==1:
-                        print("self.ini",self.ini)
+                        #print("self.ini",self.ini)
                         if self.ini==0:
                             res             = 1
                             step            = (360/(res*tmp_v.shape[0]))
                             if self.ini==step-1:
                                 start= tmp_z[0]
                                 end  = tmp_z[-1]
-                                print("start","end",start,end)
+                                #print("start","end",start,end)
-                                print(self.buffer_ini_azi[:80])
+                                ###print(self.buffer_ini_azi[:80])
                             self.ini            = self.ini+1
                         else:
                             start= tmp_z[0]
                             end  = tmp_z[-1]
-                            print("start","end",start,end)
+                            #print("start","end",start,end)
-                            print(self.buffer_ini_azi[:120])
+                            ###print(self.buffer_ini_azi[:120])
                             if step>=2:
                                 if self.flag<step-1:
                             self.buf_azi[len(tmp_z)*(self.flag):len(tmp_z)*(self.flag+1)] = tmp_z
                             self.buffer_ini=self.buf_tmp
                             self.buffer_ini_azi = self.buf_azi
-                            print("--------salida------------")
+                            ##print("--------salida------------")
                             start= tmp_z[0]
                             end  = tmp_z[-1]
-                            print("start","end",start,end)
+                            ##print("start","end",start,end)
-                            print(self.buffer_ini_azi[:120])
+                            ##print(self.buffer_ini_azi[:120])
                             self.ini= self.ini+1
                             self.flag = self.flag +1
                             if self.flag==step:
                     cbar = plt.gcf().colorbar(pm, pad=0.075)
                     caax.set_xlabel('x_range [km]')
                     caax.set_ylabel('y_range [km]')
-                    plt.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')
+                    plt.text(1.0, 1.05, 'azimuth '+str(thisDatetime), transform=caax.transAxes, va='bottom',ha='right')

schainpy/model/proc/jroproc_parameters.py +3 -2

                     self.index            = 0
-                    #print("ELVALOR DE n es:", n)
+                    ##print("ELVALOR DE n es:", n)
                     if n == None:
                         raise ValueError("n should be specified.")
                     if n != None:
-                        if n<2:
+                        if n<1:
+                            print("n should be greater than 2")
                             raise ValueError("n should be greater than 2")
                     self.n       = n

schainpy/scripts/test_sim00011.py +4 -3

                                                         walk=0)#1
             procUnitConfObjA = controllerObj.addProcUnit(datatype='ParametersProc',inputId=readUnitConfObj.getId())
             opObj11 = procUnitConfObjA.addOperation(name='Block360')
-            opObj11.addParameter(name='n', value='10', format='int')
+            opObj11.addParameter(name='n', value='40', format='int')
             opObj11= procUnitConfObjA.addOperation(name='WeatherPlot',optype='other')
-            #opObj11.addParameter(name='save', value=figpath)
+            opObj11.addParameter(name='save', value=figpath)
-            #opObj11.addParameter(name='save_period', value=1)
+            opObj11.addParameter(name='save_period', value=1)
             #opObj11 = procUnitConfObjA.addOperation(name='PowerPlot', optype='other')#PulsepairPowerPlot
             #opObj11 = procUnitConfObjA.addOperation(name='PPSignalPlot', optype='other')

schainpy/scripts/test_sim00012.py +21 -4

             import datetime
             import time
             from schainpy.controller import Project
+            #*************************************************************************
+            #**************************LECTURA config_WR.txt**************************
+            #*************************************************************************
+            from readFileconfig import ReadfileWR
+            filename= "/home/developer/Downloads/config_WR.txt"
+            dict=  ReadfileWR(filename).getDict()
+            FixRCP_IPP  = dict['ipp']*0.15 #equivalencia
+            dataBlocksPerFile= dict['b_f_adq']
+            profilesPerBlock= int(dict['n'])
+            pulsepair  = int(dict['n'])
+            #*************************************************************************
             path    = '/home/developer/Downloads/HDF5_TESTPP2V3'
             figpath = path
             desc            = "Simulator Test"
             controllerObj.setup(id='10',name='Test Simulator',description=desc)
             readUnitConfObj = controllerObj.addReadUnit(datatype='SimulatorReader',
                                                         frequency=9.345e9,
-                                                        FixRCP_IPP= 60,
+                                                        FixRCP_IPP= FixRCP_IPP,
                                                         Tau_0 = 30,
                                                         AcqH0_0=0,
                                                         samples=330,
                                                         delay=0,
                                                         online=0,
                                                         walk=0,
-                                                        profilesPerBlock=625,
+                                                        profilesPerBlock=profilesPerBlock,
-                                                        dataBlocksPerFile=100)#,#nTotalReadFiles=2)
+                                                        dataBlocksPerFile=dataBlocksPerFile)#,#nTotalReadFiles=2)
             #opObj11         = readUnitConfObj.addOperation(name='printInfo')
             procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
             opObj11 = procUnitConfObjA.addOperation(name='PulsePair')
-            opObj11.addParameter(name='n', value='625', format='int')#10
+            opObj11.addParameter(name='n', value=pulsepair, format='int')#10
             procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
             opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')
             opObj10.addParameter(name='path',value=figpath)
             #opObj10.addParameter(name='mode',value=2)
             opObj10.addParameter(name='blocksPerFile',value='100',format='int')
             opObj10.addParameter(name='metadataList',value='utctimeInit,paramInterval,heightList,profileIndex,flagDataAsBlock',format='list')
             opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'
             controllerObj.start()

schainpy/scripts/test_sim0009.py +1 0

             import datetime
             import time
             from schainpy.controller import Project
             path    = '/home/alex/Downloads/NEW_WR2/spc16removeDC'
             figpath = path
             desc            = "Simulator Test"

schainpy/scripts/wr_integradorv2.py +47 -22

@@ -54,23 +54,46 def getDatavaluefromDirFilename(path,file,value):
54	fp.close()	54	fp.close()
55	return array	55	return array
56		56
		57	opt = input ("Ingresa Modo de integracion: ")
		58
		59	if opt==1:
		60	#·········· Velocidad de Pedestal·················
		61	w = input ("Ingresa velocidad de Pedestal: ")
		62	w = 4
		63	w = float(w)
		64	#·········· Resolucion minimo en grados···········
		65	alfa = input ("Ingresa resolucion minima en grados: ")
		66	alfa = 1
		67	alfa = float(alfa)
		68	#·········· IPP del Experimento ··················
		69	IPP = input ("Ingresa el IPP del experimento: ")
		70	IPP = 0.0004
		71	IPP = float(IPP)
		72	#·········· MODE ··················
		73	mode = input ("Ingresa el MODO del experimento Time(1) or F(0): ")
		74	mode = 1
		75	mode = int(mode)
		76
		77	time_Interval_p=0.01
		78	n_perfiles_p = 100
		79	blocksPerFile = 100
		80	tiempo_file_1_adq=25
		81	tiempo_file_1_ped=1
57		82
58	#·········· Velocidad de Pedestal·················	83	else:
59	w = input ("Ingresa velocidad de Pedestal: ")	84	from readFileconfig import ReadfileWR
60	w = 4	85	filename= "/home/developer/Downloads/config_WR.txt"
61	w = float(w)	86	dict= ReadfileWR(filename).getDict()
62	#·········· Resolucion minimo en grados···········	87	mode = dict['mode']
63	alfa = input ("Ingresa resolucion minima en grados: ")	88	w = dict['vel_ped_azi']
64	alfa = 1	89	alfa= dict['resolution']
65	alfa = float(alfa)	90	IPP = dict['ipp']*1e-6
66	#·········· IPP del Experimento ··················	91	time_Interval_p = dict['t_s_ped']
67	IPP = input ("Ingresa el IPP del experimento: ")	92	n_perfiles_p = dict['len_ped']
68	IPP = 0.0004	93	blocksPerFile = int(dict['b_f_adq'])
69	IPP = float(IPP)	94
70	#·········· MODE ··················	95	tiempo_file_1_adq= dict['t_f_adq']
71	mode = input ("Ingresa el MODO del experimento T or F: ")	96	tiempo_file_1_ped= dict['t_f_ped']
72	mode = "T"
73	mode = str(mode)
74		97
75	#·········· Tiempo en generar la resolucion min···	98	#·········· Tiempo en generar la resolucion min···
76	#············ MCU ·· var_ang = w * (var_tiempo)···	99	#············ MCU ·· var_ang = w * (var_tiempo)···
@@ -82,7 +105,7 num_perfiles = int(var_tiempo/IPP)
82	#··········DATA PEDESTAL··························	105	#··········DATA PEDESTAL··························
83	dir_pedestal = "/home/developer/Downloads/Pedestal/P2021093"	106	dir_pedestal = "/home/developer/Downloads/Pedestal/P2021093"
84	#·········· DATA ADQ······························	107	#·········· DATA ADQ······························
85	if mode==~~"T"~~:	108	if mode==1:
86	dir_adq = "/home/developer/Downloads/HDF5_TESTPP2V3/d2021093" # Time domain	109	dir_adq = "/home/developer/Downloads/HDF5_TESTPP2V3/d2021093" # Time domain
87	else:	110	else:
88	dir_adq = "/home/developer/Downloads/hdf5_test/d2021053" # Frequency domain	111	dir_adq = "/home/developer/Downloads/hdf5_test/d2021053" # Frequency domain
@@ -109,8 +132,8 print("utc_adq :",utc_adq)
109	#sys.exit(0)	132	#sys.exit(0)
110		133
111	#·············Relacion: utc_adq (+/-) var_tiempo*nro_file= utc_pedestal	134	#·············Relacion: utc_adq (+/-) var_tiempo*nro_file= utc_pedestal
112	time_Interval_p = 0.01	135	time_Interval_p = time_Interval_p#0.01
113	n_perfiles_p = ~~100~~ #muestras por achivo del pedestal	136	n_perfiles_p = n_perfiles_p #100 muestras por achivo del pedestal
114		137
115		138
116	if utc_adq>utc_pedestal:	139	if utc_adq>utc_pedestal:
@@ -144,14 +167,16 list_adq = getfirstFilefromPath(path=dir_adq ,meta="D",ext=".hdf5")
144		167
145	nro_file = nro_file #10	168	nro_file = nro_file #10
146	nro_key_perfil = nro_key_p	169	nro_key_perfil = nro_key_p
147	blocksPerFile = 100##### aqui se cambia dependiendo de los blqoues por achivo en adq	170	blocksPerFile = blocksPerFile#100##### aqui se cambia dependiendo de los blqoues por achivo en adq
148	wr_path = "/home/developer/Downloads/HDF5_WR/"	171	wr_path = "/home/developer/Downloads/HDF5_WR/"
149	# Lectura de archivos de adquisicion para adicion de azimuth	172	# Lectura de archivos de adquisicion para adicion de azimuth
150	# factor de archivos	173	# factor de archivos
151	#f_a_p = tiempo_file_1_adq/tiempo_file_1_ped=25/1 = 25	174	#f_a_p = tiempo_file_1_adq/tiempo_file_1_ped=25/1 = 25
152	tiempo_file_1_adq=25	175	tiempo_file_1_adq = tiempo_file_1_adq
153	tiempo_file_1_ped=1	176	tiempo_file_1_ped = tiempo_file_1_ped
		177
154	f_a_p= int(tiempo_file_1_adq/tiempo_file_1_ped)	178	f_a_p= int(tiempo_file_1_adq/tiempo_file_1_ped)
		179
155	for thisFile in range(len(list_adq)):	180	for thisFile in range(len(list_adq)):
156	print("thisFileAdq",thisFile)	181	print("thisFileAdq",thisFile)
157	angulo_adq = numpy.zeros(blocksPerFile)	182	angulo_adq = numpy.zeros(blocksPerFile)

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