src/FDM/JSBSim/models/FGPropulsion.cpp

   1 /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   2
   3  Module:       FGPropulsion.cpp
   4  Author:       Jon S. Berndt
   5  Date started: 08/20/00
   6  Purpose:      Encapsulates the set of engines and tanks associated
   7                with this aircraft
   8
   9  ------------- Copyright (C) 2000  Jon S. Berndt (jsb@hal-pc.org) -------------
  10
  11  This program is free software; you can redistribute it and/or modify it under
  12  the terms of the GNU Lesser General Public License as published by the Free Software
  13  Foundation; either version 2 of the License, or (at your option) any later
  14  version.
  15
  16  This program is distributed in the hope that it will be useful, but WITHOUT
  17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  18  FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
  19  details.
  20
  21  You should have received a copy of the GNU Lesser General Public License along with
  22  this program; if not, write to the Free Software Foundation, Inc., 59 Temple
  23  Place - Suite 330, Boston, MA  02111-1307, USA.
  24
  25  Further information about the GNU Lesser General Public License can also be found on
  26  the world wide web at http://www.gnu.org.
  27
  28 FUNCTIONAL DESCRIPTION
  29 --------------------------------------------------------------------------------
  30 The Propulsion class is the container for the entire propulsion system, which is
  31 comprised of engines and tanks. Once the Propulsion class gets the config file,
  32 it reads in information which is specific to a type of engine. Then:
  33
  34 1) The appropriate engine type instance is created
  35 2) At least one tank object is created, and is linked to an engine.
  36
  37 At Run time each engines Calculate() method is called.
  38
  39 HISTORY
  40 --------------------------------------------------------------------------------
  41 08/20/00   JSB   Created
  42
  43 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  44 INCLUDES
  45 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
  46
  47 #include "FGPropulsion.h"
  48 #include <models/propulsion/FGRocket.h>
  49 #include <models/propulsion/FGTurbine.h>
  50 #include <models/propulsion/FGPiston.h>
  51 #include <models/propulsion/FGElectric.h>
  52 #include <models/propulsion/FGTurboProp.h>
  53 #include <input_output/FGPropertyManager.h>
  54 #include <input_output/FGXMLParse.h>
  55 #include <math/FGColumnVector3.h>
  56 #include <sstream>
  57
  58 namespace JSBSim {
  59
  60 static const char *IdSrc = "$Id$";
  61 static const char *IdHdr = ID_PROPULSION;
  62
  63 extern short debug_lvl;
  64
  65
  66 /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  67 CLASS IMPLEMENTATION
  68 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
  69
  70 FGPropulsion::FGPropulsion(FGFDMExec* exec) : FGModel(exec)
  71 {
  72   Name = "FGPropulsion";
  73
  74   InitializedEngines = 0;
  75   numSelectedFuelTanks = numSelectedOxiTanks = 0;
  76   numTanks = numEngines = 0;
  77   numOxiTanks = numFuelTanks = 0;
  78   ActiveEngine = -1; // -1: ALL, 0: Engine 1, 1: Engine 2 ...
  79   tankJ.InitMatrix();
  80   refuel = dump = false;
  81   DumpRate = 0.0;
  82   fuel_freeze = false;
  83   TotalFuelQuantity = 0.0;
  84   IsBound =
  85   HavePistonEngine =
  86   HaveTurbineEngine =
  87   HaveRocketEngine =
  88   HaveTurboPropEngine =
  89   HaveElectricEngine = false;
  90   HasInitializedEngines = false;
  91
  92   Debug(0);
  93 }
  94
  95 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  96
  97 FGPropulsion::~FGPropulsion()
  98 {
  99   for (unsigned int i=0; i<Engines.size(); i++) delete Engines[i];
 100   Engines.clear();
 101   for (unsigned int i=0; i<Tanks.size(); i++) delete Tanks[i];
 102   Tanks.clear();
 103   Debug(1);
 104 }
 105
 106 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 107
 108 bool FGPropulsion::InitModel(void)
 109 {
 110   if (!FGModel::InitModel()) return false;
 111
 112   for (unsigned int i=0; i<numTanks; i++) Tanks[i]->ResetToIC();
 113
 114   for (unsigned int i=0; i<numEngines; i++) {
 115     switch (Engines[i]->GetType()) {
 116       case FGEngine::etPiston:
 117         ((FGPiston*)Engines[i])->ResetToIC();
 118         if (HasInitializedEngines && (InitializedEngines & i)) InitRunning(i);
 119         break;
 120       case FGEngine::etTurbine:
 121         ((FGTurbine*)Engines[i])->ResetToIC();
 122         if (HasInitializedEngines && (InitializedEngines & i)) InitRunning(i);
 123         break;
 124       default:
 125         break;
 126     }
 127   }
 128
 129   return true;
 130 }
 131
 132 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 133
 134 bool FGPropulsion::Run(void)
 135 {
 136   unsigned int i;
 137
 138   if (FGModel::Run()) return true;
 139   if (FDMExec->Holding()) return false;
 140
 141   double dt = State->Getdt();
 142
 143   vForces.InitMatrix();
 144   vMoments.InitMatrix();
 145
 146   for (i=0; i<numEngines; i++) {
 147     Engines[i]->Calculate();
 148     vForces  += Engines[i]->GetBodyForces();  // sum body frame forces
 149     vMoments += Engines[i]->GetMoments();     // sum body frame moments
 150   }
 151
 152   TotalFuelQuantity = 0.0;
 153   for (i=0; i<numTanks; i++) {
 154     Tanks[i]->Calculate( dt * rate );
 155     if (Tanks[i]->GetType() == FGTank::ttFUEL) {
 156       TotalFuelQuantity += Tanks[i]->GetContents();
 157     }
 158   }
 159
 160   if (refuel) DoRefuel( dt * rate );
 161   if (dump) DumpFuel( dt * rate );
 162
 163   return false;
 164 }
 165
 166 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 167
 168 bool FGPropulsion::GetSteadyState(void)
 169 {
 170   double currentThrust = 0, lastThrust=-1;
 171   int steady_count,j=0;
 172   bool steady=false;
 173
 174   vForces.InitMatrix();
 175   vMoments.InitMatrix();
 176
 177   if (!FGModel::Run()) {
 178     for (unsigned int i=0; i<numEngines; i++) {
 179       cout << "  Finding steady state for engine " << i << endl;
 180       Engines[i]->SetTrimMode(true);
 181       steady=false;
 182       steady_count=0;
 183       j=0;
 184       while (!steady && j < 6000) {
 185         Engines[i]->Calculate();
 186         lastThrust = currentThrust;
 187         currentThrust = Engines[i]->GetThrust();
 188         if (fabs(lastThrust-currentThrust) < 0.0001) {
 189           steady_count++;
 190           if (steady_count > 120) {
 191             steady=true;
 192             cout << "    Steady state found at thrust: " << currentThrust << " lbs." << endl;
 193           }
 194         } else {
 195           steady_count=0;
 196         }
 197         j++;
 198       }
 199       if (j >= 6000) {
 200         cout << "    Could not find a steady state for this engine." << endl;
 201       }
 202       vForces  += Engines[i]->GetBodyForces();  // sum body frame forces
 203       vMoments += Engines[i]->GetMoments();     // sum body frame moments
 204       Engines[i]->SetTrimMode(false);
 205     }
 206
 207     return false;
 208   } else {
 209     return true;
 210   }
 211 }
 212
 213 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 214
 215 void FGPropulsion::InitRunning(int n)
 216 {
 217   if (n > 0) { // A specific engine is supposed to be initialized
 218
 219     if (n >= GetNumEngines() ) {
 220       cerr << "Tried to initialize a non-existent engine!" << endl;
 221       throw;
 222     }
 223     FCS->SetThrottleCmd(n,1);
 224     FCS->SetMixtureCmd(n,1);
 225     GetEngine(n)->InitRunning();
 226     GetSteadyState();
 227
 228     InitializedEngines = 1 << n;
 229     HasInitializedEngines = true;
 230
 231   } else if (n < 0) { // -1 refers to "All Engines"
 232
 233     for (unsigned int i=0; i<GetNumEngines(); i++) {
 234       FCS->SetThrottleCmd(i,1);
 235       FCS->SetMixtureCmd(i,1);
 236       GetEngine(i)->InitRunning();
 237     }
 238     GetSteadyState();
 239     InitializedEngines = -1;
 240     HasInitializedEngines = true;
 241
 242   } else if (n == 0) { // No engines are to be initialized
 243     // Do nothing
 244   }
 245 }
 246
 247 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 248
 249 bool FGPropulsion::Load(Element* el)
 250 {
 251   string type, engine_filename;
 252   bool ThrottleAdded = false;
 253
 254   Debug(2);
 255
 256   Element* engine_element = el->FindElement("engine");
 257   while (engine_element) {
 258     engine_filename = engine_element->GetAttributeValue("file");
 259
 260     if (engine_filename.empty()) {
 261       cerr << "Engine definition did not supply an engine file." << endl;
 262       return false;
 263     }
 264
 265     engine_filename = FindEngineFullPathname(engine_filename);
 266     document = LoadXMLDocument(engine_filename);
 267     document->SetParent(engine_element);
 268
 269     type = document->GetName();
 270     if (type == "piston_engine") {
 271       HavePistonEngine = true;
 272       if (!IsBound) bind();
 273       Engines.push_back(new FGPiston(FDMExec, document, numEngines));
 274     } else if (type == "turbine_engine") {
 275       HaveTurbineEngine = true;
 276       if (!IsBound) bind();
 277       Engines.push_back(new FGTurbine(FDMExec, document, numEngines));
 278     } else if (type == "turboprop_engine") {
 279       HaveTurboPropEngine = true;
 280       if (!IsBound) bind();
 281       Engines.push_back(new FGTurboProp(FDMExec, document, numEngines));
 282     } else if (type == "rocket_engine") {
 283       HaveRocketEngine = true;
 284       if (!IsBound) bind();
 285       Engines.push_back(new FGRocket(FDMExec, document, numEngines));
 286     } else if (type == "electric_engine") {
 287       HaveElectricEngine = true;
 288       if (!IsBound) bind();
 289       Engines.push_back(new FGElectric(FDMExec, document, numEngines));
 290     } else {
 291       cerr << "Unknown engine type: " << type << endl;
 292       exit(-5);
 293     }
 294
 295     FCS->AddThrottle();
 296     ThrottleAdded = true;
 297
 298     numEngines++;
 299
 300     engine_element = el->FindNextElement("engine");
 301     ResetParser();
 302   }
 303
 304   // Process tank definitions
 305
 306   Element* tank_element = el->FindElement("tank");
 307   while (tank_element) {
 308     Tanks.push_back(new FGTank(FDMExec, tank_element, numTanks));
 309     if (Tanks.back()->GetType() == FGTank::ttFUEL) numFuelTanks++;
 310     else if (Tanks.back()->GetType() == FGTank::ttOXIDIZER) numOxiTanks++;
 311     else {cerr << "Unknown tank type specified." << endl; return false;}
 312     numTanks++;
 313     tank_element = el->FindNextElement("tank");
 314   }
 315   numSelectedFuelTanks = numFuelTanks;
 316   numSelectedOxiTanks  = numOxiTanks;
 317
 318   CalculateTankInertias();
 319   if (!ThrottleAdded) FCS->AddThrottle(); // need to have at least one throttle
 320
 321   // Process fuel dump rate
 322   if (el->FindElement("dump-rate"))
 323     DumpRate = el->FindElementValueAsNumberConvertTo("dump-rate", "LBS/MIN");
 324
 325
 326   return true;
 327 }
 328
 329 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 330
 331 string FGPropulsion::FindEngineFullPathname(string engine_filename)
 332 {
 333   string fullpath, localpath;
 334   string enginePath = FDMExec->GetEnginePath();
 335   string aircraftPath = FDMExec->GetFullAircraftPath();
 336   ifstream engine_file;
 337
 338   string separator = "/";
 339 # ifdef macintosh
 340   separator = ";";
 341 # endif
 342
 343   fullpath = enginePath + separator;
 344   localpath = aircraftPath + separator + "Engines" + separator;
 345
 346   engine_file.open(string(fullpath + engine_filename + ".xml").c_str());
 347   if ( !engine_file.is_open()) {
 348     engine_file.open(string(localpath + engine_filename + ".xml").c_str());
 349       if ( !engine_file.is_open()) {
 350         cerr << " Could not open engine file: " << engine_filename << " in path "
 351              << fullpath << " or " << localpath << endl;
 352         return string("");
 353       } else {
 354         return string(localpath + engine_filename + ".xml");
 355       }
 356   }
 357   return string(fullpath + engine_filename + ".xml");
 358 }
 359
 360 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 361
 362 ifstream* FGPropulsion::FindEngineFile(string engine_filename)
 363 {
 364   string fullpath, localpath;
 365   string enginePath = FDMExec->GetEnginePath();
 366   string aircraftPath = FDMExec->GetFullAircraftPath();
 367   ifstream* engine_file = new ifstream();
 368
 369   string separator = "/";
 370 # ifdef macintosh
 371   separator = ";";
 372 # endif
 373
 374   fullpath = enginePath + separator;
 375   localpath = aircraftPath + separator + "Engines" + separator;
 376
 377   engine_file->open(string(fullpath + engine_filename + ".xml").c_str());
 378   if ( !engine_file->is_open()) {
 379     engine_file->open(string(localpath + engine_filename + ".xml").c_str());
 380       if ( !engine_file->is_open()) {
 381         cerr << " Could not open engine file: " << engine_filename << " in path "
 382              << fullpath << " or " << localpath << endl;
 383       }
 384   }
 385   return engine_file;
 386 }
 387
 388 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 389
 390 string FGPropulsion::GetPropulsionStrings(string delimeter)
 391 {
 392   unsigned int i;
 393
 394   string PropulsionStrings = "";
 395   bool firstime = true;
 396   stringstream buf;
 397
 398   for (i=0; i<Engines.size(); i++) {
 399     if (firstime)  firstime = false;
 400     else           PropulsionStrings += delimeter;
 401
 402     PropulsionStrings += Engines[i]->GetEngineLabels(delimeter);
 403   }
 404   for (i=0; i<Tanks.size(); i++) {
 405     if (Tanks[i]->GetType() == FGTank::ttFUEL) buf << delimeter << "Fuel Tank " << i;
 406     else if (Tanks[i]->GetType() == FGTank::ttOXIDIZER) buf << delimeter << "Oxidizer Tank " << i;
 407   }
 408
 409   return PropulsionStrings;
 410 }
 411
 412 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 413
 414 string FGPropulsion::GetPropulsionValues(string delimeter)
 415 {
 416   unsigned int i;
 417
 418   string PropulsionValues = "";
 419   bool firstime = true;
 420   stringstream buf;
 421
 422   for (i=0; i<Engines.size(); i++) {
 423     if (firstime)  firstime = false;
 424     else           PropulsionValues += delimeter;
 425
 426     PropulsionValues += Engines[i]->GetEngineValues(delimeter);
 427   }
 428   for (i=0; i<Tanks.size(); i++) {
 429     buf << delimeter;
 430     buf << Tanks[i]->GetContents();
 431   }
 432
 433   return PropulsionValues;
 434 }
 435
 436 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 437
 438 FGColumnVector3& FGPropulsion::GetTanksMoment(void)
 439 {
 440   vXYZtank_arm.InitMatrix();
 441   for (unsigned int i=0; i<Tanks.size(); i++) {
 442     vXYZtank_arm(eX) += Tanks[i]->GetXYZ(eX) * Tanks[i]->GetContents();
 443     vXYZtank_arm(eY) += Tanks[i]->GetXYZ(eY) * Tanks[i]->GetContents();
 444     vXYZtank_arm(eZ) += Tanks[i]->GetXYZ(eZ) * Tanks[i]->GetContents();
 445   }
 446   return vXYZtank_arm;
 447 }
 448
 449 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 450
 451 double FGPropulsion::GetTanksWeight(void)
 452 {
 453   double Tw = 0.0;
 454
 455   for (unsigned int i=0; i<Tanks.size(); i++) Tw += Tanks[i]->GetContents();
 456
 457   return Tw;
 458 }
 459
 460 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 461
 462 FGMatrix33& FGPropulsion::CalculateTankInertias(void)
 463 {
 464   unsigned int size;
 465
 466   size = Tanks.size();
 467   if (size == 0) return tankJ;
 468
 469   tankJ = FGMatrix33();
 470
 471   for (unsigned int i=0; i<size; i++)
 472     tankJ += MassBalance->GetPointmassInertia( lbtoslug * Tanks[i]->GetContents(),
 473                                                Tanks[i]->GetXYZ() );
 474
 475   return tankJ;
 476 }
 477
 478 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 479
 480 void FGPropulsion::SetMagnetos(int setting)
 481 {
 482   if (ActiveEngine < 0) {
 483     for (unsigned i=0; i<Engines.size(); i++) {
 484       // ToDo: first need to make sure the engine Type is really appropriate:
 485       //   do a check to see if it is of type Piston. This should be done for
 486       //   all of this kind of possibly across-the-board settings.
 487       ((FGPiston*)Engines[i])->SetMagnetos(setting);
 488     }
 489   } else {
 490     ((FGPiston*)Engines[ActiveEngine])->SetMagnetos(setting);
 491   }
 492 }
 493
 494 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 495
 496 void FGPropulsion::SetStarter(int setting)
 497 {
 498   if (ActiveEngine < 0) {
 499     for (unsigned i=0; i<Engines.size(); i++) {
 500       if (setting == 0)
 501         Engines[i]->SetStarter(false);
 502       else
 503         Engines[i]->SetStarter(true);
 504     }
 505   } else {
 506     if (setting == 0)
 507       Engines[ActiveEngine]->SetStarter(false);
 508     else
 509       Engines[ActiveEngine]->SetStarter(true);
 510   }
 511 }
 512
 513 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 514
 515 void FGPropulsion::SetCutoff(int setting)
 516 {
 517   if (ActiveEngine < 0) {
 518     for (unsigned i=0; i<Engines.size(); i++) {
 519       if (setting == 0)
 520         ((FGTurbine*)Engines[i])->SetCutoff(false);
 521       else
 522         ((FGTurbine*)Engines[i])->SetCutoff(true);
 523     }
 524   } else {
 525     if (setting == 0)
 526       ((FGTurbine*)Engines[ActiveEngine])->SetCutoff(false);
 527     else
 528       ((FGTurbine*)Engines[ActiveEngine])->SetCutoff(true);
 529   }
 530 }
 531
 532 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 533
 534 void FGPropulsion::SetActiveEngine(int engine)
 535 {
 536   if (engine >= (int)Engines.size() || engine < 0)
 537     ActiveEngine = -1;
 538   else
 539     ActiveEngine = engine;
 540 }
 541
 542 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 543
 544 double FGPropulsion::Transfer(int source, int target, double amount)
 545 {
 546  double shortage, overage;
 547
 548   if (source == -1) {
 549      shortage = 0.0;
 550   } else {
 551      shortage = Tanks[source]->Drain(amount);
 552   }
 553   if (target == -1) {
 554      overage = 0.0;
 555   } else {
 556      overage = Tanks[target]->Fill(amount - shortage);
 557   }
 558   return overage;
 559 }
 560
 561 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 562
 563 void FGPropulsion::DoRefuel(double time_slice)
 564 {
 565   unsigned int i;
 566
 567   double fillrate = 100 * time_slice;   // 100 lbs/sec = 6000 lbs/min
 568   int TanksNotFull = 0;
 569
 570   for (i=0; i<numTanks; i++) {
 571     if (Tanks[i]->GetPctFull() < 99.99) ++TanksNotFull;
 572   }
 573
 574   if (TanksNotFull) {
 575     for (i=0; i<numTanks; i++) {
 576       if (Tanks[i]->GetPctFull() < 99.99)
 577           Transfer(-1, i, fillrate/TanksNotFull);
 578     }
 579   }
 580 }
 581
 582 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 583
 584 void FGPropulsion::DumpFuel(double time_slice)
 585 {
 586   unsigned int i;
 587   int TanksDumping = 0;
 588
 589   for (i=0; i<numTanks; i++) {
 590     if (Tanks[i]->GetContents() > Tanks[i]->GetStandpipe()) ++TanksDumping;
 591   }
 592
 593   if (TanksDumping == 0) return;
 594
 595   double dump_rate_per_tank = DumpRate / 60.0 * time_slice / TanksDumping;
 596
 597   for (i=0; i<numTanks; i++) {
 598     if (Tanks[i]->GetContents() > Tanks[i]->GetStandpipe()) {
 599       Transfer(i, -1, dump_rate_per_tank);
 600     }
 601   }
 602 }
 603
 604 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 605
 606 void FGPropulsion::SetFuelFreeze(bool f)
 607 {
 608   fuel_freeze = f;
 609   for (unsigned int i=0; i<numEngines; i++) {
 610     Engines[i]->SetFuelFreeze(f);
 611   }
 612 }
 613
 614 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 615
 616 void FGPropulsion::bind(void)
 617 {
 618   typedef double (FGPropulsion::*PMF)(int) const;
 619   typedef int (FGPropulsion::*iPMF)(void) const;
 620
 621   IsBound = true;
 622   PropertyManager->Tie("propulsion/set-running", this, (iPMF)0, &FGPropulsion::InitRunning, true);
 623   if (HaveTurbineEngine) {
 624     PropertyManager->Tie("propulsion/starter_cmd", this, (iPMF)0, &FGPropulsion::SetStarter,  true);
 625     PropertyManager->Tie("propulsion/cutoff_cmd", this,  (iPMF)0, &FGPropulsion::SetCutoff,   true);
 626   }
 627
 628   if (HavePistonEngine) {
 629     PropertyManager->Tie("propulsion/starter_cmd", this, (iPMF)0, &FGPropulsion::SetStarter,  true);
 630     PropertyManager->Tie("propulsion/magneto_cmd", this, (iPMF)0, &FGPropulsion::SetMagnetos, true);
 631   }
 632
 633   PropertyManager->Tie("propulsion/active_engine", this, (iPMF)&FGPropulsion::GetActiveEngine,
 634                         &FGPropulsion::SetActiveEngine, true);
 635   PropertyManager->Tie("propulsion/total-fuel-lbs", this, &FGPropulsion::GetTotalFuelQuantity);
 636   PropertyManager->Tie("propulsion/refuel", this, &FGPropulsion::GetRefuel,
 637                         &FGPropulsion::SetRefuel, true);
 638   PropertyManager->Tie("propulsion/fuel_dump", this, &FGPropulsion::GetFuelDump,
 639                         &FGPropulsion::SetFuelDump, true);
 640   PropertyManager->Tie("forces/fbx-prop-lbs", this,1,
 641                        (PMF)&FGPropulsion::GetForces);
 642   PropertyManager->Tie("forces/fby-prop-lbs", this,2,
 643                        (PMF)&FGPropulsion::GetForces);
 644   PropertyManager->Tie("forces/fbz-prop-lbs", this,3,
 645                        (PMF)&FGPropulsion::GetForces);
 646   PropertyManager->Tie("moments/l-prop-lbsft", this,1,
 647                        (PMF)&FGPropulsion::GetMoments);
 648   PropertyManager->Tie("moments/m-prop-lbsft", this,2,
 649                        (PMF)&FGPropulsion::GetMoments);
 650   PropertyManager->Tie("moments/n-prop-lbsft", this,3,
 651                        (PMF)&FGPropulsion::GetMoments);
 652
 653 }
 654
 655 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 656 //    The bitmasked value choices are as follows:
 657 //    unset: In this case (the default) JSBSim would only print
 658 //       out the normally expected messages, essentially echoing
 659 //       the config files as they are read. If the environment
 660 //       variable is not set, debug_lvl is set to 1 internally
 661 //    0: This requests JSBSim not to output any messages
 662 //       whatsoever.
 663 //    1: This value explicity requests the normal JSBSim
 664 //       startup messages
 665 //    2: This value asks for a message to be printed out when
 666 //       a class is instantiated
 667 //    4: When this value is set, a message is displayed when a
 668 //       FGModel object executes its Run() method
 669 //    8: When this value is set, various runtime state variables
 670 //       are printed out periodically
 671 //    16: When set various parameters are sanity checked and
 672 //       a message is printed out when they go out of bounds
 673
 674 void FGPropulsion::Debug(int from)
 675 {
 676   if (debug_lvl <= 0) return;
 677
 678   if (debug_lvl & 1) { // Standard console startup message output
 679     if (from == 2) { // Loader
 680       cout << endl << "  Propulsion:" << endl;
 681     }
 682   }
 683   if (debug_lvl & 2 ) { // Instantiation/Destruction notification
 684     if (from == 0) cout << "Instantiated: FGPropulsion" << endl;
 685     if (from == 1) cout << "Destroyed:    FGPropulsion" << endl;
 686   }
 687   if (debug_lvl & 4 ) { // Run() method entry print for FGModel-derived objects
 688   }
 689   if (debug_lvl & 8 ) { // Runtime state variables
 690   }
 691   if (debug_lvl & 16) { // Sanity checking
 692   }
 693   if (debug_lvl & 64) {
 694     if (from == 0) { // Constructor
 695       cout << IdSrc << endl;
 696       cout << IdHdr << endl;
 697     }
 698   }
 699 }
 700 }