jhan_temp_common_coeffs

subroutine common_InitRad_coeffs(...)

...

...

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(canopyVlaiw, 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

End subroutine common_InitRad_coeffs

Previous