jhan_JaC_Redux_JAC_5.7_Radmaster_CanopyTransmitance

Canopy Transmitance

Poetic license assumed for succinctness.

! For diffuse rad, always compute canopy trasmitance mask(:) = .true.
! For beam, compute canopy trasmitance when mask(:) = sunlit (and vegetated)

DO i = 1,mp DO b = 1, nrb if( mask(i) ) then dummy(i,b) = min( ExtinctionCoeff(i,b) * reducedLAIdue2snow(i), 20. )

  CanopyTransmit(i,b) = EXP( -1.* dummy(i,b) )

endif

enddo enddo