##// END OF EJS Templates
schain-jc
RSS
Corrections to JASMET scripts
Julio Valdez -
r856:499360e025db
parent child
Show More
Show file before | Show file after | Hide comments
@@ -1,4 +1,3
1
2 import os, sys
1 import os, sys
3
2
4 path = os.path.split(os.getcwd())[0]
3 path = os.path.split(os.getcwd())[0]
@@ -14,18 +13,18 controllerObj.setup(id = '004', name='script04', description="JASMET Phase Calib
14 #-------------------------------------- Setup -----------------------------------------
13 #-------------------------------------- Setup -----------------------------------------
15 #Verificar estas variables
14 #Verificar estas variables
16
15
17 #Path para los archivos
16 #Path donde estan los archivos HDF5 de meteoros
18 path = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/meteor')
17 path = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/meteor')
19
18
20 #Path para los graficos
19 #Path para los graficos
21 pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/graphics')
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 pathfile = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/phase')
23 pathfile = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/phase')
25
24
26 #Fechas para busqueda de archivos
25 #Fechas para busqueda de archivos
27 startDate = '2016/08/24'
26 startDate = '2016/08/20'
28 endDate = '2016/08/29'
27 endDate = '2016/08/30'
29 #Horas para busqueda de archivos
28 #Horas para busqueda de archivos
30 startTime = '00:00:00'
29 startTime = '00:00:00'
31 endTime = '23:59:59'
30 endTime = '23:59:59'
Show file before | Show file after | Hide comments
@@ -1,6 +1,3
1 # DIAS 19 Y 20 FEB 2014
2 # Comprobacion de Resultados DBS con SA
3
4 import os, sys
1 import os, sys
5
2
6 path = os.path.split(os.getcwd())[0]
3 path = os.path.split(os.getcwd())[0]
@@ -10,36 +7,43 sys.path.insert(0, path)
10
7
11 from schainpy.controller import Project
8 from schainpy.controller import Project
12
9
13 desc = "JASMET Experiment Test"
14 filename = "JASMETtest.xml"
15
16 controllerObj = Project()
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
19 #Path para los graficos
21 path= os.path.join(os.environ['HOME'],'Pictures/last_campaign/meteor')
20 pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/graphics')
22 pathfile2 = os.path.join(os.environ['HOME'],'Pictures/last_campaign/winds')
23 pathfig = os.path.join(os.environ['HOME'],'Pictures/last_campaign/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 startTime = '00:00:00'
29 startTime = '00:00:00'
26 endTime = '23:59:59'
30 endTime = '23:59:59'
27 xmin ='0.0'
31
28 xmax = '24.0'
32 #Offsets optimos obtenidos con OptimumOffset.py
29
33 phaseOffsets = '-2.84, -1.77, 11.94, 9.71'
30 #------------------------------------------------------------------------------------------------
34 #------------------------------------------------------------------------------------------------
31 readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader',
35 readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader',
32 path=path,
36 path=path,
33 startDate='2016/06/02',
37 startDate=startDate,
34 endDate='2017/06/03',
38 endDate=endDate,
35 startTime=startTime,
39 startTime=startTime,
36 endTime=endTime,
40 endTime=endTime,
37 walk=1)
41 walk=1)
38 #--------------------------------------------------------------------------------------------------
42 #--------------------------------------------------------------------------------------------------
39
43
40 procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId())
44 procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId())
41 opObj10 = procUnitConfObj1.addOperation(name='CorrectMeteorPhases')
45 opObj10 = procUnitConfObj1.addOperation(name='CorrectSMPhases',optype='other')
42 opObj10.addParameter(name='phaseOffsets', value='3.4,-3.6,19.4,0.1', format='floatlist')
46 opObj10.addParameter(name='phaseOffsets', value=phaseOffsets, format='floatlist')
43
47
44 opObj13 = procUnitConfObj1.addOperation(name='SkyMapPlot', optype='other')
48 opObj13 = procUnitConfObj1.addOperation(name='SkyMapPlot', optype='other')
45 opObj13.addParameter(name='id', value='1', format='int')
49 opObj13.addParameter(name='id', value='1', format='int')
@@ -65,31 +69,22 opObj23.addParameter(name='save', value='1', format='bool')
65 opObj23.addParameter(name='figpath', value = pathfig, format='str')
69 opObj23.addParameter(name='figpath', value = pathfig, format='str')
66 opObj23.addParameter(name='zmin', value='-140', format='int')
70 opObj23.addParameter(name='zmin', value='-140', format='int')
67 opObj23.addParameter(name='zmax', value='140', format='int')
71 opObj23.addParameter(name='zmax', value='140', format='int')
68 # opObj12.addParameter(name='zmin_ver', value='-0.8', format='float')
72 opObj23.addParameter(name='xmin', value='0', format='float')
69 # opObj12.addParameter(name='zmax_ver', value='0.8', format='float')
73 opObj23.addParameter(name='xmax', value='24', 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')
75 opObj23.addParameter(name='ymin', value='70', format='float')
74 opObj23.addParameter(name='ymin', value='70', format='float')
76 opObj23.addParameter(name='ymax', value='110', format='float')
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')
77 opObj33 = procUnitConfObj1.addOperation(name='ParamWriter', optype='other')
82 # opObj24.addParameter(name='path', value=pathfile2)
78 opObj33.addParameter(name='path', value=pathfile)
83 # opObj24.addParameter(name='blocksPerFile', value='1000', format='int')
79 opObj33.addParameter(name='blocksPerFile', value='1000', format='int')
84 # opObj24.addParameter(name='metadataList',value='type,outputInterval,heightList,timeZone',format='list')
80 opObj33.addParameter(name='metadataList',value='type,outputInterval,timeZone',format='list')
85 # opObj24.addParameter(name='dataList',value='data_output,utctime,utctimeInit',format='list')
81 opObj33.addParameter(name='dataList',value='data_output,utctime',format='list')
86
87 #--------------------------------------------------------------------------------------------------
82 #--------------------------------------------------------------------------------------------------
88
83
89 print "Escribiendo el archivo XML"
84 print "Escribiendo el archivo XML"
90 controllerObj.writeXml(filename)
85 controllerObj.writeXml("JASMET05.xml")
91 print "Leyendo el archivo XML"
86 print "Leyendo el archivo XML"
92 controllerObj.readXml(filename)
87 controllerObj.readXml("JASMET05.xml")
93
88
94 controllerObj.createObjects()
89 controllerObj.createObjects()
95 controllerObj.connectObjects()
90 controllerObj.connectObjects()
General Comments 0
You need to be logged in to leave comments. Login now