plot(ctd, which="density+N2") obeys type="o"