aboutsummaryrefslogtreecommitdiffstats
path: root/mass_profile/analyze_fx.py
diff options
context:
space:
mode:
Diffstat (limited to 'mass_profile/analyze_fx.py')
-rwxr-xr-xmass_profile/analyze_fx.py32
1 files changed, 0 insertions, 32 deletions
diff --git a/mass_profile/analyze_fx.py b/mass_profile/analyze_fx.py
deleted file mode 100755
index 4f72f57..0000000
--- a/mass_profile/analyze_fx.py
+++ /dev/null
@@ -1,32 +0,0 @@
-#!/usr/bin/env python
-
-import sys
-import math
-import numpy
-
-fx1_array=[]
-fx2_array=[]
-fx3_array=[]
-for i in open('summary_fx.dat'):
- x1,x2,x3=i.split()
- x1=float(x1)
- x2=float(x2)
- x3=float(x3)
- fx1_array.append(x1)
- fx2_array.append(x2)
- fx3_array.append(x3)
-
-
-fx1_array=numpy.array(fx1_array)
-fx2_array=numpy.array(fx2_array)
-fx3_array=numpy.array(fx3_array)
-
-
-f=open('fx_result.txt','w')
-f.write("Fx(bolot)= %4.2E +/- %4.2E erg/s/cm^2\n"%(fx1_array[0],fx1_array.std()))
-print("Fx(bolot)= %4.2E +/- %4.2E erg/s/cm^2"%(fx1_array[0],fx1_array.std()))
-f.write("Fx(0.7-7)= %4.2E +/- %4.2E erg/s/cm^2\n"%(fx2_array[0],fx2_array.std()))
-print("Fx(0.7-7)= %4.2E +/- %4.2E erg/s/cm^2"%(fx2_array[0],fx2_array.std()))
-f.write("Fx(0.1-2.4)= %4.2E +/- %4.2E erg/s/cm^2\n"%(fx3_array[0],fx3_array.std()))
-print("Fx(0.1-2.4)= %4.2E +/- %4.2E erg/s/cm^2"%(fx3_array[0],fx3_array.std()))
-