schain Commit - r1429:96f17ea26f9f · Jicamarca Repository

update

avaldezp -

r1429:96f17ea26f9f

parent child

Context file:

r1429:96f17ea26f9f

Collapse all files

schainpy/model/graphics/jroplot_parameters.py +128 -68

@@ -678,22 +678,23 class WeatherRHIPlot(Plot):
678	list1=[]	678	list1=[]
679	list2=[]	679	list2=[]
680	for i in reversed(range(len(angulos))):	680	for i in reversed(range(len(angulos))):
681	diff_ = angulos[i]-angulos[i-1]	681	if not i==0:#el caso de i=0 evalula el primero de la lista con el ultimo y no es relevante
682	if diff_ >1.5:	682	diff_ = angulos[i]-angulos[i-1]
683	list1.append(i-1)	683	if abs(diff_) >1.5:
684	~~list2~~.append(~~diff_~~)	684	list1.append(i-1)
		685	list2.append(diff_)
685	return list(reversed(list1)),list(reversed(list2))	686	return list(reversed(list1)),list(reversed(list2))
686		687
687	def fixData180(self,list_,ang_):	688	def fixData90(self,list_,ang_):
688	if list_[0]==-1:	689	if list_[0]==-1:
689	vec = numpy.where(ang_<ang_[0])	690	vec = numpy.where(ang_<ang_[0])
690	ang_[vec] = ang_[vec]+180	691	ang_[vec] = ang_[vec]+90
691	return ang_	692	return ang_
692	return ang_	693	return ang_
693		694
694	def fixData180HL(self,angulos):	695	def fixData90HL(self,angulos):
695	vec = numpy.where(angulos>=180)	696	vec = numpy.where(angulos>=90)
696	angulos[vec]=angulos[vec]-180	697	angulos[vec]=angulos[vec]-90
697	return angulos	698	return angulos
698		699
699		700
@@ -704,11 +705,11 class WeatherRHIPlot(Plot):
704	i=None	705	i=None
705	return False,i	706	return False,i
706		707
707	def fixDataComp(self,ang_,list1_,list2_):	708	def fixDataComp(self,ang_,list1_,list2_,tipo_case):
708	size = len(ang_)	709	size = len(ang_)
709	size2 = 0	710	size2 = 0
710	for i in range(len(list2_)):	711	for i in range(len(list2_)):
711	size2=size2+round(list2_[i])-1	712	size2=size2+round(abs(list2_[i]))-1
712	new_size= size+size2	713	new_size= size+size2
713	ang_new = numpy.zeros(new_size)	714	ang_new = numpy.zeros(new_size)
714	ang_new2 = numpy.zeros(new_size)	715	ang_new2 = numpy.zeros(new_size)
@@ -721,23 +722,29 class WeatherRHIPlot(Plot):
721	condition , value = self.search_pos(i,list1_)	722	condition , value = self.search_pos(i,list1_)
722	if condition:	723	if condition:
723	pos = tmp + c + 1	724	pos = tmp + c + 1
724	for k in range(round(list2_[value])-1):	725	for k in range(round(abs(list2_[value]))-1):
725	ang_new[pos+k] = ang_new[pos+k-1]+1	726	if tipo_case==0 or tipo_case==3:#subida
726	ang_new2[pos+k] = ~~numpy~~.~~nan~~	727	ang_new[pos+k] = ang_new[pos+k-1]+1
		728	ang_new2[pos+k] = numpy.nan
		729	elif tipo_case==1 or tipo_case==2:#bajada
		730	ang_new[pos+k] = ang_new[pos+k-1]-1
		731	ang_new2[pos+k] = numpy.nan
		732
727	tmp = pos +k	733	tmp = pos +k
728	c = 0	734	c = 0
729	c=c+1	735	c=c+1
730	return ang_new,ang_new2	736	return ang_new,ang_new2
731		737
732	def globalCheckPED(self,angulos):	738	def globalCheckPED(self,angulos,tipo_case):
733	l1,l2 = self.get2List(angulos)	739	l1,l2 = self.get2List(angulos)
		740	print("l1",l1)
		741	print("l2",l2)
