JASMET30_Winds.py
85 lines
| 4.9 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() | ||||
|
r1383 | controllerObj.setup(id='005', name='script05', description="JASMET Wind Estimation") | ||
|
r856 | |||
#-------------------------------------- Setup ----------------------------------------- | ||||
|
r1383 | # Verificar estas variables | ||
|
r844 | |||
|
r1383 | # Path donde estan los archivos HDF5 de meteoros | ||
path = os.path.join(os.environ['HOME'], 'Pictures/JASMET30_mp/201608/meteor') | ||||
|
r844 | |||
|
r1383 | # Path para los graficos | ||
pathfig = os.path.join(os.environ['HOME'], 'Pictures/JASMET30_mp/201608/graphics') | ||||
|
r844 | |||
|
r1383 | # Path donde se almacenaran las estimaciones de vientos | ||
pathfile = os.path.join(os.environ['HOME'], 'Pictures/JASMET30_mp/201608/phase') | ||||
|
r856 | |||
|
r1383 | # Fechas para busqueda de archivos | ||
|
r993 | startDate = '2016/08/29' | ||
endDate = '2016/09/11' | ||||
|
r1383 | # Horas para busqueda de archivos | ||
|
r844 | startTime = '00:00:00' | ||
endTime = '23:59:59' | ||||
|
r856 | |||
|
r1383 | # Offsets optimos obtenidos con OptimumOffset.py | ||
|
r993 | phaseOffsets = '-2.84, -1.77, 11.94, 9.71' | ||
phaseOffsets = '-5.86, -0.93, -7.29, 23.35' | ||||
|
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()) | ||||
|
r1383 | opObj10 = procUnitConfObj1.addOperation(name='CorrectSMPhases', optype='other') | ||
|
r856 | 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') | ||||
|
r1383 | opObj23.addParameter(name='figpath', value=pathfig, format='str') | ||
|
r844 | 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') | ||||
|
r1383 | opObj33.addParameter(name='metadataList', value='type,outputInterval,timeZone', format='list') | ||
opObj33.addParameter(name='dataList', value='data_output,utctime', format='list') | ||||
|
r844 | #-------------------------------------------------------------------------------------------------- | ||
|
r1383 | controllerObj.start() | ||