Hi Rainer,
auf meiner
Homepage kannst du dir das Proggi Astronom.exe runterladen.
ansonsten hier der entsp. Quellcode in TurboPascal zum selber schreiben.
PROCEDURE STERNZEIT;
VAR k,s,gmst,lmst,gast,last,geo_laenge,delta_fp,delta_eks,
ma_sonne,ma_mond,ml_mond_mondknoten,ml_mond_sonne,ml_mondknoten,
ekliptikschiefe,ml_sonne,wahre_eks: REAL;
PROCEDURE INSTALLIEREN;
BEGIN
ma_mond:=134.962981+477198.867398*t+0.008697*t*t+1.8e-5*t*t*t;
ma_sonne:=357.527723+35999.05034*t-1.6e-4*t*t-3e-6*t*t*t;
ml_mond_mondknoten:=93.27191+483202.017538*t-0.003683*t*t+3e-6*t*t*t;
ml_mond_sonne:=297.850363+445267.11148*t-0.001914*t*t+5e-6*t*t*t;
ml_mondknoten:=125.044522-1934.136261*t+0.002071*t*t+2e-6*t*t*t;
ml_sonne:=ml_mond_mondknoten-ml_mond_sonne+ml_mondknoten;
ml_mondknoten:=RAD(ml_mondknoten);
ml_sonne:=RAD(ml_sonne);
ma_sonne:=RAD(ma_sonne);
ml_mond_mondknoten:=RAD(ml_mond_mondknoten);
ma_mond:=RAD(ma_mond);
ml_mond_sonne:=RAD(ml_mond_sonne);
ekliptikschiefe:=23.439291-0.0130042*t-1.64e-7*t*t+5.036e-7*t*t*t;
ekliptikschiefe:=RAD(ekliptikschiefe)
END; {INSTALLIEREN}
PROCEDURE AUSGABE;
BEGIN
WRITELN;
WRITELN('Julianisches Datum um 0 Uhr= ',jd:7:6);
WRITE('Mittlere Sternzeit in Greenwich: ');
WriteHMS(gmst,7);
WriteLn;
WRITE('Mittlere Ortssternzeit: ');
WriteHMS(lmst,7);
WriteLn;
WRITE('Wahre Sternzeit in Greewich: ');
WriteHMS(gast,7);
WriteLn;
WRITE('Wahre Ortssternzeit: ');
WriteHMS(last,7);
WriteLn;
WRITE('Mittlere Ekliptikschiefe: ');
WriteDMS(ekliptikschiefe,6);
WRiteLn;
WRITE('Wahre Ekliptikschiefe: ');
WriteDMS(wahre_eks,6);
WriteLn;
WRITELN('L„ngendifferenz wahrer-mittlerer Frhlingspunkt ("): ',
delta_fp:2:3,#13#10,'Differenz zwischen wahrer und mittlerer ',
'Eliptikschiefe ("): ',delta_eks:2:3)
END; {AUSGABE}
BEGIN {STERNZEIT}
CLRSCR;
DATUMEINGABE(2);
geo_laenge:=-361;
WRITELN;
WRITE('Geographische L„nge (”stlich negativ):');
WHILE(geo_laenge < -360) OR (geo_laenge > 360) DO
READREAL(geo_laenge,8,39,7);
WRITELN;
DATUM_JD;
jd:=jd-gmt/24;
t:=(jd-2451545)/36525;
gmst:=0.2790572733+100.0021391*t+1.0776e-6*t*t-7.2e-11*t*t*t
+0.0417807462*gmt;
gmst:=gmst*24;
UHRZEIT24(gmst);
lmst:=gmst*15-geo_laenge;
WINKEL360(lmst);
INSTALLIEREN;
NutationConst(t, delta_fp, delta_eks);
delta_fp:= DEG(delta_fp)*3600;
delta_eks:=DEG(delta_eks)*3600;
k:=delta_fp*COS(ekliptikschiefe);
gast:=gmst+k/54000;
last:=lmst+k/3600;
UHRZEIT24(gast);
WINKEL360(last);
wahre_eks:=ekliptikschiefe+RAD(delta_eks/3600);
gmst:=RAD(gmst*15);
lmst:=RAD(lmst);
gast:=RAD(gast*15);
last:=RAD(last);
AUSGABE;
BITTE_TASTE
END; {STERNZEIT}
Gruß Lots
