HeightVsDBH Function

public function HeightVsDBH(dbh, cra, crb, crc) result(height)

relationship between tree height and DBH (Diameter at Brest Height). Implements Chapman-Richards relationship

Reference:

Wang, C. Biomass allometric equations for 10 co-occurring tree species in Chinese temperate forests. Forest Ecology and Management, 222, 9–16, 2006

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: dbh

diameter at brest height (cm)
real(kind=float), intent(in) :: cra

Chapman-Richards asymptotic maximum height
real(kind=float), intent(in) :: crb

Chapman-Richards exponential decay parameter
real(kind=float), intent(in) :: crc

Chapman-Richards shape parameter

Return Value real(kind=float)

tree eight (m)


Source Code

FUNCTION  HeightVsDBH &
!
(dbh, cra, crb, crc) &
!
RESULT (height)


IMPLICIT NONE

! Arguments with intent(in):
REAL (KIND = float), INTENT(IN) :: dbh !!diameter at brest height (cm)
REAL (KIND = float), INTENT(IN) :: cra !!Chapman-Richards asymptotic maximum height
REAL (KIND = float), INTENT(IN) :: crb !!Chapman-Richards exponential decay parameter
REAL (KIND = float), INTENT(IN) :: crc !!Chapman-Richards shape parameter

!local declarations:
REAL (KIND = float) :: height !! tree eight (m)

!------------end of declaration------------------------------------------------ 

height = 1.3 + cra * ( 1. - EXP ( - dbh * crb) ) ** crc

RETURN
END FUNCTION HeightVsDBH