objf_diffkr(bi,bj) = 0. _d 0
& num_diffkr, objf_diffkr, maskC, myThid )
& objf_diffkr,
_RL objf_diffkr(nsx,nsy)
& + mult_diffkr * objf_diffkr(bi,bj)
f_diffkr = f_diffkr + objf_diffkr(bi,bj)