& objf_uwindm, objf_uwindsmoo,
& objf_uwind, objf_uwindm, objf_uwindsmoo,
_RL objf_uwindsmoo(nsx,nsy)
& +mult_smooth_bc*objf_uwindsmoo(bi,bj) )
f_uwindsmoo = f_uwindsmoo + objf_uwindsmoo(bi,bj)