rac = rb / (2L)
rsc = rST / (2L)
Δ = (0.00815 Ta + 0.8912)7
Rns = Rn e-CL
A = Rn - 0.2Rns
As = 0.8Rns
Ra = (Δ + γ) raa
Rs = (Δ + γ) ras + rss
Rc = (Δ + γ) rac + rsc
Cc = 1 / { 1 + (Rc Ra) / [Rs (Rc + Ra) ] }
Cs = 1 / { 1 + (Rs Ra) / [Rc (Rs + Ra) ] }
PMc = [ΔA + (ρcpD - ΔracAs) / (raa + rac)]
/ [ Δ + γ {1 + rsc/(raa + rac)}]
PMs = [ΔA + {ρcpD - Δras(A - As)} / (raa + ras)]
/ [ Δ + γ {1 + rss/(raa + ras)}]
