Python ReadClass.read_Q2_ave Examples

Example #1

File: draw_Q2_ave_diff.py Project: onuki0218/spectral-model-2D

rc('font', family='serif')

setting_name = "setting"
figure_dir = './figure/'

reader = ReadClass()
remote_data = reader.data_PID(setting_name)
reader.set_ssh(remote_data=remote_data)
reader.read_setting(setting_name)
reader.read_file_number(ssh_flag=True)

K_axis, L_axis = reader.set_axis()

c_level = 64

spec0 = reader.read_Q2_ave(0)

for time in range(1, 11):
    it = time * reader.interval_write_variable
    str_it = '{0:04d}'.format(it)
    spec = reader.read_Q2_ave(it) - spec0

    sum_spec = np.sum(spec)
    #
    print(it, sum_spec / np.sum(spec0))
    # spec_range = np.linspace(np.nanmin(spec), np.nanmax(spec), c_level)
    spec_range = np.linspace(0, np.nanmax(spec), c_level)

    fig = plt.figure(figsize=[4, 3])
    fig.subplots_adjust(left=0.15, bottom=0.15, right=0.9,
                        top=0.9, wspace=0.15, hspace=0.15)

Example #2

figure_dir = './figure/'

reader = ReadClass()
remote_data = reader.data_PID(setting_name)
reader.set_ssh(remote_data=remote_data)
reader.read_setting(setting_name)
reader.read_file_number(ssh_flag=True)

K_axis, L_axis = reader.set_axis()

c_level = 64

for time in range(0, 11):
    it = time * reader.interval_write_variable
    str_it = '{0:04d}'.format(it)
    spec = reader.read_Q2_ave(it)

    energy_all = np.sum(spec)
    spec[spec <= 0.0] = np.nan
    spec_log10 = np.log10(spec / energy_all)
    max_spec = np.nanmax(spec_log10)
    min_spec = np.nanmin(spec_log10)
    #
    print(it, max_spec, min_spec, energy_all)
    # print(spec_XY_Z_log10)
    spec_range = np.linspace(0, np.nanmax(spec), c_level)
    # spec_range_log = np.linspace(-5.0, 0.0, c_level)

    fig = plt.figure(figsize=[4, 3])
    fig.subplots_adjust(left=0.15,
                        bottom=0.15,