##// END OF EJS Templates
Fix all PlotData, add SpectraMean, CrossSpectra plots, now Parameters extends Spectra fix bugs in ParametersProc
Fix all PlotData, add SpectraMean, CrossSpectra plots, now Parameters extends Spectra fix bugs in ParametersProc

File last commit:

r856:499360e025db
r922:d680543828ae
Show More
JASMET30_Winds.py
90 lines | 5.1 KiB | text/x-python | PythonLexer
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()
controllerObj.setup(id = '005', name='script05', description="JASMET Wind Estimation")
#-------------------------------------- Setup -----------------------------------------
#Verificar estas variables
#Path donde estan los archivos HDF5 de meteoros
path = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/meteor')
#Path para los graficos
pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201608/graphics')
#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
startTime = '00:00:00'
endTime = '23:59:59'
#Offsets optimos obtenidos con OptimumOffset.py
phaseOffsets = '-2.84, -1.77, 11.94, 9.71'
#------------------------------------------------------------------------------------------------
readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader',
path=path,
startDate=startDate,
endDate=endDate,
startTime=startTime,
endTime=endTime,
walk=1)
#--------------------------------------------------------------------------------------------------
procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId())
opObj10 = procUnitConfObj1.addOperation(name='CorrectSMPhases',optype='other')
opObj10.addParameter(name='phaseOffsets', value=phaseOffsets, format='floatlist')
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')
opObj23.addParameter(name='xmin', value='0', format='float')
opObj23.addParameter(name='xmax', value='24', format='float')
opObj23.addParameter(name='ymin', value='70', format='float')
opObj23.addParameter(name='ymax', value='110', format='float')
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')
#--------------------------------------------------------------------------------------------------
print "Escribiendo el archivo XML"
controllerObj.writeXml("JASMET05.xml")
print "Leyendo el archivo XML"
controllerObj.readXml("JASMET05.xml")
controllerObj.createObjects()
controllerObj.connectObjects()
controllerObj.run()