#!/usr/bin/python3
# coding=utf8

import sys, math
from kml import *

fn = "loxodroma"

long = 0				# výchozí geodetická délka - stupně
longStep = 5
maxSteps = 9999
lat = 87				# výchozí geodetická šířka - stupně
minPoleDiff = 1			# minimální vzdálenost od pólu - stupně
R = 6378 * 1000			# zemský poloměr v metrech
azim = 105
tilt = 60
range = 4800 * 1000		# metry
speed = 1000 * 1000		# m/s		rychlost přesunu pohledu
durMin = 0.1			# secs		minimální doba kroku
durMax = 1				# secs		maximální doba kroku

def loxClong():
	azimR = math.radians(azim)
	latR = math.radians(lat)
	longStepR = math.radians(longStep)
	latNR = math.asin(math.tanh(longStepR / math.tan(azimR) + math.atanh(math.sin(latR))))
	latN = math.degrees(latNR)
	longN = canonGeodetic(long + longStep)
	move = R * math.fabs(latNR - latR) / math.cos(azimR)	# délka posunu na povrchu koule v metrech
	dur = math.fabs(move) / speed		# doba přeletu
	if(dur < durMin): dur = durMin
	elif(dur > durMax): dur = durMax
	return (longN, latN, dur, move)

def list():
	global long, lat
	k = 0
	while((90 - math.fabs(lat)) >= minPoleDiff and k < maxSteps):
		(longN, latN, dur, move) = loxClong()
		print('{}\t{}\t{}\t{}'.format(k,longN,latN, move))
		long = longN
		lat = latN
		k = k + 1

def run():
	global long, lat
	k = 0

	FlyTo.head(5, False)
	LookAt.run(long, lat, 0, azim, tilt, range)
	FlyTo.tail()

	while((90 - math.fabs(lat)) >= minPoleDiff and k < maxSteps):
		(longN, latN, dur, move) = loxClong()
		FlyTo.head(dur, True)
		comm('k={}, long={}, lat={}, longN={}, latN={}, move={}'.format(str(k), str(long), str(lat), str(longN), str(latN), str(move)))
		k = k + 1
		LookAt.run(longN, latN, 0, azim, tilt, range)
		FlyTo.tail()
		long = longN
		lat = latN

# list()
tourProlog(fn)
run()
tourEpilog()