jhan_radrhodf

Canopy Reflectance (diffuse) [CanopyRefl_dif]

rad%rhocdf

This is currently in cable_init_radiation() - pull this out as a separate F().

CanopyRefl_dif = F( rhoch, C%GAUSS_W, xk, ExtCoeff_dif )

Where ExtCoeff_dif = Extinction Coefficient (diffuse) and C%Gauss_w is a constant that needs to be imported

In terms of the formerly named variables:

rad%rhocdf = F( rhoch, C%GAUSS_W, xk, rad%extkd )