objf_swdown(bi,bj) = 0. _d 0
& objf_swdown,
& objf_swdown, objf_swdownm, objf_swdownsmoo,
_RL objf_swdown(nsx,nsy)
& + mult_swdown * ( objf_swdown(bi,bj)+objf_swdownm(bi,bj)
f_swdown = f_swdown + objf_swdown(bi,bj)