734	if len(l1)>0:	742	if len(l1)>0:
735	angulos2 = self.fixData180(list_=l1,ang_=angulos)	743	#angulos2 = self.fixData90(list_=l1,ang_=angulos)
736	l1,l2 = self.get2List(angulos2)	744	#l1,l2 = self.get2List(angulos2)
737		745	ang1_,ang2_ = self.fixDataComp(ang_=angulos,list1_=l1,list2_=l2,tipo_case=tipo_case)
738	ang1_,~~ang2~~_ = self.fixData~~Comp~~(~~ang_~~=~~angulos2~~,~~list1_~~=l1,~~list2_~~=l2)	746	#ang1_ = self.fixData90HL(ang1_)
739	~~ang1~~_ = self.fixData~~180HL~~(~~ang1~~_)	747	#ang2_ = self.fixData90HL(ang2_)
740	ang2_ = self.fixData180HL(ang2_)
741	else:	748	else:
742	ang1_= angulos	749	ang1_= angulos
743	ang2_= angulos	750	ang2_= angulos
@@ -762,52 +769,112 class WeatherRHIPlot(Plot):
762	data_T[i,:]=numpy.ones(data_T.shape[1])*numpy.nan	769	data_T[i,:]=numpy.ones(data_T.shape[1])*numpy.nan
763	return data_T	770	return data_T
764		771
765	def const_ploteo(self,data_weather,data_ele,step,res):	772	def check_case(self,data_ele,ang_max,ang_min):
		773	start = data_ele[0]
		774	end = data_ele[-1]
		775	number = (end-start)
		776	len_ang=len(data_ele)
		777
		778	if start<end and round(abs(number)+1)>=len_ang:#caso subida
		779	return 0
		780	elif start>end and round(abs(number)+1)>=len_ang:#caso bajada
		781	return 1
		782	elif round(abs(number)+1)<len_ang and data_ele[-2]>data_ele[-1]:# caso BAJADA CAMBIO ANG MAX
		783	return 2
		784	elif round(abs(number)+1)<len_ang and data_ele[-2]<data_ele[-1]:# caso SUBIDA CAMBIO ANG MIN
		785	return 3
		786
		787
		788	def const_ploteo(self,data_weather,data_ele,step,res,ang_max,ang_min):
		789	ang_max= ang_max
		790	ang_min= ang_min
766	if self.ini==0:	791	if self.ini==0:
767	#-------	792	print("**********************************************")
768	n = (180/res)-len(data_ele)	793	print("**********************************************")
		794	print("*************ini************")
		795	print("**********************************************")
		796	print("**********************************************")
		797	print("data_ele",data_ele)
		798	#----------------------------------------------------------
		799	tipo_case = self.check_case(data_ele,ang_max,ang_min)
		800	print("TIPO DE DATA",tipo_case)
769	#--------------------- new -------------------------	801	#--------------------- new -------------------------
770	data_ele_new ,data_ele_old= self.globalCheckPED(data_ele)	802	data_ele_new ,data_ele_old= self.globalCheckPED(data_ele,tipo_case)
771	#------------------------	803	print("data_ele_new",data_ele_new)
772	start = data_ele_new[-1] + res	804	print("data_ele_old",data_ele_old)
773	end = data_ele_new[0] - res	805	#-------------------------CAMBIOS RHI---------------------------------
		806	start= ang_min
		807	end = ang_min
		808	n= (ang_max-ang_min)/res
