Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from math import sin, cos, sqrt, atan2, radians
- with open('data.txt') as f:
- lines = f.readlines()
- lines = lines[3:]
- line= lines[0]
- varis = line.split()
- lat = float(varis[1])
- lon = float(varis[2])
- wind = int(varis[4])
- pres = int(varis[5])
- if(varis[6][0]=="H"):
- level = varis[6][-1]
- elif(varis[6][-1]=="N"):
- level = .5
- else:
- level = .7
- lcount =0
- lastpoint=[lat,lon]
- dist=0
- cat=0
- lines=lines[1:]
- for line in lines: #each line
- varis = line.split()
- lat = float(varis[1])
- lon = float(varis[2])
- wind = int(varis[4])
- pres = int(varis[5])
- if(varis[6][0]=='H'):
- level = int(varis[6][-1])
- elif(varis[6][-1]=='N'):
- level = .5
- else:
- level = .7
- lcount +=1
- R = 6373.0
- lat1 = radians(lastpoint[0])
- lon1 = radians(lastpoint[1])
- lat2 = radians(lat)
- lon2 = radians(lon)
- dlon = lon2 - lon1
- dlat = lat2 - lat1
- a = sin(dlat / 2)**2 + cos(lat1) * cos(lat2) * sin(dlon / 2)**2
- c = 2 * atan2(sqrt(a), sqrt(1 - a))
- distance = R * c
- dist+=distance
- lastpoint=[lat,lon]
- cat+=((1020-pres+wind)*level/5*dist/10000)
- time=(lcount-1)*6
- print cat
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement