flightgear/src/Main/fg_init.cxx

// fg_init.cxx -- Flight Gear top level initialization routines
//
// Written by Curtis Olson, started August 1997.
//
// 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 program 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.
//
// 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.
//
//
// $Id$


#ifdef HAVE_CONFIG_H
#  include <config.h>
#endif

// For BC 5.01 this must be included before OpenGL includes.
#ifdef SG_MATH_EXCEPTION_CLASH
#  include <math.h>
#endif

#include <GL/glut.h>
#include <simgear/xgl/xgl.h>

#include <stdio.h>
#include <stdlib.h>


#if defined( unix ) || defined( __CYGWIN__ )
#  include <unistd.h>		// for gethostname()
#endif

// work around a stdc++ lib bug in some versions of linux, but doesn't
// seem to hurt to have this here for all versions of Linux.
#ifdef linux
#  define _G_NO_EXTERN_TEMPLATES
#endif

#include <simgear/compiler.h>

#include STL_STRING

#include <simgear/constants.h>
#include <simgear/debug/logstream.hxx>
#include <simgear/math/point3d.hxx>
#include <simgear/math/polar3d.hxx>
#include <simgear/math/sg_geodesy.hxx>
#include <simgear/misc/fgpath.hxx>
#include <simgear/timing/sg_time.hxx>

#include <Aircraft/aircraft.hxx>
#include <FDM/UIUCModel/uiuc_aircraftdir.h>
#include <Airports/runways.hxx>
#include <Airports/simple.hxx>
#include <Autopilot/auto_gui.hxx>
#include <Autopilot/newauto.hxx>
#include <Cockpit/cockpit.hxx>
#include <Cockpit/radiostack.hxx>
#include <Cockpit/panel.hxx>
#include <Cockpit/panel_io.hxx>
#include <FDM/ADA.hxx>
#include <FDM/Balloon.h>
#include <FDM/External.hxx>
#include <FDM/JSBSim.hxx>
#include <FDM/LaRCsim.hxx>
#include <FDM/MagicCarpet.hxx>
#include <Include/general.hxx>
#include <Joystick/joystick.hxx>
#include <Objects/matlib.hxx>
#include <Navaids/fixlist.hxx>
#include <Navaids/ilslist.hxx>
#include <Navaids/mkrbeacons.hxx>
#include <Navaids/navlist.hxx>
#include <Scenery/scenery.hxx>
#include <Scenery/tilemgr.hxx>
#include <Time/event.hxx>
#include <Time/light.hxx>
#include <Time/sunpos.hxx>
#include <Time/moonpos.hxx>
#include <Time/tmp.hxx>

#ifndef FG_OLD_WEATHER
#  include <WeatherCM/FGLocalWeatherDatabase.h>
#else
#  include <Weather/weather.hxx>
#endif

#include "fg_init.hxx"
#include "fg_io.hxx"
#include "options.hxx"
#include "globals.hxx"
#include "bfi.hxx"

#if defined(FX) && defined(XMESA)
#include <GL/xmesa.h>
#endif

SG_USING_STD(string);

extern const char *default_root;


// Read in configuration (file and command line) and just set fg_root
bool fgInitFGRoot ( int argc, char **argv ) {
    string root;
    char* envp;

    // First parse command line options looking for fg-root, this will
    // override anything specified in a config file
    root = fgScanForRoot(argc, argv);

#if defined( unix ) || defined( __CYGWIN__ )
    // Next check home directory for .fgfsrc.hostname file
    if ( root == "" ) {
	envp = ::getenv( "HOME" );
	if ( envp != NULL ) {
	    FGPath config( envp );
	    config.append( ".fgfsrc" );
	    char name[256];
	    gethostname( name, 256 );
	    config.concat( "." );
	    config.concat( name );
 	    root = fgScanForRoot(config.str());
	}
    }
#endif

    // Next check home directory for .fgfsrc file
    if ( root == "" ) {
	envp = ::getenv( "HOME" );
	if ( envp != NULL ) {
	    FGPath config( envp );
	    config.append( ".fgfsrc" );
	    root = fgScanForRoot(config.str());
	}
    }
    
    // Next check if fg-root is set as an env variable
    if ( root == "" ) {
	envp = ::getenv( "FG_ROOT" );
	if ( envp != NULL ) {
	    root = envp;
	}
    }

    // Otherwise, default to a random compiled-in location if we can't
    // find fg-root any other way.
    if ( root == "" ) {
#if defined( WIN32 )
	root = "\\FlightGear";
#elif defined( macintosh )
	root = "";
#else
	root = PKGLIBDIR;
#endif
    }

    FG_LOG(FG_INPUT, FG_INFO, "fg_root = " << root );
    globals->set_fg_root(root);

    return true;
}


