;
; WINKEL II
; =========
;
; Name:           Winkel II
; Autor:          Oswald Winkel 1928
; Richtung:       Inverse Transformation
;
; Das Programm übernimmt die Koordinaten eines Punktes (x/y) und transformiert
; diese in einen Punkt (x'/y').
;
; x/y sind ebene Zielpunktkoordinaten, x'/y' geben die geogr. Breite und Länge
; der Position auf der Quell-Erdkugel, auf der der Zielpunkt gelesen werden
; kann.
;
; Maßstab und Berührungspunkt werden abgefragt.

; (C) Rolf Böhm 2004

; Benutzte Variablen
; ==================

;
; Laufende Koordinaten
;
	_name	Winkel~II
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	delta		; Poldistanz/geographisch, auch Schiefdistanz
	_var	t1		; temp
	_var	t2		; temp	
;
; Konstanten der Transformation
;
	_var	k		; 2/pi
	_var	pi^2		; pi²
	_var	lambda0		; Geogr. Länge des Bildmittelpunktes
	_var	scale		; Kartenmaßstabszahl (also 1000000, nicht 1/1000000)
;
; x, y, x', y', Cx', Cy', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten
;
; Initialisierung
; ===============
;
	tstne	initial	077$
; Dialog
	input	scale	Winkel~II~(Schnittparallel~acos(2/pi)~=~50°28')\\Maßstabszahl
	input	lambda0	Mittelpunktslänge~in~Grad
; Eingegebene Werte auf Min/Max bringen
	clip	scale	1	1E12
	clip	lambda0	-180	180
; Konstanten
	mov	pi^2	pi
	mul	pi^2	pi
	mov	k	2
	div	k	pi
; Programm ist initialisiert
	mov	initial	1
077$:

;
; SIMD-Laufbereich
; ================

;
; Maßstab, Kartenmittelpunkt etc. einrechnen
; ------------------------------------------

	sub	x	Cx'		; Bildmittelpunkt
	div	x	Rx'		; Erdradius
	mul	x	scale		; Kartenmaßstab

	sub	y	Cy'
	div	y	Ry'
	mul	y	scale
;
; Eigentlicher Entwurf, dieser invers
; -----------------------------------
;
	mov	phi	y		; Geographische Länge
	
	mov	t2	y
	power	t2	2
	mul	t2	4
	div	t2	pi^2
	mov	t1	1
	sub	t1	t2
	root	t1	2
	add	t1	k
	
	mov	t2	x
	div	t2	t1
	mul	t2	2
	mov	lambda	t2
;
; In Gradmaß umrechnen und Ausserhalbtest
; ---------------------------------------
;
	mul	phi	(°
	mul	lambda	(°
	mul	psi	(°		; Test
	cmplt	phi	-90	out
	cmpgt	phi	90	out
	cmplt	lambda	-180	out
	cmpgt	lambda	180	out

;
; Schlussarbeiten
; ---------------
;
	mov	x'	lambda
	add	x'	lambda0
	mov	y'	phi
	mov	z'	psi		; Test


	cmpgt	x'	-180	10$
	add	x'	360
10$:
	cmpgt	x'	-180	30$
	add	x'	360
30$:
	cmplt	x'	180	40$
	sub	x'	360
40$:
	cmplt	x'	180	50$
	sub	x'	360
50$:

	exit
out:
	mov	x'	-9999
	mov	y'	-9999
	exit
	_end