774	#------ new	809	#------ new
775	self.~~las~~t_data_ele = ~~end~~	810	self.start_data_ele = data_ele_new[0]
776	if start>end:	811	self.end_data_ele = data_ele_new[-1]
777	end = end + 180	812	if tipo_case==0 or tipo_case==3:
		813	n1= round(self.start_data_ele)- start
		814	n2= end - round(self.end_data_ele)
		815	if n1>0:
		816	ele1= numpy.linspace(ang_min+1,self.start_data_ele-1,n1)
		817	ele1_nan= numpy.ones(n1)*numpy.nan
		818	data_ele = numpy.hstack((ele1,data_ele_new))
		819	data_ele_old = numpy.hstack((ele1_nan,data_ele_old))
		820	if n2>0:
		821	ele2= numpy.linspace(self.end_data_ele+1,end,n2)
		822	ele2_nan= numpy.ones(n2)*numpy.nan
		823	data_ele = numpy.hstack((data_ele,ele2))
		824	data_ele_old = numpy.hstack((data_ele_old,ele2_nan))
		825	# RADAR
		826	val_mean = numpy.mean(data_weather[:,-1])
		827	self.val_mean = val_mean
		828	data_weather = self.replaceNAN(data_weather=data_weather,data_ele=data_ele_old,val=self.val_mean)
		829	'''
778	ele_vacia = numpy.linspace(start,end,int(n))	830	ele_vacia = numpy.linspace(start,end,int(n))
779	ele_vacia = numpy.where(ele_vacia>180,ele_vacia-180,ele_vacia)	831
		832
		833	ele_vacia = numpy.where(ele_vacia>ang_max,ele_vacia-ang_max,ele_vacia)
780	data_ele = numpy.hstack((data_ele_new,ele_vacia))	834	data_ele = numpy.hstack((data_ele_new,ele_vacia))
781	# RADAR	835	# RADAR
782	val_mean = numpy.mean(data_weather[:,-1])	836	val_mean = numpy.mean(data_weather[:,-1])
783	self.val_mean = val_mean	837	self.val_mean = val_mean
784	data_weather_cmp = numpy.ones([(~~180~~-data_weather.shape[0]),data_weather.shape[1]])*val_mean	838	data_weather_cmp = numpy.ones([(ang_max-data_weather.shape[0]),data_weather.shape[1]])*val_mean
785	data_weather = self.replaceNAN(data_weather=data_weather,data_ele=data_ele_old,val=self.val_mean)	839	data_weather = self.replaceNAN(data_weather=data_weather,data_ele=data_ele_old,val=self.val_mean)
786	data_weather = numpy.vstack((data_weather,data_weather_cmp))	840	data_weather = numpy.vstack((data_weather,data_weather_cmp))
		841	'''
787	else:	842	else:
788	# azimuth	843	print("**********************************************")
		844	print("**********************************************")
		845	print("**************VARIABLE********************")
		846	print("**********************************************")
		847	print("**********************************************")
		848	#-------------------------CAMBIOS RHI---------------------------------
		849	#---------------------------------------------------------------------
		850	print("INPUT data_ele",data_ele)
789	flag=0	851	flag=0
790	start_ele = self.res_ele[0]	852	start_ele = self.res_ele[0]
		853	tipo_case = self.check_case(data_ele,ang_max,ang_min)
		854	print("TIPO DE DATA",tipo_case)
791	#-----------new------------	855	#-----------new------------
792	data_ele ,data_ele_old= self.globalCheckPED(data_ele)	856	data_ele ,data_ele_old= self.globalCheckPED(data_ele,tipo_case)
		857	print("data_ele_new",data_ele)
		858	print("data_ele_old",data_ele_old)
793	data_weather = self.replaceNAN(data_weather=data_weather,data_ele=data_ele_old,val=self.val_mean)	859	data_weather = self.replaceNAN(data_weather=data_weather,data_ele=data_ele_old,val=self.val_mean)
		860	'''
794	#--------------------------	861	#--------------------------
795	start = data_ele[0]	862	start = data_ele[0]
796	end = data_ele[-1]	863	end = data_ele[-1]
797	self.last_data_ele= end	864	self.last_data_ele= end
798	if start< start_ele:	865	if start< start_ele:
799	start = start +~~180~~	866	start = start +ang_max
800	if end <start_ele:	867	if end <start_ele:
801	end = end +~~180~~	868	end = end +ang_max
802		869
803	pos_ini = int((start-start_ele)/res)	870	pos_ini = int((start-start_ele)/res)
804	len_ele = len(data_ele)	871	len_ele = len(data_ele)
805	if (~~180~~-pos_ini)<len_ele:	872	if (ang_max-pos_ini)<len_ele:
806	if pos_ini+1==~~180~~:	873	if pos_ini+1==ang_max:
807	pos_ini=0	874	pos_ini=0
808	else:	875	else:
809	flag=1	876	flag=1
810	dif= ~~180~~-pos_ini	877	dif= ang_max-pos_ini
811	comp= len_ele-dif	878	comp= len_ele-dif
812	#-----------------	879	#-----------------
813	if flag==0:	880	if flag==0:
@@ -825,6 +892,8 class WeatherRHIPlot(Plot):
825	flag=0	892	flag=0
826	data_ele = self.res_ele	893	data_ele = self.res_ele
827	data_weather = self.res_weather	894	data_weather = self.res_weather
		895	'''
		896	print("OUPUT data_ele",data_ele)
