function [DX_11_DX] = diff_scheme_prec( C,h,N_grid ) %,DZ_33_DZ,DX_12_DZ,DZ_33_DX,DZ_22_DZ,DX_33_DX,DZ_12_DX,DX_33_DZ DX_11_DX = spdiags(DPHX*c_11(:),0,N_grid^2,N ...

% (the shallower pressure), at the same values of SA and CT. This % function uses the computationally-efficient expression for specific % volume in terms of SA, CT and p (Roquet et al., 2015). The ...