Advertisement
Guest User

Untitled

a guest
Jan 19th, 2018
94
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.16 KB | None | 0 0
  1. from math import sin, cos, sqrt, atan2, radians
  2. with open('data.txt') as f:
  3. lines = f.readlines()
  4. lines = lines[3:]
  5. line= lines[0]
  6. varis = line.split()
  7. lat = float(varis[1])
  8. lon = float(varis[2])
  9. wind = int(varis[4])
  10. pres = int(varis[5])
  11. if(varis[6][0]=="H"):
  12. level = varis[6][-1]
  13. elif(varis[6][-1]=="N"):
  14. level = .5
  15. else:
  16. level = .7
  17. lcount =0
  18. lastpoint=[lat,lon]
  19. dist=0
  20. cat=0
  21. lines=lines[1:]
  22. for line in lines: #each line
  23. varis = line.split()
  24. lat = float(varis[1])
  25. lon = float(varis[2])
  26. wind = int(varis[4])
  27. pres = int(varis[5])
  28. if(varis[6][0]=='H'):
  29. level = int(varis[6][-1])
  30. elif(varis[6][-1]=='N'):
  31. level = .5
  32. else:
  33. level = .7
  34. lcount +=1
  35. R = 6373.0
  36. lat1 = radians(lastpoint[0])
  37. lon1 = radians(lastpoint[1])
  38. lat2 = radians(lat)
  39. lon2 = radians(lon)
  40. dlon = lon2 - lon1
  41. dlat = lat2 - lat1
  42. a = sin(dlat / 2)**2 + cos(lat1) * cos(lat2) * sin(dlon / 2)**2
  43. c = 2 * atan2(sqrt(a), sqrt(1 - a))
  44. distance = R * c
  45. dist+=distance
  46. lastpoint=[lat,lon]
  47. cat+=((1020-pres+wind)*level/5*dist/10000)
  48. time=(lcount-1)*6
  49. print cat
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement