def all_z(maindir='/oasis/projects/nsf/hvd115/lblecha',
subdir_arr='fid',
min_z=0.1,
max_z=0.5,
dz=0.1,
agnscale=False):
basepath_arr = cd.define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
basepath_arr = [bp.split('/')[0] for bp in basepath_arr]
print basepath_arr
for i_path, bp in enumerate(basepath_arr):
tmax = 4.2 if 'q0.2_fg0.3_sunruns' in bp else 2.95
for z in np.arange(min_z, max_z + dz, dz):
print bp + '/z%.1f' % z
print bp + '/z%.1f_agnx0' % z
cd.color_color_compare(subdirA=bp + '/z%.1f' % z,
subdirB=bp + '/z%.1f_agnx0' % z,
tmax=tmax,
pubstyle=False,
agnscale=agnscale)
def min_w2w3(maindir='/oasis/projects/nsf/hvd115/lblecha',
subdir_arr='fid',
z0=True,
min_z=0.1,
max_z=0.9,
dz=0.1,
total_only=True):
if z0:
path_arr = cd.define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
else:
basepath_arr = cd.define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
basepath_arr = [bp.split('/')[0] for bp in basepath_arr]
path_arr = []
for i_path, bp in enumerate(basepath_arr):
path_arr = path_arr + [
bp + '/z%.1f' % z for z in np.arange(min_z, max_z + dz, dz)
]
min_w2w3 = 1e6
min_w2w3_above = 1e6
for i_path, p in enumerate(path_arr):
time, snap, w1w2, w2w3 = cd.load_wise(maindir=maindir, subdir=p)
if not total_only:
print "\n" + p
print "min w2w3: %g" % w2w3.min()
if w1w2.max() > 0.5:
if not total_only:
print "min w2w3[w1w2>0.5]: %g" % w2w3[w1w2 > 0.5].min()
if w2w3[w1w2 > 0.5].min() < min_w2w3_above:
min_w2w3_above = w2w3[w1w2 > 0.5].min()
else:
if not total_only: print "no w1w2>0.5."
if w2w3.min() < min_w2w3: min_w2w3 = w2w3.min()
print "\nmin w2w3 for all sims: %g" % min_w2w3
print "min w2w3[w1w2>0.5] for all sims: %g" % min_w2w3_above
def all_z0_agnscale(maindir='/oasis/projects/nsf/hvd115/lblecha',
subdir_arr='fid',
pubstyle=True):
path_arr = cd.define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
for i_path, p in enumerate(path_arr):
cd.color_color_compare(subdirA=p,
subdirB=p + '_agnx0',
pubstyle=pubstyle)
def all_z0_fid(maindir='/oasis/projects/nsf/hvd115/lblecha',
subdir_arr='fid',
pubstyle=True):
path_arr = define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
for i_path, p in enumerate(path_arr):
print p
wa.wise_colors(path=maindir,
subdir=p,
write_llambda=True,
input_type='bb')
def onesim_z_agnx0(maindir='/oasis/projects/nsf/hvd115/lblecha',
subdir_arr='q1_fg0.3_sunruns',
min_z=1.5,
max_z=4.5,
dz=0.5):
basepath = define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
print basepath
for z in np.arange(min_z, max_z + dz, dz):
wa.wise_colors(path=maindir,
subdir=basepath[0] + '/z%.1f_agnx0' % z,
write_llambda=True,
input_type='bb')
def all_z_fid_agnx0(maindir='/oasis/projects/nsf/hvd115/lblecha',
subdir_arr='fid',
min_z=0.1,
max_z=1.0,
dz=0.1):
basepath_arr = define_sunrise_path_arr(basepath=maindir,
name=subdir_arr,
sfrhist_w_grid=False,
return_full_path=False)
basepath_arr = [bp.split('/')[0] for bp in basepath_arr]
print basepath_arr
for i_path, bp in enumerate(basepath_arr):
for z in np.arange(min_z, max_z + dz, dz):
print bp + '/z%.1f_agnx0' % z
wa.wise_colors(path=maindir,
subdir=bp + '/z%.1f_agnx0' % z,
write_llambda=True,
input_type='bb')
import column_density as cd fidres_grid_res = 0.136 ## 0.136 kpc = 0.034 h^-1 kpc / 0.7 * 2.8 (when grav smoothing kernal becomes newtonian) hires_grid_res = 0.064 ## 0.064 kpc = 0.016 h^-1 kpc / 0.7 * 2.8 (when grav smoothing kernal becomes newtonian) fid_arr = cd.define_sunrise_path_arr(name='fid', return_full_path=False) hires_arr = cd.define_sunrise_path_arr(name='hires', return_full_path=False) arr = fid_arr + hires_arr ##arr=hires_arr #print "\nmaking single-sim 'true Ngas' plots for these runs:" #print arr ##for subdir in arr: ## print "\n",subdir ## cd.merger_phases(subdir_arr=subdir,sfrhist_w_grid=True,single_sim_plots=False, ## equal_bins=False,ngas_only=True,tmax_postmrg=0.25,grid_nh_res=0) #for subdir in fid_arr: # print "\n",subdir # cd.merger_phases(subdir_arr=subdir,sfrhist_w_grid=True,single_sim_plots=False, # equal_bins=False,ngas_only=True,tmax_postmrg=0.25,grid_nh_res=fidres_grid_res) #for subdir in hires_arr: # print "\n",subdir # cd.merger_phases(subdir_arr=subdir,sfrhist_w_grid=True,single_sim_plots=False, # equal_bins=False,ngas_only=True,tmax_postmrg=0.25,grid_nh_res=fidres_grid_res) # print "\n",subdir # cd.merger_phases(subdir_arr=subdir,sfrhist_w_grid=True,single_sim_plots=False, # equal_bins=False,ngas_only=True,tmax_postmrg=0.25,grid_nh_res=hires_grid_res) name_arr = [