& objf_hfluxm, objf_hfluxsmoo,
& objf_hflux, objf_hfluxm, objf_hfluxmm, objf_hfluxsmoo,
_RL objf_hfluxsmoo (nsx,nsy)
& +mult_smooth_bc*objf_hfluxsmoo(bi,bj) )
f_hfluxsmoo = f_hfluxsmoo + objf_hfluxsmoo(bi,bj)