diff --git a/pipeline_MW.py b/pipeline_MW.py
index cfd0d58..4337522 100644
--- a/pipeline_MW.py
+++ b/pipeline_MW.py
@@ -57,6 +57,13 @@ def ring_stuff(galsec, delta_r=1):
     plt.xlabel("R [kpc]")
     plt.savefig("rotation_curve.png")
     
+    plt.figure(constrained_layout=True, figsize=(5,4))
+    velr = galsec.rings['gas']["velr"]
+    plt.plot(r_gas, velr)
+    plt.ylabel(r"$v_r$ [km/s]")
+    plt.xlabel("R [kpc]")
+    plt.savefig("vr.png")
+    
     
     plt.figure(constrained_layout=True, figsize=(5,4))
     sigma_velphi = galsec.rings['gas']["sigma_velphi"]