##// END OF EJS Templates
schain
RSS
Drifts proccesing for some heights of interest. Some parameters are needded for madrigal writting: lat, lon. NAN vectors are not written to avoid errors at madrigal view.
Drifts proccesing for some heights of interest. Some parameters are needded for madrigal writting: lat, lon. NAN vectors are not written to avoid errors at madrigal view.
jespinoza - - Load All Authors

File last commit:

r1396:f39ad5b721a3
r1690:571f81d70cd0
Show More
pedestal_client.py
87 lines | 2.8 KiB | text/x-python | PythonLexer
/ schainpy / scripts / pedestal_client.py
History | Source | Raw |Copy content |Copy permalink
avaldez
simulador pedestal client
 r1320 import numpy
import sys
import zmq
import time
import h5py
import os
jespinoza
Merge EW-Drifts
 r1396 path = "/home/alex/Downloads/pedestal"
ext = ".hdf5"
avaldez
simulador pedestal client
 r1320
jespinoza
Merge EW-Drifts
 r1396 port = "5556"
if len(sys.argv) > 1:
avaldez
simulador pedestal client
 r1320 port = sys.argv[1]
int(port)
jespinoza
Merge EW-Drifts
 r1396 if len(sys.argv) > 2:
avaldez
simulador pedestal client
 r1320 port1 = sys.argv[2]
int(port1)
jespinoza
Merge EW-Drifts
 r1396 # Socket to talk to server
avaldez
simulador pedestal client
 r1320 context = zmq.Context()
jespinoza
Merge EW-Drifts
 r1396 socket = context.socket(zmq.SUB)
avaldez
simulador pedestal client
 r1320
print("Collecting updates from weather server...")
jespinoza
Merge EW-Drifts
 r1396 socket.connect("tcp://localhost:%s" % port)
avaldez
simulador pedestal client
 r1320
jespinoza
Merge EW-Drifts
 r1396 if len(sys.argv) > 2:
socket.connect("tcp://localhost:%s" % port1)
avaldez
simulador pedestal client
 r1320
jespinoza
Merge EW-Drifts
 r1396 # Subscribe to zipcode, default is NYC,10001
avaldez
simulador pedestal client
 r1320 topicfilter = "10001"
jespinoza
Merge EW-Drifts
 r1396 socket.setsockopt_string(zmq.SUBSCRIBE, topicfilter)
# Process 5 updates
total_value = 0
count = -1
azi = []
elev = []
time0 = []
# for update_nbr in range(250):
avaldez
simulador pedestal client
 r1320 while(True):
jespinoza
Merge EW-Drifts
 r1396 string = socket.recv()
topic, ang_elev, ang_elev_dec, ang_azi, ang_azi_dec, seconds, seconds_dec = string.split()
ang_azi = float(ang_azi) + 1e-3 * float(ang_azi_dec)
ang_elev = float(ang_elev) + 1e-3 * float(ang_elev_dec)
seconds = float(seconds) + 1e-6 * float(seconds_dec)
avaldez
simulador pedestal client
 r1320 azi.append(ang_azi)
elev.append(ang_elev)
time0.append(seconds)
jespinoza
Merge EW-Drifts
 r1396 count += 1
avaldez
simulador pedestal client
 r1320 if count == 100:
jespinoza
Merge EW-Drifts
 r1396 timetuple = time.localtime()
avaldez
simulador pedestal client
 r1320 epoc = time.mktime(timetuple)
jespinoza
Merge EW-Drifts
 r1396 # print(epoc)
fullpath = path + ("/" if path[-1] != "/" else "")
avaldez
simulador pedestal client
 r1320
if not os.path.exists(fullpath):
os.mkdir(fullpath)
jespinoza
Merge EW-Drifts
 r1396 azi_array = numpy.array(azi)
avaldez
simulador pedestal client
 r1320 elev_array = numpy.array(elev)
jespinoza
Merge EW-Drifts
 r1396 time0_array = numpy.array(time0)
pedestal_array = numpy.array([azi, elev, time0])
count = 0
azi = []
elev = []
time0 = []
# print(pedestal_array[0])
# print(pedestal_array[1])
avaldez
simulador pedestal client
 r1320
jespinoza
Merge EW-Drifts
 r1396 meta = 'PE'
filex = "%s%4.4d%3.3d%10.4d%s" % (meta, timetuple.tm_year, timetuple.tm_yday, epoc, ext)
filename = os.path.join(fullpath, filex)
fp = h5py.File(filename, 'w')
# print("Escribiendo HDF5...",epoc)
# ·················· Data·....······································
avaldez
simulador pedestal client
 r1320 grp = fp.create_group("Data")
dset = grp.create_dataset("azimuth" , data=pedestal_array[0])
dset = grp.create_dataset("elevacion", data=pedestal_array[1])
dset = grp.create_dataset("utc" , data=pedestal_array[2])
jespinoza
Merge EW-Drifts
 r1396 # ·················· Metadata·······································
avaldez
simulador pedestal client
 r1320 grp = fp.create_group("Metadata")
dset = grp.create_dataset("utctimeInit", data=pedestal_array[2][0])
jespinoza
Merge EW-Drifts
 r1396 timeInterval = pedestal_array[2][1] - pedestal_array[2][0]
avaldez
simulador pedestal client
 r1320 dset = grp.create_dataset("timeInterval", data=timeInterval)
fp.close()
jespinoza
Merge EW-Drifts
 r1396 # print ("Average messagedata value for topic '%s' was %dF" % ( topicfilter,total_value / update_nbr))