objf_precip(bi,bj) = 0. _d 0
& objf_precip,
& objf_precip, objf_precipm, objf_precipsmoo,
_RL objf_precip(nsx,nsy)
& + mult_precip * ( objf_precip(bi,bj)+objf_precipm(bi,bj)
f_precip = f_precip + objf_precip(bi,bj)