Python load_res示例

示例#1

文件： validation_case05.py 项目： ThomasSchuetz/workspace

                                                  extWallIndoorSurface,
                                                  intWallIndoorSurface, 
                                                  intGainsConv, intGainsRad, ports,
                                                  model,
                                                  dt=int(3600/times_per_hour))

# Compute averaged results
T_air_c = T_air - 273.15
T_air_mean = np.array([np.mean(T_air_c[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])

T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]

# Load reference results    
(T_air_ref_1, T_air_ref_10, T_air_ref_60) = tc.load_res("inputs/case05_res.csv")
T_air_ref_1 = T_air_ref_1[:,0]
T_air_ref_10 = T_air_ref_10[:,0]
T_air_ref_60 = T_air_ref_60[:,0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(res, label="Reference", color="black", linestyle="--")
    plt.plot(ref, label="Simulation", color="blue", linestyle="-")
    plt.legend()
    plt.ylabel("Temperature in degC")
    
    plt.title(title)

示例#2

                                                         times_per_hour))

# Compute averaged results
T_air_c = T_air - 273.15
T_air_mean = np.array([
    np.mean(T_air_c[i * times_per_hour:(i + 1) * times_per_hour])
    for i in range(24 * 60)
])

T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]

# Load reference results
(T_air_ref_1, T_air_ref_10,
 T_air_ref_60) = tc.load_res("inputs/case03_res.csv")
T_air_ref_1 = T_air_ref_1[:, 0]
T_air_ref_10 = T_air_ref_10[:, 0]
T_air_ref_60 = T_air_ref_60[:, 0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(res, label="Reference", color="black", linestyle="--")
    plt.plot(ref, label="Simulation", color="blue", linestyle="-")
    plt.legend()
    plt.ylabel("Temperature in degC")

    plt.title(title)

示例#3

文件： validation_case12.py 项目： ThomasSchuetz/workspace

                                                  extWallIndoorSurface,
                                                  intWallIndoorSurface, 
                                                  intGainsConv, intGainsRad, ports,
                                                  model,
                                                  dt=int(3600/times_per_hour))

# Compute averaged results
T_air_c = T_air - 273.15
T_air_mean = np.array([np.mean(T_air_c[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])

T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]

# Load reference results    
(T_air_ref_1, T_air_ref_10, T_air_ref_60) = tc.load_res("inputs/case12_res.csv")
T_air_ref_1 = T_air_ref_1[:,0]
T_air_ref_10 = T_air_ref_10[:,0]
T_air_ref_60 = T_air_ref_60[:,0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(ref, label="Reference", color="black", linestyle="--")
    plt.plot(res, label="Simulation", color="blue", linestyle="-")
    plt.legend()
    plt.ylabel("Temperature in degC")
    
    plt.title(title)

示例#4

文件： validation_case07.py 项目： ThomasSchuetz/workspace

                                                  windowIndoorSurface, 
                                                  extWallIndoorSurface,
                                                  intWallIndoorSurface, 
                                                  intGainsConv, intGainsRad, ports,
                                                  model,
                                                  dt=int(3600/times_per_hour))

# Compute averaged results
Q_hc_mean = np.array([np.mean(Q_HC[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])

Q_hc_1 = Q_hc_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440]

# Load reference results    
(Q_hc_ref_1, Q_hc_ref_10, Q_hc_ref_60) = tc.load_res("inputs/case07_res.csv")
Q_hc_ref_1 = Q_hc_ref_1[:,0]
Q_hc_ref_10 = Q_hc_ref_10[:,0]
Q_hc_ref_60 = Q_hc_ref_60[:,0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(ref, label="Reference", color="black", linestyle="--")
    plt.plot(res, label="Simulation", color="blue", linestyle="-")
    plt.legend()
    plt.ylabel("Heat load in W")
    
    plt.title(title)

示例#5

                                                  windowIndoorSurface, 
                                                  extWallIndoorSurface,
                                                  intWallIndoorSurface, 
                                                  intGainsConv, intGainsRad, ports,
                                                  model,
                                                  dt=int(3600/times_per_hour))

# Compute averaged results
Q_hc_mean = np.array([np.mean(Q_HC[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])

Q_hc_1 = Q_hc_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440]

# Load reference results    
(Q_hc_ref_1, Q_hc_ref_10, Q_hc_ref_60) = tc.load_res("inputs/case06_res.csv")
Q_hc_ref_1 = -Q_hc_ref_1[:,0]
Q_hc_ref_10 = -Q_hc_ref_10[:,0]
Q_hc_ref_60 = -Q_hc_ref_60[:,0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(ref, label="Reference", color="black", linestyle="--")
    plt.plot(res, label="Simulation", color="blue", linestyle="-")
    plt.legend()
    plt.ylabel("Heat load in W")
    
    plt.title(title)

示例#6

文件： validation_case11.py 项目： ThomasSchuetz/workspace

Q_iw_mean = np.array([np.mean(Q_iw[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
Q_ow_mean = np.array([np.mean(Q_ow[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])

Q_hc_1 = Q_hc_mean[0:24] + Q_iw_mean[0:24] + Q_ow_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240] + Q_iw_mean[216:240] + Q_ow_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440] + Q_iw_mean[1416:1440] + Q_ow_mean[1416:1440]

T_air_c = T_air - 273.15
T_air_mean = np.array([np.mean(T_air_c[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])

T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]

# Load reference results    
(load_res_1, load_res_10, load_res_60) = tc.load_res("inputs/case11_res.csv")
Q_hc_ref_1 = load_res_1[:,1]
Q_hc_ref_10 = load_res_10[:,1]
Q_hc_ref_60 = load_res_60[:,1]

T_air_ref_1 = load_res_1[:,0]
T_air_ref_10 = load_res_10[:,0]
T_air_ref_60 = load_res_60[:,0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(ref, label="Reference", color="black", linestyle="--")
    plt.plot(res, label="Simulation", color="blue", linestyle="-")

示例#7

Q_hc_1 = Q_hc_mean[0:24] + Q_iw_mean[0:24] + Q_ow_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240] + Q_iw_mean[216:240] + Q_ow_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440] + Q_iw_mean[1416:1440] + Q_ow_mean[1416:1440]

T_air_c = T_air - 273.15
T_air_mean = np.array([
    np.mean(T_air_c[i * times_per_hour:(i + 1) * times_per_hour])
    for i in range(24 * 60)
])

T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]

# Load reference results
(load_res_1, load_res_10, load_res_60) = tc.load_res("inputs/case11_res.csv")
Q_hc_ref_1 = load_res_1[:, 1]
Q_hc_ref_10 = load_res_10[:, 1]
Q_hc_ref_60 = load_res_60[:, 1]

T_air_ref_1 = load_res_1[:, 0]
T_air_ref_10 = load_res_10[:, 0]
T_air_ref_60 = load_res_60[:, 0]


# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
    plt.figure()
    ax_top = plt.subplot(211)
    plt.plot(ref, label="Reference", color="black", linestyle="--")
    plt.plot(res, label="Simulation", color="blue", linestyle="-")