// Read in configuration (file and command line)
bool fgInitConfig ( int argc, char **argv ) {

				// First, set some sane default values
    fgSetDefaults();

    // Read global preferences from $FG_ROOT/preferences.xml
    FGPath props_path(globals->get_fg_root());
    props_path.append("preferences.xml");
    FG_LOG(FG_INPUT, FG_INFO, "Reading global preferences");
    if (!readProperties(props_path.str(), globals->get_props())) {
      FG_LOG(FG_INPUT, FG_ALERT, "Failed to read global preferences from "
	     << props_path.str());
    } else {
      FG_LOG(FG_INPUT, FG_INFO, "Finished Reading global preferences");
    }

    // Attempt to locate and parse the various config files in order
    // from least precidence to greatest precidence

    // Check for $fg_root/system.fgfsrc
    FGPath config( globals->get_fg_root() );
    config.append( "system.fgfsrc" );
    fgParseOptions(config.str());

    char name[256];
#if defined( unix ) || defined( __CYGWIN__ )
    // Check for $fg_root/system.fgfsrc.hostname
    gethostname( name, 256 );
    config.concat( "." );
    config.concat( name );
    fgParseOptions(config.str());
#endif

    // Check for ~/.fgfsrc
    char* envp = ::getenv( "HOME" );
    if ( envp != NULL ) {
	config.set( envp );
	config.append( ".fgfsrc" );
	fgParseOptions(config.str());
    }

#if defined( unix ) || defined( __CYGWIN__ )
    // Check for ~/.fgfsrc.hostname
    gethostname( name, 256 );
    config.concat( "." );
    config.concat( name );
    fgParseOptions(config.str());
#endif

    // Parse remaining command line options
    // These will override anything specified in a config file
    fgParseOptions(argc, argv);

    return true;
}


// find basic airport location info from airport database
bool fgFindAirportID( const string& id, FGAirport *a ) {
    if ( id.length() ) {
	FGPath path( globals->get_fg_root() );
	path.append( "Airports" );
	path.append( "simple.mk4" );
	FGAirports airports( path.c_str() );

	FG_LOG( FG_GENERAL, FG_INFO, "Searching for airport code = " << id );

	if ( ! airports.search( id, a ) ) {
	    FG_LOG( FG_GENERAL, FG_ALERT,
		    "Failed to find " << id << " in " << path.str() );
	    return false;
	}
    } else {
	return false;
    }

    FG_LOG( FG_GENERAL, FG_INFO,
	    "Position for " << id << " is ("
	    << a->longitude << ", "
	    << a->latitude << ")" );

    return true;
}


// Set current_options lon/lat given an airport id
bool fgSetPosFromAirportID( const string& id ) {
    FGAirport a;
    // double lon, lat;

    FG_LOG( FG_GENERAL, FG_INFO,
	    "Attempting to set starting position from airport code " << id );

    if ( fgFindAirportID( id, &a ) ) {
	fgSetDouble("/position/longitude",  a.longitude );
	fgSetDouble("/position/latitude",  a.latitude );
	FG_LOG( FG_GENERAL, FG_INFO,
		"Position for " << id << " is ("
		<< a.longitude << ", "
		<< a.latitude << ")" );

	return true;
    } else {
	return false;
    }

}


