From 9cff25b987ea25cdc96d4e2c06fb2c37cbf79a14 Mon Sep 17 00:00:00 2001 From: mccabech Date: Sat, 22 Sep 2018 11:18:17 +0100 Subject: [PATCH] Update SHMpp_hvmin.py --- SHMpp_hvmin.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/SHMpp_hvmin.py b/SHMpp_hvmin.py index 0f93204..7e0a443 100644 --- a/SHMpp_hvmin.py +++ b/SHMpp_hvmin.py @@ -157,14 +157,14 @@ def vE(Day,v0): ########################################## # Round component # Calculate the h(vmin) integral for fR(v) -# output table of vmin [km/s] and h(vmin) [s/km] +# output table of vmin [km/s] and h(vmin) [km/s] ######################################### #weight function is exp(-(x2+y2+z2)) #integral performed in the galactic frame print("Calculating hvmin for Round component") -print("vmin [km/s],","h(vmin) [s/km]") +print("vmin [km/s],","h(vmin) [km/s]") ve = vE(Day,v0) @@ -213,14 +213,14 @@ def vE(Day,v0): ########################################## # Sausage component # Calculate the h(vmin) integral for fS(v) -# output table of vmin [km/s] and h(vmin) [s/km] +# output table of vmin [km/s] and h(vmin) [km/s] ######################################### #weight function is exp(-(x2+y2+z2)) #integral performed in the galactic frame print("Calculating hvmin for Sausage component") -print("vmin [km/s],","h(vmin) [s/km]") +print("vmin [km/s],","h(vmin) [km/s]") ve = vE(Day,v0) @@ -264,7 +264,7 @@ def vE(Day,v0): ########################################## # Combine round and Sausage components for SHM++ -# output table of vmin [km/s] and h(vmin) [s/km] +# output table of vmin [km/s] and h(vmin) [km/s] ######################################### print("Calculating hvmin for SHM++ ( with eta=",eta,")") @@ -278,5 +278,5 @@ def vE(Day,v0): plot(Lvx, LhvminSHMpp,"k", label='SHM++ with eta='+str(eta)) legend(loc='upper right') xlabel("vmin [km/s]") -ylabel("h(vmin) [s/km]") +ylabel("h(vmin) [km/s]") show()