objf_theta_ini_fin(bi,bj) = 0. _d 0
& 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)