/*
* GeoTools - The Open Source Java GIS Toolkit
* http://geotools.org
*
* (C) 2000-2008, Open Source Geospatial Foundation (OSGeo)
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation;
* version 2.1 of the License.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* NOTE: permission has been given to the JScience project (http://www.jscience.org)
* to distribute this file under BSD-like license.
*/
package org.geotools.nature;
// J2SE dependencies
import java.text.DateFormat;
import java.text.ParseException;
import java.util.Date;
import java.util.TimeZone;
/**
* Approximations de quelques calculs astronomiques relatifs aux calendriers terrestres.
* Les différents cycles astronomiques (notamment le jour, le mois et l'année) ne sont pas
* constants. Par exemple, la longueur de l'année tropicale (le nombre moyen de jours entre
* deux équinoxes vernales) était d'environ 365,242196 jours en 1900 et devrait être d'environ
* 365,242184 jours en 2100, soit un changement d'environ 1 seconde. Cette classe permet de
* calculer la longueur d'une année ou d'un mois à une date spécifiée. Toutefois, il est
* important de noter que les intervalles de temps calculés par les méthodes de cette classe
* sont des <strong>moyennes</strong>. Pour une année en particulier, l'intervalle de temps
* d'un équinoxe vernale au prochain peut s'écarter de cette moyenne de plusieurs minutes.
*
* <p>Les calculs de la longueur de l'année tropicale sont basés sur les travaux de Laskar (1986).
* Les calculs de la longueur des mois synodiques sont basés sur les travaux de Chapront-Touze et
* Chapront (1988).On peut lire plus de détails au sujet des calendrier terrestre au site
* <a href="http://webexhibits.org/calendars/year-astronomy.html">http://webexhibits.org/calendars/year-astronomy.html</a> ainsi que
* <a href="http://www.treasure-troves.com/astro/TropicalYear.html">http://www.treasure-troves.com/astro/TropicalYear.html</a>.</p>
*
*
* @source $URL$
* @version $Id$
* @author Martin Desruisseaux (IRD)
*
* @since 2.1
*/
public final class Calendar {
/**
* Nombre de millisecondes dans une journée. Cette constante est
* utilisée pour convertir des intervalles de temps du Java en
* nombre de jours.
*/
private static double MILLIS_IN_DAY = 1000*60*60*24;
/**
* Jour julien correspondant à l'époch du Java (1er janvier 1970 à minuit).
* Cette constante est utilisée pour convertir des dates du Java en jour
* julien.
*
* La valeur {@link #julianDay} du 1er janvier 2000 00:00 GMT est 2451544.5 jours.
* La valeur {@link Date#getTime} du 1er janvier 2000 00:00 GMT est 10957 jours.
*/
private static double JULIAN_DAY_1970 = 2451544.5-10957;
/**
* Interdit la création de classes {@code Cycles} par l'utilisateur.
*/
private Calendar() {
}
/**
* Retourne le jour julien d'une date. Il ne s'agit pas du jour julien dans
* une année. Ce jour julien là (nommé ainsi pour <i>Julius Scaliger</i>, et
* non <i>Julius Caesar</i>) est le nombre de jours écoulés depuis midi le
* 1er janvier 4713 avant Jésus-Christ.
*/
public static double julianDay(final Date time) {
return julianDay(time.getTime());
}
/**
* Computes the {@linkplain #julianDay(Date) julian day}.
*
* @param time The date in milliseconds ellapsed since January 1st, 1970.
*/
static double julianDay(final long time) {
return (time/MILLIS_IN_DAY) + JULIAN_DAY_1970;
}
/**
* Retourne le nombre de siècles écoulés depuis le 1 janvier 2000 à midi.
* Cette information est utilisée dans les formules de Laskar (1986) pour
* calculer la longueur d'une année tropicale, ainsi que par Chapront-Touze
* et Chapront (1988) pour la longueur d'un mois synodique.
*/
static double julianCentury(final Date time) {
return ((time.getTime()/MILLIS_IN_DAY) + (JULIAN_DAY_1970-2451545.0))/36525;
}
/**
* Retourne la longueur de l'année tropicale. L'année tropicale est définie comme l'intervalle
* moyen entre deux équinoxes vernales (autour du 21 mars dans l'hémisphère nord). Il correspond
* au cycle des saisons. Cet intervalle de temps est une <strong>moyenne</strong>. Un cycle réel
* peut s'écarter de plusieurs minutes de cette moyenne. Notez aussi qu'une année tropicale
* n'est pas identique à une année sidérale, qui est le temps requis par la Terre pour compléter
* un orbite autour du Soleil. En l'an 2000, l'année tropicale avait une longueur d'environ
* 365,2422 jours tandis que l'année sidérale avait une longueur de 365,2564 jours.
*/
public static double tropicalYearLength(final Date time) {
final double T = julianCentury(time);
return 365.2421896698 + T*(-0.00000615359 + T*(-7.29E-10 + T*(2.64E-10)));
}
/**
* Retourne la longueur du mois synodique. Le mois synodique est l'intervalle de temps moyen
* entre deux conjonctions de la lune et du soleil. Il correspond au cycle des phases de la
* lune. Cet intervalle de temps est une <strong>moyenne</strong>. Un cycle réel peut s'écarter
* de plusieurs heures de cette moyenne.
*/
public static double synodicMonthLength(final Date time) {
final double T=julianCentury(time);
return 29.5305888531 + T*(0.00000021621 + T*(-3.64E-10));
}
/**
* Affiche la longueur de l'année tropicale et du mois synodique pour une date donnée.
* Cette application peut être lancée avec la syntaxe suivante:
*
* <pre>Calendar <var><date></var></pre>
*
* où <var>date</var> est un argument optionel spécifiant la date (jour, mois et année)
* d'intérêt en heure universelle (UTC). Si cet argument est omis, la date et heure
* actuelles seront utilisées.
*/
public static final void main(final String[] args) throws ParseException {
final DateFormat format=DateFormat.getDateInstance(DateFormat.SHORT);
format.setTimeZone(TimeZone.getTimeZone("UTC"));
final Date time=(args.length!=0) ? format.parse(args[0]) : new Date();
System.out.print("Date (UTC) : "); System.out.println(format.format(time));
System.out.print("Tropical year: "); System.out.println(tropicalYearLength(time));
System.out.print("Synodic month: "); System.out.println(synodicMonthLength(time));
}
}