diff --git a/schainpy/scripts/JASMET30_MetDet.py b/schainpy/scripts/JASMET30_MetDet.py new file mode 100644 index 0000000..f8f2216 --- /dev/null +++ b/schainpy/scripts/JASMET30_MetDet.py @@ -0,0 +1,95 @@ +# DIAS 19 Y 20 FEB 2014 +# Comprobacion de Resultados DBS con SA + +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 + +desc = "JASMET Experiment Test" +filename = "JASMETtest.xml" + +controllerObj = Project() + +controllerObj.setup(id = '191', name='test01', description=desc) + +pathfile1 = os.path.join(os.environ['HOME'],'Pictures/last_campaign/meteor') +pathfig = os.path.join(os.environ['HOME'],'Pictures/last_campaign/graphics') + +path = '/mnt/jars/2016_08/NOCHE' +startTime = '00:00:00' +endTime = '08:59:59' +# +# path = '/mnt/jars/2016_08/DIA' +# startTime = '12:13:00' +# endTime = '23:59:59' + +# path = '/mnt/jars/2016_08/NOCHE' +# startTime = '15:00:00' +# endTime = '23:59:59' + +#------------------------------------------------------------------------------------------------ +readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader', + path=path, + startDate='2016/08/26', + endDate='2016/08/26', + startTime=startTime, + endTime=endTime, + online=0, + delay=30, + walk=1, + getblock=1, + blocktime=100) + +opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock') + +#-------------------------------------------------------------------------------------------------- + +procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId()) + +opObj00 = procUnitConfObj0.addOperation(name='selectChannels') +opObj00.addParameter(name='channelList', value='0,1,2,3,4', format='intlist') + +opObj01 = procUnitConfObj0.addOperation(name='setRadarFrequency') +opObj01.addParameter(name='frequency', value='30.e6', format='float') + +opObj00 = procUnitConfObj0.addOperation(name='interpolateHeights') +opObj00.addParameter(name='topLim', value='73', format='int') +opObj00.addParameter(name='botLim', value='69', format='int') + +opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') + +opObj12 = procUnitConfObj0.addOperation(name='CohInt', optype='other') +opObj12.addParameter(name='n', value='2', format='int') + +#-------------------------------------------------------------------------------------------------- + +procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj0.getId()) +# +opObj10 = procUnitConfObj1.addOperation(name='MeteorDetection') +opObj10.addParameter(name='azimuth', value='45', format='float') +opObj10.addParameter(name='hmin', value='60', format='float') +opObj10.addParameter(name='hmax', value='120', format='float') + +opObj12 = procUnitConfObj1.addOperation(name='ParamWriter', optype='other') +opObj12.addParameter(name='path', value=pathfile1) +opObj12.addParameter(name='blocksPerFile', value='1000', format='int') +opObj12.addParameter(name='metadataList',value='type,heightList,paramInterval,timeZone',format='list') +opObj12.addParameter(name='dataList',value='data_param,utctime',format='list') +opObj12.addParameter(name='mode',value='2',format='int') +# # Tiene que ser de 3 dimensiones, append en lugar de aumentar una dimension +# +#-------------------------------------------------------------------------------------------------- + +print "Escribiendo el archivo XML" +controllerObj.writeXml(filename) +print "Leyendo el archivo XML" +controllerObj.readXml(filename) + +controllerObj.createObjects() +controllerObj.connectObjects() +controllerObj.run() \ No newline at end of file diff --git a/schainpy/scripts/JASMET30_PhaseCal.py b/schainpy/scripts/JASMET30_PhaseCal.py new file mode 100644 index 0000000..e0bc3b8 --- /dev/null +++ b/schainpy/scripts/JASMET30_PhaseCal.py @@ -0,0 +1,77 @@ +# DIAS 19 Y 20 FEB 2014 +# Comprobacion de Resultados DBS con SA + +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 + +desc = "JASMET Experiment Test" +filename = "JASMETtest.xml" + +controllerObj = Project() + +controllerObj.setup(id = '191', name='test01', description=desc) + +path= os.path.join(os.environ['HOME'],'Pictures/last_campaign/meteor') + +pathfile2 = os.path.join(os.environ['HOME'],'Pictures/last_campaign/phase') +pathfig = os.path.join(os.environ['HOME'],'Pictures/last_campaign/graphics') + + +startTime = '00:00:00' +endTime = '23:59:59' +# endTime = '00:01:01' +xmin ='0' +xmax = '24' + +#------------------------------------------------------------------------------------------------ +readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader', + path=path, + startDate='2016/06/02', + endDate='2017/06/03', + startTime=startTime, + endTime=endTime, + walk=1) + +#-------------------------------------------------------------------------------------------------- + +procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId()) +# # +opObj31 = procUnitConfObj1.addOperation(name='PhaseCalibration', optype='other') +opObj31.addParameter(name='nHours', value='1', format='float') +opObj31.addParameter(name='hmin', value='60', format='float') +opObj31.addParameter(name='hmax', value='120', format='float') +# opObj31.addParameter(name='channelPositions', value='(2.5,0),(0,2.5),(0,0),(0,4.5),(-2,0)', format='pairslist') + +opObj32 = procUnitConfObj1.addOperation(name='PhasePlot', optype='other') +opObj32.addParameter(name='id', value='201', format='int') +opObj32.addParameter(name='wintitle', value='PhaseCalibration', format='str') +opObj32.addParameter(name='save', value='1', format='bool') +opObj32.addParameter(name='xmin', value=xmin, format='float') +opObj32.addParameter(name='xmax', value=xmax, format='float') +opObj32.addParameter(name='ymin', value='-180', format='float') +opObj32.addParameter(name='ymax', value='180', format='float') +opObj32.addParameter(name='figpath', value=pathfig, format='str') +# # +opObj33 = procUnitConfObj1.addOperation(name='ParamWriter', optype='other') +opObj33.addParameter(name='path', value=pathfile2) +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') +# # opObj25.addParameter(name='mode',value='1,0,0',format='intlist') + +#-------------------------------------------------------------------------------------------------- + +print "Escribiendo el archivo XML" +controllerObj.writeXml(filename) +print "Leyendo el archivo XML" +controllerObj.readXml(filename) + +controllerObj.createObjects() +controllerObj.connectObjects() +controllerObj.run() \ No newline at end of file diff --git a/schainpy/scripts/JASMET30_Winds.py b/schainpy/scripts/JASMET30_Winds.py new file mode 100644 index 0000000..461f811 --- /dev/null +++ b/schainpy/scripts/JASMET30_Winds.py @@ -0,0 +1,96 @@ +# DIAS 19 Y 20 FEB 2014 +# Comprobacion de Resultados DBS con SA + +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 + +desc = "JASMET Experiment Test" +filename = "JASMETtest.xml" + +controllerObj = Project() + +controllerObj.setup(id = '191', name='test01', description=desc) + +#Verificar +path= os.path.join(os.environ['HOME'],'Pictures/last_campaign/meteor') +pathfile2 = os.path.join(os.environ['HOME'],'Pictures/last_campaign/winds') +pathfig = os.path.join(os.environ['HOME'],'Pictures/last_campaign/graphics') + +startTime = '00:00:00' +endTime = '23:59:59' +xmin ='0.0' +xmax = '24.0' + +#------------------------------------------------------------------------------------------------ +readUnitConfObj = controllerObj.addReadUnit(datatype='ParamReader', + path=path, + startDate='2016/06/02', + endDate='2017/06/03', + startTime=startTime, + endTime=endTime, + walk=1) +#-------------------------------------------------------------------------------------------------- + +procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=readUnitConfObj.getId()) +opObj10 = procUnitConfObj1.addOperation(name='CorrectMeteorPhases') +opObj10.addParameter(name='phaseOffsets', value='3.4,-3.6,19.4,0.1', 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') +# opObj12.addParameter(name='zmin_ver', value='-0.8', format='float') +# opObj12.addParameter(name='zmax_ver', value='0.8', format='float') +# opObj23.addParameter(name='SNRmin', value='-10', format='int') +# opObj23.addParameter(name='SNRmax', value='60', format='int') +# opObj23.addParameter(name='SNRthresh', value='0', format='float') +opObj23.addParameter(name='xmin', value=xmin, format='float') +opObj23.addParameter(name='xmax', value=xmax, format='float') +opObj23.addParameter(name='ymin', value='70', format='float') +opObj23.addParameter(name='ymax', value='110', format='float') +# opObj23.addParameter(name='ftp', value='1', format='int') +# opObj23.addParameter(name='exp_code', value='15', format='int') +# opObj23.addParameter(name='sub_exp_code', value='1', format='int') + +# opObj24 = procUnitConfObj1.addOperation(name='HDF5Writer', optype='other') +# opObj24.addParameter(name='path', value=pathfile2) +# opObj24.addParameter(name='blocksPerFile', value='1000', format='int') +# opObj24.addParameter(name='metadataList',value='type,outputInterval,heightList,timeZone',format='list') +# opObj24.addParameter(name='dataList',value='data_output,utctime,utctimeInit',format='list') + +#-------------------------------------------------------------------------------------------------- + +print "Escribiendo el archivo XML" +controllerObj.writeXml(filename) +print "Leyendo el archivo XML" +controllerObj.readXml(filename) + +controllerObj.createObjects() +controllerObj.connectObjects() +controllerObj.run() \ No newline at end of file