// Set current_options lon/lat given an airport id and heading (degrees)
bool fgSetPosFromAirportIDandHdg( const string& id, double tgt_hdg ) {
    FGRunway r;
    FGRunway found_r;
    double found_dir = 0.0;

    if ( id.length() ) {
	// set initial position from runway and heading

	FGPath path( globals->get_fg_root() );
	path.append( "Airports" );
	path.append( "runways.mk4" );
	FGRunways runways( path.c_str() );

	FG_LOG( FG_GENERAL, FG_INFO,
		"Attempting to set starting position from runway code "
		<< id << " heading " << tgt_hdg );

	// FGPath inpath( globals->get_fg_root() );
	// inpath.append( "Airports" );
	// inpath.append( "apt_simple" );
	// airports.load( inpath.c_str() );

	// FGPath outpath( globals->get_fg_root() );
	// outpath.append( "Airports" );
	// outpath.append( "simple.gdbm" );
	// airports.dump_gdbm( outpath.c_str() );

	if ( ! runways.search( id, &r ) ) {
	    FG_LOG( FG_GENERAL, FG_ALERT,
		    "Failed to find " << id << " in database." );
	    return false;
	}

	double diff;
	double min_diff = 360.0;

	while ( r.id == id ) {
	    // forward direction
	    diff = tgt_hdg - r.heading;
	    while ( diff < -180.0 ) { diff += 360.0; }
	    while ( diff >  180.0 ) { diff -= 360.0; }
	    diff = fabs(diff);
	    FG_LOG( FG_GENERAL, FG_INFO,
		    "Runway " << r.rwy_no << " heading = " << r.heading <<
		    " diff = " << diff );
	    if ( diff < min_diff ) {
		min_diff = diff;
		found_r = r;
		found_dir = 0;
	    }

	    // reverse direction
	    diff = tgt_hdg - r.heading - 180.0;
	    while ( diff < -180.0 ) { diff += 360.0; }
	    while ( diff >  180.0 ) { diff -= 360.0; }
	    diff = fabs(diff);
	    FG_LOG( FG_GENERAL, FG_INFO,
		    "Runway -" << r.rwy_no << " heading = " <<
		    r.heading + 180.0 <<
		    " diff = " << diff );
	    if ( diff < min_diff ) {
		min_diff = diff;
		found_r = r;
		found_dir = 180.0;
	    }

	    runways.next( &r );
	}

	FG_LOG( FG_GENERAL, FG_INFO, "closest runway = " << found_r.rwy_no
		<< " + " << found_dir );

    } else {
	return false;
    }

    double heading = found_r.heading + found_dir;
    while ( heading >= 360.0 ) { heading -= 360.0; }

    double lat2, lon2, az2;
    double azimuth = found_r.heading + found_dir + 180.0;
    while ( azimuth >= 360.0 ) { azimuth -= 360.0; }

    FG_LOG( FG_GENERAL, FG_INFO,
	    "runway =  " << found_r.lon << ", " << found_r.lat
	    << " length = " << found_r.length * SG_FEET_TO_METER * 0.5 
	    << " heading = " << azimuth );
    
    geo_direct_wgs_84 ( 0, found_r.lat, found_r.lon, 
			azimuth, found_r.length * SG_FEET_TO_METER * 0.5 - 5.0,
			&lat2, &lon2, &az2 );

    if ( fabs( fgGetDouble("/sim/startup/offset-distance") ) > SG_EPSILON ) {
	double olat, olon;
	double odist = fgGetDouble("/sim/startup/offset-distance");
	odist *= SG_NM_TO_METER;
	double oaz = azimuth;
	if ( fabs(fgGetDouble("/sim/startup/offset-azimuth")) > SG_EPSILON ) {
	    oaz = fgGetDouble("/sim/startup/offset-azimuth") + 180;
	}
	while ( oaz >= 360.0 ) { oaz -= 360.0; }
	geo_direct_wgs_84 ( 0, lat2, lon2, oaz, odist, &olat, &olon, &az2 );
	lat2=olat;
	lon2=olon;
    }
    fgSetDouble("/position/longitude",  lon2 );
    fgSetDouble("/position/latitude",  lat2 );
    fgSetDouble("/orientation/heading", heading );

    FG_LOG( FG_GENERAL, FG_INFO,
	    "Position for " << id << " is ("
	    << lon2 << ", "
	    << lat2 << ") new heading is "
	    << heading );

    return true;
}


