From d0b75259d16b6a2c1405dca149655f0d3fe070fd 2016-01-19 15:40:13 From: Julio Valdez Date: 2016-01-19 15:40:13 Subject: [PATCH] Updated JASMET scripts --- diff --git a/schainpy/scripts/Meteor_JASMET30MHz_Winds.py b/schainpy/scripts/Meteor_JASMET30MHz_Winds.py new file mode 100644 index 0000000..fad437b --- /dev/null +++ b/schainpy/scripts/Meteor_JASMET30MHz_Winds.py @@ -0,0 +1,176 @@ +# 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) + +#Experimentos + +#2014051 20 Feb 2014 +# path = '/home/soporte/Data/JASMET/JASMET_30/2014106' +# pathFigure = '/home/soporte/workspace/Graficos/JASMET/prueba1' + +#phases = 19.26890939, -7.53970335, 113.4375, 91.33928571 +remotefolder = "/home/wmaster/graficos" +path = '/media/joscanoa/84A65E64A65E5730/soporte/Data/JASMET/JASMET_30' + +pathfile1 = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201511/meteor') +pathfile2 = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201511/wind') +# opObj12 = procUnitConfObj1.addOperation(name='HDF5Writer', optype='other') +# opObj12.addParameter(name='path', value=pathfile1) +# opObj12.addParameter(name='blocksPerFile', value='1000', format='int') +# opObj12.addParameter(name='metadataList',value='type,inputUnit,heightList,paramInterval',format='list') +# opObj12.addParameter(name='dataList',value='data_param',format='list') +# opObj12.addParameter(name='mode',value='0',format='int') +#Tiene que ser de 3 dimensiones, append en lugar de aumentar una dimension + +# 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')11off/wind') + +pathfig = os.path.join(os.environ['HOME'],'Pictures/JASMET30/201511/graphics') + +startTime = '00:00:00' +endTime = '23:59:59' +xmin ='0.0' +xmax = '24.0' + +#------------------------------------------------------------------------------------------------ +readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader', + path=path, + startDate='2015/11/08', + endDate='2015/11/08', + startTime=startTime, + endTime=endTime, + online=0, + delay=5, + walk=1) + +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') +# opObj01.addParameter(name='frequency', value='50.e6', format='float') + +# opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') +# opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other') +# opObj11.addParameter(name='n', value='2', format='int') +#-------------------------------------------------------------------------------------------------- + +procUnitConfObj1 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj0.getId()) +procUnitConfObj1.addParameter(name='nSeconds', value='100', format='int') + +opObj10 = procUnitConfObj1.addOperation(name='MeteorDetection') +# opObj10.addParameter(name='predefinedPhaseShifts', value='19.27, -7.54, 0.0 ,113.44, 91.34', format='floatlist') +opObj10.addParameter(name='predefinedPhaseShifts', value='-23.8, 5.8, 0.0 ,-111.7, -87.2', format='floatlist') +# opObj10.addParameter(name='predefinedPhaseShifts', value='0, 0, 0, 0, 0', format='floatlist') +opObj10.addParameter(name='cohDetection', value='0', format='bool') +opObj10.addParameter(name='noise_multiple', value='4', format='int') +opObj10.addParameter(name='SNRThresh', value='5', format='float') +opObj10.addParameter(name='phaseThresh', value='20', format='float') +opObj10.addParameter(name='azimuth', value='45', format='float') +opObj10.addParameter(name='hmin', value='68', format='float') +opObj10.addParameter(name='hmax', value='120', format='float') +opObj10.addParameter(name='saveAll', value='1', format='bool') + +opObj12 = procUnitConfObj1.addOperation(name='HDF5Writer', optype='other') +opObj12.addParameter(name='path', value=pathfile1) +opObj12.addParameter(name='blocksPerFile', value='1000', format='int') +opObj12.addParameter(name='metadataList',value='type,inputUnit,heightList,paramInterval',format='list') +opObj12.addParameter(name='dataList',value='data_param',format='list') +opObj12.addParameter(name='mode',value='0',format='int') +#Tiene que ser de 3 dimensiones, append en lugar de aumentar una dimension + +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') + + +#-------------------------------------------------------------------------------------------------- +procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId()) + +opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other') +opObj22.addParameter(name='technique', value='Meteors', format='str') +opObj22.addParameter(name='nHours', value='1', format='float') +opObj22.addParameter(name='SNRThresh', value='12.0', format='float') +opObj22.addParameter(name='hmin', value='68', format='float') +opObj22.addParameter(name='hmax', value='120', format='float') + +opObj23 = procUnitConfObj2.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='-120', format='int') +opObj23.addParameter(name='zmax', value='120', 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='ftp', value='1', format='int') +opObj23.addParameter(name='exp_code', value='15', format='int') +opObj23.addParameter(name='sub_exp_code', value='1', format='int') +# opObj23.addParameter(name='timerange', value=str(12*60*60), format='float') + + +opObj24 = procUnitConfObj2.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,inputUnit,outputInterval',format='list') +opObj24.addParameter(name='dataList',value='data_output,utctime',format='list') + +#-------------------------------------------------------------------------------------------------- + +# procUnitConfObj4 = controllerObj.addProcUnit(name='SendToServer') +# procUnitConfObj4.addParameter(name='server', value='jro-app.igp.gob.pe', format='str') +# procUnitConfObj4.addParameter(name='username', value='wmaster', format='str') +# procUnitConfObj4.addParameter(name='password', value='mst2010vhf', format='str') +# procUnitConfObj4.addParameter(name='localfolder', value=pathfig, format='str') +# procUnitConfObj4.addParameter(name='remotefolder', value=remotefolder, format='str') +# procUnitConfObj4.addParameter(name='ext', value='.png', format='str') +# procUnitConfObj4.addParameter(name='period', value='240', format='int') +# procUnitConfObj4.addParameter(name='protocol', value='ftp', format='str') + +#-------------------------------------------------------------------------------------------------- + +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