828		897
829	return data_weather,data_ele	898	return data_weather,data_ele
830		899
@@ -837,6 +906,13 class WeatherRHIPlot(Plot):
837	r_mask = numpy.where(r>=0)[0]	906	r_mask = numpy.where(r>=0)[0]
838	r = numpy.arange(len(r_mask))*delta_height	907	r = numpy.arange(len(r_mask))*delta_height
839	self.y = 2*r	908	self.y = 2*r
		909	res = 1
		910	print("data['weather'].shape[0]",data['weather'].shape[0])
		911	ang_max = 80
		912	ang_min = 0
		913	var_ang =ang_max - ang_min
		914	step = (int(var_ang)/(res*data['weather'].shape[0]))
		915	print("step",step)
840	'''	916	'''
841	#-------------------------------------------------------------	917	#-------------------------------------------------------------
842	# RADAR	918	# RADAR
@@ -846,43 +922,26 class WeatherRHIPlot(Plot):
846	res = 1	922	res = 1
847	# STEP	923	# STEP
848	step = (180/(res*data['weather'].shape[0]))	924	step = (180/(res*data['weather'].shape[0]))
849
850
851	self.res_weather, self.res_ele = self.const_ploteo(data_weather=data['weather'][:,r_mask],data_ele=data['ele'],step=step,res=res)	925	self.res_weather, self.res_ele = self.const_ploteo(data_weather=data['weather'][:,r_mask],data_ele=data['ele'],step=step,res=res)
852	self.res_azi = numpy.mean(data['azi'])	926	self.res_azi = numpy.mean(data['azi'])
853	print("self.res_ele------------------------------:",self.res_ele)	927	print("self.res_ele------------------------------:",self.res_ele)
854	################# PLOTEO ###################
855
856	for i,ax in enumerate(self.axes):
857	if ax.firsttime:
858	plt.clf()
859	cgax, pm = wrl.vis.plot_rhi(self.res_weather,r=r,th=self.res_ele,fig=self.figures[0], proj='cg', vmin=8, vmax=35)
860	else:
861	plt.clf()
862	cgax, pm = wrl.vis.plot_rhi(self.res_weather,r=r,th=self.res_ele,fig=self.figures[0], proj='cg', vmin=8, vmax=35)
863	caax = cgax.parasites[0]
864	paax = cgax.parasites[1]
865	cbar = plt.gcf().colorbar(pm, pad=0.075)
866	caax.set_xlabel('x_range [km]')
867	caax.set_ylabel('y_range [km]')
868	plt.text(1.0, 1.05, 'Elevacion '+str(thisDatetime)+" Step "+str(self.ini)+ " Azi: "+str(round(self.res_azi,2)), transform=caax.transAxes, va='bottom',ha='right')
869
870	self.ini= self.ini+1
871
872
873	'''	928	'''
874	#--------------------------------------------------------	929	#--------------------------------------------------------
		930	print('weather',data['weather'].shape)
		931	print('ele',data['ele'].shape)
875		932
876	~~###~~self.res_weather, self.res_ele = self.const_ploteo(data_weather=data['weather'][:,r_mask],data_~~azi~~=data['ele'],step=step,res=res)	933	self.res_weather, self.res_ele = self.const_ploteo(data_weather=data['weather'][:,r_mask],data_ele=data['ele'],step=step,res=res,ang_max=ang_max,ang_min=ang_min)
877	~~###~~self.res_azi = numpy.mean(data['azi'])	934	self.res_azi = numpy.mean(data['azi'])
		935	print("self.res_ele",self.res_ele)
878	#-------------	936	#-------------
879	# 90 angulos en el axis 0	937	# 90 angulos en el axis 0
880	# 1000 step en el axis 1	938	# 1000 step en el axis 1
881	self.res_weather = numpy.ones([120,1000])	939	###self.res_weather = numpy.ones([120,1000])
882	r = numpy.linspace(0,1999,1000)	940	###r = numpy.linspace(0,1999,1000)
883	self.res_ele = numpy.arange(0,120)	941	###self.res_ele = numpy.arange(0,120)
884	self.res_azi = 240	942	###self.res_azi = 240
885	#-------------	943	#-------------
		944	'''
886	for i,ax in enumerate(self.axes):	945	for i,ax in enumerate(self.axes):
887	if ax.firsttime:	946	if ax.firsttime:
888	plt.clf()	947	plt.clf()
@@ -896,5 +955,6 class WeatherRHIPlot(Plot):
896	caax.set_xlabel('x_range [km]')	955	caax.set_xlabel('x_range [km]')
897	caax.set_ylabel('y_range [km]')	956	caax.set_ylabel('y_range [km]')
898	plt.text(1.0, 1.05, 'Elevacion '+str(thisDatetime)+" Step "+str(self.ini)+ " Azi: "+str(round(self.res_azi,2)), transform=caax.transAxes, va='bottom',ha='right')	957	plt.text(1.0, 1.05, 'Elevacion '+str(thisDatetime)+" Step "+str(self.ini)+ " Azi: "+str(round(self.res_azi,2)), transform=caax.transAxes, va='bottom',ha='right')
899		958	'''
		959	print("""""""""""""self.ini""""""""""""",self.ini)