// Set initial position and orientation
bool fgInitPosition( void ) {
    FGInterface *f = current_aircraft.fdm_state;
    string id = fgGetString("/sim/startup/airport-id");

    // set initial position from default or command line coordinates
    f->set_Longitude( fgGetDouble("/position/longitude") * SGD_DEGREES_TO_RADIANS );
    f->set_Latitude( fgGetDouble("/position/latitude") * SGD_DEGREES_TO_RADIANS );

    FG_LOG( FG_GENERAL, FG_INFO,
	    "scenery.cur_elev = " << scenery.cur_elev );
    FG_LOG( FG_GENERAL, FG_INFO,
	    "/position/altitude = " << fgGetDouble("/position/altitude") );

    // if we requested on ground startups
    if ( fgGetBool( "/sim/startup/onground" ) ) {
        fgSetDouble( "/position/altitude", (scenery.cur_elev + 1)
		     * METERS_TO_FEET );
    }

    // if requested altitude is below ground level
    if ( scenery.cur_elev >
	 fgGetDouble("/position/altitude") * METERS_TO_FEET - 1) {
	fgSetDouble("/position/altitude",
		    (scenery.cur_elev + 1) * METERS_TO_FEET );
    }

    FG_LOG( FG_GENERAL, FG_INFO,
	    "starting altitude is = " <<
	    fgGetDouble("/position/altitude") );

    f->set_Altitude( fgGetDouble("/position/altitude") );
    FG_LOG( FG_GENERAL, FG_INFO,
	    "Initial position is: ("
	    << (f->get_Longitude() * SGD_RADIANS_TO_DEGREES) << ", "
	    << (f->get_Latitude() * SGD_RADIANS_TO_DEGREES) << ", "
	    << (f->get_Altitude() * SG_FEET_TO_METER) << ")" );

    return true;
}


// General house keeping initializations
bool fgInitGeneral( void ) {
    string root;

#if defined(FX) && defined(XMESA)
    char *mesa_win_state;
#endif

    FG_LOG( FG_GENERAL, FG_INFO, "General Initialization" );
    FG_LOG( FG_GENERAL, FG_INFO, "======= ==============" );

    root = globals->get_fg_root();
    if ( ! root.length() ) {
	// No root path set? Then bail ...
	FG_LOG( FG_GENERAL, FG_ALERT,
		"Cannot continue without environment variable FG_ROOT"
		<< "being defined." );
	exit(-1);
    }
    FG_LOG( FG_GENERAL, FG_INFO, "FG_ROOT = " << '"' << root << '"' << endl );

#if defined(FX) && defined(XMESA)
    // initialize full screen flag
    global_fullscreen = false;
    if ( strstr ( general.get_glRenderer(), "Glide" ) ) {
	// Test for the MESA_GLX_FX env variable
	if ( (mesa_win_state = getenv( "MESA_GLX_FX" )) != NULL) {
	    // test if we are fullscreen mesa/glide
	    if ( (mesa_win_state[0] == 'f') ||
		 (mesa_win_state[0] == 'F') ) {
		global_fullscreen = true;
	    }
	}
    }
#endif

    return true;
}


