sw_two_stream: R implementation of ED2 two-stream radiation model
In ashiklom/edr-da: ED2 RTM testing functions

R implementation of ED2 two-stream radiation model

sw_two_stream(
  czen,
  iota_g,
  pft,
  lai,
  wai,
  cai,
  orient_factor,
  clumping_factor,
  leaf_reflect,
  leaf_trans,
  wood_reflect,
  wood_trans,
  down_sky,
  down0_sky,
  wavelengths = seq(400, 2500),
  ssa_type = "sellers"
)

`czen`	Cosine of angle of incidence ('cosaio')
`iota_g`	Ground (soil + snow) albedo (nwl)
`pft`	PFT identities of each cohort, as integer (ncoh). Shortest to tallest.
`lai`	Leaf area index of each cohort (ncoh). Shorest to tallest.
`wai`	Wood area index of each cohort (ncoh). Shorest to tallest.
`cai`	Crown area of each cohort (ncoh). Shorest to tallest.
`orient_factor`	Orient factor (npft)
`clumping_factor`	Clumping factor (npft)
`leaf_reflect`	Leaf reflectance spectra (nwl * npft)
`leaf_trans`	Leaf transmittance spectra (nwl * npft)
`wood_reflect`	Wood reflectance spectra (nwl * npft)
`wood_trans`	Wood transmittance spectra (nwl * npft)
`down_sky`	Normalized diffuse solar spectrum (nwl)
`down0_sky`	Normalized direct solar spectrum (nwl)
`wavelengths`	Numeric vector of wavelengths to use, in nm (nwl). Default is 400:2500.