FUNCTION z12dist,z1,z2

z12=(1+z2)/(1+z1)-1

return, zdist(z12)

end