// set initial aircraft speed
void
fgVelocityInit( void ) 
{
  if (!fgHasValue("/sim/startup/speed-set")) {
    current_aircraft.fdm_state->set_V_calibrated_kts(0.0);
    return;
  }

  const string speedset = fgGetString("/sim/startup/speed-set");
  if (speedset == "knots" || speedset == "KNOTS") {
    current_aircraft.fdm_state
      ->set_V_calibrated_kts(fgGetDouble("/velocities/airspeed"));
  } else if (speedset == "mach" || speedset == "MACH") {
    current_aircraft.fdm_state
      ->set_Mach_number(fgGetDouble("/velocities/mach"));
  } else if (speedset == "UVW" || speedset == "uvw") {
    current_aircraft.fdm_state
      ->set_Velocities_Wind_Body(fgGetDouble("/velocities/uBody"),
				 fgGetDouble("/velocities/vBody"),
				 fgGetDouble("/velocities/wBody"));
  } else if (speedset == "NED" || speedset == "ned") {
    current_aircraft.fdm_state
      ->set_Velocities_Local(fgGetDouble("/velocities/speed-north"),
			     fgGetDouble("/velocities/speed-east"),
			     fgGetDouble("/velocities/speed-down"));
  } else {
    FG_LOG(FG_GENERAL, FG_ALERT,
	   "Unrecognized value for /sim/startup/speed-set: " << speedset);
    current_aircraft.fdm_state->set_V_calibrated_kts(0.0);
  }
}             

        
// This is the top level init routine which calls all the other
// initialization routines.  If you are adding a subsystem to flight
// gear, its initialization call should located in this routine.
// Returns non-zero if a problem encountered.
bool fgInitSubsystems( void ) {
    fgLIGHT *l = &cur_light_params;

    FG_LOG( FG_GENERAL, FG_INFO, "Initialize Subsystems");
    FG_LOG( FG_GENERAL, FG_INFO, "========== ==========");

    // Initialize the material property lib
    FGPath mpath( globals->get_fg_root() );
    mpath.append( "materials" );
    if ( material_lib.load( mpath.str() ) ) {
    } else {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error loading material lib!" );
	exit(-1);
    }

    // Initialize the Scenery Management subsystem
    if ( fgSceneryInit() ) {
	// Material lib initialized ok.
    } else {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error in Scenery initialization!" );
	exit(-1);
    }

    if ( global_tile_mgr.init() ) {
	// Load the local scenery data
	global_tile_mgr.update( fgGetDouble("/position/longitude"),
				fgGetDouble("/position/latitude") );
    } else {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error in Tile Manager initialization!" );
	exit(-1);
    }

    FG_LOG( FG_GENERAL, FG_DEBUG,
    	    "Current terrain elevation after tile mgr init " <<
	    scenery.cur_elev );

    double dt = 1.0 / fgGetInt("/sim/model-hz");
    // cout << "dt = " << dt << endl;

    aircraft_dir = fgGetString("/sim/aircraft-dir");
    const string &model = fgGetString("/sim/flight-model");
    if (model == "larcsim") {
	cur_fdm_state = new FGLaRCsim( dt );
    } else if (model == "jsb") {
	cur_fdm_state = new FGJSBsim( dt );
    } else if (model == "ada") {
	cur_fdm_state = new FGADA( dt );
    } else if (model == "balloon") {
	cur_fdm_state = new FGBalloonSim( dt );
    } else if (model == "magic") {
	cur_fdm_state = new FGMagicCarpet( dt );
    } else if (model == "external") {
	cur_fdm_state = new FGExternal( dt );
    } else {
	FG_LOG(FG_GENERAL, FG_ALERT,
	       "Unrecognized flight model '" << model
	       << ", can't init aircraft");
	exit(-1);
    }
    cur_fdm_state->stamp();
    cur_fdm_state->set_remainder( 0 );

    // allocates structures so must happen before any of the flight
    // model or control parameters are set
    fgAircraftInit();   // In the future this might not be the case.

    fgFDMSetGroundElevation( fgGetString("/sim/flight-model"),
			     scenery.cur_elev );
    
    // set the initial position
    fgInitPosition();

    // Calculate ground elevation at starting point (we didn't have
    // tmp_abs_view_pos calculated when fgTileMgrUpdate() was called above
    //
    // calculalate a cartesian point somewhere along the line between
    // the center of the earth and our view position.  Doesn't have to
    // be the exact elevation (this is good because we don't know it
    // yet :-)

    // now handled inside of the fgTileMgrUpdate()

    // Reset our altitude if we are below ground
    FG_LOG( FG_GENERAL, FG_DEBUG, "Current altitude = "
	    << cur_fdm_state->get_Altitude() );
    FG_LOG( FG_GENERAL, FG_DEBUG, "Current runway altitude = " <<
	    cur_fdm_state->get_Runway_altitude() );

    if ( cur_fdm_state->get_Altitude() < cur_fdm_state->get_Runway_altitude() +
	 3.758099) {
	cur_fdm_state->set_Altitude( cur_fdm_state->get_Runway_altitude() +
				     3.758099 );
    }

    FG_LOG( FG_GENERAL, FG_INFO,
	    "Updated position (after elevation adj): ("
	    << (cur_fdm_state->get_Latitude() * SGD_RADIANS_TO_DEGREES) << ", "
	    << (cur_fdm_state->get_Longitude() * SGD_RADIANS_TO_DEGREES) << ", "
	    << (cur_fdm_state->get_Altitude() * SG_FEET_TO_METER) << ")" );

    // We need to calculate a few sea_level_radius here so we can pass
    // the correct value to the view class
    double sea_level_radius_meters;
    double lat_geoc;
    sgGeodToGeoc( cur_fdm_state->get_Latitude(),
		  cur_fdm_state->get_Altitude(),
		  &sea_level_radius_meters, &lat_geoc);
    cur_fdm_state->set_Sea_level_radius( sea_level_radius_meters *
					 SG_METER_TO_FEET );

    // The following section sets up the flight model EOM parameters
    // and should really be read in from one or more files.

    // Initial Velocity
    fgVelocityInit();

    // Initial Orientation
//     cur_fdm_state->
// 	set_Euler_Angles( fgGetDouble("/orientation/roll") * SGD_DEGREES_TO_RADIANS,
// 			  fgGetDouble("/orientation/pitch") * SGD_DEGREES_TO_RADIANS,
// 			  fgGetDouble("/orientation/heading") * SGD_DEGREES_TO_RADIANS );

    // Initialize the event manager
    global_events.Init();

    // Output event stats every 60 seconds
    global_events.Register( "fgEVENT_MGR::PrintStats()",
			    fgMethodCallback<fgEVENT_MGR>( &global_events,
						   &fgEVENT_MGR::PrintStats),
			    fgEVENT::FG_EVENT_READY, 60000 );

    // Initialize win_ratio parameters
    for ( int i = 0; i < globals->get_viewmgr()->size(); ++i ) {
	globals->get_viewmgr()->get_view(i)->
	    set_win_ratio( fgGetInt("/sim/startup/xsize") /
			   fgGetInt("/sim/startup/ysize") );
    }

    // Initialize pilot view
    FGViewerRPH *pilot_view =
	(FGViewerRPH *)globals->get_viewmgr()->get_view( 0 );

    pilot_view->set_geod_view_pos( cur_fdm_state->get_Longitude(), 
				   cur_fdm_state->get_Lat_geocentric(), 
				   cur_fdm_state->get_Altitude() *
				   SG_FEET_TO_METER );
    pilot_view->set_sea_level_radius( cur_fdm_state->get_Sea_level_radius() *
				      SG_FEET_TO_METER ); 
    pilot_view->set_rph( cur_fdm_state->get_Phi(),
			 cur_fdm_state->get_Theta(),
			 cur_fdm_state->get_Psi() );

    // set current view to 0 (first) which is our main pilot view
    globals->set_current_view( pilot_view );

    FG_LOG( FG_GENERAL, FG_DEBUG, "  abs_view_pos = "
	    << globals->get_current_view()->get_abs_view_pos());

    // fgUpdateSunPos() needs a few position and view parameters set
    // so it can calculate local relative sun angle and a few other
    // things for correctly orienting the sky.
    fgUpdateSunPos();
    fgUpdateMoonPos();
    global_events.Register( "fgUpdateSunPos()", fgUpdateSunPos,
			    fgEVENT::FG_EVENT_READY, 60000);
    global_events.Register( "fgUpdateMoonPos()", fgUpdateMoonPos,
			    fgEVENT::FG_EVENT_READY, 60000);

    // Initialize Lighting interpolation tables
    l->Init();

    // update the lighting parameters (based on sun angle)
    global_events.Register( "fgLight::Update()",
			    fgMethodCallback<fgLIGHT>( &cur_light_params,
						       &fgLIGHT::Update),
			    fgEVENT::FG_EVENT_READY, 30000 );
    // update the current timezone each 30 minutes
    global_events.Register( "fgUpdateLocalTime()", fgUpdateLocalTime,
			    fgEVENT::FG_EVENT_READY, 1800000);

    // Initialize the weather modeling subsystem
#ifndef FG_OLD_WEATHER
    // Initialize the WeatherDatabase
    FG_LOG(FG_GENERAL, FG_INFO, "Creating LocalWeatherDatabase");
    sgVec3 position;
    sgSetVec3( position, current_aircraft.fdm_state->get_Latitude(),
	       current_aircraft.fdm_state->get_Longitude(),
	       current_aircraft.fdm_state->get_Altitude() * SG_FEET_TO_METER );
    FGLocalWeatherDatabase::theFGLocalWeatherDatabase = 
	new FGLocalWeatherDatabase( position,
				    globals->get_fg_root() );
    // cout << theFGLocalWeatherDatabase << endl;
    // cout << "visibility = " 
    //      << theFGLocalWeatherDatabase->getWeatherVisibility() << endl;

    WeatherDatabase = FGLocalWeatherDatabase::theFGLocalWeatherDatabase;
    
    double init_vis = fgGetDouble("/environment/visibility");
    if ( init_vis > 0 ) {
	WeatherDatabase->setWeatherVisibility( init_vis );
    }

    // register the periodic update of the weather
    global_events.Register( "weather update", fgUpdateWeatherDatabase,
                            fgEVENT::FG_EVENT_READY, 30000);
#else
    current_weather.Init();
#endif

    // Initialize vor/ndb/ils/fix list management and query systems
    FG_LOG(FG_GENERAL, FG_INFO, "Loading Navaids");

    FG_LOG(FG_GENERAL, FG_INFO, "  VOR/NDB");
    current_navlist = new FGNavList;
    FGPath p_nav( globals->get_fg_root() );
    p_nav.append( "Navaids/default.nav" );
    current_navlist->init( p_nav );

    FG_LOG(FG_GENERAL, FG_INFO, "  ILS and Marker Beacons");
    current_beacons = new FGMarkerBeacons;
    current_beacons->init();
    current_ilslist = new FGILSList;
    FGPath p_ils( globals->get_fg_root() );
    p_ils.append( "Navaids/default.ils" );
    current_ilslist->init( p_ils );

    FG_LOG(FG_GENERAL, FG_INFO, "  Fixes");
    current_fixlist = new FGFixList;
    FGPath p_fix( globals->get_fg_root() );
    p_fix.append( "Navaids/default.fix" );
    current_fixlist->init( p_fix );

    // Radio stack subsystem.
    current_radiostack = new FGRadioStack;
    current_radiostack->init();
    current_radiostack->bind();

    // Initialize the Cockpit subsystem
    if( fgCockpitInit( &current_aircraft )) {
	// Cockpit initialized ok.
    } else {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error in Cockpit initialization!" );
	exit(-1);
    }

    // Initialize the flight model subsystem data structures base on
    // above values

    cur_fdm_state->init();
    cur_fdm_state->bind();
//     if ( cur_fdm_state->init( 1.0 / fgGetInt("/sim/model-hz") ) ) {
// 	// fdm init successful
//     } else {
// 	FG_LOG( FG_GENERAL, FG_ALERT, "FDM init() failed!  Cannot continue." );
// 	exit(-1);
//     }

    // *ABCD* I'm just sticking this here for now, it should probably
    // move eventually
    scenery.cur_elev = cur_fdm_state->get_Runway_altitude() * SG_FEET_TO_METER;

    if ( cur_fdm_state->get_Altitude() <
	 cur_fdm_state->get_Runway_altitude() + 3.758099)
    {
	cur_fdm_state->set_Altitude( cur_fdm_state->get_Runway_altitude() +
				     3.758099 );
    }

    FG_LOG( FG_GENERAL, FG_INFO,
	    "Updated position (after elevation adj): ("
	    << (cur_fdm_state->get_Latitude() * SGD_RADIANS_TO_DEGREES) << ", "
	    << (cur_fdm_state->get_Longitude() * SGD_RADIANS_TO_DEGREES) << ", "
	    << (cur_fdm_state->get_Altitude() * SG_FEET_TO_METER) << ")" );
    // *ABCD* end of thing that I just stuck in that I should probably
    // move

    // Joystick support
    if ( ! fgJoystickInit() ) {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error in Joystick initialization!" );
    }

    // Autopilot init
    current_autopilot = new FGAutopilot;
    current_autopilot->init();

    // initialize the gui parts of the autopilot
    NewTgtAirportInit();
    fgAPAdjustInit() ;
    NewHeadingInit();
    NewAltitudeInit();

    // Initialize I/O channels
#if ! defined( macintosh )
    fgIOInit();
#endif

    // Initialize the 2D panel.
    string panel_path = fgGetString("/sim/panel/path",
				    "Panels/Default/default.xml");
    current_panel = fgReadPanel(panel_path);
    if (current_panel == 0) {
	FG_LOG( FG_INPUT, FG_ALERT, 
		"Error reading new panel from " << panel_path );
    } else {
	FG_LOG( FG_INPUT, FG_INFO, "Loaded new panel from " << panel_path );
	current_panel->init();
	current_panel->bind();
    }

    // Initialize the BFI
    FGBFI::init();

    controls.init();
    controls.bind();

    FG_LOG( FG_GENERAL, FG_INFO, endl);

				// Save the initial state for future
				// reference.
    globals->saveInitialState();

    return true;
}


