def calculate_velocity(mom1,mom2,mom3,dens):
    """Calculates the 3-vector velocity from 3-vector momentum density and scalar density."""
    v1 = mom1/dens
    v2 = mom2/dens
    v3 = mom3/dens
    return v1,v2,v3