Sind die Intervalle zwischen den einzelnen Vollmonden immer gleich lang? Wenn ja, (was ich bezweifle-Sonne, Mond und Erde eiern sicher durchs All), hat jemand Lust, eine Formel aufzustellen, wann es in den letzten 40 Jahren Vollmond gab? Danke! Wenn nein… kennt jemand einen Jahrhundertkalender, in dem die Vollmonde verzeichnet sind?
keine Antwort, aber Buchtipp
Sind die Intervalle zwischen den
einzelnen Vollmonden immer gleich lang?
NEIN
Wenn ja, (was ich bezweifle-Sonne, Mond
und Erde eiern sicher durchs All), hat
jemand Lust, eine Formel aufzustellen,
wann es in den letzten 40 Jahren Vollmond
gab?
Nein, weil schon geschehen und mir zu mühsam.
Wenn Du es richtig genau wissen willst, dann hier ein Buchtipp:
Montebruck+Pfleger: Astronomie mit dem Personalcomputer, ISBN 3-540-57701-7 Buch anschauen
Mit dem gesamten astronomisch-physikalischen Background und vielen Programmbeispielen.
Vollmond am Freitag, dem 13.
Hier hast Du eine kleine f77-Routine, die Dir zumindest mal die Vollmonde am Freitag, den 13, zwischen 1900 und 2000 ausrechnet 
PROGRAM BADLUK
DATA IYBEG,IYEND /1900,2000/
TIMZON=-5./24
WRITE (*,’(1X,A,I5,A,I5)’) ‚Full moons on Friday the 13th from‘,
* IYBEG,’ to’,IYEND
DO 12 IYYY=IYBEG,IYEND
DO 11 IM=1,12
JDAY=JULDAY(IM,13,IYYY)
IDWK=MOD(JDAY+1,7)
IF(IDWK.EQ.5) THEN
N=12.37*(IYYY-1900+(IM-0.5)/12.)
ICON=0
1 CALL FLMOON(N,2,JD,FRAC)
IFRAC=NINT(24.*(FRAC+TIMZON))
IF(IFRAC.LT.0)THEN
JD=JD-1
IFRAC=IFRAC+24
ENDIF
IF(IFRAC.GT.12)THEN
JD=JD+1
IFRAC=IFRAC-12
ELSE
IFRAC=IFRAC+12
ENDIF
IF(JD.EQ.JDAY)THEN
WRITE (*,’(/1X,I2,A,I2,A,I4)’) IM,’/’,13,’/’,IYYY
WRITE (*,’(1X,A,I2,A)’) 'Full moon ‚,IFRAC,
* ’ hrs after midnight (EST).‘
GOTO 2
ELSE
IC=ISIGN(1,JDAY-JD)
IF(IC.EQ.-ICON) GOTO 2
ICON=IC
N=N+IC
ENDIF
GOTO 1
2 CONTINUE
ENDIF
11 CONTINUE
12 CONTINUE
END
Da fehlen julday und flmoon, um
den Spass zum Laufen zu kriegen.
Irgendwelche Bibliotheken?
MEB
Sorry, vergessen…hier:
SUBROUTINE FLMOON(N,NPH,JD,FRAC)
PARAMETER (RAD=0.017453293)
C=N+NPH/4.
T=C/1236.85
T2=T**2
AS=359.2242+29.105356*C
AM=306.0253+385.816918*C+0.010730*T2
JD=2415020+28*N+7*NPH
XTRA=0.75933+1.53058868*C+(1.178E-4-1.55E-7*T)*T2
IF(NPH.EQ.0.OR.NPH.EQ.2)THEN
XTRA=XTRA+(0.1734-3.93E-4*T)*SIN(RAD*AS)-0.4068*SIN(RAD*AM)
ELSE IF(NPH.EQ.1.OR.NPH.EQ.3)THEN
XTRA=XTRA+(0.1721-4.E-4*T)*SIN(RAD*AS)-0.6280*SIN(RAD*AM)
ELSE
PAUSE ‚NPH is unknown.‘
ENDIF
IF(XTRA.GE.0.)THEN
I=INT(XTRA)
ELSE
I=INT(XTRA-1.)
ENDIF
JD=JD+I
FRAC=XTRA-I
RETURN
END
C******************************************
FUNCTION JULDAY(MM,ID,IYYY)
PARAMETER (IGREG=15+31*(10+12*1582))
IF (IYYY.EQ.0) PAUSE ‚There is no Year Zero.‘
IF (IYYY.LT.0) IYYY=IYYY+1
IF (MM.GT.2) THEN
JY=IYYY
JM=MM+1
ELSE
JY=IYYY-1
JM=MM+13
ENDIF
JULDAY=INT(365.25*JY)+INT(30.6001*JM)+ID+1720995
IF (ID+31*(MM+12*IYYY).GE.IGREG) THEN
JA=INT(0.01*JY)
JULDAY=JULDAY+2-JA+INT(0.25*JA)
ENDIF
RETURN
END
besten Dank,
klappt nun
MEB