& mult_lwdown,
_RL mult_lwdown
& mult_lwflux, mult_lwdown, mult_evap,
mult_lwdown = 0. _d 0
& + mult_lwdown * ( objf_lwdown(bi,bj)+objf_lwdownm(bi,bj)