Example #1
0
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)
Example #2
0
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%",