& objf_sfluxm, objf_sfluxsmoo,
& objf_sflux, objf_sfluxm, objf_sfluxmm, objf_sfluxsmoo,
_RL objf_sfluxsmoo (nsx,nsy)
& +mult_smooth_bc*objf_sfluxsmoo(bi,bj) )
f_sfluxsmoo = f_sfluxsmoo + objf_sfluxsmoo(bi,bj)