def __init__(self,
scale_length=1,
bz_in=400.0,
lambda_p=5,
pert_amp='0p',
dim='2D'):
paths = hk.DataDirectoryPaths()
self.dim = dim
self.path = paths.get_path(scale_length,
bz_in,
lambda_p,
pert_amp=pert_amp,
dim=self.dim)
self.bz = bz_in
self.scale_length = scale_length
self.lambda_p = lambda_p
self.pert_amp = pert_amp
# ---> getting path ->
t_list, tc_list = cf.get_t_list(self.path, var='Te')
self.time_max = t_list[-1]
self.tag = self.get_tag()
self.save_tag = 'LT%i_%i_%s_%s' % (self.scale_length, self.bz,
self.pert_amp, self.time_max)
norm_dir = paths.norm_dir
log_file = norm_dir + 'norm.log'
[T_ref, n_ref, Z_ref, Bz_ref] = np.loadtxt(log_file)
self.norm = cf.ConversionFactors(norm_dir, Z_ref, Ar=6.51)
self.save_path = paths.save_dir
self.Te_factor = 2.0 * T_ref * 1e-3
python plotters/plot_initial_profile_lineouts.py output_folder_to_save_plot/ ''' import os import sys import matplotlib.pyplot as plt import matplotlib as mpl sys.path.extend(["./"]) import modules.house_keeping as hk import modules.chfoil_module as cf import modules.figure_latex as fprl paths = hk.DataDirectoryPaths() norm_path = paths.norm_dir cfg = cf.ConversionFactors(norm_path) q_e = 1.602e-19 m_e = 9.11e-31 m_p = 1.67e-27 k_b = 1.38e-23 epsilon0 = 8.854e-12 #--------------- # user parameters mult = 6 cfg.yi = 1 dim = '1D' bz_in = 50 # scale length (1,2,3) - (1 = 20um conduction zone sim,1.5, 2 = 40um conduction zone sim etc. ) s_len = 1
k_B = 1.38e-23
m_e = 9.11e-31
#-----
#---> file inputs
paths = hk.DataDirectoryPaths()
src_dir = paths.src_dir
data_dir = paths.data_dir_2D
if len(sys.argv) > 1:
save_path = sys.argv[1] + "/"
else:
save_path = "./"
norm_dir = paths.norm_dir
log_file = norm_dir + 'norm.log'
[T_ref, n_ref, Z_ref, Bz_ref] = np.loadtxt(log_file)
cd5 = cf.ConversionFactors(norm_dir, Z_ref, Ar=6.51)
Te_factor = 2.0 * T_ref * 1e-3
#---> getting path ->
lambda_p = 5
time_list = ['03', '05', '07']
path_list = []
bz_in = 0.0 # selected magnetic field [Tesla]
#--- generating path_list
s_list = [2]
file_name_without_path = (__file__.split('/')[-1]).split('.')[0]
save_tag = '%s_%iT_' % (file_name_without_path, bz_in)
for ip in range(len(s_list)):
path_loc = paths.get_path(s_list[ip], bz_in, lambda_p)
path_list.append(path_loc)
import modules.house_keeping as hk import modules.figure_latex as fprl from pylab import * SI_on = cd5_switches.SI_on save_on = cd5_switches.save_on hlines_on = cd5_switches.hlines_on grid_on = cd5_switches.grid_on horizontal_on = cd5_switches.horizontal_on separate_plots_on = cd5_switches.separate_plots_on marker_on = False #cd5_switches.marker_on MSIZE = 4 init_path = os.getcwd() norm_path = hk.DataDirectoryPaths().norm_dir cd5 = cfoil.ConversionFactors(norm_path) cl_index = int(cd5.cl_index) c_index = int(cd5.c_index) SI_on = cd5.SI_on tau_ei = cd5.tau_ei nu_ei = cd5.nu_ei lambda_mfp = cd5.lambda_mfp xstep_factor = cd5.xstep_factor tstep_factor = cd5.tstep_factor xlab = cd5.xlab ylab = cd5.ylab leg_title = cd5.leg_title color_lineout = cd5.color_lineout lineout_list = cd5.lineout_list