schain Commit - r1525:4d24dade1cc5 · Jicamarca Repository

1

#!/usr/bin/env python

1

#!/usr/bin/env python

2

# Ing. AVP

2

# Ing. AVP

3

# 22/06/2022

3

# 22/06/2022

4

# ARCHIVO DE LECTURA y PLOT

4

# ARCHIVO DE LECTURA y PLOT

5

import matplotlib.pyplot as pl

5

import matplotlib.pyplot as pl

6

import matplotlib

6

import matplotlib

7

import wradlib

7

import wradlib

8

import numpy

8

import numpy

9

import warnings

9

import warnings

10

import argparse

10

import argparse

11

from wradlib.io import read_generic_hdf5

11

from wradlib.io import read_generic_hdf5

12

from wradlib.util import get_wradlib_data_file

12

from wradlib.util import get_wradlib_data_file

13

from plotting_codes import sophy_cb_tables

13

from plotting_codes import sophy_cb_tables

14

import os,time

14

import os,time

15

16

for name, cb_table in sophy_cb_tables:

16

for name, cb_table in sophy_cb_tables:

17

ncmap = matplotlib.colors.ListedColormap(cb_table, name=name)

17

ncmap = matplotlib.colors.ListedColormap(cb_table, name=name)

18

matplotlib.pyplot.register_cmap(cmap=ncmap)

18

matplotlib.pyplot.register_cmap(cmap=ncmap)

19

#LINUX bash: export WRADLIB_DATA=/path/to/wradlib-data

19

#LINUX bash: export WRADLIB_DATA=/path/to/wradlib-data

20

#example

21

#export WRADLIB_DATA="/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T18-00-00/"

20

warnings.filterwarnings('ignore')

22

warnings.filterwarnings('ignore')

21

PARAM = {

23

PARAM = {

22

'S': {'var': 'power','vmin': -45, 'vmax': -15, 'cmap': 'jet', 'label': 'Power','unit': 'dBm'},

24

'S': {'var': 'power','vmin': -45, 'vmax': -15, 'cmap': 'jet', 'label': 'Power','unit': 'dBm'},

23

'V': {'var': 'velocity', 'vmin': -10, 'vmax': 10 , 'cmap': 'sophy_v', 'label': 'Velocity','unit': 'm/s'},

25

'V': {'var': 'velocity', 'vmin': -10, 'vmax': 10 , 'cmap': 'sophy_v', 'label': 'Velocity','unit': 'm/s'},

24

'Z': {'var': 'reflectivity','vmin': -30, 'vmax': 80 , 'cmap': 'sophy_r','label': 'Reflectivity','unit': 'dBZ'},

26

'Z': {'var': 'reflectivity','vmin': -30, 'vmax': 80 , 'cmap': 'sophy_r','label': 'Reflectivity','unit': 'dBZ'},

25

'W': {'var': 'spectral_width', 'vmin': 0 , 'vmax': 12 , 'cmap': 'sophy_w','label': 'Spectral Width','unit': 'hz'}

27

'W': {'var': 'spectral_width', 'vmin': 0 , 'vmax': 12 , 'cmap': 'sophy_w','label': 'Spectral Width','unit': 'hz'}

26

}

28

}

27

class Readsophy():

29

class Readsophy():

28

def __init__(self):

30

def __init__(self):

29

self.list_file = None

31

self.list_file = None

30

self.grado = None

32

self.grado = None

31

self.variable = None

33

self.variable = None

32

self.save = None

34

self.save = None

33

self.range = None

35

self.range = None

34

36

35

def read_files(self,path_file,grado=None, variable=None):

37

def read_files(self,path_file,grado=None, variable=None):

36

filter= "_E"+str(grado)+".0_"+variable

38

filter= "_E"+str(grado)+".0_"+variable

37

validFilelist = []

39

validFilelist = []

38

fileList= os.listdir(path_file)

40

fileList= os.listdir(path_file)

39

for thisFile in fileList:

41

for thisFile in fileList:

40

if (os.path.splitext(thisFile)[0][-7:] != filter):

42

if (os.path.splitext(thisFile)[0][-7:] != filter):

41

#print("s_:",os.path.splitext(thisFile)[0][-7:])

43

#print("s_:",os.path.splitext(thisFile)[0][-7:])

