;
; ECKERTSCHER ENTWURF IV
; ======================
;
; Nr. bei Wagner: 23b
; Kurzname:       Eckert IV
; Name:           Eckertscher Entwurf IV (flächentreu, elliptische Meridiane)
; Autor:          Max Eckert(-Greiffendorf)
; Quelle:         Karlheinz Wagner, Kartographische Netzentwürfe, Leipzig 1949, S. 220f.
; Richtung:       Invers
;
; 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 Mittelmeridian werden abgefragt.
;
; 12062005 RB Janos Gyoerffy hat mich auf einen Fehler aufmerksam gemacht (4pi mus an einer Stelle 4+pi sein)
;
; Literatur:
; Wagner: Kartographische Netzentwürfe, Leipzig: Bibliographisches Institut 1949

; (C) Rolf Böhm 2004

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

; Die Variablennamen entsprechen weitgehend denen von Karlheinz Wagner,

;
; Laufende Koordinaten
;
	_name	Eckerts~Entwurf~IV
	_var	phi		; Geographische Breite
	_var	lambda		; Geographische Länge
	_var	psi		; The parameter
	_var	t1		; temp
	_var	t2		; temp
	_var	t3		; temp
	_var	n		; 2/root(2+pi)
	_var	n*pi/2		; n*pi/2
	_var	4pi		; 4*pi
	_var	4+pi		; 4+pi   Köszenem!
	_var	lambda0		; Null-Länge
;
; x, y, x', y',  Cx', Cy', Rx', Ry', °(, (°, pi, pi/2 etc. sind vordefinierte globale Konstanten
;
; Initialisierung
; ===============
;
	tstne	initial	077$
; Dialog
	input	scale	Eckerts~Entwurf~IV~(flächentreu,~elliptische~Meridiane)\\Maßstabszahl
	clip	scale	1	1E12
	input	lambda0	Mittelpunktslänge~in~Grad
	clip	lambda0	-180	180
; Konstanten (Konfigurationsvariablen) berechnen
	mov 	t1	4
	mov	t2	pi
	add	t2	4
	mul	t2	pi
	root	t2	2
	div	t1	t2
	mov	n	t1
;	
	mov	n*pi/2	n
	mul	n*pi/2	pi
	div	n*pi/2	2
;
	mov	4pi	pi
	mul	4pi	4
;
	mov	4+pi	pi
	add	4+pi	4
; 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	psi	y		; psi -- the parameter
	div	psi	n*pi/2
	asin	psi
	errjump	out

	mov	t1	psi		; t1 = 1 + cos(psi) 
	cos	t1
	add	t1	1		
	
	mov	lambda	x		; lambda
	mul	lambda	2
	div	lambda	n
	div	lambda	t1

	clr	phi			; phi <-- psi
	mov	t2	psi
	mul	t2	2
	add	phi	t2
	mov	t2	psi
	sin	t2
	mul	t2	4
	add	phi	t2
	mov	t2	psi
	mul	t2	2
	sin	t2
	add	phi	t2		; the enumerator

	div	phi	4+pi		; divide through 4+pi

	asin	phi
;
; In Gradmaß umrechnen und Ausserhalbtest
; ---------------------------------------

	mul	phi	(°
	mul	lambda	(°
	mul	psi	(°
	cmplt	phi	-90	out
	cmpgt	phi	90	out
	cmplt	lambda	-180	out
	cmpgt	lambda	180	out
;
; Schlussarbeiten
; ---------------
;
	mov	x'	lambda
	mov	y'	phi
	mov	z'	psi		; test
; Lambda um lambda0 kreisen lassen
	add	x'	lambda0
	cmod	x'	-180	180
	exit
out:
	mov	x'	-9999
	mov	y'	-9999
	exit
	_end