void fgReInitSubsystems( void )
{
    FG_LOG( FG_GENERAL, FG_INFO,
	    "/position/altitude = " << fgGetDouble("/position/altitude") );

    bool freeze = globals->get_freeze();
    if( !freeze )
        globals->set_freeze( true );
    
    // Initialize the Scenery Management subsystem
    if ( ! fgSceneryInit() ) {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error in Scenery initialization!" );
	exit(-1);
    }

    if( global_tile_mgr.init() ) {
	// Load the local scenery data
	global_tile_mgr.update( fgGetDouble("/position/longitude"),
				fgGetDouble("/position/latitude") );
    } else {
    	FG_LOG( FG_GENERAL, FG_ALERT, "Error in Tile Manager initialization!" );
		exit(-1);
    }

    // cout << "current scenery elev = " << scenery.cur_elev << endl;

    fgFDMSetGroundElevation( fgGetString("/sim/flight-model"), 
			     scenery.cur_elev );
    fgInitPosition();

    // Reset our altitude if we are below ground
    FG_LOG( FG_GENERAL, FG_DEBUG, "Current altitude = "
	    << cur_fdm_state->get_Altitude() );
    FG_LOG( FG_GENERAL, FG_DEBUG, "Current runway altitude = "
	    << cur_fdm_state->get_Runway_altitude() );

    if ( cur_fdm_state->get_Altitude() <
	 cur_fdm_state->get_Runway_altitude() + 3.758099)
    {
	cur_fdm_state->set_Altitude( cur_fdm_state->get_Runway_altitude() +
				     3.758099 );
    }
    double sea_level_radius_meters;
    double lat_geoc;
    sgGeodToGeoc( cur_fdm_state->get_Latitude(), cur_fdm_state->get_Altitude(), 
		  &sea_level_radius_meters, &lat_geoc);
    cur_fdm_state->set_Sea_level_radius( sea_level_radius_meters *
					 SG_METER_TO_FEET );
	
    // The following section sets up the flight model EOM parameters
    // and should really be read in from one or more files.

    // Initial Velocity
    fgVelocityInit();

    // Initial Orientation
//     cur_fdm_state->
// 	set_Euler_Angles( fgGetDouble("/orientation/roll") * SGD_DEGREES_TO_RADIANS,
// 			  fgGetDouble("/orientation/pitch") * SGD_DEGREES_TO_RADIANS,
// 			  fgGetDouble("/orientation/heading") * SGD_DEGREES_TO_RADIANS );

    // Initialize view parameters
    FGViewerRPH *pilot_view =
	(FGViewerRPH *)globals->get_viewmgr()->get_view( 0 );

    pilot_view->set_view_offset( 0.0 );
    pilot_view->set_goal_view_offset( 0.0 );

    pilot_view->set_geod_view_pos( cur_fdm_state->get_Longitude(), 
				   cur_fdm_state->get_Lat_geocentric(), 
				   cur_fdm_state->get_Altitude() *
				   SG_FEET_TO_METER );
    pilot_view->set_sea_level_radius( cur_fdm_state->get_Sea_level_radius() *
				      SG_FEET_TO_METER ); 
    pilot_view->set_rph( cur_fdm_state->get_Phi(),
			 cur_fdm_state->get_Theta(),
			 cur_fdm_state->get_Psi() );

    // set current view to 0 (first) which is our main pilot view
    globals->set_current_view( pilot_view );

    FG_LOG( FG_GENERAL, FG_DEBUG, "  abs_view_pos = "
	    << globals->get_current_view()->get_abs_view_pos());

    cur_fdm_state->init();
//     cur_fdm_state->bind();
//     cur_fdm_state->init( 1.0 / fgGetInt("/sim/model-hz") );

    scenery.cur_elev = cur_fdm_state->get_Runway_altitude() * SG_FEET_TO_METER;

    if ( cur_fdm_state->get_Altitude() <
	 cur_fdm_state->get_Runway_altitude() + 3.758099)
    {
	cur_fdm_state->set_Altitude( cur_fdm_state->get_Runway_altitude() +
				     3.758099 );
    }

    controls.reset_all();
    current_autopilot->reset();

    if( !freeze )
        globals->set_freeze( false );
}