import sys Lib_path = "/global/home/users/kaylanb/myCode/ytScripts/Lib" sys.path.append(Lib_path) import deriv_fields import perturbation_lib as pl import init_cond as mic path = "/clusterfs/henyey/kaylanb/Test_Prob/Turb/bondiTurb/quick/driving_t0_only/HD/lev2/" p = mic.ParsNeed() # then set pars.[pars] to new vals as needed ic = mic.calcIC(p) ic.prout() # pl.perturbation(64,1,2,0.667,2.) # pl.plot_Pk_and_Hist_rmsV_all_simData(path, ic.v0, ic.v0/8., ic.tCr) pl.hist_Rho_all_data(path, 2, ic.rho0, ic.tCr)
from yt.mods import * import matplotlib.pyplot as plt from mpl_toolkits.axes_grid1 import AxesGrid import sys sys.path.append("/global/home/users/kaylanb/myCode/ytScripts/Lib") import my_fund_const as fc import init_cond as mic import deriv_fields pars= mic.ParsNeed() #then set pars.[pars] to new vals as needed ic = mic.calcIC(pars) ic.prout() ###### def onefile_multi_plot(file): pf= load(path+file) #units #myu = "r [rB]" #pf.units[myu] = 1./ic.rB #rB = ic.rB bondi radii, so 1 cm = 1/ic.rB centimeters! #plot fig = plt.figure() grid = AxesGrid(fig, (0.075,0.075,0.85,0.85), nrows_ncols = (2,2), axes_pad =1.5, label_mode = "all", share_all = False, cbar_location="right", cbar_mode="each", cbar_size="3%",