42

continue

44

continue

43

validFilelist.append(thisFile)

45

validFilelist.append(thisFile)

44

validFilelist.sort()

46

validFilelist.sort()

45

return validFilelist

47

return validFilelist

46

48

47

def setup(self, path_file,grado,range,variable,save):

49

def setup(self, path_file,grado,range,variable,save):

48

self.path_file = path_file

50

self.path_file = path_file

49

self.range = range

51

self.range = range

50

self.grado = grado

52

self.grado = grado

51

self.variable = variable

53

self.variable = variable

52

self.save = save

54

self.save = save

53

self.list_file = self.read_files(path_file=self.path_file,grado=self.grado, variable=self.variable)

55

self.list_file = self.read_files(path_file=self.path_file,grado=self.grado, variable=self.variable)

54

56

55

def selectHeights(self,heightList,minHei,maxHei):

57

def selectHeights(self,heightList,minHei,maxHei):

56

58

57

if minHei and maxHei:

59

if minHei and maxHei:

58

if (minHei < heightList[0]):

60

if (minHei < heightList[0]):

59

minHei = heightList[0]

61

minHei = heightList[0]

60

if (maxHei > heightList[-1]):

62

if (maxHei > heightList[-1]):

61

maxHei = heightList[-1]

63

maxHei = heightList[-1]

62

minIndex = 0

64

minIndex = 0

63

maxIndex = 0

65

maxIndex = 0

64

heights = heightList

66

heights = heightList

65

67

66

inda = numpy.where(heights >= minHei)

68

inda = numpy.where(heights >= minHei)

67

indb = numpy.where(heights <= maxHei)

69

indb = numpy.where(heights <= maxHei)

68

70

69

try:

71

try:

70

minIndex = inda[0][0]

72

minIndex = inda[0][0]

71

except:

73

except:

72

minIndex = 0

74

minIndex = 0

73

75

74

try:

76

try:

75

maxIndex = indb[0][-1]

77

maxIndex = indb[0][-1]

76

except:

78

except:

77

maxIndex = len(heights)

79

maxIndex = len(heights)

78

80

79

new_heightList= self.selectHeightsByIndex(heightList=heightList,minIndex=minIndex, maxIndex=maxIndex)

81

new_heightList= self.selectHeightsByIndex(heightList=heightList,minIndex=minIndex, maxIndex=maxIndex)

80

82

81

return new_heightList, minIndex,maxIndex

83

return new_heightList, minIndex,maxIndex

82

84

83

def selectHeightsByIndex(self,heightList,minIndex, maxIndex):

85

def selectHeightsByIndex(self,heightList,minIndex, maxIndex):

84

86

85

if (minIndex < 0) or (minIndex > maxIndex):

87

if (minIndex < 0) or (minIndex > maxIndex):

86

raise ValueError("Height index range (%d,%d) is not valid" % (minIndex, maxIndex))

88

raise ValueError("Height index range (%d,%d) is not valid" % (minIndex, maxIndex))

87

89

88

if (maxIndex >= len(heightList)):

90

if (maxIndex >= len(heightList)):

89

maxIndex = len(heightList)

91

maxIndex = len(heightList)

90

92

91

new_h = heightList[minIndex:maxIndex]

93

new_h = heightList[minIndex:maxIndex]

92

return new_h

94

return new_h

93

95

96

def readAttributes(self,obj,variable):

97

var = PARAM[variable]['var']

98

unit = PARAM[variable]['unit']

99

cmap = PARAM[variable]['cmap']

100

vmin = PARAM[variable]['vmin']

101

vmax = PARAM[variable]['vmax']

102

label = PARAM[variable]['label']

103

var_ = 'Data/'+var+'/H'

104

data_arr = numpy.array(obj[var_]['data']) # data

105

utc_time = numpy.array(obj['Data/time']['data'])

106

data_azi = numpy.array(obj['Metadata/azimuth']['data']) # th

107

data_ele = numpy.array(obj["Metadata/elevation"]['data'])

108

heightList = numpy.array(obj["Metadata/range"]['data']) # r

109

110

return data_arr, utc_time, data_azi,data_ele, heightList,unit,cmap,vmin,vmax,label

111

94

def run(self):

112

def run(self):

95

count= 0

113

