& objf_swfluxm, objf_swfluxsmoo,
& objf_swflux, objf_swfluxm, objf_swfluxsmoo,
_RL objf_swfluxsmoo(nsx,nsy)
& +mult_smooth_bc*objf_swfluxsmoo(bi,bj) )
f_swfluxsmoo = f_swfluxsmoo + objf_swfluxsmoo(bi,bj)