& objf_swdownm, objf_swdownsmoo,
& objf_swdown, objf_swdownm, objf_swdownsmoo,
_RL objf_swdownsmoo(nsx,nsy)
& +mult_smooth_bc*objf_swdownsmoo(bi,bj) )
f_swdownsmoo = f_swdownsmoo + objf_swdownsmoo(bi,bj)