objf_ctdsclim(bi,bj) = objf_ctdsclim(bi,bj) + fctile
& objf_ctdtclim, objf_ctdsclim,
_RL objf_ctdsclim (nsx,nsy)
& + mult_ctdsclim* objf_ctdsclim(bi,bj)
f_ctdsclim = f_ctdsclim + objf_ctdsclim(bi,bj)
objf_ctdsclim(bi,bj) = 0. _d 0