re-write depth_ntp_iter_drho:
depth_ntp_simple calls gsw_rho in nested for-loops
less evaluations of gsw_rho, much faster