##// END OF EJS Templates
25/04/18
25/04/18

File last commit:

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