Diskussion:Gamma-Wert einer Sonnenfinsternis

Quellen

Es fehlen die Quellen der Werte 1,0260 und 1,5433 und die Begründung des Satzes ... größer als die Summe aus 1,5433 und dem Kernschattenradius, .... --Analemma 22:28, 10. Nov. 2008 (CET)

Was soll das? Man kann durchaus erkennen, dass das durch die erste Quelle abgedeckt ist.?--Cactus26 07:04, 11. Nov. 2008 (CET)

Ersetze Quelle durch Begründung. 0,9972 ist begründet (erklärt), die anderen Werte (und der zitierte Satz) nicht. --Analemma 23:58, 12. Nov. 2008 (CET)

Stellungnahme von Sch

Sch, der Hauptautor dieses Artikels ist, hat mir per Mail auf die Quellen-Vorwürfe geantwortet. Er ist leider derzeit zeitlich nicht in der Lage, selbst zu diesem Thema Stellung zu nehmen. Ich versuche hier seinen Standpunkt zusammenzufassen:

In dieser Version des Artikels sind alle Zahlenwerte durch die Quelle belegt. In der Quelle selbst werden nicht alle Werte exakt, das heiß inklusive aller Nachkomastellen, hergeleitet. Teilweise handelt es sich dabei möglicherweise um empirisch ermittelte Daten, die sich aus der detaillierten Berechnung verschiedener Finsternisse ergeben. Dies kann aber nicht als Makel angesehen werden, da die Quelle (Meeus) ja allgemein anerkannt ist und es nicht Aufgabe der Wikipedia sein kann, alle aus Quellen übernommenen Zahlenwerte exakt selbst herzuleiten. Durch die Kürzung verliert der Artikel diese systematische und vollständige Aufstellung verschiedener Finsternistypen, beispielsweise wird nun nirgendwo mehr erläutert, dass es nichtzentrale totale und nichtzentrale ringförmige Finsternisse gibt.

Ich schlage also vor, den Artikel auf diese Version zurückzusetzen.--Cactus26 13:32, 21. Nov. 2008 (CET)

Ich schlage vor, dass Sch sich selbst an der Diskussion beteiligt und seinen Text (sowohl Inhalt als auch Darstellung) verteidigt. Dass es nichtzentrale totale und nichtzentrale ringförmige Finsternisse gibt, ist eines von Meeus' Morsels (Leckerbissen), an deren Entdeckung er selbst die meiste Freude hat (s. Vorwort im Band III). Die Leser unserer Enzyklopädie müssen das nicht auch noch von uns erfahren. Wer hält sich schon in diesen seltensten Fällen ausgerechnet in einem polaren Gebiet der Erde auf?
Analemma 15:01, 21. Nov. 2008 (CET)

Kurze Unterbrechung

Da Sch, wie per Mail geschrieben, nicht zeitnah anworten kann und der Artikel unter Wegfall der bisherigen Einzelnachweise gekürzt wurde, habe ich zunächst die alte Version wiederhergestellt und bitte deshalb bis zur Stellungnahme von Sch. um Geduld.--NebMaatRe 16:45, 21. Nov. 2008 (CET)

• Solch eine Behinderung hatten wir schon einmal (Meton-Zyklus), als der sich rar machende Experte Sch erst seinen Segen geben musste (er hat ihn bis heute nicht gegeben).
• unter Wegfall der bisherigen Einzelnachweise gekürzt ist in sich selbst begründet: Reduktion auf das eindeutig Grundsätzliche ohne Beilagen (quellenpflichtige Morsels) von Abweichungen <<1 (sehr viel kleiner als Eins, gegen Null gehend). Das ist auch keine einleuchtende Begründung, zum Alten zurück kehren zu müssen.

Analemma 17:30, 21. Nov. 2008 (CET)

Immerhin antwortet Sch per Mail. Cactus26 steht in direktem Kontakt. Da wird eine Antwort nicht lange auf sich warten lassen.--NebMaatRe 17:43, 21. Nov. 2008 (CET)

Programm zur Berechnung von Gamma

Habe eines in C geschrieben:

#include <math.h>
#include <stdio.h>

// number "PI" exactly = atan(1.0)*4.0

const double PIE = 3.1415926535897932384626433832795;

// prevents for crashes if sqrt of a negative number is calculated

double sqrt_sec(double s)

{
  if (s<1e-35) s=0.0;
  return(sqrt(s));
}

// sinus from degree

double sin_deg(double x)

{
  return(sin(x * PIE/180.0));
}

// cosinus from degree

double cos_deg(double x)

{
  return(cos(x * PIE/180.0));
}

/*
 * rasun   = right ascension of the sun in hour
 * decsun  = declination of the sun in degree
 * ausun   = distance of the sun in astron. units
 *
 * ramoon  = right ascension of the moon in hour
 * decmoon = declination of the moon in degree
 * aumoon  = distance of the moon in astron. units
 *
 * result  = gamma (needed for eclipses)
 *
 */

double gammaeclipse(double rasun, double decsun, double ausun, double ramoon, double decmoon, double aumoon)

