##// END OF EJS Templates
schain
RSS
jrproc_parametes.py update y script EWDriftsMP_01feb2023.py
sebastianVP -
r1743:6b1a256569db
parent child
Show More
Side by Side Unified
Context file:
r1743:6b1a256569db Loading diff...:

The requested changes are too big and content was truncated. Show full diff

Show file before | Show file after | Hide comments
@@ -0,0 +1,262
1
2 import os, sys
3 import json
4
5 #from controller import *
6 from schainpy.controller import Project
7
8 desc = "EW DRIFTS MP Experiment"
9 filename = "EWDrifts.xml"
10
11 controllerObj = Project()
12
13 controllerObj.setup(id = '191', name='test01', description=desc)
14
15 #Experimentos
16
17 #path = "/data/dia"
18 #path = '/home/pcondor/data'
19 #path = '/media/pcondor/DATA1/Database/ewdriftsschain2023prue/data'
20 #path = '/data/2024_01/MP_ISR/main_radar/rawdata/d2024023'
21 path = '/data/ISR_JULIA/d2024092'
22 #pathFigure = '/media/pcondor/DATA1/Database/ewdriftsschain2023wh5'
23 pathFile = '/media/pcondor/DATA1/Database/ewdriftsabr2024sch/EW_Drifts_01abr'
24 pathFigure = pathFile
25 pathFileavg = pathFile+'/avg'
26 pathFiledata = pathFile+'/Drifts-data'
27 #pathFileavg = '/media/pcondor/DATA1/Database/ewdriftsschain2023wh5/avg'
28 #pathFiledata = '/media/pcondor/DATA1/Database/ewdriftsschain2023wh5/Drifts-data'
29
30 xmin = 0
31 xmax = 24
32 #------------------------------------------------------------------------------------------------
33 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
34 path=path,
35 startDate='2024/04/01',
36 endDate='2024/04/01',
37 startTime='00:00:00',
38 endTime='23:59:59',
39 online=0,
40 getByBlock=1,
41 walk=0)
42
43 #--------------------------------------------------------------------------------------------------
44
45 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
46
47 #opObj11 = procUnitConfObj0.addOperation(name='selectHeights')
48 # # opObj11.addParameter(name='minHei', value='320.0', format='float')
49 # # opObj11.addParameter(name='maxHei', value='350.0', format='float')
50 #opObj11.addParameter(name='minHei', value='0.01', format='float')
51 #opObj11.addParameter(name='maxHei', value='960.0', format='float')
52
53 opObj11 = procUnitConfObj0.addOperation(name='selectChannels')
54 opObj11.addParameter(name='channelList', value='0,0,1,1', format='intlist')
55
56 #opObj11 = procUnitConfObj0.addOperation(name='Reshaper')
57 #opObj11.addParameter(name='shape', value='(500,980)', format='intlist')
58
59 opObj11 = procUnitConfObj0.addOperation(name='ProfileSelector', optype='other')
60 opObj11.addParameter(name='profileRangeList', value='0,127', format='intlist')
61
62 opObj11 = procUnitConfObj0.addOperation(name='filterByHeights')
63 opObj11.addParameter(name='window', value='10', format='int')
64
65 code=[[-1,-1,1],[1,1,-1]]
66 #code = [[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,-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,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],[-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,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,-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],[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,-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,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],[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,-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,-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],[-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,-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,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],[-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,-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,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],[-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,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,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],[-1,-1,1],[-1,-1,1],[1,1,-1],[1,1,-1]]
67 opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
68 opObj11.addParameter(name='code', value=code, format='floatlist')
69 opObj11.addParameter(name='nCode', value='2', format='int')
70 opObj11.addParameter(name='nBaud', value='3', format='int')
71
72 opObj11 = procUnitConfObj0.addOperation(name='selectHeights')
73 opObj11.addParameter(name='minHei', value='0.0', format='float')
74 opObj11.addParameter(name='maxHei', value='960', format='float')
75
76 procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId())
77 procUnitConfObj1.addParameter(name='nFFTPoints', value='128', format='int')
78 procUnitConfObj1.addParameter(name='nProfiles', value='128', format='int')
79 #procUnitConfObj1.addParameter(name='pairsList', value='(2,3),(4,5)', format='pairsList')#,(2,3)
80 procUnitConfObj1.addParameter(name='pairsList', value='(0,1),(2,3)', format='pairsList')
81
82 #opObj11 = procUnitConfObj1.addOperation(name='selectHeights')
83 # # opObj11.addParameter(name='minHei', value='320.0', format='float')
84 # # opObj11.addParameter(name='maxHei', value='350.0', format='float')
85 #opObj11.addParameter(name='minHei', value='0.0', format='float')
86 #opObj11.addParameter(name='maxHei', value='960.0', format='float')
87
88 #opObj11 = procUnitConfObj1.addOperation(name='selectChannels')
89 #opObj11.addParameter(name='channelList', value='2,3,4,5', format='intlist')
90
91 opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
92 opObj11.addParameter(name='n', value='1', format='float')
93 #opObj11.addParameter(name='timeInterval', value='300.0', format='float')
94
95 #opObj13 = procUnitConfObj1.addOperation(name='removeDC')
96
97 #opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
98 #opObj14.addParameter(name='id', value='65', format='int')
99 ## # opObj14.addParameter(name='wintitle', value='Con interf', format='str')
100 #opObj14.addParameter(name='save', value=pathFigure, format='str')
101 ##opObj14.addParameter(name='save_period', value=1, format='int')
102 #opObj14.addParameter(name='zmin', value='10', format='int')
103 #opObj14.addParameter(name='zmax', value='26', format='int')
104 #
105
106 #opObj12 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other')
107 #opObj12.addParameter(name='id', value='63', format='int')
108 #opObj12.addParameter(name='wintitle', value='RTI Plot', format='str')
109 #opObj12.addParameter(name='save', value=pathFigure, format='str')
110 #opObj12.addParameter(name='save_period', value=10, format='int')
111 ##opObj12.addParameter(name='figpath', value = pathFigure, format='str')
112 #opObj12.addParameter(name='xmin', value=xmin, format='float')
113 #opObj12.addParameter(name='xmax', value=xmax, format='float')
114 #opObj12.addParameter(name='zmin', value='20', format='int')
115 #opObj12.addParameter(name='zmax', value='36', format='int')
116
117 #--------------------------------------------------------------------------------------------------
118
119 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
120 opObj20 = procUnitConfObj2.addOperation(name='SpectralFitting', optype='other')
121 opObj20.addParameter(name='path', value='/home/pcondor/DIR_MADRIGAL/schain/schainpy/model/proc', format='str')
122 opObj20.addParameter(name='file', value='modelSpectralFitting', format='str')
123 opObj20.addParameter(name='groupList', value='(0,1),(2,3)',format='multiList')
124 opObj20.addParameter(name='taver', value='5')
125 opObj20.addParameter(name='coh_th', value='[1]',format='multiList')
126 opObj20.addParameter(name='hei_th', value='[2000]',format='multiList')
127 #opObj20.addParameter(name='filec', value='weightfit', format='str')
128
129 opObj22 = procUnitConfObj2.addOperation(name='HDFWriter', optype='other')
130 opObj22.addParameter(name='path', value=pathFiledata)
131 opObj22.addParameter(name='blocksPerFile', value='1')
132 opObj22.addParameter(name='metadataList',value='heightList,timeZone')
133 opObj22.addParameter(name='dataList',value='tmp_spectra_i,tmp_cspectra_i,tmp_spectra_c,tmp_cspectra_c,clean_num_aver,coh_num_aver,sat_spectra,sat_cspectra,index,utctime,nIncohInt,nCohInt,nProfiles,nFFTPoints,ippFactor,ippSeconds,paramInterval')
134 ##opObj22.addParameter(name='dataList',value='tmp_spectra_i,tmp_cspectra_i,tmp_spectra_c,tmp_cspectra_c,clean_num_aver,coh_num_aver,index,utctime,nIncohInt,nCohInt,nProfiles,nFFTPoints,normFactor,channelList,ippFactor,ippSeconds')
135
136 #angles :-2.41116 3.01082
137 opObj21 = procUnitConfObj2.addOperation(name='EWDriftsEstimation', optype='other')
138 opObj21.addParameter(name='zenith', value='-2.41116, 3.01082', format='floatlist')
139 opObj21.addParameter(name='zenithCorrection', value='0.0', format='float')
140 opObj21.addParameter(name='fileDrifts', value=pathFile)
141
142 # Drifts en h5
143 one = {'gdlatr': 'lat', 'gdlonr': 'lon', 'spcst':'spcst','pl':'pl','cbadn':'cbadn','inttms': 'inttms','azdir7':'azw','eldir7':'elw','azdir8':'aze','eldir8':'ele','jro14':'jro14','jro15':'jro15','jro16':'jro16','nwlos':'nwlos'}
144 two = {
145 'range': ('params', 0),
146 'gdalt': ('params', 1),
147 'VIPN': ('params', 2),
148 'dvipn': ('params', 3),
149 'vipe': ('params', 4),
150 'dvipe': ('params', 5),
151 'vi7': ('params', 6),
152 'dvi7': ('params', 7),
153 'vi8': ('params', 8),
154 'dvi8': ('params', 9),
155 'PAIWL': ('params', 10),
156 'pacwl': ('params', 11),
157 'pbiwl': ('params', 12),
158 'pbcwl': ('params', 13),
159 'pciel': ('params', 14),
160 'pccel': ('params', 15),
161 'pdiel': ('params', 16),
162 'pdcel': ('params', 17),
163 'jro10': ('params', 18),
164 'jro11': ('params', 19)
165 } #writer
166 ind = ['gdalt']
167
168 meta = {
169 'kinst': 10, #instrument code
170 'kindat': 1910, #type of data
171 'catalog': {
172 'principleInvestigator': 'Danny Scipión',
173 'expPurpose': 'Drifts'#,
174 #'sciRemarks': file_contents
175 },
176 'header': {
177 'analyst': 'Danny Scipión'
178 }
179 }
180
181 op_writer = procUnitConfObj2.addOperation(name='MADWriter')
182 op_writer.addParameter(name='path', value=pathFile)
183 op_writer.addParameter(name='format', value='hdf5')
184 op_writer.addParameter(name='oneDDict', value=json.dumps(one))
185 op_writer.addParameter(name='twoDDict', value=json.dumps(two))
186 op_writer.addParameter(name='ind2DList', value=json.dumps(ind))
187 op_writer.addParameter(name='metadata', value=json.dumps(meta))
188
189 op_writer = procUnitConfObj2.addOperation(name='setHeightDriftsavg')
190
191 # Avg Drifts
192 one_avg = {'gdlatr': 'lat', 'gdlonr': 'lon', 'spcst':'spcst','pl':'pl','cbadn':'cbadn','inttms': 'inttms'}
193 two_avg = {
194 'range': ('params_avg', 4),
195 'gdalt': ('params_avg', 5),
196 'altav': ('params_avg', 6),
197 'VIPN': ('params_avg', 0),
198 'dvipn': ('params_avg', 1),
199 'vipe': ('params_avg', 2),
200 'dvipe': ('params_avg', 3)
201 }
202 ind_avg = ['gdalt']
203 meta = {
204 'kinst': 10, #instrument code
205 'kindat': 1911, #type of data
206 'catalog': {
207 'principleInvestigator': 'Danny Scipión',
208 'expPurpose': 'Drifts'#,
209 #'sciRemarks': file_contents
210 },
211 'header': {
212 'analyst': 'Danny Scipión'
213 }
214 }
215 #dataOut.heightList = dataOut.params_avg[4]
216 op_writer = procUnitConfObj2.addOperation(name='MADWriter')
217 op_writer.addParameter(name='path', value=pathFileavg)
218 op_writer.addParameter(name='format', value='hdf5')
219 op_writer.addParameter(name='oneDDict', value=json.dumps(one_avg))
220 op_writer.addParameter(name='twoDDict', value=json.dumps(two_avg))
221 op_writer.addParameter(name='ind2DList', value=json.dumps(ind_avg))
222 op_writer.addParameter(name='metadata', value=json.dumps(meta))
223
224 op_writer = procUnitConfObj2.addOperation(name='setHeightDrifts')
225
226 opObj24 = procUnitConfObj2.addOperation(name='SpectralMomentsPlot', optype='other')
227 opObj24.addParameter(name='id', value='1', format='int')
228 ### # opObj14.addParameter(name='wintitle', value='Spectral Averaged', format='str')
229 opObj24.addParameter(name='save', value=pathFigure, format='str')
230 ###opObj24.addParameter(name='save_period', value=1, format='int')
231 opObj24.addParameter(name='zmin', value='-8', format='int')
232 opObj24.addParameter(name='zmax', value='16', format='int')
233 opObj24.addParameter(name='xaxis', value='Velocity', format='str')
234
235 #
236 titles=('SNR,Vertical Drifts,Zonal Drifts')
237 #titles=('Zonal Drifts,Vertical Drifts')
238 opObj23 = procUnitConfObj2.addOperation(name='GenericRTIPlot')
239 opObj23.addParameter(name='colormaps', value='jet,RdBu_r,RdBu_r')
240 opObj23.addParameter(name='attr_data', value='data_snr1,data_output')
241 #opObj23.addParameter(name='colormaps', value='RdBu,RdBu')
242 #opObj23.addParameter(name='attr_data', value='data_output')
243 opObj23.addParameter(name='wintitle', value='EW Drifts')
244 opObj23.addParameter(name='save', value=pathFigure)
245 opObj23.addParameter(name='titles', value=titles)
246 opObj23.addParameter(name='zfactors', value='1,1,1')
247 opObj23.addParameter(name='zlimits', value='(-5,20),(-50,50),(-150,150)')
248 opObj23.addParameter(name='cb_labels', value='dB,m/s,m/s')
249 #opObj23.addParameter(name='titles', value=titles)
250 #opObj23.addParameter(name='zfactors', value='1,1')
251 #opObj23.addParameter(name='zlimits', value='(-150,150),(-40,40)')
252 #opObj23.addParameter(name='cb_labels', value='m/s,m/s')
253 opObj23.addParameter(name='throttle', value='1')
254 opObj23.addParameter(name='xmin', value=xmin)
255 opObj23.addParameter(name='xmax', value=xmax)
256 #opObj23.addParameter(name='exp_code', value='110', format='int')
257 #opObj23.addParameter(name='server', value='10.10.110.243:4444', format='int')
258 #opObj23.addParameter(name='tag', value= 'jicamarca', format='str')
259
260 #--------------------------------------------------------------------------------------------------
261
262 controllerObj.start()
Show file before | Show file after | Hide comments
@@ -0,0 +1,186
1
2 import os, sys
3 import json
4 #from controller import *
5 from schainpy.controller import Project
6
7 desc = "EW DRIFTS MP Experiment"
8 filename = "EWDrifts.xml"
9
10 controllerObj = Project()
11
12 controllerObj.setup(id = '191', name='test01', description=desc)
13
14 #Experimentos
15
16 #path = '/media/pcondor/DATA1/Database/ewdriftsene2024sch/EW_Drifts_01ene/Drifts-data'
17 path = '/media/soporte/DATA/PERCY_SCHAIN_UPDATE/driftsschain'
18 #pathFigure = '/media/pcondor/DATA1/Database/ewdriftsschain2023proc'
19 pathFile ='/media/soporte/DATA/PERCY_SCHAIN_UPDATE/driftsschain/tmp'
20 #pathFile = '/media/pcondor/DATA1/Database/ewdriftsene2024sch/EW_Drifts_01enetmp'
21 pathFigure = pathFile
22 pathFileavg = pathFile+'/avg'
23 pathFiledata = pathFile+'/Drifts-data'
24
25 xmin = 0
26 xmax = 24
27 #------------------------------------------------------------------------------------------------
28 readUnitConfObj = controllerObj.addReadUnit(datatype='HDFReader',
29 path=path,
30 startDate='2024/01/23',
31 endDate='2024/01/23',
32 startTime='00:00:00',
33 endTime='23:59:59',
34 #online=0,
35 #getByBlock=1,
36 walk=1,
37 utcoffset='-18000')
38
39 #--------------------------------------------------------------------------------------------------
40
41 #--------------------------------------------------------------------------------------------------
42
43 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId())
44
45 opObj20 = procUnitConfObj2.addOperation(name='SpectralFitting', optype='other')
46 opObj20.addParameter(name='path', value='/home/pcondor/DIR_MADRIGAL/schain/schainpy/model/proc', format='str')
47 opObj20.addParameter(name='file', value='modelSpectralFitting', format='str')
48 opObj20.addParameter(name='groupList', value='(0,1),(2,3)',format='multiList')
49 opObj20.addParameter(name='taver', value='5')
50 opObj20.addParameter(name='coh_th', value='[1]',format='multiList')
51 opObj20.addParameter(name='hei_th', value='[2000]',format='multiList')
52 opObj20.addParameter(name='proc', value='1')
53 opObj20.addParameter(name='channelList', value='0,0,1,1')
54 opObj20.addParameter(name='filec', value='weightfit', format='str')
55
56 #opObj22 = procUnitConfObj2.addOperation(name='HDFWriter', optype='other')
57 #opObj22.addParameter(name='path', value=pathFile)
58 #opObj22.addParameter(name='blocksPerFile', value='1')
59 #opObj22.addParameter(name='metadataList',value='heightList,timeZone')
60 #opObj22.addParameter(name='dataList',value='tmp_spectra_i,tmp_cspectra_i,tmp_spectra_c,tmp_cspectra_c,clean_num_aver,coh_num_aver,index,utctime')
61 #angles :-2.41116 3.01082
62 opObj21 = procUnitConfObj2.addOperation(name='EWDriftsEstimation', optype='other')
63 opObj21.addParameter(name='zenith', value='-2.41116, 3.01082', format='floatlist')
64 opObj21.addParameter(name='zenithCorrection', value='0.0', format='float')
65 opObj21.addParameter(name='fileDrifts', value=pathFile)
66
67 # Drifts en h5
68 one = {'gdlatr': 'lat', 'gdlonr': 'lon', 'spcst':'spcst','pl':'pl','cbadn':'cbadn','inttms': 'inttms','azdir7':'azw','eldir7':'elw','azdir8':'aze','eldir8':'ele','jro14':'jro14','jro15':'jro15','jro16':'jro16','nwlos':'nwlos'}
69 two = {
70 'range': ('params', 0),
71 'gdalt': ('params', 1),
72 'VIPN': ('params', 2),
73 'dvipn': ('params', 3),
74 'vipe': ('params', 4),
75 'dvipe': ('params', 5),
76 'vi7': ('params', 6),
77 'dvi7': ('params', 7),
78 'vi8': ('params', 8),
79 'dvi8': ('params', 9),
80 'PAIWL': ('params', 10),
81 'pacwl': ('params', 11),
82 'pbiwl': ('params', 12),
83 'pbcwl': ('params', 13),
84 'pciel': ('params', 14),
85 'pccel': ('params', 15),
86 'pdiel': ('params', 16),
87 'pdcel': ('params', 17),
88 'jro10': ('params', 18),
89 'jro11': ('params', 19)
90 } #writer
91 ind = ['gdalt']
92
93 #f=open('/home/roberto/moder_test.txt','r')
94 #file_contents=f.read()
95
96 meta = {
97 'kinst': 10, #instrument code
98 'kindat': 1910, #type of data
99 'catalog': {
100 'principleInvestigator': 'Danny Scipión',
101 'expPurpose': 'Drifts'#,
102 #'sciRemarks': file_contents
103 },
104 'header': {
105 'analyst': 'Danny Scipión'
106 }
107 }
108 #f.close()
109
110 op_writer = procUnitConfObj2.addOperation(name='MADWriter')
111 op_writer.addParameter(name='path', value=pathFile)
112 op_writer.addParameter(name='format', value='hdf5')
113 op_writer.addParameter(name='oneDDict', value=json.dumps(one))
114 op_writer.addParameter(name='twoDDict', value=json.dumps(two))
115 op_writer.addParameter(name='ind2DList', value=json.dumps(ind))
116 op_writer.addParameter(name='metadata', value=json.dumps(meta))
117
118 op_writer = procUnitConfObj2.addOperation(name='setHeightDriftsavg')
119
120 # Avg Drifts
121 one_avg = {'gdlatr': 'lat', 'gdlonr': 'lon', 'spcst':'spcst','pl':'pl','cbadn':'cbadn','inttms': 'inttms'}
122 two_avg = {
123 'range': ('params_avg', 4),
124 'gdalt': ('params_avg', 5),
125 'altav': ('params_avg', 6),
126 'VIPN': ('params_avg', 0),
127 'dvipn': ('params_avg', 1),
128 'vipe': ('params_avg', 2),
129 'dvipe': ('params_avg', 3)
130 }
131 ind_avg = ['gdalt']
132 meta = {
133 'kinst': 10, #instrument code
134 'kindat': 1911, #type of data
135 'catalog': {
136 'principleInvestigator': 'Danny Scipión',
137 'expPurpose': 'Drifts'#,
138 #'sciRemarks': file_contents
139 },
140 'header': {
141 'analyst': 'Danny Scipión'
142 }
143 }
144
145 op_writer = procUnitConfObj2.addOperation(name='MADWriter')
146 op_writer.addParameter(name='path', value=pathFileavg)
147 op_writer.addParameter(name='format', value='hdf5')
148 op_writer.addParameter(name='oneDDict', value=json.dumps(one_avg))
149 op_writer.addParameter(name='twoDDict', value=json.dumps(two_avg))
150 op_writer.addParameter(name='ind2DList', value=json.dumps(ind_avg))
151 op_writer.addParameter(name='metadata', value=json.dumps(meta))
152
153 op_writer = procUnitConfObj2.addOperation(name='setHeightDrifts')
154
155 opObj24 = procUnitConfObj2.addOperation(name='SpectralMomentsPlot', optype='other')
156 opObj24.addParameter(name='id', value='1', format='int')
157 ### # opObj14.addParameter(name='wintitle', value='Spectral Averaged', format='str')
158 opObj24.addParameter(name='save', value=pathFigure, format='str')
159 ###opObj24.addParameter(name='save_period', value=1, format='int')
160 opObj24.addParameter(name='zmin', value='-8', format='int')
161 opObj24.addParameter(name='zmax', value='16', format='int')
162 opObj24.addParameter(name='xaxis', value='Velocity', format='str')
163
164 #
165 titles=('SNR,Vertical Drifts,Zonal Drifts')
166 opObj23 = procUnitConfObj2.addOperation(name='GenericRTIPlot')
167 #opObj23.addParameter(name='colormaps', value='jet,RdBu_r,RdBu_r')
168 opObj23.addParameter(name='colormaps', value='jro,seismic,seismic')
169 #opObj23.addParameter(name='colormaps', value='jro,bwr,bwr')
170 opObj23.addParameter(name='attr_data', value='data_snr1,data_output')
171 opObj23.addParameter(name='wintitle', value='EW Drifts')
172 opObj23.addParameter(name='save', value=pathFigure)
173 opObj23.addParameter(name='titles', value=titles)
174 opObj23.addParameter(name='zfactors', value='1,1,1')
175 opObj23.addParameter(name='zlimits', value='(0,13),(-50,50),(-150,150)')
176 opObj23.addParameter(name='cb_labels', value='dB,m/s,m/s')
177 opObj23.addParameter(name='throttle', value='1')
178 opObj23.addParameter(name='xmin', value=xmin)
179 opObj23.addParameter(name='xmax', value=xmax)
180 #opObj23.addParameter(name='exp_code', value='110', format='int')
181 #opObj23.addParameter(name='server', value='10.10.110.243:4444', format='int')
182 #opObj23.addParameter(name='tag', value= 'jicamarca', format='str')
183
184 #--------------------------------------------------------------------------------------------------
185
186 controllerObj.start()
Show file before | Show file after | Hide comments
1 NO CONTENT: modified file
The requested commit or file is too big and content was truncated. Show full diff
General Comments 0
You need to be logged in to leave comments. Login now