count= 0

96

len_files = len(self.list_file)

114

len_files = len(self.list_file)

97

115

98

for thisFile in self.list_file:

116

for thisFile in self.list_file:

99

count = count +1

117

count = count +1

100

fullpathfile = self.path_file + thisFile

118

fullpathfile = self.path_file + thisFile

101

filename = get_wradlib_data_file(fullpathfile)

119

filename = get_wradlib_data_file(fullpathfile)

102

test_hdf5 = read_generic_hdf5(filename)

120

test_hdf5 = read_generic_hdf5(filename)

103

121

104

var = PARAM[self.variable]['var']

122

# LECTURA

105

unit = PARAM[self.variable]['unit']

123

data_arr, utc_time, data_azi,data_ele, heightList,unit,cmap,vmin,vmax,label = self.readAttributes(obj= test_hdf5,variable=self.variable)

106

cmap = PARAM[self.variable]['cmap']

107

vmin = PARAM[self.variable]['vmin']

108

vmax = PARAM[self.variable]['vmax']

109

label = PARAM[self.variable]['label']

110

var_ = 'Data/'+var+'/H'

111

data_arr = numpy.array(test_hdf5[var_]['data']) # data

112

utc_time = numpy.array(test_hdf5['Data/time']['data'])

113

data_azi = numpy.array(test_hdf5['Metadata/azimuth']['data']) # th

114

data_ele = numpy.array(test_hdf5["Metadata/elevation"]['data'])

115

heightList = numpy.array(test_hdf5["Metadata/range"]['data']) # r

116

124

117

if self.range==0:

125

if self.range==0:

118

self.range == heightList[-1]

126

self.range == heightList[-1]

119

new_heightList,minIndex,maxIndex = self.selectHeights(heightList,0.06,self.range)

127

new_heightList,minIndex,maxIndex = self.selectHeights(heightList,0.06,self.range)

120

128

121

129

# TIEMPO

122

my_time = time.strftime('%Y-%m-%d %H:%M:%S', time.localtime(utc_time[0]))

130

my_time = time.strftime('%Y-%m-%d %H:%M:%S', time.localtime(utc_time[0]))

123

time_save = time.strftime('%Y%m%d_%H%M%S',time.localtime(utc_time[0]))

131

time_save = time.strftime('%Y%m%d_%H%M%S',time.localtime(utc_time[0]))

124

132

125

# PLOT DATA WITH ANNOTATION

133

# PLOT DATA WITH ANNOTATION

126

if count ==1:

134

if count ==1:

127

fig = pl.figure(figsize=(10,8))

135

fig = pl.figure(figsize=(10,8))

128

cgax, pm = wradlib.vis.plot_ppi(data_arr[:,minIndex:maxIndex],r=new_heightList,az=data_azi,rf=1,fig=fig, ax=111,proj='cg',cmap=cmap,vmin=vmin, vmax=vmax)

136

cgax, pm = wradlib.vis.plot_ppi(data_arr[:,minIndex:maxIndex],r=new_heightList,az=data_azi,rf=1,fig=fig, ax=111,proj='cg',cmap=cmap,vmin=vmin, vmax=vmax)

129

caax = cgax.parasites[0]

137

caax = cgax.parasites[0]

130

title = 'Simple PPI'+"-"+ my_time+" E."+self.grado

138

title = 'Simple PPI'+"-"+ my_time+" E."+self.grado

131

t = pl.title(title, fontsize=12,y=1.05)

139

t = pl.title(title, fontsize=12,y=1.05)

132

cbar = pl.gcf().colorbar(pm, pad=0.075)

140

cbar = pl.gcf().colorbar(pm, pad=0.075)

133

pl.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')

141

pl.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')

134

cbar.set_label(label+'[' + unit + ']')

142

cbar.set_label(label+'[' + unit + ']')

135

gh = cgax.get_grid_helper()

143

gh = cgax.get_grid_helper()

136

else:

144

else:

137

cgax, pm = wradlib.vis.plot_ppi(data_arr[:,minIndex:maxIndex],r=new_heightList,az=data_azi,rf=1,fig=fig, ax=111,proj='cg',cmap=cmap,vmin=vmin, vmax=vmax)

145

