Rev 909 | Rev 975 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 909 | Rev 949 | ||
---|---|---|---|
Line 879... | Line 879... | ||
879 | {
|
879 | {
|
880 | double hd_low, hd_high; |
880 | double hd_low, hd_high; |
881 | mSpecies->hdRange(mDbh, hd_low, hd_high); |
881 | mSpecies->hdRange(mDbh, hd_low, hd_high); |
882 | 882 | ||
883 | DBG_IF_X(hd_low>hd_high, "Tree::relative_height_growth", "hd low higher dann hd_high for ", dump()); |
883 | DBG_IF_X(hd_low>hd_high, "Tree::relative_height_growth", "hd low higher dann hd_high for ", dump()); |
884 | DBG_IF_X(hd_low < |
884 | DBG_IF_X(hd_low < 10 || hd_high>250, "Tree::relative_height_growth", "hd out of range ", dump() + QString(" hd-low: %1 hd-high: %2").arg(hd_low).arg(hd_high)); |
885 | 885 | ||
886 | // scale according to LRI: if receiving much light (LRI=1), the result is hd_low (for open grown trees)
|
886 | // scale according to LRI: if receiving much light (LRI=1), the result is hd_low (for open grown trees)
|
887 | // use the corrected LRI (see tracker#11)
|
887 | // use the corrected LRI (see tracker#11)
|
888 | double lri = limit(mLRI * mRU->LRImodifier(),0.,1.); |
888 | double lri = limit(mLRI * mRU->LRImodifier(),0.,1.); |
889 | double hd_ratio = hd_high - (hd_high-hd_low)*lri; |
889 | double hd_ratio = hd_high - (hd_high-hd_low)*lri; |