#!/bin/env python

import pyart
import glob
from netCDF4 import Dataset
import signal


def plotdata(i):
	radar = pyart.io.read_cfradial(i)
	fh = Dataset(i)
	nyquist = radar.get_nyquist_vel(0)
	dist = round(fh.variables['range'][-1]/1000) # get the last cell from the range data, convert to KM and round it
	nm_dist = round(dist/1.852) # Convert to NM
	vel = max(fh.variables['VEL'][0]) # get the max value from VEL data
	#print(i, ": Max Range:", nm_dist,"NM,", dist, "KM,", "Nyquist Velocity:",nyquist, "M/S", "MaxVelinData", vel, "M/S")
	info = i, ": Max Range:", nm_dist,"NM,", dist, "KM,", "Nyquist Vel:",nyquist, "M/S", "MaxVelData", vel, "M/S"
	try:
		fh = open('data.txt', 'a')
		fh.write(str(info)+'\n')
		print(str(info))
		fh.close()
	except IOError as e:
		print(e)

mylist = [f for f in glob.glob("HWX*.nc")]

try:
	for i in mylist:
		try:
			plotdata(i)
		except KeyboardInterrupt:
			exit()
		except:
			pass
except IndexError as e:
	#print(e)
	exit()