##// END OF EJS Templates
schain
RSS
Corrections to JASMET scripts
Julio Valdez -
r856:499360e025db
parent child
Show More
Show file before | Show file after | Hide comments
@@ -1,78 +1,77
1
2 1 import os, sys
3 2
4 3 path = os.path.split(os.getcwd())[0]
5 4 path = os.path.split(path)[0]
6 5
7 6 sys.path.insert(0, path)
8 7
9 8 from schainpy.controller import Project
10 9
11 10 controllerObj = Project()
12 11 controllerObj.setup(id = '004', name='script04', description="JASMET Phase Calibration")
13 12
14 13 #-------------------------------------- Setup -----------------------------------------
15 14 #Verificar estas variables
16 15
17 #Path para los archivos
16 #Path donde estan los archivos HDF5 de meteoros
18 17 path = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/meteor')
19 18
20 19 #Path para los graficos
21 20 pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/graphics')
22 21
23 #Path para los archivos HDF5 de meteoros
22 #Path donde se almacenaran las fases calculadas
24 23 pathfile = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/phase')
25 24
26 25 #Fechas para busqueda de archivos
27 startDate = '2016/08/24'
28 endDate = '2016/08/29'
26 startDate = '2016/08/20'
27 endDate = '2016/08/30'
29 28 #Horas para busqueda de archivos
30 29 startTime = '00:00:00'
31 30 endTime = '23:59:59'
32 31
33 32 #------------------------------------------------------------------------------------------------
34 33 readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader',
35 34 path=path,
36 35 startDate='2016/06/02',
37 36 endDate='2017/06/03',
38 37 startTime=startTime,
39 38 endTime=endTime,
40 39 walk=1)
41 40
42 41 #--------------------------------------------------------------------------------------------------
43 42
44 43 procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId())
45 44 # #
46 45 opObj31 = procUnitConfObj1.addOperation(name='SMPhaseCalibration', optype='other')
47 46 opObj31.addParameter(name='nHours', value='1', format='float')
48 47 opObj31.addParameter(name='hmin', value='60', format='float')
49 48 opObj31.addParameter(name='hmax', value='120', format='float')
50 49 # opObj31.addParameter(name='channelPositions', value='(2.5,0),(0,2.5),(0,0),(0,4.5),(-2,0)', format='pairslist')
51 50
52 51 opObj32 = procUnitConfObj1.addOperation(name='PhasePlot', optype='other')
53 52 opObj32.addParameter(name='id', value='201', format='int')
54 53 opObj32.addParameter(name='wintitle', value='PhaseCalibration', format='str')
55 54 opObj32.addParameter(name='save', value='1', format='bool')
56 55 opObj32.addParameter(name='xmin', value='0', format='float')
57 56 opObj32.addParameter(name='xmax', value='24', format='float')
58 57 opObj32.addParameter(name='ymin', value='-180', format='float')
59 58 opObj32.addParameter(name='ymax', value='180', format='float')
60 59 opObj32.addParameter(name='figpath', value=pathfig, format='str')
61 60 # #
62 61 opObj33 = procUnitConfObj1.addOperation(name='ParamWriter', optype='other')
63 62 opObj33.addParameter(name='path', value=pathfile)
64 63 opObj33.addParameter(name='blocksPerFile', value='1000', format='int')
65 64 opObj33.addParameter(name='metadataList',value='type,outputInterval,timeZone',format='list')
66 65 opObj33.addParameter(name='dataList',value='data_output,utctime',format='list')
67 66 # # opObj25.addParameter(name='mode',value='1,0,0',format='intlist')
68 67
69 68 #--------------------------------------------------------------------------------------------------
70 69
71 70 print "Escribiendo el archivo XML"
72 71 controllerObj.writeXml("JASMET04.xml")
73 72 print "Leyendo el archivo XML"
74 73 controllerObj.readXml("JASMET04.xml")
75 74
76 75 controllerObj.createObjects()
77 76 controllerObj.connectObjects()
78 77 controllerObj.run() No newline at end of file
Show file before | Show file after | Hide comments
@@ -1,96 +1,91
1 # DIAS 19 Y 20 FEB 2014
2 # Comprobacion de Resultados DBS con SA
3
4 1 import os, sys
5 2
6 3 path = os.path.split(os.getcwd())[0]
7 4 path = os.path.split(path)[0]
8 5
9 6 sys.path.insert(0, path)
10 7
11 8 from schainpy.controller import Project
12 9
13 desc = "JASMET Experiment Test"
14 filename = "JASMETtest.xml"
15
16 10 controllerObj = Project()
11 controllerObj.setup(id = '005', name='script05', description="JASMET Wind Estimation")
12
13 #-------------------------------------- Setup -----------------------------------------
14 #Verificar estas variables
17 15
18 controllerObj.setup(id = '191', name='test01', description=desc)
16 #Path donde estan los archivos HDF5 de meteoros
17 path = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/meteor')
19 18
20 #Verificar
21 path= os.path.join(os.environ['HOME'],'Pictures/last_campaign/meteor')
22 pathfile2 = os.path.join(os.environ['HOME'],'Pictures/last_campaign/winds')
23 pathfig = os.path.join(os.environ['HOME'],'Pictures/last_campaign/graphics')
19 #Path para los graficos
20 pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/graphics')
24 21
22 #Path donde se almacenaran las estimaciones de vientos
23 pathfile = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/phase')
24
25 #Fechas para busqueda de archivos
26 startDate = '2016/08/20'
27 endDate = '2016/08/30'
28 #Horas para busqueda de archivos
25 29 startTime = '00:00:00'
26 30 endTime = '23:59:59'
27 xmin ='0.0'
28 xmax = '24.0'
29
31
32 #Offsets optimos obtenidos con OptimumOffset.py
33 phaseOffsets = '-2.84, -1.77, 11.94, 9.71'
30 34 #------------------------------------------------------------------------------------------------
31 35 readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader',
32 36 path=path,
33 startDate='2016/06/02',
34 endDate='2017/06/03',
37 startDate=startDate,
38 endDate=endDate,
35 39 startTime=startTime,
36 40 endTime=endTime,
37 41 walk=1)
38 42 #--------------------------------------------------------------------------------------------------
39 43
40 44 procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId())
41 opObj10 = procUnitConfObj1.addOperation(name='CorrectMeteorPhases')
42 opObj10.addParameter(name='phaseOffsets', value='3.4,-3.6,19.4,0.1', format='floatlist')
45 opObj10 = procUnitConfObj1.addOperation(name='CorrectSMPhases',optype='other')
46 opObj10.addParameter(name='phaseOffsets', value=phaseOffsets, format='floatlist')
43 47
44 48 opObj13 = procUnitConfObj1.addOperation(name='SkyMapPlot', optype='other')
45 49 opObj13.addParameter(name='id', value='1', format='int')
46 50 opObj13.addParameter(name='wintitle', value='Sky Map', format='str')
47 51 opObj13.addParameter(name='save', value='1', format='bool')
48 52 opObj13.addParameter(name='figpath', value=pathfig, format='str')
49 53 opObj13.addParameter(name='ftp', value='1', format='int')
50 54 opObj13.addParameter(name='exp_code', value='15', format='int')
51 55 opObj13.addParameter(name='sub_exp_code', value='1', format='int')
52 56 opObj13.addParameter(name='tmin', value='0', format='int')
53 57 opObj13.addParameter(name='tmax', value='24', format='int')
54 58
55 59 opObj22 = procUnitConfObj1.addOperation(name='WindProfiler', optype='other')
56 60 opObj22.addParameter(name='technique', value='Meteors', format='str')
57 61 opObj22.addParameter(name='nHours', value='1', format='float')
58 62 opObj22.addParameter(name='hmin', value='70', format='float')
59 63 opObj22.addParameter(name='hmax', value='120', format='float')
60 64
61 65 opObj23 = procUnitConfObj1.addOperation(name='WindProfilerPlot', optype='other')
62 66 opObj23.addParameter(name='id', value='2', format='int')
63 67 opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str')
64 68 opObj23.addParameter(name='save', value='1', format='bool')
65 69 opObj23.addParameter(name='figpath', value = pathfig, format='str')
66 70 opObj23.addParameter(name='zmin', value='-140', format='int')
67 71 opObj23.addParameter(name='zmax', value='140', format='int')
68 # opObj12.addParameter(name='zmin_ver', value='-0.8', format='float')
69 # opObj12.addParameter(name='zmax_ver', value='0.8', format='float')
70 # opObj23.addParameter(name='SNRmin', value='-10', format='int')
71 # opObj23.addParameter(name='SNRmax', value='60', format='int')
72 # opObj23.addParameter(name='SNRthresh', value='0', format='float')
73 opObj23.addParameter(name='xmin', value=xmin, format='float')
74 opObj23.addParameter(name='xmax', value=xmax, format='float')
72 opObj23.addParameter(name='xmin', value='0', format='float')
73 opObj23.addParameter(name='xmax', value='24', format='float')
75 74 opObj23.addParameter(name='ymin', value='70', format='float')
76 75 opObj23.addParameter(name='ymax', value='110', format='float')
77 # opObj23.addParameter(name='ftp', value='1', format='int')
78 # opObj23.addParameter(name='exp_code', value='15', format='int')
79 # opObj23.addParameter(name='sub_exp_code', value='1', format='int')
80 76
81 # opObj24 = procUnitConfObj1.addOperation(name='HDF5Writer', optype='other')
82 # opObj24.addParameter(name='path', value=pathfile2)
83 # opObj24.addParameter(name='blocksPerFile', value='1000', format='int')
84 # opObj24.addParameter(name='metadataList',value='type,outputInterval,heightList,timeZone',format='list')
85 # opObj24.addParameter(name='dataList',value='data_output,utctime,utctimeInit',format='list')
86
77 opObj33 = procUnitConfObj1.addOperation(name='ParamWriter', optype='other')
78 opObj33.addParameter(name='path', value=pathfile)
79 opObj33.addParameter(name='blocksPerFile', value='1000', format='int')
80 opObj33.addParameter(name='metadataList',value='type,outputInterval,timeZone',format='list')
81 opObj33.addParameter(name='dataList',value='data_output,utctime',format='list')
87 82 #--------------------------------------------------------------------------------------------------
88 83
89 84 print "Escribiendo el archivo XML"
90 controllerObj.writeXml(filename)
85 controllerObj.writeXml("JASMET05.xml")
91 86 print "Leyendo el archivo XML"
92 controllerObj.readXml(filename)
87 controllerObj.readXml("JASMET05.xml")
93 88
94 89 controllerObj.createObjects()
95 90 controllerObj.connectObjects()
96 91 controllerObj.run() No newline at end of file
General Comments 0
You need to be logged in to leave comments. Login now