{

/*

const double maxT = 365.; // maximal simulation time

const SeedSpecificParameter* rs = this->param(); // rename

Vector3d iheading(0,0,-1);

// Taproot

Root* taproot = new Root(plant, 1, iheading ,0, nullptr, 0, 0); // tap root has root type 1

taproot->addNode(rs->seedPos,0);

this->addChild(taproot);

auto pmap = plant->getOrganRandomParameter(Organism::ot_root);

// Basal roots

int bt = getParamSubType(Organism::ot_root, "basal");

if (bt>0) {

basalType = bt;

} // otherwise stick with default

if (rs->maxB>0) {

try {

plant->getOrganRandomParameter(Organism::ot_root, basalType); // if the type is not defined an exception is thrown

} catch (...) {

std::cout << "Seed::initialize: Basal root type #" << basalType << " was not defined, using tap root parameters instead\n" << std::flush;

RootRandomParameter* brtp = (RootRandomParameter*)plant->getOrganRandomParameter(Organism::ot_root, 1)->copy(plant);

brtp->subType = basalType;

plant->setOrganRandomParameter(brtp);

}

int maxB = rs->maxB;

if (rs->delayB>0) { // limit if possible

maxB = std::min(maxB,int(ceil((maxT-rs->firstB)/rs->delayB))); // maximal for simtime maxT

}

double delay = rs->firstB;

for (int i=0; i<maxB; i++) {

Root* basalroot = new Root(plant, basalType, iheading, delay, nullptr, 0, 0);

basalroot->addNode(taproot->getNode(0), taproot->getNodeId(0), delay);

this->addChild(basalroot);

delay += rs->delayB;

}

}

// Shoot borne roots

int st = getParamSubType(Organism::ot_root, "shootborne");

if (st>0) {

shootborneType = st;

} // otherwise stick with default

if ((rs->nC>0) && (rs->firstSB+rs->delaySB<maxT)) { // only if there are any shootborne roots

try {

plant->getOrganRandomParameter(Organism::ot_root, shootborneType); // if the type is not defined an exception is thrown

} catch (...) {

std::cout << "Seed::initialize:Shootborne root type #" << shootborneType << " was not defined, using tap root parameters instead\n";

RootRandomParameter* srtp = (RootRandomParameter*)plant->getOrganRandomParameter(Organism::ot_root, 1)->copy(plant);

srtp->subType = shootborneType;

plant->setOrganRandomParameter(srtp);

}

Vector3d sbpos = rs->seedPos;

sbpos.z=sbpos.z/2.; // half way up the mesocotyl

numberOfRootCrowns = ceil((maxT-rs->firstSB)/rs->delayRC); // maximal number of root crowns

double delay = rs->firstSB;

for (int i=0; i<numberOfRootCrowns; i++) {

Root* shootborne0 = new Root(plant, shootborneType, iheading ,delay, nullptr, 0, 0);

// TODO fix the initial radial heading

shootborne0->addNode(sbpos,delay);

this->addChild(shootborne0);

delay += rs->delaySB;

for (int j=1; j<rs->nC; j++) {

Root* shootborne = new Root(plant, shootborneType, iheading ,delay, nullptr, 0, 0);

// TODO fix the initial radial heading

shootborne->addNode(shootborne0->getNode(0), shootborne0->getNodeId(0),delay);

this->addChild(shootborne);

delay += rs->delaySB;

}

sbpos.z+=rs->nz; // move up, for next root crown

delay = rs->firstSB + i*rs->delayRC; // reset age

}

} else {

numberOfRootCrowns = 0;

}

/*

{

std::vector<Organ*> organs;

for (auto& o : children) {

organs.push_back(o->copy(plant));

}

return organs;

{

auto orp = plant->getOrganRandomParameter(organtype);

for (auto& o:orp) {

if (o->name == str) {

return o->subType;

}

}

return -1;

{

std::stringstream str;

str << "Seed #"<< id <<": type "<< param_->subType << ", length: "<< length << ", age: " << age;

return str.str();