& objf_tauvm, objf_tauvsmoo,
& objf_tauv, objf_tauvm, objf_tauvsmoo,
_RL objf_tauvsmoo (nsx,nsy)
& +mult_smooth_bc*objf_tauvsmoo(bi,bj) )
f_tauvsmoo = f_tauvsmoo + objf_tauvsmoo(bi,bj)