]> git.mxchange.org Git - simgear.git/blobdiff - simgear/timing/sg_time.cxx
Patch from Melchior Franz:
[simgear.git] / simgear / timing / sg_time.cxx
index 3fca4a8b0ef94ac893b2591e8b0cbd42ef0c3fcf..d2823b7794a9738efe26349a900fc11d154dc2e7 100644 (file)
@@ -4,31 +4,33 @@
 //
 // Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
 //
-// This program is free software; you can redistribute it and/or
-// modify it under the terms of the GNU General Public License as
-// published by the Free Software Foundation; either version 2 of the
-// License, or (at your option) any later version.
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Library General Public
+// License as published by the Free Software Foundation; either
+// version 2 of the License, or (at your option) any later version.
 //
-// This program is distributed in the hope that it will be useful, but
-// WITHOUT ANY WARRANTY; without even the implied warranty of
+// 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
-// General Public License for more details.
+// Library General Public License for more details.
 //
-// You should have received a copy of the GNU General Public License
-// along with this program; if not, write to the Free Software
-// Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+// You should have received a copy of the GNU Library General Public
+// License along with this library; if not, write to the
+// Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+// Boston, MA  02111-1307, USA.
 //
 // $Id$
 
 
 #ifdef HAVE_CONFIG_H
-#  include <config.h>
+#  include <simgear_config.h>
 #endif
 
 #include <simgear/compiler.h>
 
-#ifdef FG_HAVE_STD_INCLUDES
-#  include <cmath>
+#include <errno.h>             // for errno
+
+#ifdef SG_HAVE_STD_INCLUDES
 #  include <cstdio>
 #  include <cstdlib>
 #  include <ctime>
 #  include <math.h>
 #  include <stdio.h>
 #  include <stdlib.h>
-#  include <time.h>
 #endif
 
+#ifdef HAVE_SYS_TIME_H
+#  include <sys/time.h>  // for get/setitimer, gettimeofday, struct timeval
+#endif
 #ifdef HAVE_SYS_TIMEB_H
 #  include <sys/timeb.h> // for ftime() and struct timeb
 #endif
 #ifdef HAVE_UNISTD_H
 #  include <unistd.h>    // for gettimeofday()
 #endif
-#ifdef HAVE_SYS_TIME_H
-#  include <sys/time.h>  // for get/setitimer, gettimeofday, struct timeval
-#endif
+
+#include <math.h>        // for NAN
 
 #include <simgear/constants.h>
 #include <simgear/debug/logstream.hxx>
-#include <simgear/misc/fgpath.hxx>
+#include <simgear/misc/sg_path.hxx>
 
 #include "sg_time.hxx"
 #include "timezone.h"
 
 
 #define DEGHR(x)        ((x)/15.)
-#define RADHR(x)        DEGHR(x*RAD_TO_DEG)
-
+#define RADHR(x)        DEGHR(x*SGD_RADIANS_TO_DEGREES)
 
