Wrong unit shown in `workflow` plots

@belletva @marinebcht While revisiting this today, I realized that `crs.linear_units` is the horizontal units, not the vertical. The function cannot allow to differentiate vertical/horizontal axes. We need to use `axis.unit_name` after isolating the vertical axis (if it exists).

I'm fixing this in #656 with the 3D CRS override.