cgax, pm = wradlib.vis.plot_ppi(data_arr[:,minIndex:maxIndex],r=new_heightList,az=data_azi,rf=1,fig=fig, ax=111,proj='cg',cmap=cmap,vmin=vmin, vmax=vmax)

138

caax = cgax.parasites[0]

146

caax = cgax.parasites[0]

139

title = 'Simple PPI'+"-"+my_time+" E."+self.grado

147

title = 'Simple PPI'+"-"+my_time+" E."+self.grado

140

t = pl.title(title, fontsize=12,y=1.05)

148

t = pl.title(title, fontsize=12,y=1.05)

141

cbar = pl.gcf().colorbar(pm, pad=0.075)

149

cbar = pl.gcf().colorbar(pm, pad=0.075)

142

pl.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')

150

pl.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')

143

cbar.set_label(label+'[' + unit + ']')

151

cbar.set_label(label+'[' + unit + ']')

144

gh = cgax.get_grid_helper()

152

gh = cgax.get_grid_helper()

145

if self.save == 1:

153

if self.save == 1:

146

if count ==1:

154

if count ==1:

147

filename = "SOPHY"+"_"+time_save+"_"+"E."+self.grado+"_"+self.variable+".png"

155

filename = "SOPHY"+"_"+time_save+"_"+"E."+self.grado+"_"+self.variable+".png"

148

dir =self.variable+"_"+"E."+self.grado+"CH0/"

156

dir =self.variable+"_"+"E."+self.grado+"CH0/"

149

filesavepath = os.path.join(self.path_file,dir)

157

filesavepath = os.path.join(self.path_file,dir)

150

try:

158

try:

151

os.mkdir(filesavepath)

159

os.mkdir(filesavepath)

152

except:

160

except:

153

pass

161

pass

154

else:

162

else:

155

filename = "SOPHY"+"_"+time_save+"_"+"E."+self.grado+"_"+self.variable+".png"

163

filename = "SOPHY"+"_"+time_save+"_"+"E."+self.grado+"_"+self.variable+".png"

156

pl.savefig(filesavepath+filename)

164

pl.savefig(filesavepath+filename)

157

165

158

pl.pause(1)

166

pl.pause(1)

159

pl.clf()

167

pl.clf()

160

if count==len_files:

168

if count==len_files:

161

pl.close()

169

pl.close()

162

pl.show()

170

pl.show()

163

171

164

PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T18-00-00/"

172

PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T18-00-00/"

165

#PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T19-00-00/"

173

#PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T19-00-00/"

166

174

167

#PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-05-31T12-00-17/paramC0N36.0/2022-05-31T16-00-00/"

175

#PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-05-31T12-00-17/paramC0N36.0/2022-05-31T16-00-00/"

168

176

169

def main(args):

177

def main(args):

170

grado = args.grado

178

grado = args.grado

171

parameters = args.parameters

179

parameters = args.parameters

172

save = args.save

180

save = args.save

173

range = args.range

181

range = args.range

174

obj = Readsophy()

182

obj = Readsophy()

175

for param in parameters:

183

for param in parameters:

176

obj.setup(path_file = PATH,grado = grado,range=range, variable=param,save=int(save))

184

obj.setup(path_file = PATH,grado = grado,range=range, variable=param,save=int(save))

177

obj.run()

185

obj.run()

178

186

179

if __name__ == '__main__':

187

if __name__ == '__main__':

180

188

181

parser = argparse.ArgumentParser(description='Script to process SOPHy data.')

189

parser = argparse.ArgumentParser(description='Script to process SOPHy data.')

182

parser.add_argument('--parameters', nargs='*', default=['S'],

190

parser.add_argument('--parameters', nargs='*', default=['S'],

183

help='Variables to process: P, Z, V ,W')

191

help='Variables to process: P, Z, V ,W')

184

parser.add_argument('--grado', default=2,

192

parser.add_argument('--grado', default=2,

185

help='Angle in Elev to plot')

193

help='Angle in Elev to plot')

186

parser.add_argument('--save', default=0,

194

parser.add_argument('--save', default=0,

187

help='Save plot')

195

help='Save plot')

188

parser.add_argument('--range', default=0, type=float,

196

parser.add_argument('--range', default=0, type=float,

189

help='Max range to plot')

197

help='Max range to plot')

190

args = parser.parse_args()

