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