Subversion Repositories public iLand

Rev

Rev 909 | Rev 975 | Go to most recent revision | Show entire file | Ignore 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 < 20 || 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));


            -
            
 

          

          

            -
            
 

            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;