# file-dependent variables.
istart = 33
istop = 34
OUTDIR = "../run/"
RM = 4
mp = 1.67e-24
kb = 1.3801e-16
USECENTERING = 0
for i in xrange(istart, istop):
number = repr(i).zfill(5)
filename = "cell" + number + ".dat"
f = open(filename, 'r')
header = TK.getRamsesInfo(i, OUTDIR + "./output_")
ncpu, ndim, levelmin, levelmax, ngridmax, nstep_coarse, boxlen, time, aexp, H0, omega_m, omega_l, omega_k, omega_b, unit_l, unit_d, unit_t = header
norm = unit_d / mp
dd = []
tk = []
h2 = []
l = 0
ddMax = 0.
offSet = 0.
for line in f:
values = line.split()
dd.append(float(values[6]))
tk.append(float(values[10]))
# first guess
xc = 0.547045824143
yc = 0.512265650452
zc = 0.458181241497
xc = 0.429377013232
yc = 0.396090123252
zc = 0.420656749876
xc = 0.534073102248
yc = 0.506075144996
zc = 0.467231232454
pc = 3.08e18
for i in xrange(fstart, fstop):
number = repr(i).zfill(5)
header = TK.getRamsesInfo(i, "../output_")
ncpu, ndim, levelmin, levelmax, ngridmax, nstep_coarse, boxlen, time, aexp, H0, omega_m, omega_l, omega_k, omega_b, unit_l, unit_d, unit_t = header
infile = DIR + "star" + number + ".dat"
f = open(infile, "r")
outfile = DIR + "star" + number + extension + ".png"
print "#open file ", f
mass = []
age = []
x = []
y = []
z = []
metal = []
xc_star = 0.
yc_star = 0.
zc_star = 0.
i = len(sys.argv)
if i < 5 or (i > 5 and i < 8):
print "Usage: get_profile.py filename input max_dist(kpc) value [x0 y0 z0]"
sys.exit()
f = open(sys.argv[1], "r")
index = 5 + int(sys.argv[4])
array = []
x = []
y = []
z = []
dx = []
header = TK.getRamsesInfo(int(sys.argv[2]), OUTDIR + "./output_")
ncpu, ndim, levelmin, levelmax, ngridmax, nstep_coarse, boxlen, time, aexp, H0, omega_m, omega_l, omega_k, omega_b, unit_l, unit_d, unit_t = header
array_conv = 1.
if index == 6:
array_conv = unit_d / mH
amax = 0.
i = 0
for line in f:
values = line.split()
x.append(float(values[0]) * unit_l / kpc)
y.append(float(values[1]) * unit_l / kpc)
z.append(float(values[2]) * unit_l / kpc)
dx.append(float(values[3]) * unit_l / kpc)
array.append(float(values[index]) * array_conv)
# file-dependent variables. istart=32 istop =33 RM=4 mp=1.67e-24 kb=1.3801e-16 USECENTERING=0 OUTDIR="/scratch/psteger/sim/" DIR="./" for i in xrange(istart,istop): number=repr(i).zfill(5) filename=DIR+"cell"+number+".dat" f = open(filename,'r') header=TK.getRamsesInfo(i,OUTDIR+"./output_") ncpu,ndim,levelmin,levelmax,ngridmax,nstep_coarse,boxlen,time,aexp,H0,omega_m,omega_l,omega_k,omega_b,unit_l,unit_d,unit_t=header norm=unit_d/mp polytropic_constant=( 4.*0.5**float(levelmax)*unit_l/aexp)**2/(math.pi)*6.67e-8*unit_d*(unit_t/unit_l)**2 polytropic_constant=0. dd=[] tk=[] metal=[] l=0 ddMax=0. offSet=0. xx=0.;yy=0.;zz=0. for line in f:
i = len(sys.argv)
if i < 5 or (i > 5 and i < 8):
print "Usage: get_profile.py filename input max_dist(kpc) value [x0 y0 z0]"
sys.exit()
f = open(sys.argv[1], "r")
index = 5 + int(sys.argv[4])
array = []
x = []
y = []
z = []
dx = []
header = TK.getRamsesInfo(int(sys.argv[2]), OUTDIR + "./output_")
ncpu, ndim, levelmin, levelmax, ngridmax, nstep_coarse, boxlen, time, aexp, H0, omega_m, omega_l, omega_k, omega_b, unit_l, unit_d, unit_t = (
header
)
array_conv = 1.0
if index == 6:
array_conv = unit_d / mH
amax = 0.0
i = 0
for line in f:
values = line.split()
x.append(float(values[0]) * unit_l / kpc)
y.append(float(values[1]) * unit_l / kpc)
z.append(float(values[2]) * unit_l / kpc)
ZMAX=0
ZMIN=-7.
DZ=.1
BINN=(ZMAX-ZMIN)/DZ
NBINN=int(BINN)
# initial guess for location.
# search length
# offset for plotting
extent=0.001
pc=3.08e18
for i in xrange(fstart,fstop):
number=repr(i).zfill(5)
header=TK.getRamsesInfo(i,"../run/output_")
ncpu,ndim,levelmin,levelmax,ngridmax,nstep_coarse,boxlen,time,aexp,H0,omega_m,omega_l,omega_k,omega_b,unit_l,unit_d,unit_t=header
infile=DIR+prefix+number+".dat"
f=open(infile,"r")
outfile=DIR+prefix+number+extension+".png"
print "#open file ",f
for j in xrange(8):
line=f.readline()
if line.find("time")>0:
key,val=line.split()
simtime=float(val)
if line.find("simu")>0:
key,val=line.split("=")
simage=float(val)*1e3 # 1e3 is to convert to Myr
print "#found simage of ", simage
mass=[]