198

args = parser.parse_args()

191

199

192

main(args)

200

main(args)

193

201

194

#python sophy_proc_rev006.py --parameters Z --grado 8 --save 1 --range 28

202

#python sophy_proc_rev006.py --parameters Z --grado 8 --save 1 --range 28

195

'''

203

'''

196

def read_and_overview(filename):

204

def read_and_overview(filename):

197

"""Read HDF5 using read_generic_hdf5 and print upper level dictionary keys

205

"""Read HDF5 using read_generic_hdf5 and print upper level dictionary keys

198

"""

206

"""

199

test = read_generic_hdf5(filename)

207

test = read_generic_hdf5(filename)

200

print("\nPrint keys for file %s" % os.path.basename(filename))

208

print("\nPrint keys for file %s" % os.path.basename(filename))

201

for key in test.keys():

209

for key in test.keys():

202

print("\t%s" % key)

210

print("\t%s" % key)

203

return test

211

return test

204

212

205

file__ = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0_FD_PL_R15.0km/2022-06-09T18-00-00/SOPHY_20220609_180229_E2.0_Z.hdf5"

213

file__ = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0_FD_PL_R15.0km/2022-06-09T18-00-00/SOPHY_20220609_180229_E2.0_Z.hdf5"

206

214

207

215

208

filename = get_wradlib_data_file(file__)

216

filename = get_wradlib_data_file(file__)

209

217

210

print("filename:\n",filename )

218

print("filename:\n",filename )

211

219

212

test= read_and_overview(filename)

220

test= read_and_overview(filename)

213

# ANADIR INFORMACION

221

# ANADIR INFORMACION

214

# informacion de los pulsos de TX

222

# informacion de los pulsos de TX

215

# informacion de los ruidos

223

# informacion de los ruidos

216

# informacion de los SNR ¿?

224

# informacion de los SNR ¿?

217

# Aumentar la amplitud de la USRP

225

# Aumentar la amplitud de la USRP

218

LAST_UPDATE

226

LAST_UPDATE

219

---- Noise

227

---- Noise

220

---- Mapas

228

---- Mapas

221

'''

229

'''

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

             #!/usr/bin/env python
             # Ing. AVP
             # 22/06/2022
             # ARCHIVO DE LECTURA y PLOT
             import matplotlib.pyplot as pl
             import matplotlib
             import wradlib
             import numpy
             import warnings
             import argparse
             from wradlib.io import read_generic_hdf5
             from wradlib.util import get_wradlib_data_file
             from plotting_codes import sophy_cb_tables
             import os,time
             for name, cb_table in sophy_cb_tables:
                 ncmap = matplotlib.colors.ListedColormap(cb_table, name=name)
                 matplotlib.pyplot.register_cmap(cmap=ncmap)
             #LINUX bash: export WRADLIB_DATA=/path/to/wradlib-data
