Rev 165 | Rev 169 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 165 | Rev 167 | ||
|---|---|---|---|
| Line 524... | Line 524... | ||
| 524 | if (mFoliageMass<0.) mFoliageMass=0.; // limit to zero |
524 | if (mFoliageMass<0.) mFoliageMass=0.; // limit to zero |
| 525 | 525 | ||
| 526 | mLeafArea = mFoliageMass * species()->specificLeafArea(); // update leaf area |
526 | mLeafArea = mFoliageMass * species()->specificLeafArea(); // update leaf area |
| 527 | 527 | ||
| 528 | // stress index
|
528 | // stress index
|
| 529 | d.stress_index =qMax(1. - npp |
529 | d.stress_index =qMax(1. - (npp-reserve_size) / to_fol*foliage_mass_allo, 0.); |
| 530 | 530 | ||
| 531 | // Woody compartments
|
531 | // Woody compartments
|
| 532 | // (1) transfer to reserve pool
|
532 | // (1) transfer to reserve pool
|
| 533 | double gross_woody = apct_wood * npp; |
533 | double gross_woody = apct_wood * npp; |
| 534 | double to_reserve = qMin(reserve_size, gross_woody); |
534 | double to_reserve = qMin(reserve_size, gross_woody); |
| Line 654... | Line 654... | ||
| 654 | // death if leaf area is 0
|
654 | // death if leaf area is 0
|
| 655 | if (mFoliageMass<0.00001) |
655 | if (mFoliageMass<0.00001) |
| 656 | die(); |
656 | die(); |
| 657 | 657 | ||
| 658 | double p_death, p_stress; |
658 | double p_death, p_stress; |
| 659 | p_stress |
659 | //p_stress = d.stress_index * species()->deathProb_stress();
|
| - | 660 | if (d.stress_index>0) |
|
| - | 661 | p_stress = species()->deathProb_stress(); |
|
| 660 | p_death = species()->deathProb_intrinsic() + p_stress; |
662 | p_death = species()->deathProb_intrinsic() + p_stress; |
| 661 | double p = random(); //0..1 |
663 | double p = random(); //0..1 |
| 662 | if (p<p_death) { |
664 | if (p<p_death) { |
| 663 | // die...
|
665 | // die...
|
| 664 | die(); |
666 | die(); |