Lai july2022

[rcnijzink]

This branch had my latest edits for the LAI in the VOM (just before EGU). I think it would be best to merge this into a branch on git@github.com:schymans/VOM.git, then you could start from there. I tried below to give an overview on what I did.

This branch includes a first attempt to separate PAR into PAR shade and PAR sunlit, based on the approach in Campbell & Norman. I followed the example on p260 (example 15.2). Specifically:

• estimate average PAR on shaded leaves (transpmode.f90, lines 1685 and 1713, first equation p261 Campbell & Norman):

                    par_shade_avg = (pardiff_h(th_) * ( 1.d0 - p_E ** ( -sqrt(i_alpha_abs) * kappa_d * lai_lt(ii) )   ) ) / &
                                    ( sqrt(i_alpha_abs) * kappa_d* lai_lt(ii)    ) 

• determine scattered PAR (transpmode.f90, lines 1689 and 1717)::

                    par_sct = 0.5* ( pardir_h(th_) * p_E ** (-sqrt(i_alpha_abs) * kappa_t * lai_lt(ii) )  -  &
                                    pardir_h(th_) * p_E ** (-kappa_t * lai_lt(ii) ) )

the above is the same as the calculation on p261 of Campbell & Norman, with scattered PAR as Qsc = 0.5 (Qbt - Qb), with Qbt the top beam radiation and Qb bottom beam radiation.

• Eventually, PARshade and PARsun are determined like this:

PARshade = (par_sct + par_shade_avg) * i_alpha_abs

PARsun = i_alpha_abs * (kappa_t * pardir_h(th_) + PARshade )

On p261 of Campbell & Norman, these are the equations for Qsh and Qsl.

[schymans]

This has been merged in the LAI branch and a new PR has been created: #69
Thanks a lot, @rcnijzink!