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