objf_swflux(bi,bj) = 0. _d 0
& objf_swflux,
& objf_swflux, objf_swfluxm, objf_swfluxsmoo,
_RL objf_swflux(nsx,nsy)
& + mult_swflux * ( objf_swflux(bi,bj)+objf_swfluxm(bi,bj)
f_swflux = f_swflux + objf_swflux(bi,bj)