900	self.ini= self.ini+1	960	self.ini= self.ini+1

schainpy/model/proc/jroproc_parameters.py +4 -4

                         #print(i)# OJO IDENTIFICADOR DE SINCRONISMO
                         self.utc_ped_list.append(self.gettimeutcfromDirFilename(path=self.path_ped,file=self.list_pedestal[i]))
                     dataOut.wr_exp     = wr_exp
-                    print("SETUP READY")
+                    #print("SETUP READY")
                 def setNextFileP(self,dataOut):
                     if self.__dataReady:
                         dataOut.data_360         = data_360 # S
                         ##print("---------------------------------------------------------------------------------")
-                        print("---------------------------DATAREADY---------------------------------------------")
+                        ###print("---------------------------DATAREADY---------------------------------------------")
                         ##print("---------------------------------------------------------------------------------")
                         ##print("data_360",dataOut.data_360.shape)
                         dataOut.data_azi         = data_p
                         dataOut.data_ele         = data_e
-                        print("azi:    ",dataOut.data_azi)
+                        ###print("azi:    ",dataOut.data_azi)
-                        print("ele:    ",dataOut.data_ele)
+                        ###print("ele:    ",dataOut.data_ele)
                         #print("jroproc_parameters",data_p[0],data_p[-1])#,data_360.shape,avgdatatime)
                         dataOut.utctime         = avgdatatime
                         dataOut.flagNoData      = False

schainpy/scripts/USRP_PLOT_THOR.py +1 -1

             #######################################################################
             procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
-            procUnitConfObjB.addParameter(name='nFFTPoints', value='250', format='int')
+            procUnitConfObjB.addParameter(name='nFFTPoints', value='0', format='int')
             procUnitConfObjB.addParameter(name='nProfiles', value='250', format='int')
             #procUnitConfObjC = controllerObj.addProcUnit(datatype='SpectraHeisProc', inputId=procUnitConfObjA.getId())

schainpy/scripts/check_pedestal.py +3 -2

             #--------------------------------
             path_ped= "/DATA_RM/TEST_PEDESTAL/P20220401-172744"
             # Metodo para verificar numero
             def isNumber(str):
                 try:
             # CONDICION POR LEER EN TIEMPO REAL NO OFFLINE
             for filename in LIST:
-                tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_pos")
+                #tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_pos")
                 tmp_ele_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_pos")
-                #tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_vel")
+                tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_vel")
                 #tmp_ele_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_vel")
                 # CONDICION POR LEER EN TIEMPO REAL NO OFFLINE

schainpy/scripts/test_wradlib_RHI.py +14 0

             warnings.filterwarnings('ignore')
             # lectura de gaMIC hdf5 file
             filename        = wrl.util.get_wradlib_data_file("/home/soporte/Downloads/2014-06-09--185000.rhi.mvol")
+            #filename        = wrl.util.get_wradlib_data_file("2014-06-09--185000.rhi.mvol")
             data1, metadata = wrl.io.read_gamic_hdf5(filename)
             print(data1)
             data1 = data1['SCAN0']['ZH']['data']
             print("Longitud,Latitud,Altura",site)
             ma1 =  np.array(data1)
+            for i in range(3):
+                print("dark",ma1[i])
             '''
             mask_ind = np.where(data1 <= np.nanmin(data1))
             data1[mask_ind] = np.nan
             ma1 = np.ma.array(data1, mask=np.isnan(data1))
             '''
