JASMET30_Winds.py
84 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() | |||
|
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 | |
|
r993 | path = os.path.join(os.environ['HOME'],'Pictures/JASMET30_mp/201608/meteor') | |
|
r844 | ||
|
r856 | #Path para los graficos | |
|
r993 | pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30_mp/201608/graphics') | |
|
r844 | ||
|
r856 | #Path donde se almacenaran las estimaciones de vientos | |
|
r993 | pathfile = os.path.join(os.environ['HOME'],'Pictures/JASMET30_mp/201608/phase') | |
|
r856 | ||
#Fechas para busqueda de archivos | |||
|
r993 | startDate = '2016/08/29' | |
endDate = '2016/09/11' | |||
|
r856 | #Horas para busqueda de archivos | |
|
r844 | startTime = '00:00:00' | |
endTime = '23:59:59' | |||
|
r856 | ||
#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()) | |||
|
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 | #-------------------------------------------------------------------------------------------------- | |
|
r993 | controllerObj.start() |