jhan_JaC_Redux_JAC_5.7_Radmaster_CanopyReflectance

canopy reflectance

Computed separately for direct beam/diffuse components

! Canopy reflection (6.21) beam:

DO i = 1,mp DO b = 1, 2 IF( sunlit_veg_mask(i) ) & CanopyRefl_beam(i,b) = 2. * ExtCoeff_beam(i) / & ( ExtCoeff_beam(i) + ExtCoeff_dif(i) ) * rhoch(i,b) END DO END DO

! Canopy REFLection of diffuse radiation for black leaves:

DO ictr=1,nrb

CanopyRefl_dif(:,ictr) = rhoch(:,ictr) * 2. * ( CGAUSS_W(1) * xk(:,1) / ( xk(:,1) + ExtCoeff_dif(:) )& + CGAUSS_W(2) * xk(:,2) / ( xk(:,2) + ExtCoeff_dif(:) )& + CGAUSS_W(3) * xk(:,3) / ( xk(:,3) + ExtCoeff_dif(:) ) )

ENDDO