src/Systems/electrical.cxx

   1 // electrical.cxx - a flexible, generic electrical system model.
   2 //
   3 // Written by Curtis Olson, started September 2002.
   4 //
   5 // Copyright (C) 2002  Curtis L. Olson  - http://www.flightgear.org/~curt
   6 //
   7 // This program is free software; you can redistribute it and/or
   8 // modify it under the terms of the GNU General Public License as
   9 // published by the Free Software Foundation; either version 2 of the
  10 // License, or (at your option) any later version.
  11 //
  12 // This program is distributed in the hope that it will be useful, but
  13 // WITHOUT ANY WARRANTY; without even the implied warranty of
  14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 // General Public License for more details.
  16 //
  17 // You should have received a copy of the GNU General Public License
  18 // along with this program; if not, write to the Free Software
  19 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  20 //
  21 // $Id$
  22
  23 #ifdef HAVE_CONFIG_H
  24 #  include <config.h>
  25 #endif
  26
  27 #include <cstring>
  28
  29 #include <simgear/structure/exception.hxx>
  30 #include <simgear/misc/sg_path.hxx>
  31 #include <simgear/debug/logstream.hxx>
  32 #include <simgear/props/props_io.hxx>
  33
  34 #include <Main/fg_props.hxx>
  35 #include <Main/globals.hxx>
  36
  37 #include "electrical.hxx"
  38
  39
  40 FGElectricalComponent::FGElectricalComponent() :
  41     kind(-1),
  42     name(""),
  43     volts(0.0),
  44     load_amps(0.0)
  45 {
  46 }
  47
  48
  49 FGElectricalSupplier::FGElectricalSupplier ( SGPropertyNode *node ) {
  50     kind = FG_SUPPLIER;
  51
  52     // cout << "Creating a supplier" << endl;
  53     name = node->getStringValue("name");
  54     string _model = node->getStringValue("kind");
  55     // cout << "_model = " << _model << endl;
  56     if ( _model == "battery" ) {
  57         model = FG_BATTERY;
  58         amp_hours = node->getFloatValue("amp-hours", 40.0);
  59         percent_remaining = node->getFloatValue("percent-remaining", 1.0);
  60         charge_amps = node->getFloatValue("charge-amps", 7.0);
  61     } else if ( _model == "alternator" ) {
  62         model = FG_ALTERNATOR;
  63         rpm_src = node->getStringValue("rpm-source");
  64         rpm_threshold = node->getFloatValue("rpm-threshold", 600.0);
  65         ideal_amps = node->getFloatValue("amps", 60.0);
  66     } else if ( _model == "external" ) {
  67         model = FG_EXTERNAL;
  68         ideal_amps = node->getFloatValue("amps", 60.0);
  69     } else {
  70         model = FG_UNKNOWN;
  71     }
  72     ideal_volts = node->getFloatValue("volts");
  73
  74     int i;
  75     for ( i = 0; i < node->nChildren(); ++i ) {
  76         SGPropertyNode *child = node->getChild(i);
  77         // cout << " scanning: " << child->getName() << endl;
  78         if ( !strcmp(child->getName(), "prop") ) {
  79             string prop = child->getStringValue();
  80             // cout << "  Adding prop = " << prop << endl;
  81             add_prop( prop );
  82             fgSetFloat( prop.c_str(), ideal_amps );
  83         }
  84     }
  85
  86     _rpm_node = fgGetNode( rpm_src.c_str(), true);
  87 }
  88
  89
  90 float FGElectricalSupplier::apply_load( float amps, float dt ) {
  91     if ( model == FG_BATTERY ) {
  92         // calculate amp hours used
  93         float amphrs_used = amps * dt / 3600.0;
  94
  95         // calculate percent of total available capacity
  96         float percent_used = amphrs_used / amp_hours;
  97         percent_remaining -= percent_used;
  98         if ( percent_remaining < 0.0 ) {
  99             percent_remaining = 0.0;
 100         } else if ( percent_remaining > 1.0 ) {
 101             percent_remaining = 1.0;
 102         }
 103         // cout << "battery percent = " << percent_remaining << endl;
 104         return amp_hours * percent_remaining;
 105     } else if ( model == FG_ALTERNATOR ) {
 106         // scale alternator output for rpms < 600.  For rpms >= 600
 107         // give full output.  This is just a WAG, and probably not how
 108         // it really works but I'm keeping things "simple" to start.
 109         float rpm = _rpm_node->getFloatValue();
 110         float factor = rpm / rpm_threshold;
 111         if ( factor > 1.0 ) {
 112             factor = 1.0;
 113         }
 114         // cout << "alternator amps = " << amps * factor << endl;
 115         float available_amps = ideal_amps * factor;
 116         return available_amps - amps;
 117     } else if ( model == FG_EXTERNAL ) {
 118         // cout << "external amps = " << 0.0 << endl;
 119         float available_amps = ideal_amps;
 120         return available_amps - amps;
 121     } else {
 122         SG_LOG( SG_ALL, SG_ALERT, "unknown supplier type" );
 123     }
 124
 125     return 0.0;
 126 }
 127
 128
 129 float FGElectricalSupplier::get_output_volts() {
 130     if ( model == FG_BATTERY ) {
 131         // cout << "battery amps = " << amps << endl;
 132         float x = 1.0 - percent_remaining;
 133         float tmp = -(3.0 * x - 1.0);
 134         float factor = (tmp*tmp*tmp*tmp*tmp + 32) / 32;
 135         // cout << "battery % = " << percent_remaining <<
 136         //         " factor = " << factor << endl;
 137         // percent_remaining -= 0.001;
 138         return ideal_volts * factor;
 139     } else if ( model == FG_ALTERNATOR ) {
 140         // scale alternator output for rpms < 600.  For rpms >= 600
 141         // give full output.  This is just a WAG, and probably not how
 142         // it really works but I'm keeping things "simple" to start.
 143         float rpm = _rpm_node->getFloatValue();
 144         float factor = rpm / rpm_threshold;
 145         if ( factor > 1.0 ) {
 146             factor = 1.0;
 147         }
 148         // cout << "alternator amps = " << amps * factor << endl;
 149         return ideal_volts * factor;
 150     } else if ( model == FG_EXTERNAL ) {
 151         // cout << "external amps = " << 0.0 << endl;
 152         return ideal_volts;
 153     } else {
 154         SG_LOG( SG_ALL, SG_ALERT, "unknown supplier type" );
 155     }
 156
 157     return 0.0;
 158 }
 159
 160
 161 float FGElectricalSupplier::get_output_amps() {
 162     if ( model == FG_BATTERY ) {
 163         // cout << "battery amp_hours = " << amp_hours << endl;
 164
 165         // This is a WAG, but produce enough amps to burn the entire
 166         // battery in one minute.
 167         return amp_hours * 60.0;
 168     } else if ( model == FG_ALTERNATOR ) {
 169         // scale alternator output for rpms < 600.  For rpms >= 600
 170         // give full output.  This is just a WAG, and probably not how
 171         // it really works but I'm keeping things "simple" to start.
 172         float rpm = _rpm_node->getFloatValue();
 173         float factor = rpm / rpm_threshold;
 174         if ( factor > 1.0 ) {
 175             factor = 1.0;
 176         }
 177         // cout << "alternator amps = " << ideal_amps * factor << endl;
 178         return ideal_amps * factor;
 179     } else if ( model == FG_EXTERNAL ) {
 180         // cout << "external amps = " << 0.0 << endl;
 181         return ideal_amps;
 182     } else {
 183         SG_LOG( SG_ALL, SG_ALERT, "unknown supplier type" );
 184     }
 185
 186     return 0.0;
 187 }
 188
 189
 190 FGElectricalBus::FGElectricalBus ( SGPropertyNode *node ) {
 191     kind = FG_BUS;
 192
 193     name = node->getStringValue("name");
 194     int i;
 195     for ( i = 0; i < node->nChildren(); ++i ) {
 196         SGPropertyNode *child = node->getChild(i);
 197         if ( !strcmp(child->getName(), "prop") ) {
 198             string prop = child->getStringValue();
 199             add_prop( prop );
 200         }
 201     }
 202 }
 203
 204
 205 FGElectricalOutput::FGElectricalOutput ( SGPropertyNode *node ) {
 206     kind = FG_OUTPUT;
 207     load_amps = 0.1;            // arbitrary default value
 208
 209     name = node->getStringValue("name");
 210     SGPropertyNode *draw = node->getNode("rated-draw");
 211     if ( draw != NULL ) {
 212         load_amps = draw->getFloatValue();
 213     }
 214     // cout << "rated draw = " << output_amps << endl;
 215
 216     int i;
 217     for ( i = 0; i < node->nChildren(); ++i ) {
 218         SGPropertyNode *child = node->getChild(i);
 219         if ( !strcmp(child->getName(), "prop") ) {
 220             string prop = child->getStringValue();
 221             add_prop( prop );
 222         }
 223     }
 224 }
 225
 226
 227 FGElectricalSwitch::FGElectricalSwitch( SGPropertyNode *node ) :
 228     switch_node( NULL ),
 229     rating_amps( 0.0f ),
 230     circuit_breaker( false )
 231 {
 232     bool initial_state = true;
 233     int i;
 234     for ( i = 0; i < node->nChildren(); ++i ) {
 235         SGPropertyNode *child = node->getChild(i);
 236         string cname = child->getName();
 237         string cval = child->getStringValue();
 238         if ( cname == "prop" ) {
 239             switch_node = fgGetNode( cval.c_str(), true );
 240             // cout << "switch node = " << cval << endl;
 241         } else if ( cname == "initial-state" ) {
 242             if ( cval == "off" || cval == "false" ) {
 243                 initial_state = false;
 244             }
 245             // cout << "initial state = " << initial_state << endl;
 246         } else if ( cname == "rating-amps" ) {
 247             rating_amps = atof( cval.c_str() );
 248             circuit_breaker = true;
 249             // cout << "initial state = " << initial_state << endl;
 250         }
 251     }
 252
 253     switch_node->setBoolValue( initial_state );
 254     // cout << "  value = " << switch_node->getBoolValue() << endl;
 255 }
 256
 257
 258 FGElectricalConnector::FGElectricalConnector ( SGPropertyNode *node,
 259                                                FGElectricalSystem *es ) {
 260     kind = FG_CONNECTOR;
 261     name = "connector";
 262     int i;
 263     for ( i = 0; i < node->nChildren(); ++i ) {
 264         SGPropertyNode *child = node->getChild(i);
 265         string cname = child->getName();
 266         string cval = child->getStringValue();
 267         // cout << "  " << cname << " = " << cval << endl;
 268         if ( cname == "input" ) {
 269             FGElectricalComponent *s = es->find( child->getStringValue() );
 270             if ( s != NULL ) {
 271                 add_input( s );
 272                 if ( s->get_kind() == FG_SUPPLIER ) {
 273                     s->add_output( this );
 274                 } else if ( s->get_kind() == FG_BUS ) {
 275                     s->add_output( this );
 276                 } else {
 277                     SG_LOG( SG_ALL, SG_ALERT,
 278                             "Attempt to connect to something that can't provide an output: "
 279                             << child->getStringValue() );
 280                 }
 281             } else {
 282                 SG_LOG( SG_ALL, SG_ALERT, "Can't find named source: "
 283                         << child->getStringValue() );
 284             }
 285         } else if ( cname == "output" ) {
 286             FGElectricalComponent *s = es->find( child->getStringValue() );
 287             if ( s != NULL ) {
 288                 add_output( s );
 289                 if ( s->get_kind() == FG_BUS ) {
 290                     s->add_input( this );
 291                 } else if ( s->get_kind() == FG_OUTPUT ) {
 292                     s->add_input( this );
 293                 } else if ( s->get_kind() == FG_SUPPLIER &&
 294                             ((FGElectricalSupplier *)s)->get_model()
 295                             == FGElectricalSupplier::FG_BATTERY ) {
 296                     s->add_output( this );
 297                 } else {
 298                     SG_LOG( SG_ALL, SG_ALERT,
 299                             "Attempt to connect to something that can't provide an input: "
 300                             << child->getStringValue() );
 301                 }
 302             } else {
 303                 SG_LOG( SG_ALL, SG_ALERT, "Can't find named source: "
 304                         << child->getStringValue() );
 305             }
 306         } else if ( cname == "switch" ) {
 307              // cout << "Switch = " << child->getStringValue() << endl;
 308             FGElectricalSwitch s( child );
 309             add_switch( s );
 310         }
 311     }
 312 }
 313
 314
 315 // set all switches to the specified state
 316 void FGElectricalConnector::set_switches( bool state ) {
 317     // cout << "setting switch state to " << state << endl;
 318     for ( unsigned int i = 0; i < switches.size(); ++i ) {
 319         switches[i].set_state( state );
 320     }
 321 }
 322
 323
 324 // return true if all switches are true, false otherwise.  A connector
 325 // could have multiple switches, but they all need to be true(closed)
 326 // for current to get through.
 327 bool FGElectricalConnector::get_state() {
 328     unsigned int i;
 329     for ( i = 0; i < switches.size(); ++i ) {
 330         if ( ! switches[i].get_state() ) {
 331             return false;
 332         }
 333     }
 334
 335     return true;
 336 }
 337
 338
 339 FGElectricalSystem::FGElectricalSystem ( SGPropertyNode *node ) :
 340     name(node->getStringValue("name", "electrical")),
 341     num(node->getIntValue("number", 0)),
 342     path(node->getStringValue("path")),
 343     enabled(false)
 344 {
 345 }
 346
 347
 348 FGElectricalSystem::~FGElectricalSystem () {
 349 }
 350
 351
 352 void FGElectricalSystem::init () {
 353     SGPropertyNode_ptr config_props = new SGPropertyNode;
 354
 355     _volts_out = fgGetNode( "/systems/electrical/volts", true );
 356     _amps_out = fgGetNode( "/systems/electrical/amps", true );
 357
 358     // allow the electrical system to be specified via the
 359     // aircraft-set.xml file (for backwards compatibility) or through
 360     // the aircraft-systems.xml file.  If a -set.xml entry is
 361     // specified, that overrides the system entry.
 362     SGPropertyNode *path_n = fgGetNode("/sim/systems/electrical/path");
 363     if ( path_n ) {
 364         if ( path.length() ) {
 365             SG_LOG( SG_ALL, SG_INFO,
 366                     "NOTICE: System manager configuration specifies an " <<
 367                     "electrical system: " << path << " but it is " <<
 368                     "being overridden by the one specified in the -set.xml " <<
 369                     "file: " << path_n->getStringValue() );
 370         }
 371
 372         path = path_n->getStringValue();
 373     }
 374
 375     if ( path.length() ) {
 376         SGPath config( globals->get_fg_root() );
 377         config.append( path );
 378
 379         // load an obsolete xml configuration
 380         SG_LOG( SG_ALL, SG_WARN,
 381                 "Reading deprecated xml electrical system model from\n    "
 382                 << config.str() );
 383         try {
 384             readProperties( config.str(), config_props );
 385
 386             if ( build(config_props) ) {
 387                 enabled = true;
 388             } else {
 389                 SG_LOG( SG_ALL, SG_ALERT,
 390                         "Detected a logic error in the electrical system ");
 391                 SG_LOG( SG_ALL, SG_ALERT,
 392                         "specification file.  See earlier errors for " );
 393                 SG_LOG( SG_ALL, SG_ALERT,
 394                         "details.");
 395                 exit(-1);
 396             }
 397         } catch (const sg_exception&) {
 398             SG_LOG( SG_ALL, SG_ALERT,
 399                     "Failed to load electrical system model: "
 400                     << config.str() );
 401         }
 402     } else {
 403         SG_LOG( SG_ALL, SG_INFO,
 404                 "No xml-based electrical model specified for this model!");
 405     }
 406
 407     if ( !enabled ) {
 408         _amps_out->setDoubleValue(0);
 409     }
 410
 411 }
 412
 413
 414 void FGElectricalSystem::bind () {
 415 }
 416
 417
 418 void FGElectricalSystem::unbind () {
 419 }
 420
 421
 422 void FGElectricalSystem::update (double dt) {
 423     if ( !enabled ) {
 424         return;
 425     }
 426
 427     // cout << "Updating electrical system, dt = " << dt << endl;
 428
 429     unsigned int i;
 430
 431     // zero out the voltage before we start, but don't clear the
 432     // requested load values.
 433     for ( i = 0; i < suppliers.size(); ++i ) {
 434         suppliers[i]->set_volts( 0.0 );
 435     }
 436     for ( i = 0; i < buses.size(); ++i ) {
 437         buses[i]->set_volts( 0.0 );
 438     }
 439     for ( i = 0; i < outputs.size(); ++i ) {
 440         outputs[i]->set_volts( 0.0 );
 441     }
 442     for ( i = 0; i < connectors.size(); ++i ) {
 443         connectors[i]->set_volts( 0.0 );
 444     }
 445
 446     // for each "external" supplier, propagate the electrical current
 447     for ( i = 0; i < suppliers.size(); ++i ) {
 448         FGElectricalSupplier *node = (FGElectricalSupplier *)suppliers[i];
 449         if ( node->get_model() == FGElectricalSupplier::FG_EXTERNAL ) {
 450             float load;
 451             // cout << "Starting propagation: " << suppliers[i]->get_name()
 452             //      << endl;
 453             load = propagate( suppliers[i], dt,
 454                               node->get_output_volts(),
 455                               node->get_output_amps(),
 456                               " " );
 457
 458             if ( node->apply_load( load, dt ) < 0.0 ) {
 459                 SG_LOG(SG_ALL, SG_ALERT,
 460                        "Error drawing more current than available!");
 461             }
 462         }
 463     }
 464
 465     // for each "alternator" supplier, propagate the electrical
 466     // current
 467     for ( i = 0; i < suppliers.size(); ++i ) {
 468         FGElectricalSupplier *node = (FGElectricalSupplier *)suppliers[i];
 469         if ( node->get_model() == FGElectricalSupplier::FG_ALTERNATOR) {
 470             float load;
 471             // cout << "Starting propagation: " << suppliers[i]->get_name()
 472             //      << endl;
 473             load = propagate( suppliers[i], dt,
 474                               node->get_output_volts(),
 475                               node->get_output_amps(),
 476                               " " );
 477
 478             if ( node->apply_load( load, dt ) < 0.0 ) {
 479                 SG_LOG(SG_ALL, SG_ALERT,
 480                        "Error drawing more current than available!");
 481             }
 482         }
 483     }
 484
 485     // for each "battery" supplier, propagate the electrical
 486     // current
 487     for ( i = 0; i < suppliers.size(); ++i ) {
 488         FGElectricalSupplier *node = (FGElectricalSupplier *)suppliers[i];
 489         if ( node->get_model() == FGElectricalSupplier::FG_BATTERY ) {
 490             float load;
 491             // cout << "Starting propagation: " << suppliers[i]->get_name()
 492             //      << endl;
 493             load = propagate( suppliers[i], dt,
 494                               node->get_output_volts(),
 495                               node->get_output_amps(),
 496                               " " );
 497             // cout << "battery load = " << load << endl;
 498
 499             if ( node->apply_load( load, dt ) < 0.0 ) {
 500                 SG_LOG(SG_ALL, SG_ALERT,
 501                        "Error drawing more current than available!");
 502             }
 503         }
 504     }
 505
 506     float alt_norm
 507         = fgGetFloat("/systems/electrical/suppliers/alternator") / 60.0;
 508
 509     // impliment an extremely simplistic voltage model (assumes
 510     // certain naming conventions in electrical system config)
 511     // FIXME: we probably want to be able to feed power from all
 512     // engines if they are running and the master-alt is switched on
 513     float volts = 0.0;
 514     if ( fgGetBool("/controls/engines/engine[0]/master-bat") ) {
 515         volts = 24.0;
 516     }
 517     if ( fgGetBool("/controls/engines/engine[0]/master-alt") ) {
 518         if ( fgGetFloat("/engines/engine[0]/rpm") > 800 ) {
 519             float alt_contrib = 28.0;
 520             if ( alt_contrib > volts ) {
 521                 volts = alt_contrib;
 522             }
 523         } else if ( fgGetFloat("/engines/engine[0]/rpm") > 200 ) {
 524             float alt_contrib = 20.0;
 525             if ( alt_contrib > volts ) {
 526                 volts = alt_contrib;
 527             }
 528         }
 529     }
 530     _volts_out->setFloatValue( volts );
 531
 532     // impliment an extremely simplistic amps model (assumes certain
 533     // naming conventions in the electrical system config) ... FIXME:
 534     // make this more generic
 535     float amps = 0.0;
 536     if ( fgGetBool("/controls/engines/engine[0]/master-bat") ) {
 537         if ( fgGetBool("/controls/engines/engine[0]/master-alt") &&
 538              fgGetFloat("/engines/engine[0]/rpm") > 800 )
 539         {
 540             amps += 40.0 * alt_norm;
 541         }
 542         amps -= 15.0;            // normal load
 543         if ( fgGetBool("/controls/switches/flashing-beacon") ) {
 544             amps -= 7.5;
 545         }
 546         if ( fgGetBool("/controls/switches/nav-lights") ) {
 547             amps -= 7.5;
 548         }
 549         if ( amps > 7.0 ) {
 550             amps = 7.0;
 551         }
 552     }
 553     _amps_out->setFloatValue( amps );
 554 }
 555
 556
 557 bool FGElectricalSystem::build (SGPropertyNode* config_props) {
 558     SGPropertyNode *node;
 559     int i;
 560
 561     int count = config_props->nChildren();
 562     for ( i = 0; i < count; ++i ) {
 563         node = config_props->getChild(i);
 564         string name = node->getName();
 565         // cout << name << endl;
 566         if ( name == "supplier" ) {
 567             FGElectricalSupplier *s =
 568                 new FGElectricalSupplier( node );
 569             suppliers.push_back( s );
 570         } else if ( name == "bus" ) {
 571             FGElectricalBus *b =
 572                 new FGElectricalBus( node );
 573             buses.push_back( b );
 574         } else if ( name == "output" ) {
 575             FGElectricalOutput *o =
 576                 new FGElectricalOutput( node );
 577             outputs.push_back( o );
 578         } else if ( name == "connector" ) {
 579             FGElectricalConnector *c =
 580                 new FGElectricalConnector( node, this );
 581             connectors.push_back( c );
 582         } else {
 583             SG_LOG( SG_ALL, SG_ALERT, "Unknown component type specified: "
 584                     << name );
 585             return false;
 586         }
 587     }
 588
 589     return true;
 590 }
 591
 592
 593 // propagate the electrical current through the network, returns the
 594 // total current drawn by the children of this node.
 595 float FGElectricalSystem::propagate( FGElectricalComponent *node, double dt,
 596                                      float input_volts, float input_amps,
 597                                      string s ) {
 598     s += " ";
 599
 600     float total_load = 0.0;
 601
 602     // determine the current to carry forward
 603     float volts = 0.0;
 604     if ( !fgGetBool("/systems/electrical/serviceable") ) {
 605         volts = 0;
 606     } else if ( node->get_kind() == FGElectricalComponent::FG_SUPPLIER ) {
 607         // cout << s << "is a supplier (" << node->get_name() << ")" << endl;
 608         FGElectricalSupplier *supplier = (FGElectricalSupplier *)node;
 609         if ( supplier->get_model() == FGElectricalSupplier::FG_BATTERY ) {
 610             // cout << s << " (and is a battery)" << endl;
 611             float battery_volts = supplier->get_output_volts();
 612             if ( battery_volts < (input_volts - 0.1) ) {
 613                 // special handling of a battery charge condition
 614                 // cout << s << "  (and is being charged) in v = "
 615                 //      << input_volts << " current v = " << battery_volts
 616                 //      << endl;
 617                 supplier->apply_load( -supplier->get_charge_amps(), dt );
 618                 return supplier->get_charge_amps();
 619             }
 620         }
 621         volts = input_volts;
 622     } else if ( node->get_kind() == FGElectricalComponent::FG_BUS ) {
 623         // cout << s << "is a bus (" << node->get_name() << ")" << endl;
 624         volts = input_volts;
 625     } else if ( node->get_kind() == FGElectricalComponent::FG_OUTPUT ) {
 626         // cout << s << "is an output (" << node->get_name() << ")" << endl;
 627         volts = input_volts;
 628         if ( volts > 1.0 ) {
 629             // draw current if we have voltage
 630             total_load = node->get_load_amps();
 631         }
 632     } else if ( node->get_kind() == FGElectricalComponent::FG_CONNECTOR ) {
 633         // cout << s << "is a connector (" << node->get_name() << ")" << endl;
 634         if ( ((FGElectricalConnector *)node)->get_state() ) {
 635             volts = input_volts;
 636         } else {
 637             volts = 0.0;
 638         }
 639         // cout << s << "  input_volts = " << volts << endl;
 640     } else {
 641         SG_LOG( SG_ALL, SG_ALERT, "unkown node type" );
 642     }
 643
 644     int i;
 645
 646     // if this node has found a stronger power source, update the
 647     // value and propagate to all children
 648     if ( volts > node->get_volts() ) {
 649         node->set_volts( volts );
 650         for ( i = 0; i < node->get_num_outputs(); ++i ) {
 651             FGElectricalComponent *child = node->get_output(i);
 652             // send current equal to load
 653             total_load += propagate( child, dt,
 654                                      volts, child->get_load_amps(),
 655                                      s );
 656         }
 657
 658         // if not an output node, register the downstream current draw
 659         // (sum of all children) with this node.  If volts are zero,
 660         // current draw should be zero.
 661         if ( node->get_kind() != FGElectricalComponent::FG_OUTPUT ) {
 662             node->set_load_amps( total_load );
 663         }
 664
 665         node->set_available_amps( input_amps - total_load );
 666
 667         // publish values to specified properties
 668         for ( i = 0; i < node->get_num_props(); ++i ) {
 669             fgSetFloat( node->get_prop(i).c_str(), node->get_volts() );
 670         }
 671
 672         /*
 673         cout << s << node->get_name() << " -> (volts) " << node->get_volts()
 674              << endl;
 675         cout << s << node->get_name() << " -> (load amps) " << total_load
 676              << endl;
 677         cout << s << node->get_name() << " -> (input amps) " << input_amps
 678              << endl;
 679         cout << s << node->get_name() << " -> (extra amps) "
 680              << node->get_available_amps() << endl;
 681         */
 682
 683         return total_load;
 684     } else {
 685         // cout << s << "no further propagation" << endl;
 686         return 0.0;
 687     }
 688 }
 689
 690
 691 // search for the named component and return a pointer to it, NULL otherwise
 692 FGElectricalComponent *FGElectricalSystem::find ( const string &name ) {
 693     unsigned int i;
 694     string s;
 695
 696     // search suppliers
 697     for ( i = 0; i < suppliers.size(); ++i ) {
 698         s = suppliers[i]->get_name();
 699         // cout <<  "    " << s << endl;
 700         if ( s == name ) {
 701             return suppliers[i];
 702         }
 703     }
 704
 705     // then search buses
 706     for ( i = 0; i < buses.size(); ++i ) {
 707         s = buses[i]->get_name();
 708         // cout <<  "    " << s << endl;
 709         if ( s == name ) {
 710             return buses[i];
 711         }
 712     }
 713
 714     // then search outputs
 715     for ( i = 0; i < outputs.size(); ++i ) {
 716         s = outputs[i]->get_name();
 717         // cout <<  "    " << s << endl;
 718         if ( s == name ) {
 719             return outputs[i];
 720         }
 721     }
 722
 723     // nothing found
 724     return NULL;
 725 }