& objf_lwfluxm, objf_lwfluxsmoo,
& objf_lwflux, objf_lwfluxm, objf_lwfluxsmoo,
_RL objf_lwfluxsmoo(nsx,nsy)
& +mult_smooth_bc*objf_lwfluxsmoo(bi,bj) )
f_lwfluxsmoo = f_lwfluxsmoo + objf_lwfluxsmoo(bi,bj)