+            ####################### test ####################s
+            th=(np.arange(450)/10.0)+5
+            #th= np.roll(th,-2)
+            #th=np.where(a<7,np.nan,a)
+            ma1=np.roll(ma1,-2,axis=0)
+            for i in range(3):
+                print("green",ma1[i])
+            print("a",th)
+            #th = [i for i in reversed(a)]
+            ######################### test
             #cgax, pm = wrl.vis.plot_rhi(ma1,r=r,th=th,rf=1e3)
             fig = plt.figure(figsize=(10,8))
             cgax, pm = wrl.vis.plot_rhi(ma1,r=r,th=th,rf=1e3,fig=fig, ax=111,proj='cg')

schainpy/scripts/wr_marzo_v1_001.py +26 -10

             #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-171722"
             #path = "/DATA_RM/DRONE01ABRIL1701"
             path = "/DATA_RM/DATA/Torre_con_bola_1649092242/rawdata"
+            path="/DATA_RM/DRONE01ABRIL1727"
-            path_ped = "/DATA_RM/DRONE01ABRIL1450"
+            #path_ped = "/DATA_RM/DRONE01ABRIL1450"
+            path_ped="/DATA_RM/TEST_PEDESTAL/P20220401-172744"
+            #path_ped = "/DATA_RM/DATA/Torre_con_bola_1649092242/position/2022-04-04T17-00-00"
             #-------------------------------------------------------------------------------
             figpath_pp     = "/home/soporte/Pictures/Torre_con_bola_1649092242"
             #figpath_pp     = "/home/soporte/Pictures/MARTES_22_PP_1M_1us"
             figpath_spec   = "/home/soporte/Pictures/MARTES_22_1M_1us"
             figpath_pp_ppi = "/home/soporte/Pictures/MARTES_22_1M_1us_PPI"
+            figpath_pp_rhi  = "/DATA_RM/LUNES04ABRIL_1200_RHI"
             #--------------------------OPCIONES---------------------------------------------
-            plot        = 1
             plot_ppi    = 0
-            integration = 0
+            plot        = 0
+            plot_rhi    = 1
+            integration = 1
             save        = 0
             plot_spec   = 0
             #---------------------------SAVE HDF5 PROCESADO/--------------------------------
             #------------------------ UNIDAD DE LECTURA-------------------------------------
             readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
                                                         path=path,
-                                                        startDate="2022/04/04",#today,
+                                                        startDate="2022/04/01",#today,
-                                                        endDate="2022/04/04",#today,
+                                                        endDate="2022/04/01",#today,
-                                                        startTime='12:11:05',#'17:39:25',
+                                                        startTime='00:10:05',#'17:39:25',
                                                         endTime='23:59:59',#23:59:59',
                                                         delay=0,
                                                         #set=0,
              # CUARTA PARTE de 60 Km POR ESO ENTRE 20  - 3 Km
+            '''
             procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId())
             procUnitConfObjB.addParameter(name='nFFTPoints', value='32', format='int')
             procUnitConfObjB.addParameter(name='nProfiles', value='32', format='int')
             opObj11.addParameter(name='showprofile', value='1', format='int')
             #opObj11.addParameter(name='save', value=figpath, format='str')
             opObj11.addParameter(name='save_period', value=10, format='int')
+            '''
             if mode_proc ==0:
                 ####################### METODO PULSE PAIR ######################################################################
                     opObj11 = procUnitConfObjB.addOperation(name='PedestalInformation')
                     opObj11.addParameter(name='path_ped', value=path_ped)
                     opObj11.addParameter(name='t_Interval_p', value='0.01', format='float')
-                    opObj11.addParameter(name='wr_exp', value='PPI')
+                    #opObj11.addParameter(name='wr_exp', value='PPI')
+                    opObj11.addParameter(name='wr_exp', value='RHI')
                 if plot_ppi==1:
                     opObj11 = procUnitConfObjB.addOperation(name='Block360')
                     opObj11.addParameter(name='n', value='10', format='int')
                     opObj11= procUnitConfObjB.addOperation(name='WeatherPlot',optype='other')
                     opObj11.addParameter(name='save', value=figpath_pp_ppi)
                     opObj11.addParameter(name='save_period', value=1)
+                if plot_rhi==1:
+                    opObj11 = procUnitConfObjB.addOperation(name='Block360')
+                    opObj11.addParameter(name='n', value='10', format='int')
+                    opObj11.addParameter(name='mode', value=mode_proc, 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='WeatherRHIPlot',optype='other')
+                    opObj11.addParameter(name='save', value=figpath_pp_rhi)
+                    opObj11.addParameter(name='save_period', value=1)
             controllerObj.start()

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