src/AIModel/submodel.cxx

   1 //// submodel.cxx - models a releasable submodel.
   2 // Written by Dave Culp, started Aug 2004
   3 // With major additions by Vivian Meaaza 2004 - 2007
   4 //
   5 // This file is in the Public Domain and comes with no warranty.
   6
   7 #ifdef HAVE_CONFIG_H
   8 #  include "config.h"
   9 #endif
  10
  11 #include "submodel.hxx"
  12
  13 #include <simgear/structure/exception.hxx>
  14 #include <simgear/misc/sg_path.hxx>
  15 #include <simgear/math/sg_geodesy.hxx>
  16 #include <simgear/props/props_io.hxx>
  17
  18 #include <Main/fg_props.hxx>
  19 #include <Main/util.hxx>
  20
  21
  22 #include "AIBase.hxx"
  23 #include "AIManager.hxx"
  24 #include "AIBallistic.hxx"
  25
  26 using std::cout;
  27 using std::endl;
  28 using std::string;
  29 using std::vector;
  30
  31 const double FGSubmodelMgr::lbs_to_slugs = 0.031080950172;
  32
  33 FGSubmodelMgr::FGSubmodelMgr()
  34 {
  35     x_offset = y_offset = z_offset = 0.0;
  36     pitch_offset = 0.0;
  37     yaw_offset = 0.0;
  38
  39     //out[0] = out[1] = out[2] = 0;
  40     //string contents_node;
  41     contrail_altitude = 30000;
  42     _count = 0;
  43     _found_sub = true;
  44 }
  45
  46 FGSubmodelMgr::~FGSubmodelMgr()
  47 {
  48 }
  49
  50 FGAIManager* FGSubmodelMgr::aiManager()
  51 {
  52    return (FGAIManager*)globals->get_subsystem("ai-model");
  53 }
  54
  55 void FGSubmodelMgr::init()
  56 {
  57     index = 0;
  58
  59     _serviceable_node = fgGetNode("/sim/submodels/serviceable", true);
  60     _serviceable_node->setBoolValue(true);
  61
  62     _user_lat_node = fgGetNode("/position/latitude-deg", true);
  63     _user_lon_node = fgGetNode("/position/longitude-deg", true);
  64     _user_alt_node = fgGetNode("/position/altitude-ft", true);
  65
  66     _user_heading_node = fgGetNode("/orientation/heading-deg", true);
  67     _user_pitch_node   = fgGetNode("/orientation/pitch-deg", true);
  68     _user_roll_node    = fgGetNode("/orientation/roll-deg", true);
  69     _user_yaw_node     = fgGetNode("/orientation/yaw-deg", true);
  70     _user_alpha_node   = fgGetNode("/orientation/alpha-deg", true);
  71
  72     _user_speed_node = fgGetNode("/velocities/uBody-fps", true);
  73
  74     _user_wind_from_east_node  = fgGetNode("/environment/wind-from-east-fps", true);
  75     _user_wind_from_north_node = fgGetNode("/environment/wind-from-north-fps", true);
  76
  77     _user_speed_down_fps_node   = fgGetNode("/velocities/speed-down-fps", true);
  78     _user_speed_east_fps_node   = fgGetNode("/velocities/speed-east-fps", true);
  79     _user_speed_north_fps_node  = fgGetNode("/velocities/speed-north-fps", true);
  80
  81     _contrail_altitude_node = fgGetNode("/environment/params/contrail-altitude", true);
  82     contrail_altitude       = _contrail_altitude_node->getDoubleValue();
  83     _contrail_trigger       = fgGetNode("ai/submodels/contrails", true);
  84     _contrail_trigger->setBoolValue(false);
  85
  86     load();
  87 }
  88
  89 void FGSubmodelMgr::postinit()
  90 {
  91     // postinit, so that the AI list is populated
  92     loadAI();
  93
  94     while (_found_sub)
  95         loadSubmodels();
  96
  97     //TODO reload submodels if an MP ac joins
  98
  99     //_model_added_node = fgGetNode("ai/models/model-added", true);
 100     //_model_added_node->addChangeListener(this, false);
 101 }
 102
 103 void FGSubmodelMgr::bind()
 104 {
 105 }
 106
 107 void FGSubmodelMgr::unbind()
 108 {
 109     submodel_iterator = submodels.begin();
 110     while (submodel_iterator != submodels.end()) {
 111         (*submodel_iterator)->prop->untie("count");
 112         ++submodel_iterator;
 113     }
 114 }
 115
 116 void FGSubmodelMgr::update(double dt)
 117 {
 118     if (!_serviceable_node->getBoolValue())
 119         return;
 120
 121     _impact = false;
 122     _hit = false;
 123     _expiry = false;
 124
 125     // Check if the submodel hit an object or terrain
 126     FGAIManager::ai_list_type sm_list(aiManager()->get_ai_list());
 127     FGAIManager::ai_list_iterator sm_list_itr = sm_list.begin();
 128     FGAIManager::ai_list_iterator end         = sm_list.end();
 129
 130     for (; sm_list_itr != end; ++sm_list_itr) {
 131         FGAIBase::object_type object_type =(*sm_list_itr)->getType();
 132
 133         // Continue if object is not ballistic
 134         if (object_type != FGAIBase::otBallistic) {
 135             continue;
 136         }
 137
 138         int parent_subID = (*sm_list_itr)->_getSubID();
 139         int id = (*sm_list_itr)->getID();
 140
 141         if (parent_subID == 0 || id == -1) // this entry in the list has no associated submodel
 142             continue;                       // or is invalid so we can continue
 143
 144         //SG_LOG(SG_AI, SG_DEBUG, "Submodel: Impact " << _impact << " hit! "
 145         //        << _hit <<" parent_subID " << parent_subID);
 146
 147         _hit = (*sm_list_itr)->_getCollisionData();
 148         _impact = (*sm_list_itr)->_getImpactData();
 149         _expiry = (*sm_list_itr)->_getExpiryData();
 150
 151         //SG_LOG(SG_AI, SG_ALERT, "Submodel: " << (*sm_list_itr)->_getName()
 152         //    << " Impact " << _impact << " hit! " << _hit
 153         //    << " exipiry :-( " << _expiry );
 154
 155         if (_impact || _hit || _expiry) {
 156     //        SG_LOG(SG_AI, SG_ALERT, "Submodel: Impact " << _impact << " hit! " << _hit
 157                 //<< " exipiry :-( " << _expiry );
 158
 159             submodel_iterator = submodels.begin();
 160
 161             while (submodel_iterator != submodels.end()) {
 162                 int child_ID = (*submodel_iterator)->id;
 163                 //cout << "Impact: parent SubID " << parent_subID << " child_ID " << child_ID << endl;
 164
 165                 if ( parent_subID == child_ID ) {
 166                     _parent_lat = (*sm_list_itr)->_getImpactLat();
 167                     _parent_lon = (*sm_list_itr)->_getImpactLon();
 168                     _parent_elev = (*sm_list_itr)->_getImpactElevFt();
 169                     _parent_hdg = (*sm_list_itr)->_getImpactHdg();
 170                     _parent_pitch = (*sm_list_itr)->_getImpactPitch();
 171                     _parent_roll = (*sm_list_itr)->_getImpactRoll();
 172                     _parent_speed = (*sm_list_itr)->_getImpactSpeed();
 173                     (*submodel_iterator)->first_time = true;
 174                     //cout << "Impact: parent SubID = child_ID elev " << _parent_elev << endl;
 175
 176                     if (release(*submodel_iterator, dt)) {
 177                         (*sm_list_itr)->setDie(true);
 178                         //cout << "Impact: set die" << (*sm_list_itr)->_getName() << endl;
 179                     }
 180
 181                 }
 182
 183                 ++submodel_iterator;
 184             }
 185         }
 186     }
 187
 188     _contrail_trigger->setBoolValue(_user_alt_node->getDoubleValue() > contrail_altitude);
 189
 190     bool trigger = false;
 191     int i = -1;
 192
 193     submodel_iterator = submodels.begin();
 194     while (submodel_iterator != submodels.end())  {
 195         i++;
 196
 197         /*SG_LOG(SG_AI, SG_DEBUG,
 198                 "Submodels:  " << (*submodel_iterator)->id
 199                 << " name " << (*submodel_iterator)->name
 200                 );*/
 201
 202         if ((*submodel_iterator)->trigger_node != 0) {
 203             _trigger_node = (*submodel_iterator)->trigger_node;
 204             trigger = _trigger_node->getBoolValue();
 205             //cout << (*submodel_iterator)->name << "trigger node found " <<  trigger << endl;
 206         }
 207         else {
 208             trigger = false;
 209             //cout << (*submodel_iterator)->name << " trigger node not found " << trigger << endl;
 210         }
 211
 212         if (trigger && (*submodel_iterator)->count != 0) {
 213             //int id = (*submodel_iterator)->id;
 214             //const string& name = (*submodel_iterator)->name;
 215
 216             SG_LOG(SG_AI, SG_DEBUG,
 217             "Submodels release:  " << (*submodel_iterator)->id
 218             << " name " << (*submodel_iterator)->name
 219             << " count " << (*submodel_iterator)->count
 220             << " slaved " << (*submodel_iterator)->slaved
 221             );
 222
 223             release(*submodel_iterator, dt);
 224         } else
 225             (*submodel_iterator)->first_time = true;
 226
 227         ++submodel_iterator;
 228     }
 229 }
 230
 231 bool FGSubmodelMgr::release(submodel *sm, double dt)
 232 {
 233     // Only run if first time or repeat is set to true
 234     if (!sm->first_time && !sm->repeat) {
 235         return false;
 236     }
 237
 238     sm->timer += dt;
 239
 240     if (sm->timer < sm->delay) {
 241         //cout << "not yet: timer " << sm->timer << " delay " << sm->delay << endl;
 242         return false;
 243     }
 244
 245     //cout << "released timer: " << sm->timer << " delay " << sm->delay << endl;
 246
 247     sm->timer = 0.0;
 248
 249     if (sm->first_time) {
 250         dt = 0.0;
 251         sm->first_time = false;
 252     }
 253
 254     // Calculate submodel's initial conditions in world-coordinates
 255     transform(sm);
 256
 257     FGAIBallistic* ballist = new FGAIBallistic;
 258     ballist->setPath(sm->model.c_str());
 259     ballist->setName(sm->name);
 260     ballist->setSlaved(sm->slaved);
 261     ballist->setRandom(sm->random);
 262     ballist->setLifeRandomness(sm->life_randomness->get_value());
 263     ballist->setLatitude(offsetpos.getLatitudeDeg());
 264     ballist->setLongitude(offsetpos.getLongitudeDeg());
 265     ballist->setAltitude(offsetpos.getElevationFt());
 266     ballist->setAzimuthRandomError(sm->azimuth_error->get_value());
 267     ballist->setAzimuth(IC.azimuth);
 268     ballist->setElevationRandomError(sm->elevation_error->get_value());
 269     ballist->setElevation(IC.elevation);
 270     ballist->setRoll(IC.roll);
 271     ballist->setSpeed(IC.speed / SG_KT_TO_FPS);
 272     ballist->setWind_from_east(IC.wind_from_east);
 273     ballist->setWind_from_north(IC.wind_from_north);
 274     ballist->setMass(IC.mass);
 275     ballist->setDragArea(sm->drag_area);
 276     ballist->setLife(sm->life);
 277     ballist->setBuoyancy(sm->buoyancy);
 278     ballist->setWind(sm->wind);
 279     ballist->setCdRandomness(sm->cd_randomness->get_value());
 280     ballist->setCd(sm->cd);
 281     ballist->setStabilisation(sm->aero_stabilised);
 282     ballist->setNoRoll(sm->no_roll);
 283     ballist->setCollision(sm->collision);
 284     ballist->setExpiry(sm->expiry);
 285     ballist->setImpact(sm->impact);
 286     ballist->setImpactReportNode(sm->impact_report);
 287     ballist->setFuseRange(sm->fuse_range);
 288     ballist->setSubmodel(sm->submodel.c_str());
 289     ballist->setSubID(sm->sub_id);
 290     ballist->setForceStabilisation(sm->force_stabilised);
 291     ballist->setExternalForce(sm->ext_force);
 292     ballist->setForcePath(sm->force_path.c_str());
 293     ballist->setXoffset(sm->x_offset);
 294     ballist->setYoffset(sm->y_offset);
 295     ballist->setZoffset(sm->z_offset);
 296     ballist->setPitchoffset(sm->pitch_offset->get_value());
 297     ballist->setYawoffset(sm->yaw_offset->get_value());
 298     ballist->setParentNodes(_selected_ac);
 299     ballist->setContentsNode(sm->contents_node);
 300     ballist->setWeight(sm->weight);
 301
 302     aiManager()->attach(ballist);
 303
 304     if (sm->count > 0)
 305         sm->count--;
 306     return true;
 307 }
 308
 309 void FGSubmodelMgr::load()
 310 {
 311     SGPropertyNode_ptr path_node = fgGetNode("/sim/submodels/path");
 312
 313     if (path_node) {
 314         const int id = 0;
 315         const string& path = path_node->getStringValue();
 316         bool serviceable =_serviceable_node->getBoolValue();
 317         setData(id, path, serviceable, "/ai/submodels/submodel", submodels);
 318     }
 319 }
 320
 321 void FGSubmodelMgr::transform(submodel *sm)
 322 {
 323     // Set initial conditions
 324     if (sm->contents_node != 0 && !sm->slaved) {
 325         // Get the weight of the contents (lbs) and convert to mass (slugs)
 326         sm->contents = sm->contents_node->getChild("level-lbs",0,1)->getDoubleValue();
 327         //cout << "transform: contents " << sm->contents << endl;
 328         IC.mass = (sm->weight + sm->contents) * lbs_to_slugs;
 329         //cout << "mass inc contents"  << IC.mass << endl;
 330
 331         // Set contents to 0 in the parent
 332         sm->contents_node->getChild("level-gal_us",0,1)->setDoubleValue(0);
 333         /*cout << "contents " << sm->contents_node->getChild("level-gal_us")->getDoubleValue()
 334         << " " << sm->contents_node->getChild("level-lbs",0,1)->getDoubleValue()
 335         << endl;*/
 336     }
 337     else {
 338         IC.mass = sm->weight * lbs_to_slugs;
 339     }
 340
 341     int id = sm->id;
 342     //int sub_id = sm->sub_id;
 343     //const string& name = sm->name;
 344
 345     if (sm->speed_node != 0)
 346         sm->speed = sm->speed_node->getDoubleValue();
 347
 348     double yaw_offset = sm->yaw_offset->get_value();
 349     double pitch_offset = sm->pitch_offset->get_value();
 350
 351     //cout << " name " << name << " id " << id << " sub id" << sub_id << endl;
 352
 353     // set the Initial Conditions for the types of submodel parent
 354
 355     if (_impact || _hit || _expiry) {
 356         _count++;
 357
 358         // Set the data for a submodel tied to a submodel
 359         IC.lat             = _parent_lat;
 360         IC.lon             = _parent_lon;
 361         IC.alt             = _parent_elev;
 362         IC.roll            = _parent_roll;    // rotation about x axis
 363         IC.elevation       = _parent_pitch;   // rotation about y axis
 364         IC.azimuth         = _parent_hdg;     // rotation about z axis
 365         IC.speed           = _parent_speed;
 366         IC.speed_down_fps  = 0;
 367         IC.speed_east_fps  = 0;
 368         IC.speed_north_fps = 0;
 369     }
 370     else if (id == 0) {
 371         // Set the data for a submodel tied to the main model
 372         IC.lat             = _user_lat_node->getDoubleValue();
 373         IC.lon             = _user_lon_node->getDoubleValue();
 374         IC.alt             = _user_alt_node->getDoubleValue();
 375         IC.roll            = _user_roll_node->getDoubleValue();    // rotation about x axis
 376         IC.elevation       = _user_pitch_node->getDoubleValue();   // rotation about y axis
 377         IC.azimuth         = _user_heading_node->getDoubleValue(); // rotation about z axis
 378         IC.speed           = _user_speed_node->getDoubleValue();
 379         IC.speed_down_fps  = _user_speed_down_fps_node->getDoubleValue();
 380         IC.speed_east_fps  = _user_speed_east_fps_node->getDoubleValue();
 381         IC.speed_north_fps = _user_speed_north_fps_node->getDoubleValue();
 382     }
 383     else {
 384         // Set the data for a submodel tied to an AI Object
 385         setParentNode(id);
 386     }
 387
 388     //cout << "Submodel: setting IC "<< name << endl;
 389     //cout << "heading " << IC.azimuth << endl ;
 390     //cout << "speed down " << IC.speed_down_fps << endl ;
 391     //cout << "speed east " << IC.speed_east_fps << endl ;
 392     //cout << "speed north " << IC.speed_north_fps << endl ;
 393     //cout << "parent speed fps in " << IC.speed << "sm speed in " << sm->speed << endl ;
 394     //cout << "lat " << IC.lat;
 395     //cout << "alt " << IC.alt <<  endl ;
 396
 397     // Set the Initial Conditions that are common to all types of parent
 398     IC.wind_from_east =  _user_wind_from_east_node->getDoubleValue();
 399     IC.wind_from_north = _user_wind_from_north_node->getDoubleValue();
 400
 401     //cout << "wind e " << IC.wind_from_east << " n " << IC.wind_from_north << endl;
 402
 403     userpos.setLatitudeDeg(IC.lat);
 404     userpos.setLongitudeDeg(IC.lon);
 405     userpos.setElevationFt(IC.alt);
 406
 407     _x_offset = sm->x_offset;
 408     _y_offset = sm->y_offset;
 409     _z_offset = sm->z_offset;
 410
 411     setOffsetPos();
 412
 413     // Pre-process the trig functions
 414     cosRx = cos(-IC.roll * SG_DEGREES_TO_RADIANS);
 415     sinRx = sin(-IC.roll * SG_DEGREES_TO_RADIANS);
 416     cosRy = cos(-IC.elevation * SG_DEGREES_TO_RADIANS);
 417     sinRy = sin(-IC.elevation * SG_DEGREES_TO_RADIANS);
 418     cosRz = cos(IC.azimuth * SG_DEGREES_TO_RADIANS);
 419     sinRz = sin(IC.azimuth * SG_DEGREES_TO_RADIANS);
 420
 421     // Get submodel initial velocity vector angles in XZ and XY planes.
 422     // This vector should be added to aircraft's vector.
 423     IC.elevation += (yaw_offset * sinRx) + (pitch_offset * cosRx);
 424     IC.azimuth   += (yaw_offset * cosRx) - (pitch_offset * sinRx);
 425
 426     // Calculate the total speed north
 427     IC.total_speed_north = sm->speed * cos(IC.elevation * SG_DEGREES_TO_RADIANS)
 428             * cos(IC.azimuth * SG_DEGREES_TO_RADIANS) + IC.speed_north_fps;
 429
 430     // Calculate the total speed east
 431     IC.total_speed_east = sm->speed * cos(IC.elevation * SG_DEGREES_TO_RADIANS)
 432             * sin(IC.azimuth * SG_DEGREES_TO_RADIANS) + IC.speed_east_fps;
 433
 434     // Calculate the total speed down
 435     IC.total_speed_down = sm->speed * -sin(IC.elevation * SG_DEGREES_TO_RADIANS)
 436             + IC.speed_down_fps;
 437
 438     // Re-calculate speed, elevation and azimuth
 439     IC.speed = sqrt(IC.total_speed_north * IC.total_speed_north
 440             + IC.total_speed_east * IC.total_speed_east
 441             + IC.total_speed_down * IC.total_speed_down);
 442
 443     // If speeds are low this calculation can become unreliable
 444     if (IC.speed > 1) {
 445         IC.azimuth = atan2(IC.total_speed_east, IC.total_speed_north) * SG_RADIANS_TO_DEGREES;
 446         //        cout << "azimuth1 " << IC.azimuth<<endl;
 447
 448         // Rationalize the output
 449         if (IC.azimuth < 0)
 450             IC.azimuth += 360;
 451         else if (IC.azimuth >= 360)
 452             IC.azimuth -= 360;
 453         // cout << "azimuth2 " << IC.azimuth<<endl;
 454
 455         IC.elevation = -atan(IC.total_speed_down / sqrt(IC.total_speed_north
 456             * IC.total_speed_north + IC.total_speed_east * IC.total_speed_east))
 457             * SG_RADIANS_TO_DEGREES;
 458     }
 459     //cout << "IC.speed " << IC.speed / SG_KT_TO_FPS << endl;
 460 }
 461
 462 void FGSubmodelMgr::updatelat(double lat)
 463 {
 464     ft_per_deg_latitude = 366468.96 - 3717.12 * cos(lat / SG_RADIANS_TO_DEGREES);
 465     ft_per_deg_longitude = 365228.16 * cos(lat / SG_RADIANS_TO_DEGREES);
 466 }
 467
 468 void FGSubmodelMgr::loadAI()
 469 {
 470     SG_LOG(SG_AI, SG_DEBUG, "Submodels: Loading AI submodels");
 471
 472     FGAIManager::ai_list_type sm_list(aiManager()->get_ai_list());
 473
 474     if (sm_list.empty()) {
 475         SG_LOG(SG_AI, SG_ALERT, "Submodels: Unable to read AI submodel list");
 476         return;
 477     }
 478
 479     FGAIManager::ai_list_iterator sm_list_itr = sm_list.begin();
 480     FGAIManager::ai_list_iterator end         = sm_list.end();
 481
 482     while (sm_list_itr != end) {
 483         string path = (*sm_list_itr)->_getSMPath();
 484
 485         if (path.empty()) {
 486             ++sm_list_itr;
 487             continue;
 488         }
 489
 490         int id = (*sm_list_itr)->getID();
 491         bool serviceable = (*sm_list_itr)->_getServiceable();
 492
 493         //string type = (*sm_list_itr)->getTypeString();
 494         //cout << "loadAI: type " << type << " path "<< path << " serviceable " << serviceable << endl;
 495
 496         setData(id, path, serviceable, "/ai/submodels/submodel", submodels);
 497         ++sm_list_itr;
 498     }
 499 }
 500
 501 void FGSubmodelMgr::setData(int id, const string& path, bool serviceable, const string& property_path, submodel_vector_type& models)
 502 {
 503     SGPropertyNode root;
 504
 505     SGPath config = globals->resolve_aircraft_path(path);
 506     try {
 507         SG_LOG(SG_AI, SG_DEBUG,
 508                 "Submodels: Trying to read AI submodels file: " << config.str());
 509         readProperties(config.str(), &root);
 510     }
 511     catch (const sg_exception &) {
 512         SG_LOG(SG_AI, SG_ALERT,
 513                 "Submodels: Unable to read AI submodels file: " << config.str());
 514         return;
 515     }
 516
 517     vector<SGPropertyNode_ptr> children = root.getChildren("submodel");
 518     vector<SGPropertyNode_ptr>::iterator it = children.begin();
 519     vector<SGPropertyNode_ptr>::iterator end = children.end();
 520
 521     for (int i = 0; it != end; ++it, i++) {
 522         //cout << "Reading AI submodel " << (*it)->getPath() << endl;
 523         submodel *sm = new submodel;
 524         SGPropertyNode_ptr entry_node = *it;
 525         sm->name             = entry_node->getStringValue("name", "none_defined");
 526         sm->model            = entry_node->getStringValue("model", "Models/Geometry/rocket.ac");
 527         sm->speed            = entry_node->getDoubleValue("speed", 2329.4);
 528         sm->repeat           = entry_node->getBoolValue("repeat", false);
 529         sm->delay            = entry_node->getDoubleValue("delay", 0.25);
 530         sm->count            = entry_node->getIntValue("count", 1);
 531         sm->slaved           = entry_node->getBoolValue("slaved", false);
 532         sm->x_offset         = entry_node->getDoubleValue("x-offset", 0.0);
 533         sm->y_offset         = entry_node->getDoubleValue("y-offset", 0.0);
 534         sm->z_offset         = entry_node->getDoubleValue("z-offset", 0.0);
 535         sm->drag_area        = entry_node->getDoubleValue("eda", 0.034);
 536         sm->life             = entry_node->getDoubleValue("life", 900.0);
 537         sm->buoyancy         = entry_node->getDoubleValue("buoyancy", 0);
 538         sm->wind             = entry_node->getBoolValue("wind", false);
 539         sm->cd               = entry_node->getDoubleValue("cd", 0.193);
 540         sm->weight           = entry_node->getDoubleValue("weight", 0.25);
 541         sm->aero_stabilised  = entry_node->getBoolValue("aero-stabilised", true);
 542         sm->no_roll          = entry_node->getBoolValue("no-roll", false);
 543         sm->collision        = entry_node->getBoolValue("collision", false);
 544         sm->expiry           = entry_node->getBoolValue("expiry", false);
 545         sm->impact           = entry_node->getBoolValue("impact", false);
 546         sm->impact_report    = entry_node->getStringValue("impact-reports");
 547         sm->fuse_range       = entry_node->getDoubleValue("fuse-range", 0.0);
 548         sm->contents_node    = fgGetNode(entry_node->getStringValue("contents", "none"), false);
 549         sm->speed_node       = fgGetNode(entry_node->getStringValue("speed-prop", "none"), false);
 550         sm->submodel         = entry_node->getStringValue("submodel-path", "");
 551         sm->force_stabilised = entry_node->getBoolValue("force-stabilised", false);
 552         sm->ext_force        = entry_node->getBoolValue("external-force", false);
 553         sm->force_path       = entry_node->getStringValue("force-path", "");
 554         sm->random           = entry_node->getBoolValue("random", false);
 555
 556         SGPropertyNode_ptr prop_root = fgGetNode("/", true);
 557         SGPropertyNode n;
 558
 559         SGPropertyNode_ptr a = entry_node->getNode("yaw-offset");
 560         sm->yaw_offset   = new FGXMLAutopilot::InputValue(*prop_root, a ? *a : n);
 561
 562         a = entry_node->getNode("pitch-offset");
 563         sm->pitch_offset = new FGXMLAutopilot::InputValue(*prop_root, a ? *a : n);
 564
 565         // Randomness
 566         a = entry_node->getNode("randomness", true);
 567         SGPropertyNode_ptr b;
 568
 569         // Maximum azimuth randomness error in degrees
 570         b = a->getNode("azimuth");
 571         sm->azimuth_error   = new FGXMLAutopilot::InputValue(*prop_root, b ? *b : n);
 572
 573         // Maximum elevation randomness error in degrees
 574         b = a->getNode("elevation");
 575         sm->elevation_error = new FGXMLAutopilot::InputValue(*prop_root, b ? *b : n);
 576
 577         // Randomness of Cd (plus or minus)
 578         b = a->getNode("cd");
 579         if (!b) {
 580             b = a->getNode("cd", true);
 581             b->setDoubleValue(0.1);
 582         }
 583         sm->cd_randomness   = new FGXMLAutopilot::InputValue(*prop_root, *b);
 584
 585         // Randomness of life (plus or minus)
 586         b = a->getNode("life");
 587         if (!b) {
 588             b = a->getNode("life", true);
 589             b->setDoubleValue(0.5);
 590         }
 591         sm->life_randomness = new FGXMLAutopilot::InputValue(*prop_root, *b);
 592
 593         if (sm->contents_node != 0)
 594             sm->contents = sm->contents_node->getDoubleValue();
 595
 596         const char *trigger_path = entry_node->getStringValue("trigger", 0);
 597         if (trigger_path) {
 598             sm->trigger_node = fgGetNode(trigger_path, true);
 599             sm->trigger_node->setBoolValue(sm->trigger_node->getBoolValue());
 600         }
 601         else {
 602             sm->trigger_node = 0;
 603         }
 604
 605         if (sm->speed_node != 0)
 606             sm->speed = sm->speed_node->getDoubleValue();
 607
 608         sm->timer = sm->delay;
 609         sm->id = id;
 610         sm->first_time = false;
 611         sm->serviceable = serviceable;
 612         sm->sub_id = 0;
 613
 614         sm->prop = fgGetNode(property_path, index, true);
 615         sm->prop->tie("delay", SGRawValuePointer<double>(&(sm->delay)));
 616         sm->prop->tie("count", SGRawValuePointer<int>(&(sm->count)));
 617         sm->prop->tie("repeat", SGRawValuePointer<bool>(&(sm->repeat)));
 618         sm->prop->tie("id", SGRawValuePointer<int>(&(sm->id)));
 619         sm->prop->tie("sub-id", SGRawValuePointer<int>(&(sm->sub_id)));
 620         sm->prop->tie("serviceable", SGRawValuePointer<bool>(&(sm->serviceable)));
 621         sm->prop->tie("random", SGRawValuePointer<bool>(&(sm->random)));
 622         sm->prop->tie("slaved", SGRawValuePointer<bool>(&(sm->slaved)));
 623         const string& name = sm->name;
 624         sm->prop->setStringValue("name", name.c_str());
 625
 626         const string& submodel = sm->submodel;
 627         sm->prop->setStringValue("submodel", submodel.c_str());
 628
 629         const string& force_path = sm->force_path;
 630         sm->prop->setStringValue("force_path", force_path.c_str());
 631
 632         if (sm->contents_node != 0)
 633             sm->prop->tie("contents-lbs", SGRawValuePointer<double>(&(sm->contents)));
 634
 635         index++;
 636         models.push_back(sm);
 637     }
 638 }
 639
 640 void FGSubmodelMgr::loadSubmodels()
 641 {
 642     SG_LOG(SG_AI, SG_DEBUG, "Submodels: Loading sub submodels");
 643
 644     _found_sub = false;
 645
 646     submodel_iterator = submodels.begin();
 647
 648     while (submodel_iterator != submodels.end()) {
 649         const string& submodel  = (*submodel_iterator)->submodel;
 650         if (!submodel.empty()) {
 651             //int id = (*submodel_iterator)->id;
 652             bool serviceable = true;
 653             SG_LOG(SG_AI, SG_DEBUG, "found path sub sub "
 654                     << submodel
 655                     << " index " << index
 656                     << " name " << (*submodel_iterator)->name);
 657
 658             if ((*submodel_iterator)->sub_id == 0) {
 659                 (*submodel_iterator)->sub_id = index;
 660                 _found_sub = true;
 661                 setData(index, submodel, serviceable, "/ai/submodels/subsubmodel", subsubmodels);
 662             }
 663         }
 664
 665         ++submodel_iterator;
 666     }
 667
 668     submodels.reserve(submodels.size() + subsubmodels.size());
 669
 670     // Add all elements from subsubmodels to submodels
 671     subsubmodel_iterator = subsubmodels.begin();
 672     while (subsubmodel_iterator != subsubmodels.end()) {
 673         submodels.push_back(*subsubmodel_iterator);
 674         ++subsubmodel_iterator;
 675     }
 676
 677     subsubmodels.clear();
 678 }
 679
 680 SGVec3d FGSubmodelMgr::getCartOffsetPos() const
 681 {
 682     // Convert geodetic positions to geocentered
 683     SGVec3d cartuserPos = SGVec3d::fromGeod(userpos);
 684
 685     // Transform to the right coordinate frame, configuration is done in
 686     // the x-forward, y-right, z-up coordinates (feet), computation
 687     // in the simulation usual body x-forward, y-right, z-down coordinates
 688     // (meters) )
 689     SGVec3d _off(_x_offset * SG_FEET_TO_METER,
 690         _y_offset * SG_FEET_TO_METER,
 691         -_z_offset * SG_FEET_TO_METER);
 692
 693     // Transform the user position to the horizontal local coordinate system.
 694     SGQuatd hlTrans = SGQuatd::fromLonLat(userpos);
 695
 696     // And postrotate the orientation of the user model wrt the horizontal
 697     // local frame
 698     hlTrans *= SGQuatd::fromYawPitchRollDeg(
 699        IC.azimuth,
 700        IC.elevation,
 701        IC.roll);
 702
 703     // The offset converted to the usual body fixed coordinate system
 704     // rotated to the earth-fixed coordinates axis
 705     SGVec3d off = hlTrans.backTransform(_off);
 706
 707     // Add the position offset of the user model to get the geocentered position
 708     SGVec3d offsetPos = cartuserPos + off;
 709     return offsetPos;
 710 }
 711
 712 void FGSubmodelMgr::setOffsetPos()
 713 {
 714     // Convert the offset geocentered position to geodetic
 715     SGVec3d cartoffsetPos = getCartOffsetPos();
 716     SGGeodesy::SGCartToGeod(cartoffsetPos, offsetpos);
 717 }
 718
 719 void FGSubmodelMgr::valueChanged(SGPropertyNode *prop)
 720 {
 721     return; // this isn't working atm
 722
 723     const char* _model_added = _model_added_node->getStringValue();
 724
 725     std::basic_string <char>::size_type indexCh2b;
 726
 727     string str2 = _model_added;
 728     const char *cstr2b = "multiplayer";
 729     indexCh2b = str2.find(cstr2b, 0);
 730
 731     // Ignoring Ballistic Objects; there are potentially too many
 732     if (indexCh2b != string::npos ) {
 733         //cout << "Submodels: model added - " << str2 <<" read path "<< endl;
 734         //return;
 735         SGPropertyNode *a_node = fgGetNode(_model_added, true);
 736         SGPropertyNode *sub_node = a_node->getChild("sim", 0, true);
 737         SGPropertyNode_ptr path_node = sub_node->getChild("path", 0, true);
 738         SGPropertyNode_ptr callsign_node = a_node->getChild("callsign", 0, true);
 739
 740         //const string& callsign = callsign_node->getStringValue();
 741         //cout << "Submodels: model added - " << callsign <<" read callsign "<< endl;
 742     }
 743     else {
 744         cout << "model added - " << str2 <<" returning " << endl;
 745     }
 746 }
 747
 748 void FGSubmodelMgr::setParentNode(int id)
 749 {
 750     const SGPropertyNode_ptr ai = fgGetNode("/ai/models", true);
 751
 752     for (int i = ai->nChildren() - 1; i >= -1; i--) {
 753         SGPropertyNode_ptr model;
 754
 755         if (i < 0) {
 756             // Last iteration: selected model
 757             model = _selected_ac;
 758         }
 759         else {
 760             model = ai->getChild(i);
 761             int parent_id = model->getIntValue("id");
 762             if (!model->nChildren()) {
 763                 continue;
 764             }
 765             if (parent_id == id) {
 766                 // Save selected model for last iteration
 767                 _selected_ac = model;
 768                 break;
 769             }
 770         }
 771         if (!model)
 772             continue;
 773     }
 774
 775     if (_selected_ac != 0) {
 776         const string name  = _selected_ac->getStringValue("name");
 777         IC.lat             = _selected_ac->getDoubleValue("position/latitude-deg");
 778         IC.lon             = _selected_ac->getDoubleValue("position/longitude-deg");
 779         IC.alt             = _selected_ac->getDoubleValue("position/altitude-ft");
 780         IC.roll            = _selected_ac->getDoubleValue("orientation/roll-deg");
 781         IC.elevation       = _selected_ac->getDoubleValue("orientation/pitch-deg");
 782         IC.azimuth         = _selected_ac->getDoubleValue("orientation/true-heading-deg");
 783         IC.speed           = _selected_ac->getDoubleValue("velocities/true-airspeed-kt") * SG_KT_TO_FPS;
 784         IC.speed_down_fps  = -_selected_ac->getDoubleValue("velocities/vertical-speed-fps");
 785         IC.speed_east_fps  = _selected_ac->getDoubleValue("velocities/speed-east-fps");
 786         IC.speed_north_fps = _selected_ac->getDoubleValue("velocities/speed-north-fps");
 787     }
 788     else {
 789         SG_LOG(SG_AI, SG_ALERT, "AISubmodel: parent node not found ");
 790     }
 791 }