JASMET30_Winds.py
90 lines
| 5.1 KiB
| text/x-python
|
PythonLexer
|
r844 | import os, sys | |
path = os.path.split(os.getcwd())[0] | |||
path = os.path.split(path)[0] | |||
sys.path.insert(0, path) | |||
from schainpy.controller import Project | |||
controllerObj = Project() | |||
|
r856 | controllerObj.setup(id = '005', name='script05', description="JASMET Wind Estimation") | |
#-------------------------------------- Setup ----------------------------------------- | |||
#Verificar estas variables | |||
|
r844 | ||
|
r856 | #Path donde estan los archivos HDF5 de meteoros | |
path = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/meteor') | |||
|
r844 | ||
|
r856 | #Path para los graficos | |
pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/graphics') | |||
|
r844 | ||
|
r856 | #Path donde se almacenaran las estimaciones de vientos | |
pathfile = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/phase') | |||
#Fechas para busqueda de archivos | |||
startDate = '2016/08/20' | |||
endDate = '2016/08/30' | |||
#Horas para busqueda de archivos | |||
|
r844 | startTime = '00:00:00' | |
endTime = '23:59:59' | |||
|
r856 | ||
#Offsets optimos obtenidos con OptimumOffset.py | |||
phaseOffsets = '-2.84, -1.77, 11.94, 9.71' | |||
|
r844 | #------------------------------------------------------------------------------------------------ | |
readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader', | |||
path=path, | |||
|
r856 | startDate=startDate, | |
endDate=endDate, | |||
|
r844 | startTime=startTime, | |
endTime=endTime, | |||
walk=1) | |||
#-------------------------------------------------------------------------------------------------- | |||
procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId()) | |||
|
r856 | opObj10 = procUnitConfObj1.addOperation(name='CorrectSMPhases',optype='other') | |
opObj10.addParameter(name='phaseOffsets', value=phaseOffsets, format='floatlist') | |||
|
r844 | ||
opObj13 = procUnitConfObj1.addOperation(name='SkyMapPlot', optype='other') | |||
opObj13.addParameter(name='id', value='1', format='int') | |||
opObj13.addParameter(name='wintitle', value='Sky Map', format='str') | |||
opObj13.addParameter(name='save', value='1', format='bool') | |||
opObj13.addParameter(name='figpath', value=pathfig, format='str') | |||
opObj13.addParameter(name='ftp', value='1', format='int') | |||
opObj13.addParameter(name='exp_code', value='15', format='int') | |||
opObj13.addParameter(name='sub_exp_code', value='1', format='int') | |||
opObj13.addParameter(name='tmin', value='0', format='int') | |||
opObj13.addParameter(name='tmax', value='24', format='int') | |||
opObj22 = procUnitConfObj1.addOperation(name='WindProfiler', optype='other') | |||
opObj22.addParameter(name='technique', value='Meteors', format='str') | |||
opObj22.addParameter(name='nHours', value='1', format='float') | |||
opObj22.addParameter(name='hmin', value='70', format='float') | |||
opObj22.addParameter(name='hmax', value='120', format='float') | |||
opObj23 = procUnitConfObj1.addOperation(name='WindProfilerPlot', optype='other') | |||
opObj23.addParameter(name='id', value='2', format='int') | |||
opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str') | |||
opObj23.addParameter(name='save', value='1', format='bool') | |||
opObj23.addParameter(name='figpath', value = pathfig, format='str') | |||
opObj23.addParameter(name='zmin', value='-140', format='int') | |||
opObj23.addParameter(name='zmax', value='140', format='int') | |||
|
r856 | opObj23.addParameter(name='xmin', value='0', format='float') | |
opObj23.addParameter(name='xmax', value='24', format='float') | |||
|
r844 | opObj23.addParameter(name='ymin', value='70', format='float') | |
opObj23.addParameter(name='ymax', value='110', format='float') | |||
|
r856 | opObj33 = procUnitConfObj1.addOperation(name='ParamWriter', optype='other') | |
opObj33.addParameter(name='path', value=pathfile) | |||
opObj33.addParameter(name='blocksPerFile', value='1000', format='int') | |||
opObj33.addParameter(name='metadataList',value='type,outputInterval,timeZone',format='list') | |||
opObj33.addParameter(name='dataList',value='data_output,utctime',format='list') | |||
|
r844 | #-------------------------------------------------------------------------------------------------- | |
print "Escribiendo el archivo XML" | |||
|
r856 | controllerObj.writeXml("JASMET05.xml") | |
|
r844 | print "Leyendo el archivo XML" | |
|
r856 | controllerObj.readXml("JASMET05.xml") | |
|
r844 | ||
controllerObj.createObjects() | |||
controllerObj.connectObjects() | |||
controllerObj.run() |