objf_lwdown(bi,bj) = 0. _d 0
& objf_lwdown,
& objf_lwdown, objf_lwdownm, objf_lwdownsmoo,
_RL objf_lwdown(nsx,nsy)
& + mult_lwdown * ( objf_lwdown(bi,bj)+objf_lwdownm(bi,bj)
f_lwdown = f_lwdown + objf_lwdown(bi,bj)