objf_lwflux(bi,bj) = 0. _d 0
& objf_lwflux,
& objf_lwflux, objf_lwfluxm, objf_lwfluxsmoo,
_RL objf_lwflux(nsx,nsy)
& + mult_lwflux * ( objf_lwflux(bi,bj)+objf_lwfluxm(bi,bj)
f_lwflux = f_lwflux + objf_lwflux(bi,bj)