& objf_lwdownm, objf_lwdownsmoo,
& objf_lwdown, objf_lwdownm, objf_lwdownsmoo,
_RL objf_lwdownsmoo(nsx,nsy)
& +mult_smooth_bc*objf_lwdownsmoo(bi,bj) )
f_lwdownsmoo = f_lwdownsmoo + objf_lwdownsmoo(bi,bj)