The neutral density ( ) or empirical neutral density is a density variable used in oceanography, introduced in 1997 by David R. Jackett and Trevor McDougall.[1] It is a function of the three state variables (salinity, temperature, and pressure) and the geographical location (longitude and latitude). It has the typical units of density (M/V). Isosurfaces of form “neutral density surfaces”, which are closely aligned with the "neutral tangent plane". It is widely believed, although this has yet to be rigorously proven, that the flow in the deep ocean is almost entirely aligned with the neutral tangent plane, and strong lateral mixing occurs along this plane ("epineutral mixing") vs weak mixing across this plane ("dianeutral mixing"). These surfaces are widely used in water mass analyses. Neutral density is a density variable that depends on the particular state of the ocean, and hence is also a function of time, though this is often ignored. In practice, its construction from a given hydrographic dataset is achieved by means of a computational code (available for Matlab and Fortran), that contains the computational algorithm developed by Jackett and McDougall. Use of this code is currently restricted to the present day ocean.