{
  double sQValue,msQValue,smProduct,help,gamma,zplane;
  double xsun,ysun,zsun,xmoon,ymoon,zmoon,dmsx,dmsy,dmsz;

  // help=ASTROUNIT*1000000.0/EARTHRAD;
  // help=149597870.691/6378.137;
  help=23454.7910606184887;

  ausun  *= help;
  aumoon *= help;

  rasun *= 15.0;
  help = cos_deg(decsun)*ausun;
  xsun = sin_deg(rasun)*help;
  ysun = cos_deg(rasun)*help;
  zsun = sin_deg(decsun)*ausun;

  ramoon *= 15.0;
  help  = cos_deg(decmoon)*aumoon;
  xmoon = sin_deg(ramoon)*help;
  ymoon = cos_deg(ramoon)*help;
  zmoon = sin_deg(decmoon)*aumoon;

  dmsx = xmoon - xsun;
  dmsy = ymoon - ysun;
  dmsz = zmoon - zsun;

  sQValue   = xsun*xsun + ysun*ysun + zsun*zsun;
  msQValue  = dmsx*dmsx + dmsy*dmsy + dmsz*dmsz;
  smProduct = xsun*dmsx + ysun*dmsy + zsun*dmsz;

  // prevents for divisions through zero
  if (msQValue<1e-20) return(0.0);

  // calculate gamma
  gamma = sqrt_sec(sQValue - (smProduct*smProduct)/msQValue);

  // first version for zplane
  // t = xsun/(xsun - xmoon);
  // zplane = zsun + (zmoon - zsun)*t;

  // second version for zplane
  // t = -xsun/dmsx;
  // zplane = zsun + dmsz*t;

  // shortest version for zplane
  zplane = zsun - dmsz*xsun/dmsx;

  printf("\n Z-Plane: %.7f ", zplane);

  // zplane negative -> gamma also
  return ((zplane<0.0) ? -gamma : gamma);
}

int main(void)

{
  double gamma,gammaex;
  double rasun,decsun,ausun;
  double ramoon,decmoon,aumoon;

  printf("\n Total Eclipse of the Sun August 11th, 1999, 11:03:08 UTC \n");

  rasun   =  9.385640074365;
  decsun  = 15.3277038455988652;
  ausun   =  1.0135790889737213;

  ramoon  =  9.392955931008;
  decmoon = 15.8105919224704010;
  aumoon  =  0.0024954509874679;

  // http://en.wikipedia.org/wiki/Solar_eclipse_of_August_11,_1999
  // gamma 0.5062

  gammaex = 0.5062;

  gamma = gammaeclipse(rasun, decsun, ausun, ramoon, decmoon, aumoon);

  printf("\n Gamma:   %.7f ", gamma);
  printf("\n (Info Gamma Wikipedia: %.4f) \n",gammaex);

  printf("\n Total Eclipse of the Sun May 29th, 1919, 13:08:55 UTC \n");

  rasun   =  4.352035473895;
  decsun  = 21.5044408467224315;
  ausun   =  1.0137195673879480;

  ramoon  =  4.353762507108;
  decmoon = 21.2054246082255915;
  aumoon  =  0.0024004724728975;

  // http://en.wikipedia.org/wiki/Solar_eclipse_of_May_29,_1919
  // gamma -0.2955

  gammaex = -0.2955;

  gamma = gammaeclipse(rasun, decsun, ausun, ramoon, decmoon, aumoon);

  printf("\n Gamma:   %.7f ", gamma);
  printf("\n (Info Gamma Wikipedia: %.4f) \n",gammaex);

  return(0);
}

und hier das Ergebnis:

 Total Eclipse of the Sun August 11th, 1999, 11:03:08 UTC 

 Z-Plane: 0.5490825 
 Gamma:   0.5062142 
 (Info Gamma Wikipedia: 0.5062) 

 Total Eclipse of the Sun May 29th, 1919, 13:08:55 UTC 

 Z-Plane: -0.3208799 
 Gamma:   -0.2954862 
 (Info Gamma Wikipedia: -0.2955) 

-- Fraba 09:17, 8. Dez. 2010 (CET)

Danke, Fraba. Werde das mal ausprobieren. Viele Grüße --Cactus26 11:49, 8. Dez. 2010 (CET)

Zum Glück hat mir jemand in der alt.astronomy Newsgroup geholfen:

If you can cast your Sun and Moon positions into Earth-centered rectangular coordinates then you can apply the point-line distance formula to find the distance from the center of the Earth to the Sun-Moon central line. Even better if your coordinate units are Earth equitorial radii; Then the gamma ratio is the calculated distance.

3D Point-Line distance:

http://mathworld.wolfram.com/Point-LineDistance3-Dimensional.html

Of course, accuracy may be affected by atmospheric refraction, which will "bend" the apparent Sun-Moon line's direction. If you correct the apparent positions of the Sun and Moon for refraction ahead of time then you need only plug their adjusted coordinates into the distance formula.

-- Fraba 20:34, 9. Dez. 2010 (CET)