From 62b0ced60d81bf214fdf9d03a593a17765f67f81 Mon Sep 17 00:00:00 2001
From: Noe Brucy <noe.brucy@cea.fr>
Date: Thu, 28 Nov 2024 11:44:28 +0100
Subject: [PATCH] add velr plot

---
 pipeline_MW.py | 7 +++++++
 1 file changed, 7 insertions(+)

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"]