# 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=[]