jhan_JaC_Redux_JAC_5.7_Radmaster_Common_InitRad_Scalings

Common_InitRad_Scalings

Very minor liberties have been taken for display purposes. Otherwise this like every other page is the exact code running.

cos3 = COS(CPI180 * (/ 15.0, 45.0, 75.0 /))

! See Sellers 1985, eq.13 (leaf angle parameters):

WHERE ( veg_mask ) xphi1 = 0.5 - VegXfang * (0.633 + 0.33 * VegXfang) xphi2 = 0.877 * (1.0 - 2.0 * xphi1) END WHERE

! 2 dimensional LAI

xvlai2 = SPREAD(reducedLAIdue2snow, 2, 3)

! Extinction coefficient for beam radiation and black leaves;
! eq. B6, Wang and Leuning, 1998

WHERE (xvlai2 > cLAI_THRESH) ! vegetated xk = SPREAD(xphi1, 2, 3) / SPREAD(cos3, 1, mp) + SPREAD(xphi2, 2, 3) ELSEWHERE ! i.e. bare soil xk = 0.0 END WHERE

CALL calc_rhoch( .... )