& objf_tauum, objf_tauusmoo,
& objf_tauu, objf_tauum, objf_tauusmoo,
_RL objf_tauusmoo (nsx,nsy)
& +mult_smooth_bc*objf_tauusmoo(bi,bj) )
f_tauusmoo = f_tauusmoo + objf_tauusmoo(bi,bj)