objf_theta_ini_fin(bi,bj) =
& objf_theta_ini_fin(bi,bj) + fctile
& objf_theta_ini_fin, objf_salt_ini_fin,
_RL objf_theta_ini_fin(nsx,nsy)
& + mult_ini_fin *(objf_theta_ini_fin(bi,bj) +
& objf_theta_ini_fin(bi,bj) + objf_salt_ini_fin(bi,bj)
objf_theta_ini_fin(bi,bj) = 0. _d 0