+            #example
+            #export WRADLIB_DATA="/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T18-00-00/"
             warnings.filterwarnings('ignore')
             PARAM = {
                 'S': {'var': 'power','vmin': -45, 'vmax': -15, 'cmap': 'jet', 'label': 'Power','unit': 'dBm'},
                 'V': {'var': 'velocity', 'vmin': -10, 'vmax': 10 , 'cmap': 'sophy_v', 'label': 'Velocity','unit': 'm/s'},
                 'Z': {'var': 'reflectivity','vmin': -30, 'vmax': 80 , 'cmap': 'sophy_r','label': 'Reflectivity','unit': 'dBZ'},
                 'W': {'var': 'spectral_width',    'vmin': 0  , 'vmax': 12 , 'cmap': 'sophy_w','label': 'Spectral Width','unit': 'hz'}
                 }
             class Readsophy():
                 def __init__(self):
                     self.list_file = None
                     self.grado     = None
                     self.variable  = None
                     self.save      = None
                     self.range     = None
                 def read_files(self,path_file,grado=None, variable=None):
                     filter= "_E"+str(grado)+".0_"+variable
                     validFilelist = []
                     fileList= os.listdir(path_file)
                     for thisFile in fileList:
                         if (os.path.splitext(thisFile)[0][-7:] != filter):
                             #print("s_:",os.path.splitext(thisFile)[0][-7:])
                             continue
                         validFilelist.append(thisFile)
                         validFilelist.sort()
                     return validFilelist
                 def setup(self, path_file,grado,range,variable,save):
                     self.path_file = path_file
                     self.range     = range
                     self.grado     = grado
                     self.variable  = variable
                     self.save      = save
                     self.list_file = self.read_files(path_file=self.path_file,grado=self.grado, variable=self.variable)
                 def selectHeights(self,heightList,minHei,maxHei):
                     if minHei and maxHei:
                         if (minHei < heightList[0]):
                             minHei = heightList[0]
                         if (maxHei > heightList[-1]):
                             maxHei = heightList[-1]
                         minIndex = 0
                         maxIndex = 0
                         heights = heightList
                         inda = numpy.where(heights >= minHei)
                         indb = numpy.where(heights <= maxHei)
                         try:
                             minIndex = inda[0][0]
                         except:
                             minIndex = 0
                         try:
                             maxIndex = indb[0][-1]
                         except:
                             maxIndex = len(heights)
                         new_heightList= self.selectHeightsByIndex(heightList=heightList,minIndex=minIndex, maxIndex=maxIndex)
                     return new_heightList, minIndex,maxIndex
                 def selectHeightsByIndex(self,heightList,minIndex, maxIndex):
                     if (minIndex < 0) or (minIndex > maxIndex):
                         raise ValueError("Height index range (%d,%d) is not valid" % (minIndex, maxIndex))
                     if (maxIndex >= len(heightList)):
                         maxIndex = len(heightList)
                     new_h = heightList[minIndex:maxIndex]
                     return new_h
+                def readAttributes(self,obj,variable):
+                    var    = PARAM[variable]['var']
+                    unit   = PARAM[variable]['unit']
+                    cmap   = PARAM[variable]['cmap']
+                    vmin   = PARAM[variable]['vmin']
+                    vmax   = PARAM[variable]['vmax']
+                    label  = PARAM[variable]['label']
+                    var_    = 'Data/'+var+'/H'
+                    data_arr   = numpy.array(obj[var_]['data']) # data
+                    utc_time   = numpy.array(obj['Data/time']['data'])
+                    data_azi   = numpy.array(obj['Metadata/azimuth']['data']) # th
+                    data_ele   = numpy.array(obj["Metadata/elevation"]['data'])
+                    heightList = numpy.array(obj["Metadata/range"]['data']) # r
+                    return data_arr, utc_time, data_azi,data_ele, heightList,unit,cmap,vmin,vmax,label
                 def run(self):
                     count= 0
                     len_files = len(self.list_file)
                     for thisFile in self.list_file:
                         count = count +1
                         fullpathfile = self.path_file + thisFile
                         filename   = get_wradlib_data_file(fullpathfile)
                         test_hdf5  = read_generic_hdf5(filename)
-                        var    = PARAM[self.variable]['var']
+                        # LECTURA
-                        unit   = PARAM[self.variable]['unit']
+                        data_arr, utc_time, data_azi,data_ele, heightList,unit,cmap,vmin,vmax,label = self.readAttributes(obj= test_hdf5,variable=self.variable)
-                        cmap   = PARAM[self.variable]['cmap']
-                        vmin   = PARAM[self.variable]['vmin']
-                        vmax   = PARAM[self.variable]['vmax']
-                        label  = PARAM[self.variable]['label']
-                        var_    = 'Data/'+var+'/H'
-                        data_arr   = numpy.array(test_hdf5[var_]['data']) # data
-                        utc_time   = numpy.array(test_hdf5['Data/time']['data'])
-                        data_azi   = numpy.array(test_hdf5['Metadata/azimuth']['data']) # th
-                        data_ele   = numpy.array(test_hdf5["Metadata/elevation"]['data'])
-                        heightList = numpy.array(test_hdf5["Metadata/range"]['data']) # r
                         if self.range==0:
                             self.range == heightList[-1]
                         new_heightList,minIndex,maxIndex = self.selectHeights(heightList,0.06,self.range)
