jhan_JaC_Redux_JAC_5.7_Radmaster_albedo

cbl_albedo

CanopyTransmit_dif(:,:) = 0.0

CanopyTransmit_beam(:,:) = 0.0

CanopyRefl_dif(:,:) = 0.0

CanopyRefl_beam(:,:) = 0.0

AlbSnow(:,:) = 0.0

!Modify parametrised soil albedo based on snow coverage 

call surface_albedosn( .... )

! Define canopy Reflectance for diffuse/direct radiation
! Formerly rad%rhocbm, rad%rhocdf

call CanopyReflectance( .... )

! Define canopy diffuse transmittance 
! Formerly rad%cexpkbm, rad%cexpkdm

call CanopyTransmitance( .... )

!---1 # visible, 2 nir radiaition
! Finally compute Effective 4-band albedo for diffuse/direct radiation- 
! In the UM this is the required variable to be passed back on the rad call
! Formerly rad%reffbm, rad%reffdf

! Even when there is no vegetation, albedo is at least snow modified soil albedo

EffSurfRefl_dif = AlbSnow

EffSurfRefl_beam = AlbSnow

call EffectiveSurfaceReflectance( .... )

! Compute total albedo to SW given the Effective Surface Reflectance 
! (considering Canopy/Soil/Snow contributions) 
! we dont need to do this on rad call AND may not haveappropriate RadFbeam

RadAlbedo = AlbSnow if(.NOT. jls_radiation) & call FbeamRadAlbedo( .... )