##// END OF EJS Templates
DBS and SA new scripts
Julio Valdez -
r855:48fed4a7e9fc
parent child
Show More
@@ -30,8 +30,8 controllerObj.setup(id = '191', name='test01', description=desc)
30 #2014051 20 Feb 2014
30 #2014051 20 Feb 2014
31 path = '/media/joscanoa/84A65E64A65E5730/soporte/Data/MST/DBS/d2014051'
31 path = '/media/joscanoa/84A65E64A65E5730/soporte/Data/MST/DBS/d2014051'
32 # path = '/media/joscanoa/disco4/Data/2014/DBS_SA JAN 2014/DBS_SA/250/d2014050'
32 # path = '/media/joscanoa/disco4/Data/2014/DBS_SA JAN 2014/DBS_SA/250/d2014050'
33 pathFigure = '/home/joscanoa/Pictures/testHDF5/figure'
33 pathfile1 = os.path.join(os.environ['HOME'],'Pictures/testHDF5/moments')
34 xmax = '4'
34 xmax = '1'
35 xmin = '0'
35 xmin = '0'
36 startTime = '00:00:00'
36 startTime = '00:00:00'
37 filehdf5 = "DBS_2014051.hdf5"
37 filehdf5 = "DBS_2014051.hdf5"
@@ -52,7 +52,7 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
52 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
52 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
53
53
54
54
55 #--------------------------------------------------------------------------------------------------
55 #------------------------------ Voltage Processing Unit -------------------------------------
56
56
57 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
57 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
58
58
@@ -66,29 +66,18 opObj11 = procUnitConfObj0.addOperation(name='selectHeightsByIndex')
66 opObj11.addParameter(name='minIndex', value='10', format='float')
66 opObj11.addParameter(name='minIndex', value='10', format='float')
67 opObj11.addParameter(name='maxIndex', value='60', format='float')
67 opObj11.addParameter(name='maxIndex', value='60', format='float')
68
68
69 #---------------------------------------------------------------------------------------------------
69 # opObj12 = procUnitConfObj0.addOperation(name='selectChannels')
70 # opObj12.addParameter(name='channelList', value='0,1', format='intlist')
71
72 #------------------------------ Spectra Processing Unit -------------------------------------
70
73
71 procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId())
74 procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId())
72 procUnitConfObj1.addParameter(name='nFFTPoints', value='64', format='int')
75 procUnitConfObj1.addParameter(name='nFFTPoints', value='64', format='int')
73 procUnitConfObj1.addParameter(name='nProfiles', value='64', format='int')
74 # procUnitConfObj1.addParameter(name='ippFactor', value='2', format='int')
76 # procUnitConfObj1.addParameter(name='ippFactor', value='2', format='int')
75 procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(0,1),(2,1)', format='pairsList')
76
77
77 opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
78 opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
78 opObj11.addParameter(name='n', value='5', format='int')
79 opObj11.addParameter(name='n', value='5', format='int')
79
80
80 # opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
81 # opObj14.addParameter(name='id', value='1', format='int')
82 # opObj14.addParameter(name='wintitle', value='Con interf', format='str')
83 # opObj14.addParameter(name='save', value='0', format='bool')
84 # opObj14.addParameter(name='figpath', value=pathFigure, format='str')
85 # opObj14.addParameter(name='zmin', value='5', format='int')
86 # opObj14.addParameter(name='zmax', value='90', format='int')
87
88 opObj12 = procUnitConfObj1.addOperation(name='removeInterference')
89 opObj13 = procUnitConfObj1.addOperation(name='removeDC')
90 opObj13.addParameter(name='mode', value='1', format='int')
91 #
92 opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
81 opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
93 opObj14.addParameter(name='id', value='1', format='int')
82 opObj14.addParameter(name='id', value='1', format='int')
94 opObj14.addParameter(name='wintitle', value='Con interf', format='str')
83 opObj14.addParameter(name='wintitle', value='Con interf', format='str')
@@ -96,8 +85,11 opObj14.addParameter(name='save', value='0', format='bool')
96 opObj14.addParameter(name='figpath', value=pathFigure, format='str')
85 opObj14.addParameter(name='figpath', value=pathFigure, format='str')
97 opObj14.addParameter(name='zmin', value='5', format='int')
86 opObj14.addParameter(name='zmin', value='5', format='int')
98 opObj14.addParameter(name='zmax', value='90', format='int')
87 opObj14.addParameter(name='zmax', value='90', format='int')
99
88
100
89 opObj12 = procUnitConfObj1.addOperation(name='removeInterference')
90 opObj13 = procUnitConfObj1.addOperation(name='removeDC')
91 opObj13.addParameter(name='mode', value='1', format='int')
92
101 opObj12 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other')
93 opObj12 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other')
102 opObj12.addParameter(name='id', value='2', format='int')
94 opObj12.addParameter(name='id', value='2', format='int')
103 opObj12.addParameter(name='wintitle', value='RTI Plot', format='str')
95 opObj12.addParameter(name='wintitle', value='RTI Plot', format='str')
@@ -107,20 +99,29 opObj12.addParameter(name='xmin', value=xmin, format='float')
107 opObj12.addParameter(name='xmax', value=xmax, format='float')
99 opObj12.addParameter(name='xmax', value=xmax, format='float')
108 opObj12.addParameter(name='zmin', value='5', format='int')
100 opObj12.addParameter(name='zmin', value='5', format='int')
109 opObj12.addParameter(name='zmax', value='90', format='int')
101 opObj12.addParameter(name='zmax', value='90', format='int')
102
110
103
111 #--------------------------------------------------------------------------------------------------
104 #------------------------------ Parameters Processing Unit -------------------------------------
112
105
113 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
106 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
114 opObj20 = procUnitConfObj2.addOperation(name='GetMoments')
107
115
108 opObj11 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
109
110 # opObj12 = procUnitConfObj2.addOperation(name='HDF5Writer', optype='other')
111 # opObj12.addParameter(name='path', value=pathfile1)
112 # opObj12.addParameter(name='blocksPerFile', value='10', format='int')
113 # opObj12.addParameter(name='metadataList',value='type,inputUnit,heightList,paramInterval,timeZone',format='list')
114 # opObj12.addParameter(name='dataList',value='data_param,data_SNR,noise,utctime',format='list')
115 # opObj12.addParameter(name='mode',value='1',format='int')
116
116 # opObj21 = procUnitConfObj2.addOperation(name='MomentsPlot', optype='other')
117 # opObj21 = procUnitConfObj2.addOperation(name='MomentsPlot', optype='other')
117 # opObj21.addParameter(name='id', value='3', format='int')
118 # opObj21.addParameter(name='id', value='3', format='int')
118 # opObj21.addParameter(name='wintitle', value='Moments Plot', format='str')
119 # opObj21.addParameter(name='wintitle', value='Moments Plot', format='str')
119 # opObj21.addParameter(name='save', value='1', format='bool')
120 # opObj21.addParameter(name='save', value='0', format='bool')
120 # opObj21.addParameter(name='figpath', value=pathFigure, format='str')
121 # # opObj21.addParameter(name='figpath', value=pathFigure, format='str')
121 # opObj21.addParameter(name='zmin', value='5', format='int')
122 # opObj21.addParameter(name='zmin', value='5', format='int')
122 # opObj21.addParameter(name='zmax', value='90', format='int')
123 # opObj21.addParameter(name='zmax', value='90', format='int')
123
124 #
124 # opObj21 = procUnitConfObj2.addOperation(name='ParametersPlot', optype='other')
125 # opObj21 = procUnitConfObj2.addOperation(name='ParametersPlot', optype='other')
125 # opObj21.addParameter(name='id', value='5', format='int')
126 # opObj21.addParameter(name='id', value='5', format='int')
126 # opObj21.addParameter(name='wintitle', value='Radial Velocity Plot', format='str')
127 # opObj21.addParameter(name='wintitle', value='Radial Velocity Plot', format='str')
@@ -133,21 +134,19 opObj20 = procUnitConfObj2.addOperation(name='GetMoments')
133 # opObj21.addParameter(name='SNRthresh', value='0', format='float')
134 # opObj21.addParameter(name='SNRthresh', value='0', format='float')
134 # opObj21.addParameter(name='xmin', value=xmin, format='float')
135 # opObj21.addParameter(name='xmin', value=xmin, format='float')
135 # opObj21.addParameter(name='xmax', value=xmax, format='float')
136 # opObj21.addParameter(name='xmax', value=xmax, format='float')
136
137
137 opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
138 opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
138 opObj22.addParameter(name='technique', value='DBS', format='str')
139 opObj22.addParameter(name='technique', value='DBS', format='str')
139 opObj22.addParameter(name='correctAzimuth', value='51.06', format='float')
140 opObj22.addParameter(name='correctAzimuth', value='51.06', format='float')
140 opObj22.addParameter(name='correctFactor', value='-1', format='float')
141 opObj22.addParameter(name='correctFactor', value='-1', format='float')
141 opObj22.addParameter(name='dirCosx', value='0.041016, 0, -0.054688', format='floatlist')
142 opObj22.addParameter(name='dirCosx', value='0.041016, 0, -0.054688', format='floatlist')
142 opObj22.addParameter(name='dirCosy', value='-0.041016, 0.025391, -0.023438', format='floatlist')
143 opObj22.addParameter(name='dirCosy', value='-0.041016, 0.025391, -0.023438', format='floatlist')
143 # opObj22.addParameter(name='horizontalOnly', value='0', format='bool')
144
144 # opObj22.addParameter(name='channelList', value='1,2,3', format='intlist')
145
146 opObj23 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
145 opObj23 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
147 opObj23.addParameter(name='id', value='3', format='int')
146 opObj23.addParameter(name='id', value='4', format='int')
148 opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str')
147 opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str')
149 opObj23.addParameter(name='save', value='1', format='bool')
148 opObj23.addParameter(name='save', value='0', format='bool')
150 opObj23.addParameter(name='figpath', value = pathFigure, format='str')
149 # opObj23.addParameter(name='figpath', value = pathFigure, format='str')
151 opObj23.addParameter(name='zmin', value='-10', format='int')
150 opObj23.addParameter(name='zmin', value='-10', format='int')
152 opObj23.addParameter(name='zmax', value='10', format='int')
151 opObj23.addParameter(name='zmax', value='10', format='int')
153 opObj23.addParameter(name='zmin_ver', value='-80', format='float')
152 opObj23.addParameter(name='zmin_ver', value='-80', format='float')
@@ -44,7 +44,9 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
44 endTime='23:59:59',
44 endTime='23:59:59',
45 online=0,
45 online=0,
46 delay=5,
46 delay=5,
47 walk=0)
47 walk=0,
48 getblock=1,
49 blocksize=32768)
48
50
49 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
51 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
50
52
@@ -56,78 +58,64 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=rea
56 opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
58 opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
57
59
58 opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other')
60 opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other')
59 opObj11.addParameter(name='n', value='600', format='int')
61 # opObj11.addParameter(name='n', value='600', format='int')
60 # opObj11.addParameter(name='n', value='10', format='int')
62 opObj11.addParameter(name='n', value='256', format='int')
61
63
62 opObj11 = procUnitConfObj0.addOperation(name='selectHeightsByIndex')
64 opObj11 = procUnitConfObj0.addOperation(name='selectHeightsByIndex')
63 opObj11.addParameter(name='minIndex', value='10', format='float')
65 opObj11.addParameter(name='minIndex', value='10', format='float')
64 opObj11.addParameter(name='maxIndex', value='60', format='float')
66 opObj11.addParameter(name='maxIndex', value='60', format='float')
65 #---------------------------------------------------------------------------------------------------
67 #---------------------------------------------------------------------------------------------------
66 procUnitConfObj1 = controllerObj.addProcUnit(datatype='CorrelationProc', inputId=procUnitConfObj0.getId())
68 procUnitConfObj1 = controllerObj.addProcUnit(datatype='CorrelationProc', inputId=procUnitConfObj0.getId())
67 # procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(1,0),(2,3)', format='pairsList')
69 procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(1,0),(2,3)', format='pairsList')
68 procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(0,3),(0,2),(1,3),(1,2),(0,1),(2,3)', format='pairsList')
70 # procUnitConfObj1.addParameter(name='removeDC', value='1', format='bool')
69 procUnitConfObj1.addParameter(name='fullT', value='1', format='bool')
71 # #procUnitConfObj1.addParameter(name='lagT', value='0,1,2,3', format='intlist')
70 procUnitConfObj1.addParameter(name='removeDC', value='1', format='bool')
72 #
71 #procUnitConfObj1.addParameter(name='lagT', value='0,1,2,3', format='intlist')
73 # opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other')
72
74 # opObj12.addParameter(name='id', value='1', format='int')
73 opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other')
75 # opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str')
74 opObj12.addParameter(name='id', value='1', format='int')
76 # opObj12.addParameter(name='save', value='1', format='bool')
75 opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str')
77 # opObj12.addParameter(name='zmin', value='0', format='int')
76 opObj12.addParameter(name='save', value='1', format='bool')
78 # opObj12.addParameter(name='zmax', value='1', format='int')
77 opObj12.addParameter(name='zmin', value='0', format='int')
79 # opObj12.addParameter(name='figpath', value = pathFigure, format='str')
78 opObj12.addParameter(name='zmax', value='1', format='int')
80 #
79 opObj12.addParameter(name='figpath', value = pathFigure, format='str')
81 # opObj12 = procUnitConfObj1.addOperation(name='removeNoise')
80
82 # opObj12.addParameter(name='mode', value='2', format='int')
81 opObj12 = procUnitConfObj1.addOperation(name='removeNoise')
83 # opObj12 = procUnitConfObj1.addOperation(name='calculateNormFactor')
82 opObj12.addParameter(name='mode', value='2', format='int')
84 #
83 opObj12 = procUnitConfObj1.addOperation(name='calculateNormFactor')
85 # opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other')
84
86 # opObj12.addParameter(name='id', value='2', format='int')
85 opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other')
87 # opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str')
86 opObj12.addParameter(name='id', value='2', format='int')
88 # opObj12.addParameter(name='save', value='1', format='bool')
87 opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str')
89 # opObj12.addParameter(name='zmin', value='0', format='int')
88 opObj12.addParameter(name='save', value='1', format='bool')
90 # opObj12.addParameter(name='zmax', value='1', format='int')
89 opObj12.addParameter(name='zmin', value='0', format='int')
91 # opObj12.addParameter(name='figpath', value = pathFigure, format='str')
90 opObj12.addParameter(name='zmax', value='1', format='int')
92 #
91 opObj12.addParameter(name='figpath', value = pathFigure, format='str')
93 # #---------------------------------------------------------------------------------------------------
92
93 #---------------------------------------------------------------------------------------------------
94 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
94 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
95 opObj20 = procUnitConfObj2.addOperation(name='GetLags')
96
95
96 opObj20 = procUnitConfObj2.addOperation(name='SALags', optype='other')
97 #
97 opObj21 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
98 opObj21 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
98 opObj21.addParameter(name='technique', value='SA', format='str')
99 opObj21.addParameter(name='technique', value='SA', format='str')
99 # opObj21.addParameter(name='correctFactor', value='-1', format='float')
100 # # opObj21.addParameter(name='correctFactor', value='-1', format='float')
100 opObj21.addParameter(name='positionX', value='36,0,36,0', format='floatlist')
101 opObj21.addParameter(name='positionX', value='36,0,36,0', format='floatlist')
101 opObj21.addParameter(name='positionY', value='36,0,0,36', format='floatlist')
102 opObj21.addParameter(name='positionY', value='36,0,0,36', format='floatlist')
102 opObj21.addParameter(name='azimuth', value='51.06', format='float')
103 opObj21.addParameter(name='azimuth', value='51.06', format='float')
103 opObj21.addParameter(name='crosspairsList', value='(0,3),(0,2),(1,3),(1,2),(0,1),(2,3)', format='pairsList')#COrregir
104
104 #
105 # opObj22 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
105 opObj22 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
106 # opObj22.addParameter(name='id', value='4', format='int')
106 opObj22.addParameter(name='id', value='4', format='int')
107 # opObj22.addParameter(name='wintitle', value='Wind Profiler', format='str')
107 opObj22.addParameter(name='wintitle', value='Wind Profiler', format='str')
108 opObj22.addParameter(name='save', value='1', format='bool')
109 opObj22.addParameter(name='figpath', value = pathFigure, format='str')
110 opObj22.addParameter(name='zmin', value='-15', format='int')
111 opObj22.addParameter(name='zmax', value='15', format='int')
112 opObj22.addParameter(name='zmin_ver', value='-80', format='float')
113 opObj22.addParameter(name='zmax_ver', value='80', format='float')
114 opObj22.addParameter(name='SNRmin', value='-20', format='int')
115 opObj22.addParameter(name='SNRmax', value='40', format='int')
116 opObj22.addParameter(name='SNRthresh', value='-3.5', format='float')
117 opObj22.addParameter(name='xmin', value=xmin, format='float')
118 opObj22.addParameter(name='xmax', value=xmax, format='float')
119 # #-----------------------------------------------------------------------------------
120 #
121 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId())
122 # procUnitConfObj2.addParameter(name='nFFTPoints', value='128', format='int')
123 # procUnitConfObj2.addParameter(name='nProfiles', value='128', format='int')
124 # procUnitConfObj2.addParameter(name='pairsList', value='(0,0),(0,1),(2,1)', format='pairsList')
125 #
126 # opObj22 = procUnitConfObj2.addOperation(name='SpectraPlot', optype='other')
127 # opObj22.addParameter(name='id', value='5', format='int')
128 # opObj22.addParameter(name='wintitle', value='Spectra Plot', format='str')
129 # opObj22.addParameter(name='save', value='1', format='bool')
108 # opObj22.addParameter(name='save', value='1', format='bool')
130 # opObj22.addParameter(name='figpath', value = pathFigure, format='str')
109 # opObj22.addParameter(name='figpath', value = pathFigure, format='str')
110 # opObj22.addParameter(name='zmin', value='-15', format='int')
111 # opObj22.addParameter(name='zmax', value='15', format='int')
112 # opObj22.addParameter(name='zmin_ver', value='-80', format='float')
113 # opObj22.addParameter(name='zmax_ver', value='80', format='float')
114 # opObj22.addParameter(name='SNRmin', value='-20', format='int')
115 # opObj22.addParameter(name='SNRmax', value='40', format='int')
116 # opObj22.addParameter(name='SNRthresh', value='-3.5', format='float')
117 # opObj22.addParameter(name='xmin', value=xmin, format='float')
118 # opObj22.addParameter(name='xmax', value=xmax, format='float')
131
119
132 #-----------------------------------------------------------------------------------
120 #-----------------------------------------------------------------------------------
133
121
General Comments 0
You need to be logged in to leave comments. Login now