-SGTime::SGTime( const string& root )
-{
-    if (cur_time_params) {
-       FG_LOG( FG_GENERAL, FG_ALERT, 
-               "Error: only one instance of SGTime allowed" );
-       exit(-1);
-    }
 
-    cur_time_params = this;
+static const double MJD0    = 2415020.0;
+static const double J2000   = 2451545.0 - MJD0;
+static const double SIDRATE = 0.9972695677;
 
-    FGPath zone( root );
-    zone.append( "Timezone" );
-    zone.append( "zone.tab" );
 
-    FG_LOG( FG_EVENT, FG_DEBUG, "Reading timezone info from: " << zone.str() );
-    tzContainer = new TimezoneContainer( zone.c_str() );
-
-    zonename = NULL;
-}
+void SGTime::init( double lon, double lat,
+                   const string& root, time_t init_time )
+{
+    SG_LOG( SG_EVENT, SG_INFO, "Initializing Time" );
 
+    gst_diff = -9999.0;
 
-SGTime::~SGTime()
-{
-    delete tzContainer;
+    if ( init_time ) {
+       cur_time = init_time;
+    } else {
+       cur_time = time(NULL); 
+    }
 
-    if ( zonename != NULL ) {
-       delete zonename;
+    cout << "Current greenwich mean time = " << asctime(gmtime(&cur_time))
+         << endl;
+    cout << "Current local time          = " 
+         << asctime(localtime(&cur_time)) << endl;
+
+    if ( root != (string)"" ) {
+        SGPath zone( root );
+        zone.append( "zone.tab" );
+        SG_LOG( SG_EVENT, SG_INFO, "Reading timezone info from: "
+                << zone.str() );
+        tzContainer = new TimezoneContainer( zone.c_str() );
+
+        GeoCoord location( SGD_RADIANS_TO_DEGREES * lat, SGD_RADIANS_TO_DEGREES * lon );
+        GeoCoord* nearestTz = tzContainer->getNearest(location);
+
+        SGPath name( root );
+        name.append( nearestTz->getDescription() );
+        zonename = strdup( name.c_str() );
+        SG_LOG( SG_EVENT, SG_INFO, "Using zonename = " << zonename );
+    } else {
+        SG_LOG( SG_EVENT, SG_INFO, "*** NO TIME ZONE NAME ***" );
+        tzContainer = NULL;
+        zonename = NULL;
     }
 }
 
-void SGTime::updateLocal( double lon, double lat, const string& root )
+SGTime::SGTime( double lon, double lat, const string& root, time_t init_time )
 {
-  time_t currGMT;
-  time_t aircraftLocalTime;
-  GeoCoord location( RAD_TO_DEG * lat, RAD_TO_DEG * lon );
-  GeoCoord* nearestTz = tzContainer->getNearest(location);
-  FGPath zone( root );
-  zone.append ("Timezone" );
-  zone.append ( nearestTz->getDescription() );
-  if ( zonename ) {
-      delete zonename;
-  }
-  zonename = strdup( zone.c_str() );
-  currGMT = get_gmt( gmtime(&cur_time) );
-  aircraftLocalTime = get_gmt( (fgLocaltime(&cur_time, zone.c_str())) );
-  local_offset = aircraftLocalTime - currGMT;
-  // cout << "Using " << local_offset << " as local time offset Timezone is " 
-  //      << zonename << endl;
+    init( lon, lat, root, init_time );
 }
 
-// Initialize the time dependent variables (maybe I'll put this in the
-// constructor later)
-void SGTime::init( double lon, double lat, const string& root )
-// time_t timeOffset, sgTimingOffsetType offsetType )
-{
-    FG_LOG( FG_EVENT, FG_INFO, "Initializing Time" );
-    gst_diff = -9999.0;
-
-    // time_t currGMT;
-    // time_t systemLocalTime;
-    // time_t aircraftLocalTime;
-
-    // would it be better to put these sanity checks in the options
-    // parsing code? (CLO)
 
-    cur_time = time(NULL); 
-
-    // printf ("Current greenwich mean time = %24s", asctime(gmtime(&cur_time)));
-    // printf ("Current local time          = %24s", asctime(localtime(&cur_time)));
-    // time_t tmp = cur_time;
-    GeoCoord location( RAD_TO_DEG * lat, RAD_TO_DEG * lon );
-
-    GeoCoord* nearestTz = tzContainer->getNearest(location);
-
-    FGPath zone( root );
-    zone.append( "Timezone" );
-    zone.append( nearestTz->getDescription() );
-
-    // cout << "Using " << zone.str() << " for timezone information" << endl;
-    zonename = strdup( zone.c_str() );
-    // cout << "zonename = " << zonename << endl;
-    //show( buffer.c_str(), cur_time, 1); 
-    //printf ("Current greenwich mean time = %24s", asctime(gmtime(&cur_time)));
-    //printf ("Current local time          = %24s", asctime(localtime(&cur_time)));
-    // currGMT = get_gmt( gmtime(&cur_time) );
-    // cout << "currGMT = " << currGMT << endl;
-
-    // systemLocalTime = get_gmt( localtime(&cur_time) );
-    // cout << "systemLocalTime = " << systemLocalTime << endl;
-
-    // aircraftLocalTime = get_gmt( fgLocaltime(&cur_time, zone.c_str()) ); 
-    //printf ("Current greenwich mean time = %24s", asctime(gmtime(&cur_time)));
-    //printf ("Current local time          = %24s", asctime(localtime(&cur_time)));
+SGTime::SGTime( const string& root ) {
+    init( 0.0, 0.0, root, 0 );
 }
 
 
-// given a date in months, mn, days, dy, years, yr, return the
-// modified Julian date (number of days elapsed since 1900 jan 0.5),
-// mjd.  Adapted from Xephem.
-
-void SGTime::cal_mjd (int mn, double dy, int yr) 
-{
-    //static double last_mjd, last_dy;
-    //double mjd;
-    //static int last_mn, last_yr;
-    int b, d, m, y;
-    long c;
-  
-    if (mn == last_mn && yr == last_yr && dy == last_dy) {
-       mjd = last_mjd;
-    }
-  
-    m = mn;
-    y = (yr < 0) ? yr + 1 : yr;
-    if (mn < 3) {
-       m += 12;
-       y -= 1;
-    }
-  
-    if (yr < 1582 || (yr == 1582 && (mn < 10 || (mn == 10 && dy < 15)))) {
-       b = 0;
-    } else {
-       int a;
-       a = y/100;
-       b = 2 - a + a/4;
-    }
-  
-    if (y < 0) {
-       c = (long)((365.25*y) - 0.75) - 694025L;
-    } else {
-       c = (long)(365.25*y) - 694025L;
-    }
-  
-    d = (int)(30.6001*(m+1));
-  
-    mjd = b + c + d + dy - 0.5;
-  
-    last_mn = mn;
-    last_dy = dy;
-    last_yr = yr;
-    last_mjd = mjd;
+SGTime::SGTime() {
+    init( 0.0, 0.0, "", 0 );
 }
 
 
-// given an mjd, calculate greenwich mean sidereal time, gst
-void SGTime::utc_gst () 
+SGTime::~SGTime()
 {
-    double day = floor(mjd-0.5)+0.5;
-    double hr = (mjd-day)*24.0;
-    double T, x;
-
-    T = ((int)(mjd - 0.5) + 0.5 - J2000)/36525.0;
-    x = 24110.54841 + (8640184.812866 + (0.093104 - 6.2e-6 * T) * T) * T;
-    x /= 3600.0;
-    gst = (1.0/SIDRATE)*hr + x;
+    if ( tzContainer != NULL ) {
+        TimezoneContainer *tmp = tzContainer;
+        tzContainer = NULL;
+       delete tmp;
+    }
 
-    FG_LOG( FG_EVENT, FG_DEBUG, "  gst => " << gst );
+    if ( zonename != NULL ) {
+        char *tmp = zonename;
+        zonename = NULL;
+       free(tmp);
+    }
 }
 
 
@@ -226,48 +145,41 @@ void SGTime::utc_gst ()
 // Local Sidereal Time, in decimal hours.
 //
 // Provided courtesy of ecdowney@noao.edu (Elwood Downey) 
-
-double SGTime::sidereal_precise (double lng) 
+static double sidereal_precise( double mjd, double lng )
 {
-    double lstTmp;
-
     /* printf ("Current Lst on JD %13.5f at %8.4f degrees West: ", 
        mjd + MJD0, lng); */
 
     // convert to required internal units
-    lng *= DEG_TO_RAD;
+    lng *= SGD_DEGREES_TO_RADIANS;
 
     // compute LST and print
-    utc_gst();
-    lstTmp = gst - RADHR (lng);
-    lstTmp -= 24.0*floor(lstTmp/24.0);
+    double gst = sgTimeCalcGST( mjd );
+    double lst = gst - RADHR( lng );
+    lst -= 24.0 * floor( lst / 24.0 );
     // printf ("%7.4f\n", lstTmp);
 
-    // that's all
-    return (lstTmp);
+    return lst;
 }
 
 
 // return a courser but cheaper estimate of sidereal time
-double SGTime::sidereal_course(double lng) 
+static double sidereal_course( time_t cur_time, const struct tm *gmt, double lng )
 {
-    //struct tm *gmt;
-    //double lstTmp;
     time_t start_gmt, now;
     double diff, part, days, hours, lstTmp;
     char tbuf[64];
   
-    //gmt = t->gmt;
-    //now = t->cur_time;
     now = cur_time;
-    start_gmt = get_gmt(gmt->tm_year, 2, 21, 12, 0, 0);
+    start_gmt = sgTimeGetGMT(gmt->tm_year, 2, 21, 12, 0, 0);
   
-    FG_LOG( FG_EVENT, FG_DEBUG, "  COURSE: GMT = " << format_time(gmt, tbuf) );
-    FG_LOG( FG_EVENT, FG_DEBUG, "  March 21 noon (GMT) = " << start_gmt );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  COURSE: GMT = "
+           << sgTimeFormatTime(gmt, tbuf) );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  March 21 noon (GMT) = " << start_gmt );
   
     diff = (now - start_gmt) / (3600.0 * 24.0);
   
-    FG_LOG( FG_EVENT, FG_DEBUG, 
+    SG_LOG( SG_EVENT, SG_DEBUG, 
            "  Time since 3/21/" << gmt->tm_year << " GMT = " << diff );
   
     part = fmod(diff, 1.0);
@@ -280,40 +192,49 @@ double SGTime::sidereal_course(double lng)
        lstTmp += 24.0;
     }
   
-    FG_LOG( FG_EVENT, FG_DEBUG,
+    SG_LOG( SG_EVENT, SG_DEBUG,
            "  days = " << days << "  hours = " << hours << "  lon = " 
            << lng << "  lst = " << lstTmp );
   
-    return(lstTmp);
+    return lstTmp;
 }
 
 
-// Update time variables such as gmt, julian date, and sidereal time
-void SGTime::update( double lon, double lat, double alt_m, long int warp ) {
+// Update the time related variables
+void SGTime::update( double lon, double lat, time_t ct, long int warp ) {
     double gst_precise, gst_course;
 
-    FG_LOG( FG_EVENT, FG_DEBUG, "Updating time" );
+#if defined(_MSC_VER) || defined(__MINGW32__)
+    tm * gmt = &m_gmt;
+#endif
+
+    SG_LOG( SG_EVENT, SG_DEBUG, "Updating time" );
 
     // get current Unix calendar time (in seconds)
     // warp += warp_delta;
-    cur_time = time(NULL) + warp;
-    FG_LOG( FG_EVENT, FG_DEBUG, 
+    if ( ct ) {
+       cur_time = ct + warp;
+    } else {
+       cur_time = time(NULL) + warp;
+    }
+    SG_LOG( SG_EVENT, SG_DEBUG, 
            "  Current Unix calendar time = " << cur_time 
            << "  warp = " << warp );
 
     // get GMT break down for current time
+#if defined(_MSC_VER) || defined(__MINGW32__)
+    memcpy( gmt, gmtime(&cur_time), sizeof(tm) );
+#else
     gmt = gmtime(&cur_time);
-    FG_LOG( FG_EVENT, FG_DEBUG, 
+#endif
+    SG_LOG( SG_EVENT, SG_DEBUG, 
            "  Current GMT = " << gmt->tm_mon+1 << "/" 
-           << gmt->tm_mday << "/" << gmt->tm_year << " "
+           << gmt->tm_mday << "/" << (1900 + gmt->tm_year) << " "
            << gmt->tm_hour << ":" << gmt->tm_min << ":" 
            << gmt->tm_sec );
 
-    // calculate modified Julian date
-    // t->mjd = cal_mjd ((int)(t->gmt->tm_mon+1), (double)t->gmt->tm_mday, 
-    //     (int)(t->gmt->tm_year + 1900));
-    cal_mjd ((int)(gmt->tm_mon+1), (double)gmt->tm_mday, 
-            (int)(gmt->tm_year + 1900));
+    // calculate modified Julian date starting with current
+    mjd = sgTimeCurrentMJD( ct, warp );
 
     // add in partial day
     mjd += (gmt->tm_hour / 24.0) + (gmt->tm_min / (24.0 * 60.0)) +
@@ -321,36 +242,220 @@ void SGTime::update( double lon, double lat, double alt_m, long int warp ) {
 
     // convert "back" to Julian date + partial day (as a fraction of one)
     jd = mjd + MJD0;
-    FG_LOG( FG_EVENT, FG_DEBUG, "  Current Julian Date = " << jd );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  Current Julian Date = " << jd );
 
-    // printf("  Current Longitude = %.3f\n", FG_Longitude * RAD_TO_DEG);
+    // printf("  Current Longitude = %.3f\n", FG_Longitude * SGD_RADIANS_TO_DEGREES);
 
     // Calculate local side real time
     if ( gst_diff < -100.0 ) {
        // first time through do the expensive calculation & cheap
         // calculation to get the difference.
-       FG_LOG( FG_EVENT, FG_INFO, "  First time, doing precise gst" );
-       gst_precise = gst = sidereal_precise(0.00);
-       gst_course = sidereal_course(0.00);
+       SG_LOG( SG_EVENT, SG_INFO, "  First time, doing precise gst" );
+       gst_precise = gst = sidereal_precise( mjd, 0.00 );
+       gst_course = sidereal_course( cur_time, gmt, 0.00 );
       
        gst_diff = gst_precise - gst_course;
 
-       lst = sidereal_course(-(lon * RAD_TO_DEG)) + gst_diff;
+       lst = sidereal_course( cur_time, gmt, -(lon * SGD_RADIANS_TO_DEGREES) ) + gst_diff;
     } else {
        // course + difference should drift off very slowly
-       gst = sidereal_course( 0.00 ) + gst_diff;
-       lst = sidereal_course( -(lon * RAD_TO_DEG)) + gst_diff;
+       gst = sidereal_course( cur_time, gmt, 0.00 ) + gst_diff;
+       lst = sidereal_course( cur_time, gmt, -(lon * SGD_RADIANS_TO_DEGREES) ) + gst_diff;
     }
 
-    FG_LOG( FG_EVENT, FG_DEBUG,
+    SG_LOG( SG_EVENT, SG_DEBUG,
            "  Current lon=0.00 Sidereal Time = " << gst );
-    FG_LOG( FG_EVENT, FG_DEBUG,
+    SG_LOG( SG_EVENT, SG_DEBUG,
            "  Current LOCAL Sidereal Time = " << lst << " (" 
-           << sidereal_precise(-(lon * RAD_TO_DEG)
+           << sidereal_precise( mjd, -(lon * SGD_RADIANS_TO_DEGREES) 
            << ") (diff = " << gst_diff << ")" );
 }
 
 
+// Given lon/lat, update timezone information and local_offset
+void SGTime::updateLocal( double lon, double lat, const string& root ) {
+    // sanity checking
+    if ( lon < -SGD_PI || lon > SGD_PI ) {
+        // not within -180 ... 180
+        lon = 0.0;
+    }
+    if ( lat < -SGD_PI * 0.5 || lat > SGD_PI * 0.5 ) {
+        // not within -90 ... 90
+        lat = 0.0;
+    }
+    if ( lon != lon ) {
+        // only true if lon == nan
+        SG_LOG( SG_EVENT, SG_ALERT, "  Detected lon == nan, resetting to 0.0" );
+        lon = 0.0;
+    }
+    if ( lat != lat ) {
+        // only true if lat == nan
+        SG_LOG( SG_EVENT, SG_ALERT, "  Detected lat == nan, resetting to 0.0" );
+        lat = 0.0;
+    }
+    time_t currGMT;
+    time_t aircraftLocalTime;
+    GeoCoord location( SGD_RADIANS_TO_DEGREES * lat, SGD_RADIANS_TO_DEGREES * lon );
+    GeoCoord* nearestTz = tzContainer->getNearest(location);
+    SGPath zone( root );
+    zone.append ( nearestTz->getDescription() );
+    if ( zonename ) {
+        char *ptr = zonename;
+        zonename = NULL;
+        free(ptr);
+    }
+    zonename = strdup( zone.c_str() );
+
+    //Avoid troubles when zone.tab hasn't got the right line endings
+    if (zonename[strlen(zonename)-1] == '\r')
+    {
+      zonename[strlen(zonename)-1]=0;
+      zone.set( zonename );
+    }
+
+    currGMT = sgTimeGetGMT( gmtime(&cur_time) );
+    aircraftLocalTime = sgTimeGetGMT( (fgLocaltime(&cur_time, zone.c_str())) );
+    local_offset = aircraftLocalTime - currGMT;
+    // cout << "Using " << local_offset << " as local time offset Timezone is " 
+    //      << zonename << endl;
+}
+
+
+// given a date in months, mn, days, dy, years, yr, return the
+// modified Julian date (number of days elapsed since 1900 jan 0.5),
+// mjd.  Adapted from Xephem.
+double sgTimeCalcMJD(int mn, double dy, int yr) {
+    double mjd;
+
+    // internal book keeping data
+    static double last_mjd, last_dy;
+    static int last_mn, last_yr;
+
+    int b, d, m, y;
+    long c;
+  
+    if (mn == last_mn && yr == last_yr && dy == last_dy) {
+       mjd = last_mjd;
+    }
+  
+    m = mn;
+    y = (yr < 0) ? yr + 1 : yr;
+    if (mn < 3) {
+       m += 12;
+       y -= 1;
+    }
+  
+    if (yr < 1582 || (yr == 1582 && (mn < 10 || (mn == 10 && dy < 15)))) {
+       b = 0;
+    } else {
+       int a;
+       a = y/100;
+       b = 2 - a + a/4;
+    }
+  
+    if (y < 0) {
+       c = (long)((365.25*y) - 0.75) - 694025L;
+    } else {
+       c = (long)(365.25*y) - 694025L;
+    }
+  
+    d = (int)(30.6001*(m+1));
+  
+    mjd = b + c + d + dy - 0.5;
+  
+    last_mn = mn;
+    last_dy = dy;
+    last_yr = yr;
+    last_mjd = mjd;
+
+    return mjd;
+}
+
+
+// return the current modified Julian date (number of days elapsed
+// since 1900 jan 0.5), mjd.
+double sgTimeCurrentMJD( time_t ct, long int warp ) {
+
+#if defined(_MSC_VER) || defined(__MINGW32__)
+    struct tm m_gmt;    // copy of system gmtime(&time_t) structure
+    struct tm *gmt = &m_gmt;
+#else
+    struct tm *gmt;
+#endif
+
+    // get current Unix calendar time (in seconds)
+    // warp += warp_delta;
+    time_t cur_time;
+    if ( ct ) {
+        cur_time = ct + warp;
+    } else {
+        cur_time = time(NULL) + warp;
+    }
+    SG_LOG( SG_EVENT, SG_DEBUG, 
+           "  Current Unix calendar time = " << cur_time 
+           << "  warp = " << warp );
+
+    // get GMT break down for current time
+#if defined(_MSC_VER) || defined(__MINGW32__)
+    memcpy( gmt, gmtime(&cur_time), sizeof(tm) );
+#else
+    gmt = gmtime(&cur_time);
+#endif
+    SG_LOG( SG_EVENT, SG_DEBUG, 
+           "  Current GMT = " << gmt->tm_mon+1 << "/" 
+           << gmt->tm_mday << "/" << (1900 + gmt->tm_year) << " "
+           << gmt->tm_hour << ":" << gmt->tm_min << ":" 
+           << gmt->tm_sec );
+
+    // calculate modified Julian date
+    // t->mjd = cal_mjd ((int)(t->gmt->tm_mon+1), (double)t->gmt->tm_mday, 
+    //     (int)(t->gmt->tm_year + 1900));
+    double mjd = sgTimeCalcMJD( (int)(gmt->tm_mon+1), (double)gmt->tm_mday, 
+                                (int)(gmt->tm_year + 1900) );
+
+    return mjd;
+}
+
+
+// given an mjd, calculate greenwich mean sidereal time, gst
+double sgTimeCalcGST( double mjd ) {
+    double gst;
+
+    double day = floor(mjd-0.5)+0.5;
+    double hr = (mjd-day)*24.0;
+    double T, x;
+
+    T = ((int)(mjd - 0.5) + 0.5 - J2000)/36525.0;
+    x = 24110.54841 + (8640184.812866 + (0.093104 - 6.2e-6 * T) * T) * T;
+    x /= 3600.0;
+    gst = (1.0/SIDRATE)*hr + x;
+
+    SG_LOG( SG_EVENT, SG_DEBUG, "  gst => " << gst );
+
+    return gst;
+}
+
+
+#if defined( HAVE_TIMEGM ) 
+    // ignore this function
+#elif defined( MK_TIME_IS_GMT )
+    // ignore this function
+#else // ! defined ( MK_TIME_IS_GMT )
+
+    // Fix up timezone if using ftime()
+    static long int fix_up_timezone( long int timezone_orig ) {
+#   if !defined( HAVE_GETTIMEOFDAY ) && defined( HAVE_FTIME )
+       // ftime() needs a little extra help finding the current timezone
+       struct timeb current;
+       ftime(&current);
+       return( current.timezone * 60 );
+#   else
+        return( timezone_orig );
+#   endif
+    }
+#endif
+
+
 /******************************************************************
  * The following are some functions that were included as SGTime
  * members, although they currently don't make use of any of the
@@ -382,7 +487,7 @@ void SGTime::update( double lon, double lat, double alt_m, long int warp ) {
 // If you are having problems with incorrectly positioned astronomical
 // bodies, this is a really good place to start looking.
 
-time_t SGTime::get_gmt(int year, int month, int day, int hour, int min, int sec)
+time_t sgTimeGetGMT(int year, int month, int day, int hour, int min, int sec)
 {
     struct tm mt;
 
@@ -402,77 +507,86 @@ time_t SGTime::get_gmt(int year, int month, int day, int hour, int min, int sec)
 #  define MK_TIME_IS_GMT 1
 #endif
 
-#if defined( HAVE_TIMEGM ) 
+#if defined( HAVE_TIMEGM )
     return ( timegm(&mt) );
 #elif defined( MK_TIME_IS_GMT )
-    return ( mktime(&mt) );
+    time_t ret = mktime(&mt);
+
+#ifdef __CYGWIN__
+       ret -= _timezone;
+#endif
+
+    // This is necessary as some mktime() calls may
+    // try to access the system timezone files
+    // if this open fails errno is set to 2
+    // CYGWIN for one does this
+    // if ( errno ) {
+    //     perror( "sgTimeGetGMT()" );
+    //     errno = 0;
+    // }
+
+    // reset errno in any event.
+    errno = 0;
+
+    return ret;
 #else // ! defined ( MK_TIME_IS_GMT )
 
     // timezone seems to work as a proper offset for Linux & Solaris
-#   if defined( __linux__ ) || defined( __sun__ ) 
+#   if defined( __linux__ ) || defined( __sun__ ) ||defined(__CYGWIN__)
 #       define TIMEZONE_OFFSET_WORKS 1
 #   endif
 
-    long int start = mktime(&mt);
+#if defined(__CYGWIN__)
+#define TIMEZONE _timezone
+#else
+#define TIMEZONE timezone
+#endif
+       
+    time_t start = mktime(&mt);
 
-    FG_LOG( FG_EVENT, FG_DEBUG, "start1 = " << start );
+    SG_LOG( SG_EVENT, SG_DEBUG, "start1 = " << start );
     // the ctime() call can screw up time progression on some versions
     // of Linux
-    // fgPrintf( FG_EVENT, FG_DEBUG, "start2 = %s", ctime(&start));
-    FG_LOG( FG_EVENT, FG_DEBUG, "(tm_isdst = " << mt.tm_isdst << ")" );
+    // fgPrintf( SG_EVENT, SG_DEBUG, "start2 = %s", ctime(&start));
+    SG_LOG( SG_EVENT, SG_DEBUG, "(tm_isdst = " << mt.tm_isdst << ")" );
 
-    timezone = fix_up_timezone( timezone );
+    TIMEZONE = fix_up_timezone( TIMEZONE );
 
 #  if defined( TIMEZONE_OFFSET_WORKS )
-    FG_LOG( FG_EVENT, FG_DEBUG, 
-           "start = " << start << ", timezone = " << timezone );
-    return( start - timezone );
+    SG_LOG( SG_EVENT, SG_DEBUG,
+           "start = " << start << ", timezone = " << TIMEZONE );
+    return( start - TIMEZONE );
 #  else // ! defined( TIMEZONE_OFFSET_WORKS )
 
     daylight = mt.tm_isdst;
     if ( daylight > 0 ) {
        daylight = 1;
     } else if ( daylight < 0 ) {
-       FG_LOG( FG_EVENT, FG_WARN, 
+       SG_LOG( SG_EVENT, FG_WARN, 
                "OOOPS, problem in sg_time.cxx, no daylight savings info." );
     }
 
-    long int offset = -(timezone / 3600 - daylight);
+    long int offset = -(TIMEZONE / 3600 - daylight);
 
-    FG_LOG( FG_EVENT, FG_DEBUG, "  Raw time zone offset = " << timezone );
-    FG_LOG( FG_EVENT, FG_DEBUG, "  Daylight Savings = " << daylight );
-    FG_LOG( FG_EVENT, FG_DEBUG, "  Local hours from GMT = " << offset );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  Raw time zone offset = " << TIMEZONE );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  Daylight Savings = " << daylight );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  Local hours from GMT = " << offset );
     
-    long int start_gmt = start - timezone + (daylight * 3600);
+    long int start_gmt = start - TIMEZONE + (daylight * 3600);
     
-    FG_LOG( FG_EVENT, FG_DEBUG, "  March 21 noon (CST) = " << start );
+    SG_LOG( SG_EVENT, SG_DEBUG, "  March 21 noon (CST) = " << start );
 
     return ( start_gmt );
 #  endif // ! defined( TIMEZONE_OFFSET_WORKS )
 #endif // ! defined ( MK_TIME_IS_GMT )
 }
 
-// Fix up timezone if using ftime()
-long int SGTime::fix_up_timezone( long int timezone_orig ) 
-{
-#if !defined( HAVE_GETTIMEOFDAY ) && defined( HAVE_FTIME )
-    // ftime() needs a little extra help finding the current timezone
-    struct timeb current;
-    ftime(&current);
-    return( current.timezone * 60 );
-#else
-    return( timezone_orig );
-#endif
-}
 
-
-char* SGTime::format_time( const struct tm* p, char* buf )
+// format time
+char* sgTimeFormatTime( const struct tm* p, char* buf )
 {
     sprintf( buf, "%d/%d/%2d %d:%02d:%02d", 
             p->tm_mon, p->tm_mday, p->tm_year,
             p->tm_hour, p->tm_min, p->tm_sec);
     return buf;
 }
-
-
-SGTime* SGTime::cur_time_params = 0;