Niet alleen schaduwpunten moeten berekend worden, maar ook de stijlverheffing en de hoek voor de substijl.
Ook hiervoor is de hoofdprocedure te gebruiken.
Om deze constanten van de zonnewijzer te berekenen wordt als invoer voor de 5 routines de waarden
decl = 90° and t = 0° gebruikt.
Pas nu echter geen conrole toe op de 2 beslissingspunten in de procedure.
De stijlverheffing v is dan: v
= arcsin(z3). (Testen op kwadrant is niet nodig)
Als v = 0 dan ligt de poolstijl evenwijdig aan het zonnewijzervlak en is er geen doordringingspunt van de stijl met het vlak.
Anders: de coördinaten van het doordringingspunt zijn x,y.
De substijl kan getekend worden vanaf het voetpunt van de gnomon ( point 0, 0 ) naar dit doordringingspunt.
De waarde van v is ook gelijk aan de breedtegraad waar de zonnewijzer een horizontale zonnewijzer wordt.
Merk ook op dat als v > 0 het eindpunt van de poolstijl naar de noordpool wijst en dat de uurlijnen met de klok mee lopen, maar als v < 0 wijst de stijl naar de zuidpool en en lopen de uurlijnen tegen de klok in.
Hoewel het niet noodzakelijk is om de zonnewijzer te construeren wordt ook een formule gegven om de uurhoek van de substijl ts te kunnnen berekenen.
tan ts = sin i * sin d / (cos phi * cos i + sin phi * sin i * cos d)
Controleer ts voor ligging in het juiste kwadrant.
Fomules om de tijdvereffening en de declinatie van de zon te bereken uit een gegeven dagnummer:
L = DAGNR*360/365.2422 - 80.535132 GRADEN
TIJDVEREFFENIG = - 107.0605*SIN( L) - 428.6697*COS( L)
+ 596.1009*SIN(2L) - 2.0898*COS(2L)
+ 4.4173*SIN(3L) + 19.2776*COS(3L)
- 12.7338*SIN(4L) SECONDEN (in tijd)
LAMBDA = L + 0.4277*SIN( L) + 1.8664*COS( L)
- 0.0180*SIN(2L) + 0.0087*COS(2L) GRADEN
EPSILON = 23.43954 GRADEN
DECLINATIE = ARCSIN(SIN LAMBDA * SIN EPSILON) GRADEN
Strikt genomen gelden deze formules voor 1998 en voor 12.00 UT, maar voor zonnewijzers kunt u ze gedurende
langere tijd gerust toepassen.
See ook 'Cousins', pagina 236, waar zo'n formule voor 1931 is vermeld.
Als u zelf deze formules voor een bepaald jaar in de periode van 1901 - 2099 wil berekenen kunt u het programma
equadecl.zip binnen halen. 26kB
Of download een
tabel
met de gemiddelde waarden over de periode 2000-2099.
12kB
