##// 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,4 +1,3
1
2 1 import os, sys
3 2
4 3 path = os.path.split(os.getcwd())[0]
@@ -14,18 +13,18 controllerObj.setup(id = '004', name='script04', description="JASMET Phase Calib
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'
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 1 import os, sys
5 2
6 3 path = os.path.split(os.getcwd())[0]
@@ -10,36 +7,43 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')
@@ -65,31 +69,22 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()
General Comments 0
You need to be logged in to leave comments. Login now