+                        # TIEMPO
                         my_time    = time.strftime('%Y-%m-%d %H:%M:%S', time.localtime(utc_time[0]))
                         time_save  = time.strftime('%Y%m%d_%H%M%S',time.localtime(utc_time[0]))
                         # PLOT DATA WITH ANNOTATION
                         if count ==1:
                             fig = pl.figure(figsize=(10,8))
                             cgax, pm = wradlib.vis.plot_ppi(data_arr[:,minIndex:maxIndex],r=new_heightList,az=data_azi,rf=1,fig=fig, ax=111,proj='cg',cmap=cmap,vmin=vmin, vmax=vmax)
                             caax  = cgax.parasites[0]
                             title = 'Simple PPI'+"-"+ my_time+"  E."+self.grado
                             t = pl.title(title, fontsize=12,y=1.05)
                             cbar = pl.gcf().colorbar(pm, pad=0.075)
                             pl.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')
                             cbar.set_label(label+'[' + unit + ']')
                             gh = cgax.get_grid_helper()
                         else:
                             cgax, pm = wradlib.vis.plot_ppi(data_arr[:,minIndex:maxIndex],r=new_heightList,az=data_azi,rf=1,fig=fig, ax=111,proj='cg',cmap=cmap,vmin=vmin, vmax=vmax)
                             caax  = cgax.parasites[0]
                             title = 'Simple PPI'+"-"+my_time+"  E."+self.grado
                             t = pl.title(title, fontsize=12,y=1.05)
                             cbar = pl.gcf().colorbar(pm, pad=0.075)
                             pl.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')
                             cbar.set_label(label+'[' + unit + ']')
                             gh = cgax.get_grid_helper()
                         if self.save == 1:
                             if count ==1:
                                 filename     = "SOPHY"+"_"+time_save+"_"+"E."+self.grado+"_"+self.variable+".png"
                                 dir =self.variable+"_"+"E."+self.grado+"CH0/"
                                 filesavepath = os.path.join(self.path_file,dir)
                                 try:
                                     os.mkdir(filesavepath)
                                 except:
                                     pass
                             else:
                                 filename     = "SOPHY"+"_"+time_save+"_"+"E."+self.grado+"_"+self.variable+".png"
                             pl.savefig(filesavepath+filename)
                         pl.pause(1)
                         pl.clf()
                     if count==len_files:
                         pl.close()
                     pl.show()
             PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T18-00-00/"
             #PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0N36.0/2022-06-09T19-00-00/"
             #PATH = "/home/soporte/Documents/EVENTO/HYO_PM@2022-05-31T12-00-17/paramC0N36.0/2022-05-31T16-00-00/"
             def main(args):
                 grado      = args.grado
                 parameters = args.parameters
                 save       = args.save
                 range      = args.range
                 obj        = Readsophy()
                 for param in parameters:
                     obj.setup(path_file = PATH,grado = grado,range=range, variable=param,save=int(save))
                     obj.run()
             if __name__ == '__main__':
                 parser = argparse.ArgumentParser(description='Script to process SOPHy data.')
                 parser.add_argument('--parameters', nargs='*', default=['S'],
                                     help='Variables to process: P, Z, V ,W')
                 parser.add_argument('--grado', default=2,
                                     help='Angle in Elev to plot')
                 parser.add_argument('--save', default=0,
                                     help='Save plot')
                 parser.add_argument('--range', default=0, type=float,
                                     help='Max range to plot')
                 args = parser.parse_args()
                 main(args)
             #python sophy_proc_rev006.py  --parameters Z --grado 8 --save 1 --range 28
             '''
             def read_and_overview(filename):
                 """Read HDF5 using read_generic_hdf5 and print upper level dictionary keys
                 """
                 test = read_generic_hdf5(filename)
                 print("\nPrint keys for file %s" % os.path.basename(filename))
                 for key in test.keys():
                     print("\t%s" % key)
                 return test
             file__ = "/home/soporte/Documents/EVENTO/HYO_PM@2022-06-09T15-05-12/paramC0_FD_PL_R15.0km/2022-06-09T18-00-00/SOPHY_20220609_180229_E2.0_Z.hdf5"
             filename = get_wradlib_data_file(file__)
             print("filename:\n",filename )
             test= read_and_overview(filename)
             # ANADIR INFORMACION
             # informacion de los pulsos de TX
             # informacion de los ruidos
             # informacion de los SNR ¿?
             # Aumentar la amplitud de la USRP
             LAST_UPDATE
             ---- Noise
             ---- Mapas
             '''