Public Cell_type As String
Public Cell_efficiency As Variant
Public Prod_scale As Variant
Public system_type As String
Public grid As String
Public capacity As Variant
Public Irradiation As Variant
Public PR As Variant
Public r As Variant
Public om As Variant
Public Power As Variant
Public Power_MJ As Variant

Public D_module As Variant


Public ELC_sys_sog As Variant
Public OIL_sys_sog As Variant
Public COAL_sys_sog As Variant
Public PER_sys_sog As Variant
Public CO2_sys_sog As Variant
Public COST_sys_sog As Variant
Public ELC_sys_wafer As Variant
Public OIL_sys_wafer As Variant
Public COAL_sys_wafer As Variant
Public PER_sys_wafer As Variant
Public CO2_sys_wafer As Variant
Public COST_sys_wafer As Variant
Public ELC_sys_cell As Variant
Public OIL_sys_cell As Variant
Public COAL_sys_cell As Variant
Public PER_sys_cell As Variant
Public CO2_sys_cell As Variant
Public COST_sys_cell As Variant
Public ELC_sys_module As Variant
Public OIL_sys_module As Variant
Public COAL_sys_module As Variant
Public PER_sys_module As Variant
Public CO2_sys_module As Variant
Public COST_sys_module As Variant


Public AC_cell As Variant
Public eata_cell As Variant
Public eata_module As Variant
Public eata_elc As Variant
Public Tech_type As String
Public Tech As Integer
Public filename1 As String
Public filename2 As String
Dim dbs1(252, 7) As Variant
Dim dbs2(100, 7) As Variant
Public eata_r As Variant
Public Y_module As Variant
Public AC_module As Variant
Public n_cell As Integer
Public l_module As Variant
Public l_cell As Variant
Public w_module As Variant
Public w_cell As Variant
Public s_cell As Variant
Public AN_module As Variant
Public AS_module As Variant
Public Y_cell As Variant
Public AN_cell As Variant
Public AS_cell As Variant
Public AN_wafer As Variant
Public AS_wafer As Variant
Public AW_ingot As Variant
Public t_kerf As Variant
Public t_wafer As Variant
Public Y_wafer As Variant
Public Y_v As Variant
Public nn_ingot As Integer
Public w_ingot As Variant
Public l_ingot As Variant
Public s_ingot As Variant
Public AW_sog As Variant
Public Y_cast As Variant
Public AW_silicon As Variant
Public Y_dec As Variant
Public Y_ds As Variant
Public n_ds As Variant
Public Y_r As Variant
Public AW_quartz As Variant
Public Y_i As Variant
Public Y_m As Variant
Public Y_p As Variant
Public AN_cargo As Variant
Public CAP_cargo As Variant
Dim E_cargo(3) As Variant
Public UW_cargo As Variant
Public f_w As Variant
Public L_cargo As Integer
Public UC_cargo As Variant
Public f_m_cargo As Variant
Public f_b_cargo As Variant
Public r_c_cargo As Variant
Public C_cargo As Variant
Dim E_p_cargo(3) As Variant
Public UE_cargo As Variant
Public ADAYS_cargo As Integer
Public UL_cargo As Variant
Public c_l As Variant
Public AW_wg As Variant
Public CAP_ac As Variant
Public h_ac As Variant
Public ucap_ac As Variant
Public N_ac As Variant
Public LINE_ac As Variant
Public N0_ac As Integer
Public UW_ac As Variant
Dim E_ac(3) As Variant
Public L_equip As Variant
Public C_ac As Variant
Public UC_ac As Variant
Public f_c As Variant
Public r_c_equip As Variant
Dim E_p_ac(3) As Variant
Public UE_ac As Variant
Public C_p_ac As Variant
Dim E_na2co3(3) As Variant
Public C_l_ac As Variant
Public UL_ac As Variant
Public f_l As Variant
Public AREA_ac As Variant
Public uarea_ac As Variant
Public uarea_mf As Variant
Public f_a As Variant
Public CO2_ac As Variant
Public N_p As Variant
Public UCAP_p As Variant
Public N0_p As Variant
Public LINE_P As Variant
Public W_p As Variant
Public UW_p As Variant
Dim E_pu(3) As Variant
Public C_pu As Variant
Public UC_p As Variant
Public E_p_pu(3) As Variant
Public UE_p As Variant
Public C_p_pu As Variant
Public c_e As Variant
Public c_oil As Variant
Public AW_ma As Variant
Public mr_ma As Variant
Dim E_ma(3) As Variant
Public C_ma As Variant
Public DW_pu As Variant
Public UDW_p As Variant
Dim E_dw_pu(3) As Variant
Public C_dw_pu As Variant
Public AW_lime As Variant
Public mr_lime As Variant
Public C_l_pu As Variant
Public UL_P As Variant
Public AREA_pu As Variant
Public uarea_p As Variant
Dim E_lime(3) As Variant
Public C_lime As Variant
Dim E_build_sio2(3) As Variant
Dim ue_build(3) As Variant
Public w_build_steal As Variant
Public w_build_cement As Variant
Public L_build As Integer
Public C_build_sio2 As Variant
Public c_build As Variant
Public r_c_build As Variant
Public CO2_build_cement_sio2 As Variant
Public E_oh_sio2_elc As Variant
Public e_oh As Variant
Public C_oh_sio2 As Variant
Dim E_sio2(3) As Variant
Public C_sio2 As Variant
Public co2_elc As Variant
Public co2_OIL As Variant
Public co2_coal As Variant
Public uco2_sio2 As Variant
Public w_cp_ab As Variant
Public Y_cp_h As Variant
Public w_cp_re As Variant
Public AW_cp_c As Variant
Public AW_cp As Variant
Public UCAP_cp_n As Variant
Public CAP_cp_n As Variant
Public N_cp_n As Integer
Public MINCAP_cp_n As Variant
Public W_cp_n As Variant
Public UW_cp_n As Variant
Dim E_cp_n(3) As Variant
Public C_cp_n As Variant
Public UC_cp_n As Variant
Public UE_cp_n As Variant
Public UNITE_cp_n As Variant
Public E_p_cp_n_elc As Variant
Public C_p_cp_n As Variant
Public AREA_cp_n As Variant
Public uarea_cp_n As Variant
Public UCAP_cp_r As Variant
Public CAP_cp_r As Variant
Public N_cp_r As Integer
Public W_cp_r As Variant
Public UW_cp_r As Variant
Dim E_cp_r(3) As Variant
Public C_cp_r As Variant
Public UC_cp_r As Variant
Public UNITE_cp_r As Variant
Public UE_cp_r As Variant
Public E_p_cp_r_elc As Variant
Public C_p_cp_r As Variant
Public AREA_cp_r As Variant
Public uarea_cp_r As Variant
Public N_cp_f As Integer
Public UCAP_cp_f As Variant
Public W_cp_f As Variant
Public UW_cp_f As Variant
Dim E_cp_f(3) As Variant
Public C_cp_f As Variant
Public UC_cp_f As Variant
Public E_p_cp_f_elc As Variant
Public UE_cp_f As Variant
Public C_p_cp_f As Variant
Public AREA_cp_f As Variant
Public uarea_cp_f As Variant
Public AW_cp_ab As Variant
Dim E_cp_ab(3) As Variant
Public C_cp_ab As Variant
Public AW_cp_re As Variant
Dim E_cp_re(3) As Variant
Public C_cp_re As Variant
Public AW_cp_cl2 As Variant
Public w_cp_cl2 As Variant
Dim E_cp_cl2(3) As Variant
Public C_cp_cl2 As Variant
Public AW_cp_naoh As Variant
Public w_cp_naoh As Variant
Dim E_cp_naoh(3) As Variant
Public C_cp_naoh As Variant
Public DW_cp_f As Variant
Public UDW_cp_f As Variant
Dim E_dw_cp_f(3) As Variant
Public C_dw_cp_f As Variant
Public LBR_cp As Variant
Public UL_cp_f As Variant
Public UL_cp_n As Variant
Public UL_cp_r As Variant
Public C_l_cp As Variant
Public AREA_cp As Variant
Dim E_build_cp(3) As Variant
Public C_build_cp As Variant
Public CO2_build_cement_cp As Variant
Public E_oh_cp_elc As Variant
Public C_oh_cp As Variant
Dim E_cp(3) As Variant
Public C_cp As Variant
Public CO2_cp As Variant
Public uco2_cp As Variant
Public N_af As Variant
Public CAP_af As Variant
Public W_af_g As Variant
Public uw_af_g As Variant
Public W_af_brick As Variant
Public uw_af_brick As Variant
Public W_af As Variant
Dim E_af(3) As Variant
Public C_af As Variant
Public UC_af As Variant
Public E_p_af_elc As Variant
Public UE_af As Variant
Public C_p_af As Variant
Public AW_af_ar As Variant
Public v_af_ar As Variant
Dim E_af_ar(3) As Variant
Public C_af_ar As Variant
Public AW_af_n As Variant
Public w_af_n As Variant
Dim E_af_n(3) As Variant
Public C_af_n As Variant
Public AW_af_e As Variant
Public w_af_e As Variant
Dim E_af_e(3) As Variant
Public C_af_e As Variant
Public C_l_af As Variant
Public UL_af As Variant
Public AREA_af As Variant
Public uarea_af As Variant
Public CO2_af As Variant
Public E_co_af As Variant
Public N_dec As Integer
Public CAP_dec As Variant
Public LINE_dec As Integer
Public N0_dec As Integer
Public W_dec_sio2 As Variant
Public uw_dec_sio2 As Variant
Public W_dec As Variant
Dim E_dec(3) As Variant
Public C_dec As Variant
Public UC_dec As Variant
Public E_p_dec_elc As Variant
Public UE_dec As Variant
Public C_l_dec As Variant
Public UL_dec As Variant
Public AREA_dec As Variant
Public uarea_dec As Variant
Public AW_si_ds As Variant
Public NN_ds As Integer
Public CAP_ds As Variant
Public LINE_ds As Integer
Public N0_ds As Integer
Public W_ds As Variant
Public w_ds_brick As Variant
Dim E_ds(3) As Variant
Public C_ds As Variant
Public UC_ds As Variant
Public E_p_ds_elc As Variant
Public UE_ds As Variant
Public C_p_ds As Variant
Public C_l_ds As Variant
Public UL_ds As Variant
Public AREA_ds As Variant
Public uarea_ds As Variant
Public AREA_sog As Variant
Dim E_build_sog(3) As Variant
Public C_build_sog As Variant
Public CO2_build_cement_sog As Variant
Public E_oh_sog_elc As Variant
Public C_oh_sog As Variant
Public W_cp As Variant
Dim E2_cp(3) As Variant
Public C2_cp As Variant
Public CO22_cp As Variant
Dim E_sog(3) As Variant
Public C_sog As Variant
Public CO2_sog As Variant
Public uco2_sog As Variant
Public UCAP_em As Variant
Public UCAP_em_0 As Variant
Public r_em As Variant
Public r_em_0 As Variant
Public CF_em_max As Variant
Public CF_em_max_0 As Variant
Public N_em As Integer
Public UE_em As Variant
Public UE_em_0 As Variant
Public s_ingot_0 As Variant
Public w_ingot_0 As Variant
Public l_ingot_0 As Variant
Public UW_em As Variant
Public UW_em_0 As Variant
Public UC_em As Variant
Public UC_em_0 As Variant
Dim E_em(3) As Variant
Public f_w_wafer As Variant
Public C_em As Variant
Public CF_em As Variant
Public E_p_em_elc As Variant
Public f_e_wafer As Variant
Public C_p_em As Variant
Public UL_em As Variant
Public UN_wafer As Variant
Public h_ingot As Variant
Public n_ingot As Variant
Public UCAP_mws As Variant
Public CF_mws_max As Variant
Public T_mws_0 As Variant
Public N_mws As Variant
Public UE_mws As Variant
Public s_cell_0 As Variant
Public w_cell_0 As Variant
Public l_cell_0 As Variant
Public r_mws As Variant
Public r_mws_0 As Variant
Public n_ingot_0 As Variant
Public UE_mws_0 As Variant
Public UW_mws As Variant
Public UW_mws_0 As Variant
Public UC_mws As Variant
Public UC_mws_0 As Variant
Dim E_mws(3) As Variant
Public C_mws As Variant
Public T_mws As Variant
Public E_p_mws_elc As Variant
Public C_p_mws As Variant
Dim E_mws_sic(3) As Variant
Public w_sic As Variant
Dim E_mws_wire(3) As Variant
Public w_wire As Variant
Public C_mws_wire As Variant
Dim E_mws_om(4) As Variant
Public C_mws_om As Variant
Public C_l_mws As Variant
Public UL_mws As Variant
Public AREA_wafer As Variant
Public uarea_wafer As Variant
Dim E_build_wafer(3) As Variant
Public C_build_wafer As Variant
Public CO2_build_cement_wafer As Variant
Public E_oh_wafer_elc As Variant
Public C_oh_wafer As Variant
Dim E_wafer(3) As Variant
Public C_wafer As Variant
Public CO2_wafer As Variant
Dim ec_wafer(4) As Variant
Public uco2_wafer As Variant
Dim CAP_cell(8) As Variant
Dim UCAP_cell_0(8) As Variant
Public r_cell_s As Variant
Public T_cell_max As Variant
Public T_cell_max_0 As Variant
Dim NN_cell(8) As Integer
Dim UW_cell(8) As Variant
Dim UW_cell_0(8) As Variant
Dim UC_cell(8) As Variant
Dim UC_cell_0(8) As Variant
Dim E_cell_eq(3) As Variant
Public C_cell_eq As Variant
Dim UE_cell_elc(8) As Variant
Dim UE_cell_elc_0(8) As Variant
Dim T_cell(8) As Variant
Public E_p_cell_elc As Variant
Public C_p_cell_elc As Variant
Public AW_tef As Variant
Public uw_tef As Variant
Dim E_cell_tef(3) As Variant
Public C_cell_tef As Variant
Public uw_KOH As Variant
Public uw_KOH_0 As Variant
Dim E_cell_KOH(3) As Variant
Public C_cell_KOH As Variant
Public uw_IPA As Variant
Public uw_IPA_0 As Variant
Dim E_cell_IPA(3) As Variant
Public C_cell_IPA As Variant
Dim uw_dw(2) As Variant
Dim uw_dw_0(2) As Variant
Dim E_cell_dw(3) As Variant
Public C_cell_dw As Variant
Public uw_POCl3 As Variant
Public uw_POCl3_0 As Variant
Dim E_cell_POCl3(3) As Variant
Public C_cell_POCl3 As Variant
Dim uw_N2O2(3) As Variant
Dim uw_N2O2_0(3) As Variant
Dim E_cell_N2O2(3) As Variant
Public C_cell_N2O2 As Variant
Public uw_tube As Variant
Public uw_tube_0 As Variant
Public un_tube As Variant
Dim E_cell_tube(3) As Variant
Public uw_folder As Variant
Public uw_folder_0 As Variant
Public un_folder As Variant
Dim E_cell_folder(3) As Variant
Public C_cell_folder As Variant
Public uw_TPT As Variant
Public uw_TPT_0 As Variant
Dim E_cell_TPT(3) As Variant
Public C_cell_TPT As Variant
Public uw_AgAl As Variant
Public uw_AgAl_0 As Variant
Dim E_cell_AgAl(3) As Variant
Public C_cell_AgAl As Variant
Public uw_Ag As Variant
Public uw_Ag_0 As Variant
Dim E_cell_Ag(3) As Variant
Public C_cell_Ag As Variant
Public uw_screen As Variant
Public uw_screen_0 As Variant
Public un_screen As Variant
Dim E_cell_screen(3) As Variant
Public C_cell_screen As Variant
Public un_Xe As Variant
Public un_Xe_0 As Variant
Dim E_cell_Xe(3) As Variant
Public C_cell_Xe As Variant
Public un_case As Variant
Public un_case_0 As Variant
Dim E_cell_case(3) As Variant
Public C_cell_case As Variant
Public C_l_cell As Variant
Dim UL_cell(5) As Variant
Public AREA_cell As Variant
Dim uarea_cell(8) As Variant
Dim E_build_cell(3) As Variant
Public C_build_cell As Variant
Public CO2_build_cement_cell As Variant
Public E_oh_cell_elc As Variant
Public C_oh_cell As Variant
Dim E_cell(3) As Variant
Public C_cell As Variant
Public CO2_cell As Variant
Public LINE_module As Integer
Public UCAP_module As Variant
Public CF_module_max As Variant
Public CF_module As Variant
Dim E_module_eq(3) As Variant
Public UW_module As Variant
Public C_module_eq As Variant
Public UC_module As Variant
Public E_p_module_elc As Variant
Public C_p_module As Variant
Public UE_module As Variant
Dim E_glass(3) As Variant
Public t_glass As Variant
Public C_glass As Variant
Dim E_Al(3) As Variant
Public C_Al As Variant
Public w_Al_l As Variant
Public w_Al_w As Variant
Dim E_wire(3) As Variant
Public C_wire As Variant
Public w_Cu As Variant
Dim E_EVA(3) As Variant
Public C_EVA As Variant
Public w_EVA As Variant
Dim E_BS(3) As Variant
Public C_BS As Variant
Public w_Tedlar As Variant
Public w_AS As Variant
Dim E_rubber(3) As Variant
Public C_rubber As Variant
Public w_rubber As Variant
Public C_l_module As Variant
Public UL_module As Variant
Public AREA_module As Variant
Public uarea_module As Variant
Dim E_build_module(3) As Variant
Public C_build_module As Variant
Public CO2_build_cement_module As Variant
Public E_oh_module_elc As Variant
Public C_oh_module As Variant
Dim E_module(3) As Variant
Public C_module As Variant
Public CO2_module As Variant
Dim ec_cell(4) As Variant
Public uco2_cell As Variant
Dim ec_module(4) As Variant
Public uco2_module As Variant
Public ingots As Variant
Dim Total_ec_module(4) As Variant
Public Total_co2_module As Variant
Dim ec_perW_module(4) As Variant
Public uco2_perW_module As Variant
Public alpha As Variant

Public sf1  As Variant
Public sf2  As Variant
Public sf3  As Variant
Public sf4  As Variant

Dim ec_steal1(4) As Variant
Dim ec_steal2(4) As Variant
Dim ec_ss(4) As Variant
Dim ec_na2co3(4) As Variant
Dim ec_ma(4) As Variant
Dim ec_dw(4) As Variant
Dim ec_lime(4) As Variant
Dim ec_cement(4) As Variant
Dim ec_sio2(4) As Variant
Dim ec_ab(4) As Variant
Dim ec_re(4) As Variant
Dim ec_cl2(4) As Variant
Dim ec_naoh(4) As Variant
Public ec_w As Variant
Dim ec_cp(4) As Variant
Dim ec_brick(4) As Variant
Dim ec_ar(4) As Variant
Dim ec_sog(4) As Variant
Dim ec_sic(4) As Variant
Dim ec_wire(4) As Variant
Dim ec_tef(4) As Variant
Dim ec_KOH(4) As Variant
Dim ec_IPA(4) As Variant
Dim ec_POCl3(4) As Variant
Dim ec_N2O2(4) As Variant
Dim ec_tube(4) As Variant
Dim ec_folder(4) As Variant
Dim ec_TPT(4) As Variant
Dim ec_AgAl(4) As Variant
Dim ec_screen(4) As Variant
Dim ec_Xe(4) As Variant
Dim ec_case(4) As Variant
Dim ec_glass(4) As Variant
Dim ec_Al(4) As Variant
Dim ec_Cu(4) As Variant
Dim ec_EVA(4) As Variant
Dim ec_Tedlar(4) As Variant
Dim ec_BS(4) As Variant
Dim ec_rubber(4) As Variant
Dim ec_in(4) As Variant


'物理定数などの定義
Public Const sw_Si As Variant = 2.34 '[t/m3]
Public Const mw_SiO2 As Variant = 60.08 '[g/mol]
Public Const mw_na2co3 As Variant = 105.99 '[g/mol]
Public Const aw_c As Variant = 12.01 '[g/mol]
Public Const aw_Si As Variant = 28.09 '[g/mol]
Public Const e_oil As Variant = 10500 '[Mcal/t]
Public Const mw_wg As Variant = 602.54 '[g/mol]
Public Const c_cement_p As Variant = 0.214 '[t-C/t-cement]
Public Const r_cp_c As Variant = 0.6 '[-]
Public Const w_c_re As Variant = 0.7 '[-]
Public Const sw_ar As Variant = 0.00178 '[t/m3]
Public Const h_co As Variant = 2420 '[Mcal/t-CO]
Public Const mw_co As Variant = 28.0104 '[g/mol]
Public Const sw_glass As Variant = 2.5 '[t/m3]
Public Const JCal As Variant = 4.1861 '[J/cal]
Public Const kWhMcal As Variant = 0.86 '[kWh/Mcal]
Public Const kWhMJ As Variant = 3.6 '[MJ/kWh]
Public Const CO2_oil_power As Variant = 200 '[g-C/kWh}


Public filename3 As String
Dim dbs3(5, 9) As Variant
Public AS_steal As Variant
Public AS_Al As Variant
Public AS_Cu As Variant
Public AS_in As Variant
Public AS_elc As Variant
Public AS_OIL As Variant
Public AS_cost As Variant
Public PC_steal As Variant
Public PC_Al As Variant
Public PC_Cu As Variant
Public PC_in As Variant
Public PC_elc As Variant
Public PC_OIL As Variant
Public PC_cost As Variant
Public WM_steal As Variant
Public WM_Al As Variant
Public WM_Cu As Variant
Public WM_in As Variant
Public WM_elc As Variant
Public WM_OIL As Variant
Public WM_cost As Variant
Public LC_steal As Variant
Public LC_Al As Variant
Public LC_Cu As Variant
Public LC_in As Variant
Public LC_elc As Variant
Public LC_OIL As Variant
Public LC_cost As Variant

Public AREA_as As Variant
Public Life_module As Variant
Public Life_pc As Variant
Public N_pc As Integer


Public Steal_pc As Variant
Public Al_pc As Variant
Public Cu_pc As Variant
Public In_pc As Variant
Public Steal_as As Variant

Dim EC_pc_steal(3) As Variant
Dim EC_pc_Al(3) As Variant
Dim EC_pc_Cu(3) As Variant
Dim EC_pc_in(3) As Variant
Dim EC_as_steal(3) As Variant
Dim EC_wm_cu(3) As Variant

Dim Req_sys_module(5) As Variant
Dim Req_sys_pc(5) As Variant
Dim Req_sys_as(5) As Variant
Dim Req_sys_wm(5) As Variant
Public Req_sys_LC As Variant

Public Total_elc As Variant
Public Total_OIL As Variant
Public Total_coal As Variant
Public Total_cost As Variant
Public Total_co2 As Variant


Sub Input_poly_assumption()

'技術前提条件ファイルの読み込み
filename1 = "poly-assumption.csv"
Open filename1 For Input As #1
n = 0
Do While Not EOF(1)
    n = n + 1
    For m = 1 To 7
        Input #1, dbs1(n, m)
        'Debug.Print dbs(n, m)
    Next m
    'Debug.Print n, dbs1(n, 7)
Loop

Close #1

Tech_type = PV_LCA8.Label7.Caption
'製造技術オプションケースの格納
    If Tech_type = "現在技術" Then
        Tech = 1
    ElseIf Tech_type = "近未来技術" Then
        Tech = 2
    Else
        Tech = 3
    End If

'Debug.Print Tech

'*****データ配列から変数名への格納*****
eata_r = dbs1(32, Tech + 3)
Y_module = dbs1(29, Tech + 3)
w_cell = dbs1(3, Tech + 3)
w_cell_0 = dbs1(3, 4)
l_cell = dbs1(4, Tech + 3)
l_cell_0 = dbs1(4, 4)
w_module = dbs1(30, Tech + 3)
l_module = dbs1(31, Tech + 3)
Y_cell = dbs1(49, Tech + 3)
t_kerf = dbs1(127, Tech + 3)
t_wafer = dbs1(129, Tech + 3)
Y_wafer = dbs1(122, Tech + 3)
nn_ingot = dbs1(128, Tech + 3)
w_ingot = dbs1(147, Tech + 3)
l_ingot = dbs1(148, Tech + 3)
Y_cast = dbs1(144, Tech + 3)
Y_dec = dbs1(235, Tech + 3)
Y_ds = dbs1(244, Tech + 3)
n_ds = dbs1(245, Tech + 3)
Y_r = dbs1(223, Tech + 3)
Y_i = dbs1(160, Tech + 3)
Y_m = dbs1(172, Tech + 3)
Y_p = dbs1(210, Tech + 3)
CAP_cargo = dbs1(161, Tech + 3)
UW_cargo = dbs1(162, Tech + 3)
f_w = dbs1(12, Tech + 3)
L_cargo = dbs1(163, Tech + 3)
UC_cargo = dbs1(164, Tech + 3)
f_m_cargo = dbs1(166, Tech + 3)
f_b_cargo = dbs1(167, Tech + 3)
r_c_cargo = dbs1(165, Tech + 3)
UE_cargo = dbs1(168, Tech + 3)
ADAYS_cargo = dbs1(169, Tech + 3)
UL_cargo = dbs1(170, Tech + 3)
c_l = dbs1(9, Tech + 3)
h_ac = dbs1(173, Tech + 3)
ucap_ac = dbs1(174, Tech + 3)
N0_ac = dbs1(175, Tech + 3)
UW_ac = dbs1(176, Tech + 3)
sf2 = dbs1(25, Tech + 3)
L_equip = dbs1(15, Tech + 3)
UC_ac = dbs1(177, Tech + 3)
f_c = dbs1(13, Tech + 3)
r_c_equip = dbs1(16, Tech + 3)
UE_ac = dbs1(178, Tech + 3)
UL_ac = dbs1(179, Tech + 3)
f_l = dbs1(14, Tech + 3)
uarea_ac = dbs1(180, Tech + 3)
uarea_mf = dbs1(181, Tech + 3)
f_a = dbs1(17, Tech + 3)
UCAP_p = dbs1(212, Tech + 3)
N0_p = dbs1(214, Tech + 3)
UW_p = dbs1(216, Tech + 3)
sf3 = dbs1(26, Tech + 3)
sf1 = dbs1(24, Tech + 3)
sf4 = dbs1(27, Tech + 3)
UC_p = dbs1(217, Tech + 3)
UE_p = dbs1(220, Tech + 3)
c_e = dbs1(5, Tech + 3)
c_oil = dbs1(11, Tech + 3)
mr_ma = dbs1(218, Tech + 3)
UDW_p = dbs1(221, Tech + 3)
mr_lime = dbs1(219, Tech + 3)
UL_P = dbs1(215, Tech + 3)
uarea_p = dbs1(213, Tech + 3)
w_build_cement = dbs1(18, Tech + 3)
w_build_steal = dbs1(19, Tech + 3)
L_build = dbs1(22, Tech + 3)
c_build = dbs1(20, Tech + 3)
r_c_build = dbs1(23, Tech + 3)
e_oh = dbs1(21, Tech + 3)
Y_cp_h = dbs1(183, Tech + 3)
UCAP_cp_n = dbs1(188, Tech + 3)
MINCAP_cp_n = dbs1(187, Tech + 3)
UW_cp_n = dbs1(190, Tech + 3)
UC_cp_n = dbs1(189, Tech + 3)
UE_cp_n = dbs1(191, Tech + 3)
uarea_cp_n = dbs1(193, Tech + 3)
UCAP_cp_r = dbs1(195, Tech + 3)
UW_cp_r = dbs1(197, Tech + 3)
UC_cp_r = dbs1(196, Tech + 3)
UE_cp_r = dbs1(198, Tech + 3)
uarea_cp_r = dbs1(200, Tech + 3)
UCAP_cp_f = dbs1(202, Tech + 3)
UW_cp_f = dbs1(204, Tech + 3)
UC_cp_f = dbs1(203, Tech + 3)
UE_cp_f = dbs1(205, Tech + 3)
uarea_cp_f = dbs1(207, Tech + 3)
w_cp_cl2 = dbs1(184, Tech + 3)
w_cp_naoh = dbs1(185, Tech + 3)
UDW_cp_f = dbs1(208, Tech + 3)
ec_w = dbs1(8, Tech + 3)
UL_cp_f = dbs1(206, Tech + 3)
UL_cp_n = dbs1(192, Tech + 3)
UL_cp_r = dbs1(199, Tech + 3)
co2_OIL = dbs1(6, Tech + 3)
co2_coal = dbs1(7, Tech + 3)
CAP_af = dbs1(224, Tech + 3)
uw_af_g = dbs1(228, Tech + 3)
uw_af_brick = dbs1(229, Tech + 3)
UC_af = dbs1(225, Tech + 3)
UE_af = dbs1(233, Tech + 3)
v_af_ar = dbs1(230, Tech + 3)
w_af_n = dbs1(231, Tech + 3)
w_af_e = dbs1(232, Tech + 3)
UL_af = dbs1(226, Tech + 3)
uarea_af = dbs1(227, Tech + 3)
CAP_dec = dbs1(236, Tech + 3)
N0_dec = dbs1(237, Tech + 3)
uw_dec_sio2 = dbs1(241, Tech + 3)
UC_dec = dbs1(238, Tech + 3)
UE_dec = dbs1(242, Tech + 3)
UL_dec = dbs1(239, Tech + 3)
uarea_dec = dbs1(240, Tech + 3)
CAP_ds = dbs1(246, Tech + 3)
N0_ds = dbs1(247, Tech + 3)
w_ds_brick = dbs1(251, Tech + 3)
UC_ds = dbs1(248, Tech + 3)
UE_ds = dbs1(252, Tech + 3)
UL_ds = dbs1(249, Tech + 3)
uarea_ds = dbs1(250, Tech + 3)
r_em = dbs1(149, Tech + 3)
r_em_0 = dbs1(149, 4)
CF_em_max = dbs1(150, Tech + 3)
CF_em_max_0 = dbs1(150, 4)
UE_em_0 = dbs1(154, 4)
w_ingot_0 = dbs1(147, 4)
l_ingot_0 = dbs1(148, 4)
UW_em_0 = dbs1(155, 4)
UC_em_0 = dbs1(156, 4)
f_w_wafer = dbs1(123, Tech + 3)
f_e_wafer = dbs1(124, Tech + 3)
UL_em = dbs1(157, Tech + 3)
h_ingot = dbs1(132, Tech + 3)
n_ingot = dbs1(130, Tech + 3)
n_ingot_0 = dbs1(130, 4)
CF_mws_max = dbs1(125, Tech + 3)
T_mws_0 = dbs1(126, 4)
r_mws = dbs1(131, Tech + 3)
r_mws_0 = dbs1(131, 4)
UE_mws_0 = dbs1(137, 4)
UW_mws_0 = dbs1(139, 4)
UC_mws_0 = dbs1(140, 4)
w_sic = dbs1(145, Tech + 3)
w_wire = dbs1(146, Tech + 3)
UL_mws = dbs1(141, Tech + 3)
uarea_wafer = dbs1(142, Tech + 3)
UCAP_cell_0(1) = dbs1(55, Tech + 3)
UW_cell_0(1) = dbs1(56, Tech + 3)
UC_cell_0(1) = dbs1(57, Tech + 3)
UE_cell_elc_0(1) = dbs1(59, Tech + 3)
UCAP_cell_0(2) = dbs1(65, Tech + 3)
UW_cell_0(2) = dbs1(66, Tech + 3)
UC_cell_0(2) = dbs1(67, Tech + 3)
UE_cell_elc_0(2) = dbs1(69, Tech + 3)
UCAP_cell_0(3) = dbs1(79, Tech + 3)
UW_cell_0(3) = dbs1(80, Tech + 3)
UC_cell_0(3) = dbs1(81, Tech + 3)
UE_cell_elc_0(3) = dbs1(83, Tech + 3)
UCAP_cell_0(5) = dbs1(94, Tech + 3)
UW_cell_0(5) = dbs1(95, Tech + 3)
UC_cell_0(5) = dbs1(96, Tech + 3)
UE_cell_elc_0(5) = dbs1(98, Tech + 3)
UCAP_cell_0(6) = dbs1(101, Tech + 3)
UW_cell_0(6) = dbs1(102, Tech + 3)
UC_cell_0(6) = dbs1(103, Tech + 3)
UE_cell_elc_0(6) = dbs1(105, Tech + 3)
UCAP_cell_0(7) = dbs1(107, Tech + 3)
UW_cell_0(7) = dbs1(108, Tech + 3)
UC_cell_0(7) = dbs1(109, Tech + 3)
UE_cell_elc_0(7) = dbs1(111, Tech + 3)
UCAP_cell_0(8) = dbs1(113, Tech + 3)
UW_cell_0(8) = dbs1(114, Tech + 3)
UC_cell_0(8) = dbs1(115, Tech + 3)
UE_cell_elc_0(8) = dbs1(117, Tech + 3)
r_cell_s = dbs1(52, Tech + 3)
T_cell_max_0 = dbs1(50, 4)
T_cell_max = dbs1(50, Tech + 3)
uw_tef = dbs1(51, Tech + 3)
uw_KOH_0 = dbs1(61, Tech + 3)
uw_IPA_0 = dbs1(62, Tech + 3)
uw_dw_0(1) = dbs1(63, Tech + 3)
uw_dw_0(2) = dbs1(77, Tech + 3)
uw_POCl3_0 = dbs1(71, Tech + 3)
uw_N2O2_0(1) = dbs1(72, Tech + 3)
uw_N2O2_0(2) = dbs1(89, Tech + 3)
uw_N2O2_0(3) = dbs1(90, Tech + 3)
uw_tube_0 = dbs1(73, Tech + 3)
un_tube = dbs1(74, Tech + 3)
uw_folder_0 = dbs1(75, Tech + 3)
un_folder = dbs1(76, Tech + 3)
uw_TPT_0 = dbs1(85, Tech + 3)
uw_AgAl_0 = dbs1(87, Tech + 3)
uw_Ag_0 = dbs1(88, Tech + 3)
uw_screen_0 = dbs1(91, Tech + 3)
un_screen = dbs1(92, Tech + 3)
un_Xe_0 = dbs1(119, Tech + 3)
un_case_0 = dbs1(120, Tech + 3)
UL_cell(1) = dbs1(60, Tech + 3)
UL_cell(2) = dbs1(70, Tech + 3)
UL_cell(3) = dbs1(84, Tech + 3)
UL_cell(4) = dbs1(99, Tech + 3)
UL_cell(5) = dbs1(118, Tech + 3)
uarea_cell(1) = dbs1(58, Tech + 3)
uarea_cell(2) = dbs1(68, Tech + 3)
uarea_cell(3) = dbs1(82, Tech + 3)
uarea_cell(5) = dbs1(97, Tech + 3)
uarea_cell(6) = dbs1(104, Tech + 3)
uarea_cell(7) = dbs1(110, Tech + 3)
uarea_cell(8) = dbs1(116, Tech + 3)
UCAP_module = dbs1(41, Tech + 3)
CF_module_max = dbs1(42, Tech + 3)
UW_module = dbs1(43, Tech + 3)
UC_module = dbs1(44, Tech + 3)
UE_module = dbs1(45, Tech + 3)
t_glass = dbs1(33, Tech + 3)
l_module = dbs1(31, Tech + 3)
w_module = dbs1(30, Tech + 3)
w_Al_l = dbs1(35, Tech + 3)
w_Al_w = dbs1(34, Tech + 3)
w_Cu = dbs1(36, Tech + 3)
w_EVA = dbs1(37, Tech + 3)
w_Tedlar = dbs1(38, Tech + 3)
w_AS = dbs1(39, Tech + 3)
w_rubber = dbs1(40, Tech + 3)
UL_module = dbs1(46, Tech + 3)
uarea_module = dbs1(47, Tech + 3)

'材料原単位ファイルの読み込み
filename2 = "materials.csv"
Open filename2 For Input As #2
n = 0
Do While Not EOF(2)
    n = n + 1
    For m = 1 To 7
        Input #2, dbs2(n, m)
        'Debug.Print dbs(n, m)
    Next m
    'Debug.Print n, dbs1(n, 7)
Loop
Close #2

'*****データ配列から変数名への格納*****
For n = 1 To 4
    ec_steal1(n) = dbs2(2, n + 3)
    ec_steal2(n) = dbs2(3, n + 3)
    ec_ss(n) = dbs2(4, n + 3)
    ec_na2co3(n) = dbs2(5, n + 3)
    ec_ma(n) = dbs2(6, n + 3)
    ec_dw(n) = dbs2(7, n + 3)
    ec_lime(n) = dbs2(8, n + 3)
    ec_cement(n) = dbs2(9, n + 3)
    ec_ab(n) = dbs2(10, n + 3)
    ec_re(n) = dbs2(11, n + 3)
    ec_cl2(n) = dbs2(12, n + 3)
    ec_naoh(n) = dbs2(13, n + 3)
    ec_brick(n) = dbs2(14, n + 3)
    ec_ar(n) = dbs2(15, n + 3)
    ec_sic(n) = dbs2(16, n + 3)
    ec_wire(n) = dbs2(17, n + 3)
    ec_tef(n) = dbs2(18, n + 3)
    ec_KOH(n) = dbs2(19, n + 3)
    ec_IPA(n) = dbs2(20, n + 3)
    ec_POCl3(n) = dbs2(21, n + 3)
    ec_N2O2(n) = dbs2(22, n + 3)
    ec_tube(n) = dbs2(23, n + 3)
    ec_folder(n) = dbs2(24, n + 3)
    ec_TPT(n) = dbs2(25, n + 3)
    ec_AgAl(n) = dbs2(26, n + 3)
    ec_screen(n) = dbs2(27, n + 3)
    ec_Xe(n) = dbs2(28, n + 3)
    ec_case(n) = dbs2(29, n + 3)
    ec_glass(n) = dbs2(30, n + 3)
    ec_Al(n) = dbs2(31, n + 3)
    ec_Cu(n) = dbs2(32, n + 3)
    ec_EVA(n) = dbs2(33, n + 3)
    ec_Tedlar(n) = dbs2(34, n + 3)
    ec_BS(n) = dbs2(35, n + 3)
    ec_rubber(n) = dbs2(36, n + 3)
    ec_in(n) = dbs2(37, n + 3)
Next n

End Sub


Sub poly_LCA() '==========多結晶SiモジュールのLCA計算==========

AC_cell = Val(PV_LCA8.Label5.Caption) * 10 ^ 6  '[W/y]

eata_cell = Val(PV_LCA8.Label18.Caption) / 100  '[-]
co2_elc = Val(PV_LCA8.Label19.Caption) '[g-C/kWh]
eata_elc = Val(PV_LCA8.Label16.Caption) / 100 '[-]

alpha = PV_LCA8.Label20.Caption

n_ds = Val(PV_LCA8.Label35.Caption)

'*****0:共通パラメータ算出************************************************
'単位面積あたりの建屋建設エネルギー
For n = 1 To 3
    ue_build(n) = ec_steal1(n) * w_build_steal + ec_cement(n) * w_build_cement
Next n

'*****B-1:Si収支計算*****************************************************

'(1)太陽電池モジュール年間生産量
AC_module = eata_r * Y_module * AC_cell '[W/y]
n_cell = Int(l_module / l_cell) * Int(w_module / w_cell)
s_cell = w_cell * l_cell    '[m2/cell]
AN_module = Y_module / (eata_cell * n_cell * s_cell * 10 ^ 3) * AC_cell '[枚/y]
AS_module = AN_module * (l_module * w_module)   '[m2/y]

'(2)太陽電池セル年間生産量
AN_cell = AC_cell / (eata_cell * 10 ^ 3 * s_cell)   '[枚/y]
AS_cell = AN_cell * s_cell '[m2/y]

'(3)ウェハ年間生産量
AN_wafer = AN_cell / Y_cell '[枚/y]
AS_wafer = AS_cell / Y_cell '[m2/y]

'Debug.Print "AS_wafer="; AS_wafer

'(4)Siインゴット(切り出し後)年間生産量
AW_ingot = AS_wafer * (t_kerf + t_wafer) / 10 ^ 6 * sw_Si / Y_wafer '[t/y]

'Debug.Print "AW_ingot="; AW_ingot

'(5)SOG-Si年間生産量
s_ingot = w_ingot * l_ingot
ingots = Int((w_ingot / w_cell) * (l_ingot / l_cell))
Y_v = ingots * (s_cell / s_ingot)
AW_sog = AW_ingot / (Y_v * Y_cast)

'Debug.Print w_ingot; l_ingot; ingots; s_cell
'Debug.Print "AW_sog="; AW_sog

'(6)粗精製Si年間生産量
AW_silicon = AW_sog / (Y_dec * Y_ds ^ n_ds)

'(7)高純度SiO2年間生産量
AW_sio2 = (AW_silicon / Y_r) * mw_SiO2 / aw_Si

'(8)原料珪石の年間必要量
AW_quartz = AW_sio2 / (Y_i * Y_m * Y_p)

'Debug.Print "AW_quartz="; AW_quartz

'*****B-2:高純度SiO2製造原単位計算*****************************************

'B-2-1 原料珪石輸送工程---------------------------------------------------

'(1)装置製造

AN_cargo = AW_quartz / CAP_cargo
If AN_cargo < 1 Then
    AN_cargo = AN_cargo
Else
    If AN_cargo = Int(AN_cargo) Then
        AN_cargo = Int(AN_cargo) '輸送船年間必要隻数
    Else
        AN_cargo = Int(AN_cargo) + 1 '輸送船年間必要隻数
    End If
End If

For n = 1 To 3
    E_cargo(n) = (AN_cargo * UW_cargo * f_w * ec_steal1(n)) / L_cargo '年間輸送船製造エネルギー
Next n

If CAP_cargo >= AW_quartz Then
    f_b_cargo = 1.5
End If
C_cargo = AN_cargo * UC_cargo * f_m_cargo * f_b_cargo * r_c_cargo  '輸送船年経費

'(2)直接投入エネルギー
E_p_cargo(2) = AN_cargo * UE_cargo * ADAYS_cargo * e_oil
C_p_cargo = AN_cargo * UE_cargo * ADAYS_cargo * c_oil * f_m_cargo * f_b_cargo

'(3)人件費
C_l_cargo = AN_cargo * UL_cargo * c_l * f_m_cargo * f_b_cargo

'Debug.Print E_cargo(1) + E_p_cargo(1); E_p_cargo(2) * JCal; E_cargo(3) * JCal; C_cargo + C_p_cargo + C_l_cargo

'B-2-2 水ガラス化工程-----------------------------------------------------

AW_wg = AW_quartz * mw_wg / (3 * mw_SiO2) '水ガラス年間生産量
CAP_ac = (8760 / h_ac) * ucap_ac    'オートクレーブ年間最大処理量

If AW_wg / CAP_ac = Int(AW_wg / CAP_ac) Then
    N_ac = Int(AW_wg / CAP_ac) 'オートクレーブ必要台数
Else
    N_ac = Int(AW_wg / CAP_ac) + 1 'オートクレーブ必要台数
End If

If N_ac > N0_ac Then
    If AW_wg / (CAP_ac * N0_ac) = Int(AW_wg / (CAP_ac * N0_ac)) Then
        LINE_ac = Int(AW_wg / (CAP_ac * N0_ac)) '必要ライン数
    Else
        LINE_ac = Int(AW_wg / (CAP_ac * N0_ac)) + 1 '必要ライン数
    End If
Else
    LINE_ac = 1
End If

'(1)装置製造

If N_ac > N0_ac Then    '装置総重量
    W_ac = UW_ac * LINE_ac * f_w
Else
    W_ac = UW_ac * LINE_ac * (N_ac / N0_ac) ^ sf2 * f_w
End If

For n = 1 To 3  '年間装置製造エネルギー
    E_ac(n) = (W_ac * ec_ss(n)) / L_equip
    'Debug.Print E_ac(n)
Next n

If N_ac > N0_ac Then    '装置年経費
    C_ac = UC_ac * LINE_ac * f_c * r_c_equip
Else
    C_ac = UC_ac * LINE_ac * (N_ac / N0_ac) ^ sf2 * f_c * r_c_equip
End If
'Debug.Print C_ac

'(2)直接投入エネルギー

E_p_ac(2) = AW_wg * UE_ac
C_p_ac = (E_p_ac(2) / e_oil) * c_oil
'Debug.Print AW_wg, E_p_ac(2), C_p_ac

'(3)材料投入

    'Na2CO3
    AW_na2co3 = AW_quartz * mw_na2co3 / (3 * mw_SiO2)
    For n = 1 To 3  '年間投入エネルギー
        E_na2co3(n) = AW_na2co3 * ec_na2co3(n)
        'Debug.Print E_na2co3(n)
    Next n
    C_na2co3 = AW_na2co3 * ec_na2co3(4)
    'Debug.Print C_na2co3

'(4)人件費
C_l_ac = LINE_ac * UL_ac * c_l * f_l
'Debug.Print C_l_ac

'(5)装置占有面積
AREA_ac = (uarea_ac * N_ac + uarea_mf * LINE_ac) * f_a

'(6)反応からのCO2直接排出 [g-C/y]
CO2_ac = AW_na2co3 * aw_c / mw_na2co3 * 10 ^ 6

'Debug.Print E_ac(1) + E_na2co3(1); (E_ac(2) + E_p_ac(2) + E_na2co3(2)) * JCal; (E_ac(3) + E_na2co3(3)) * JCal; C_ac + C_p_ac + C_na2co3 + C_l_ac

'B-2-3 析出・精製工程---------------------------------------------------

'(1)装置製造

If AW_sio2 / UCAP_p = Int(AW_sio2 / UCAP_p) Then
    N_p = Int(AW_sio2 / UCAP_p) '洗浄槽必要台数
Else
    N_p = Int(AW_sio2 / UCAP_p) + 1 '洗浄槽必要台数
End If

If N_p >= N0_p Then
    If N_p / N0_p = Int(N_p / N0_p) Then
        LINE_P = Int(N_p / N0_p)
    Else
        LINE_P = Int(N_p / N0_p) + 1
    End If
    W_p = UW_p * LINE_P * f_w
    For n = 1 To 3
        E_pu(n) = (W_p * ec_ss(n)) / L_equip
    Next n
    C_pu = UC_p * LINE_P * f_c * r_c_equip
Else
    LINE_P = 1
    W_p = UW_p * LINE_P * f_w * (N_p / N0_p) ^ sf3
    For n = 1 To 3
        E_pu(n) = (W_p * ec_ss(n)) / L_equip
    Next n
    C_pu = UC_p * LINE_P * f_c * r_c_equip * (N_p / N0_p) ^ sf3
End If

'(2)直接投入エネルギー
E_p_pu(1) = AW_sio2 * UE_p
C_p_pu = E_p_pu(1) * c_e


'(3)材料投入

    '鉱酸
    AW_ma = AW_sio2 * mr_ma
    For n = 1 To 3
        E_ma(n) = AW_ma * ec_ma(n)
    Next n
    C_ma = AW_ma * ec_ma(4)
    
    '純水
    DW_pu = AW_sio2 * UDW_p
    For n = 1 To 3
        E_dw_pu(n) = DW_pu * ec_dw(n)
    Next n
    C_dw_pu = DW_pu * ec_w
    
    '石灰石
    AW_lime = AW_sio2 * mr_lime
    For n = 1 To 3
        E_lime(n) = AW_lime * ec_lime(n)
    Next n
    C_lime = AW_lime * ec_lime(4)
        
'(4)人件費
C_l_pu = LINE_P * UL_P * f_l * c_l

'(5)装置占有面積
AREA_pu = N_p * uarea_p * f_a

'B-2-4 高純度SiO2製造工程の建屋建設および空調・照明--------------------------

'(1)建屋建設の年間エネルギー
For n = 1 To 3
    E_build_sio2(n) = ue_build(n) * (AREA_ac + AREA_pu) / L_build
Next n

'(2)建屋建設の年経費
C_build_sio2 = c_build * (AREA_ac + AREA_pu) * r_c_build

'(3)セメント由来の年間CO2排出量 [g-C/y]
CO2_build_cement_sio2 = c_cement_p * w_build_cement * (AREA_ac + AREA_pu) / L_build * 10 ^ 6

'(4)照明・空調の年間エネルギー
E_oh_sio2_elc = e_oh * 8760 * (AREA_ac + AREA_pu)

'(5)照明・空調の年経費
C_oh_sio2 = c_e * E_oh_sio2_elc

'B-2-5 高純度SiO2の製造原単位----------------------------------------------

'(1)高純度SiO2製造の年間投入エネルギー
    '<電力> [kWh/y]
    E_sio2(1) = E_cargo(1) + E_p_cargo(1) + E_ac(1) + E_p_ac(1) + E_na2co3(1) + E_pu(1) + E_p_pu(1) + E_ma(1) + E_dw_pu(1) + E_lime(1) + E_build_sio2(1) + E_oh_sio2_elc
    '<石油> [Mcal/y]
    E_sio2(2) = E_cargo(2) + E_p_cargo(2) + E_ac(2) + E_p_ac(2) + E_na2co3(2) + E_pu(2) + E_p_pu(2) + E_ma(2) + E_dw_pu(2) + E_lime(2) + E_build_sio2(2)
    '<石炭> [Mcal/y]
    E_sio2(3) = E_cargo(3) + E_p_cargo(3) + E_ac(3) + E_p_ac(3) + E_na2co3(3) + E_pu(3) + E_p_pu(3) + E_ma(3) + E_dw_pu(3) + E_lime(3) + E_build_sio2(3)
'(2)高純度SiO2製造の年経費 [円/y]
    C_sio2 = C_cargo + C_p_cargo + C_l_cargo + C_ac + C_p_ac + C_na2co3 + C_l_ac + C_pu + C_p_pu + C_ma + C_dw_pu + C_lime + C_l_pu + C_build_sio2 + C_oh_sio2
'(3)高純度SiO2製造の年間CO2排出量 [g-C/y]
    CO2_sio2 = co2_elc * E_sio2(1) + co2_OIL * 10 ^ 6 * E_sio2(2) + co2_coal * 10 ^ 6 * E_sio2(3) + CO2_ac + CO2_build_cement_sio2

'(4)高純度SiO2のエネルギー原単位
    For n = 1 To 3
        ec_sio2(n) = E_sio2(n) / AW_sio2
    Next n
'(5)高純度SiO2の単価
    ec_sio2(4) = C_sio2 / AW_sio2
'(6)高純度SiO2のCO2排出原単位
    uco2_sio2 = CO2_sio2 / AW_sio2
    
'*****B-3:高純度カーボン製造原単位計算*****************************************

'B-3-1 原料の年間投入量---------------------------------------------------

'(1)高純度カーボン年間必要量
AW_cp = AW_sio2 * (2 * aw_c / mw_SiO2)
'(2)単位重量あたりのアセチレンブラック必要量
w_cp_ab = 1 / ((1 + r_cp_c * w_c_re) * Y_cp_h)
'(3)単位重量あたりのレゾール必要量
w_cp_re = r_cp_c / ((1 + r_cp_c * w_c_re) * Y_cp_h)
'(4)年間投入原料
AW_cp_c = (w_cp_ab + w_cp_re) * AW_cp

'B-3-2 ニーダ関連---------------------------------------------------

'(0)ニーダ一台あたりの年間原料処理量および必要台数
If AW_cp_c >= UCAP_cp_n Then
    CAP_cp_n = UCAP_cp_n
Else
    CAP_cp_n = AW_cp_c
End If
If AW_cp_c / CAP_cp_n = Int(AW_cp_c / CAP_cp_n) Then
    N_cp_n = Int(AW_cp_c / CAP_cp_n)
Else
    N_cp_n = Int(AW_cp_c / CAP_cp_n) + 1
End If
'(1)装置製造関連
    '(1-1)ニーダ総重量
    If CAP_cp_n >= MINCAP_cp_n Then
        W_cp_n = UW_cp_n * (CAP_cp_n / MINCAP_cp_n) ^ sf3 * N_cp_n * f_w
    Else
        W_cp_n = UW_cp_n * N_cp_n * f_w
    End If
    '(1-2)ニーダ製造の年間投入エネルギー
    For n = 1 To 3
        E_cp_n(n) = W_cp_n * ec_ss(n) / L_equip
    Next n
    '(1-3)ニーダ製造の年経費
    If CAP_cp_n >= MINCAP_cp_n Then
        C_cp_n = UC_cp_n * (CAP_cp_n / MINCAP_cp_n) ^ sf3 * N_cp_n * f_c * r_c_equip
    Else
        C_cp_n = UC_cp_n * N_cp_n * f_c * r_c_equip
    End If

'(2)直接投入エネルギー
    '(2-1)ニーダの電力消費原単位
    If CAP_cp_n >= MINCAP_cp_n Then
        UNITE_cp_n = (UE_cp_n / (CAP_cp_n / MINCAP_cp_n)) * (CAP_cp_n / MINCAP_cp_n) ^ sf1
    Else
        UNITE_cp_n = UE_cp_n
    End If
    '(2-2)ニーダの年間電力消費量
    E_p_cp_n_elc = AW_cp_c * UNITE_cp_n
    '(2-3)ニーダの消費電力の年経費
    C_p_cp_n = E_p_cp_n_elc * c_e

'(3)装置占有面積
    If CAP_cp_n >= MINCAP_cp_n Then
        AREA_cp_n = uarea_cp_n * (CAP_cp_n / MINCAP_cp_n) ^ sf4 * N_cp_n * f_a
    Else
        AREA_cp_n = uarea_cp_n * N_cp_n * f_a
    End If

'B-3-3 ロータリーキルン関連---------------------------------------------------

'(0)キルン一台あたりの年間原料処理量および必要台数
If AW_cp_c >= UCAP_cp_r Then
    CAP_cp_r = UCAP_cp_r
Else
    CAP_cp_r = AW_cp_c
End If
If AW_cp_c / CAP_cp_r = Int(AW_cp_c / CAP_cp_r) Then
    N_cp_r = Int(AW_cp_c / CAP_cp_r)
Else
    N_cp_r = Int(AW_cp_c / CAP_cp_r) + 1
End If

'(1)装置製造関連
    '(1-2)キルン総重量
    If CAP_cp_r >= UCAP_cp_r Then
        W_cp_r = UW_cp_r * N_cp_r * f_w
    Else
        W_cp_r = UW_cp_r * (CAP_cp_r / UCAP_cp_r) ^ sf3 * N_cp_r * f_w
    End If
    '(1-2)キルン製造の年間投入エネルギー
    For n = 1 To 3
        E_cp_r(n) = W_cp_r * ec_ss(n) / L_equip
    Next n
    '(1-3)キルン製造の年経費
    If CAP_cp_r >= UCAP_cp_r Then
        C_cp_r = UC_cp_r * N_cp_r * f_c * r_c_equip
    Else
        C_cp_r = UC_cp_r * (CAP_cp_r / UCAP_cp_r) ^ sf3 * N_cp_r * f_c * r_c_equip
    End If

'(2)直接投入エネルギー
    '(2-1)キルンの電力消費原単位
    If CAP_cp_r >= UCAP_cp_r Then
        UNITE_cp_r = UE_cp_r
    Else
        UNITE_cp_r = (UE_cp_r / (CAP_cp_r / UCAP_cp_r)) * (CAP_cp_r / UCAP_cp_r) ^ sf1
    End If
    '(2-2)キルンの年間電力消費量
    E_p_cp_r_elc = AW_cp_c * UNITE_cp_r
    '(2-3)キルンの消費電力の年経費
    C_p_cp_r = E_p_cp_r_elc * c_e

'(3)装置占有面積
    If CAP_cp_r >= UCAP_cp_r Then
        AREA_cp_r = uarea_cp_r * N_cp_r * f_a
    Else
        AREA_cp_r = uarea_cp_r * (CAP_cp_r / UCAP_cp_r) ^ sf4 * N_cp_r * f_a
    End If

'B-3-4 誘導加熱炉関連---------------------------------------------------

'誘導加熱炉必要台数
If AW_cp_c / UCAP_cp_f = Int(AW_cp_c / UCAP_cp_f) Then
    N_cp_f = Int(AW_cp_c / UCAP_cp_f)
Else
    N_cp_f = Int(AW_cp_c / UCAP_cp_f) + 1
End If

'(1)装置製造関連
    '(1-1)装置総重量
    W_cp_f = UW_cp_f * N_cp_f * f_w
    '(1-2)装置製造年間エネルギー
    For n = 1 To 3
        E_cp_f(n) = W_cp_f * ec_ss(n) / L_equip
    Next n
    '(1-3)装置年経費
    C_cp_f = UC_cp_f * N_cp_f * f_c * r_c_equip

'(2)直接投入エネルギー
E_p_cp_f_elc = AW_cp_c * UE_cp_f
C_p_cp_f = E_p_cp_f_elc * c_e

'(3)装置占有面積
AREA_cp_f = uarea_cp_f * N_cp_f * f_a

'B-3-5 投入材料関連---------------------------------------------------

'(1)アセチレンブラック
    '(1-1)年間消費量
    AW_cp_ab = w_cp_ab * AW_cp
    '(1-2)年間投入エネルギー
    For n = 1 To 3
        E_cp_ab(n) = AW_cp_ab * ec_ab(n)
    Next n
    '(1-3)年経費
    C_cp_ab = AW_cp_ab * ec_ab(4)
    
'(2)レゾール
    '(2-1)年間消費量
    AW_cp_re = w_cp_re * AW_cp
    '(2-2)年間投入エネルギー
    For n = 1 To 3
        E_cp_re(n) = AW_cp_re * ec_re(n)
    Next n
    '(2-3)年経費
    C_cp_re = AW_cp_re * ec_re(4)

'(3)塩素ガス
    '(3-1)年間消費量
    AW_cp_cl2 = w_cp_cl2 * AW_cp
    '(3-2)年間投入エネルギー
    For n = 1 To 3
        E_cp_cl2(n) = AW_cp_cl2 * ec_cl2(n)
    Next n
    '(3-3)年経費
    C_cp_cl2 = AW_cp_cl2 * ec_cl2(4)

'(4)水酸化ナトリウム
    '(4-1)年間消費量
    AW_cp_naoh = w_cp_naoh * AW_cp
    '(4-2)年間投入エネルギー
    For n = 1 To 3
        E_cp_naoh(n) = AW_cp_naoh * ec_naoh(n)
    Next n
    '(4-3)年経費
    C_cp_naoh = AW_cp_naoh * ec_naoh(4)

'(5)冷却水
    '(5-1)年間消費量
    DW_cp_f = AW_cp_c * UDW_cp_f
    '(5-2)年間投入エネルギー
    For n = 1 To 3
        E_dw_cp_f(n) = DW_cp_f * ec_dw(n)
    Next n
    '(5-3)年経費
    C_dw_cp_f = DW_cp_f * ec_w

'B-3-6 人件費、建屋、照明・空調-----------------------------------------------

'(1)人件費
If UCAP_cp_n <= MINCAP_cp_n Then
    LBR_cp = 1
Else
    If Int(N_cp_f / UL_cp_f) = N_cp_f / UL_cp_f Then
        LBR_cp = UL_cp_n * N_cp_n + UL_cp_r * N_cp_r + Int(N_cp_f / UL_cp_f)
    Else
        LBR_cp = UL_cp_n * N_cp_n + UL_cp_r * N_cp_r + (Int(N_cp_f / UL_cp_f) + 1)
    End If
End If
C_l_cp = c_l * LBR_cp * f_l

'(2)建屋建設
    '全占有面積
    AREA_cp = AREA_cp_n + AREA_cp_r + AREA_cp_f
    '年間投入エネルギー
    For n = 1 To 3
        E_build_cp(n) = ue_build(n) * AREA_cp / L_build
    Next n
    '年経費
    C_build_cp = c_build * AREA_cp * r_c_build
    'セメント由来の年間CO2排出量 [g-C/y]
    CO2_build_cement_cp = c_cement_p * w_build_cement * AREA_cp / L_build * 10 ^ 6
'(3)照明・空調
E_oh_cp_elc = e_oh * 8760 * AREA_cp
C_oh_cp = c_e * E_oh_cp_elc
    
'B-3-7 高純度カーボンの製造原単位-----------------------------------------------
    
'(1)高純度カーボン製造の年間投入エネルギー
    '<電力> [kWh/y]
    E_cp(1) = E_cp_n(1) + E_p_cp_n_elc + E_cp_r(1) + E_p_cp_r_elc + E_cp_f(1) + E_p_cp_f_elc + E_cp_ab(1) + E_cp_re(1) + E_cp_cl2(1) + E_cp_naoh(1) + E_dw_cp_f(1) + E_build_cp(1) + E_oh_cp_elc
    '<石油> [Mcal/y]
    E_cp(2) = E_cp_n(2) + E_cp_r(2) + E_cp_f(2) + E_cp_ab(2) + E_cp_re(2) + E_cp_cl2(2) + E_cp_naoh(2) + E_dw_cp_f(2) + E_build_cp(2)
    '<石炭> [Mcal/y]
    E_cp(3) = E_cp_n(3) + E_cp_r(3) + E_cp_f(3) + E_cp_ab(3) + E_cp_re(3) + E_cp_cl2(3) + E_cp_naoh(3) + E_dw_cp_f(3) + E_build_cp(3)
'(2)高純度カーボンの年経費 [円/y]
    C_cp = C_cp_n + C_p_cp_n + C_cp_r + C_p_cp_r + C_cp_f + C_p_cp_f + C_cp_ab + C_cp_re + C_cp_cl2 + C_cp_naoh + C_dw_cp_f + C_l_cp + C_build_cp + C_oh_cp
'(3)高純度カーボン製造の年間CO2排出量 [g-C/y]
    CO2_cp = co2_elc * E_cp(1) + co2_OIL * 10 ^ 6 * E_cp(2) + co2_coal * 10 ^ 6 * E_cp(3) + CO2_build_cement_cp
    
'(4)高純度カーボンのエネルギー原単位
    For n = 1 To 3
        ec_cp(n) = E_cp(n) / AW_cp
    Next n
'(5)高純度カーボンの単価
    ec_cp(4) = C_cp / AW_cp
'(6)高純度カーボンのCO2排出原単位
    uco2_cp = CO2_cp / AW_cp


'*****B-4:SOG-Si製造原単位計算***********************************************

'B-4-1 高純度SiO2還元工程-----------------------------------------------

'アーク炉必要台数
N_af = Int(AW_silicon / CAP_af)
If N_af <> AW_silicon / CAP_af Then
    N_af = N_af + 1
End If

'(1)装置製造関連
    '(1-1)グラファイトおよびレンガ重量
    If AW_silicon > CAP_af Then
        W_af_g = N_af * uw_af_g
        W_af_brick = N_af * uw_af_brick
    Else
        W_af_g = N_af * uw_af_g * (AW_silicon / CAP_af) ^ sf3
        W_af_brick = N_af * uw_af_brick * (AW_silicon / CAP_af) ^ sf3
    End If
    '(1-2)総重量
    W_af = (W_af_g + W_af_brick) * f_w
    '(1-3)年間投入エネルギー
    For n = 1 To 3
        E_af(n) = (W_af_g * ec_cp(n) + W_af_brick * ec_brick(n) + (W_af_g + W_af_brick) * ec_ss(n)) / L_equip
    Next n
    '(1-4)年経費
    If AW_silicon > CAP_af Then
        C_af = N_af * UC_af * f_c * r_c_equip
    Else
        C_af = N_af * UC_af * (AW_silicon / CAP_af) ^ sf3 * f_c * r_c_equip
    End If
    
'(2)直接投入エネルギー
    E_p_af_elc = AW_silicon * UE_af
    C_p_af = E_p_af_elc * c_e
    
'(3)投入材料
    '(3-1)Arガス
    AW_af_ar = AW_silicon * v_af_ar * sw_ar
    For n = 1 To 3
        E_af_ar(n) = AW_af_ar * ec_ar(n)
    Next n
    C_af_ar = AW_af_ar * ec_ar(4)
    '(3-2)黒鉛ノズル
    AW_af_n = AW_silicon * w_af_n
    For n = 1 To 3
        E_af_n(n) = AW_af_n * ec_cp(n)
    Next n
    C_af_n = AW_af_n * ec_cp(4)
    '(3-3)黒鉛電極
    AW_af_e = AW_silicon * w_af_e
    For n = 1 To 3
        E_af_e(n) = AW_af_e * ec_cp(n)
    Next n
    C_af_e = AW_af_e * ec_cp(4)
    
'(4)人件費
    C_l_af = N_af * UL_af * f_l * c_l

'(5)装置占有面積
    AREA_af = N_af * uarea_af * f_a
    
'(6)反応からのCO2直接排出
    CO2_af = AW_silicon * (2 * aw_c / aw_Si) * 10 ^ 6
    
'(7)CO燃焼熱
    E_co_af = AW_silicon * (2 * mw_co / aw_Si) * h_co

'B-4-2 Si精製工程-----------------------------------------------

'シリカ坩堝必要台数およびライン数
N_dec = Int(AW_silicon / CAP_dec)
If N_dec <> AW_silicon / CAP_dec Then
    N_dec = N_dec + 1
End If
LINE_dec = Int(N_dec / N0_dec)
If LINE_dec <> N_dec / N0_dec Then
    LINE_dec = LINE_dec + 1
End If

'(1)装置製造関連
    W_dec_sio2 = N_dec * uw_dec_sio2
    W_dec = 2 * W_dec_sio2
    For n = 1 To 3
        E_dec(n) = W_dec / 2 * (ec_sio2(n) + ec_ss(n)) / L_equip
    Next n
    If N_dec >= N0_dec Then
        C_dec = LINE_dec * UC_dec * f_c * r_c_equip
    Else
        C_dec = UC_dec * (N_dec / N0_dec) ^ sf2 * f_c * r_c_equip
    End If

'(2)直接投入エネルギー
    E_p_dec_elc = AW_silicon * UE_dec
    C_p_dec = E_p_dec_elc * c_e

'(3)人件費
    C_l_dec = LINE_dec * UL_dec * f_l * c_l

'(4)装置占有面積
    AREA_dec = N_dec * uarea_dec * f_a

'B-4-3 方向性凝固工程-----------------------------------------------

'累積Si処理量
AW_si_ds = AW_silicon * Y_dec

For n = 1 To n_ds - 1
    AW_si_ds = AW_si_ds + AW_si_ds * Y_ds
Next n
'シリカ坩堝必要台数およびライン数
NN_ds = Int(AW_si_ds / CAP_ds)
If NN_ds <> AW_si_ds / CAP_ds Then
    NN_ds = NN_ds + 1
End If
LINE_ds = Int(NN_ds / N0_ds)
If LINE_ds <> NN_ds / N0_ds Then
    LINE_ds = LINE_ds + 1
End If

'(1)装置製造関連
    W_ds = NN_ds * w_ds_brick
    For n = 1 To 3
        E_ds(n) = W_ds * (ec_brick(n) + ec_ss(n)) / L_equip
    Next n
    If NN_ds >= N0_ds Then
        C_ds = LINE_ds * UC_ds * f_c * r_c_equip
    Else
        C_ds = UC_ds * (NN_ds / N0_ds) ^ sf2 * f_c * r_c_equip
    End If
'(2)直接投入エネルギー
    E_p_ds_elc = AW_si_ds * UE_ds
    C_p_ds = E_p_ds_elc * c_e
'(3)人件費
    C_l_ds = LINE_ds * UL_ds * f_l * c_l
'(4)装置占有面積
    AREA_ds = NN_ds * uarea_ds * f_a

'B-4-4 SOG-Si製造工程の建屋および照明・空調---------------------------------------

'占有面積
AREA_sog = AREA_af + AREA_dec + AREA_ds
'建屋の年間投入エネルギー
For n = 1 To 3
    E_build_sog(n) = ue_build(n) * AREA_sog / L_build
Next n
'年経費
C_build_sog = c_build * AREA_sog * r_c_build
'セメント由来の年間CO2排出量 [g-C/y]
CO2_build_cement_sog = c_cement_p * w_build_cement * AREA_sog / L_build * 10 ^ 6
'照明・空調の年間エネルギー
E_oh_sog_elc = e_oh * 8760 * AREA_sog
'照明・空調の年経費
C_oh_sog = c_e * E_oh_sog_elc

'B-4-5 SOG-Si製造原単位---------------------------------------

'(1)高純度カーボン必要量
    W_cp = AW_sio2 * (2 * aw_c / mw_SiO2)
'(2)高純度カーボン由来のエネルギー・コスト・CO2
    For n = 1 To 3
        E2_cp(n) = ec_cp(n) * W_cp
    Next n
    C2_cp = ec_cp(4) * W_cp
    CO22_cp = uco2_cp * W_cp
'(3)SOG-Si製造の年間投入エネルギー・コスト・CO2
    '<電気(kWh/y)>
    E_sog(1) = E_sio2(1) + E2_cp(1) + E_af(1) + E_p_af_elc + E_af_ar(1) + E_af_n(1) + E_af_e(1) + E_dec(1) + E_p_dec_elc + E_ds(1) + E_p_ds_elc + E_build_sog(1) + E_oh_sog_elc
    '<石油(Mcal/y)>
    E_sog(2) = E_sio2(2) + E2_cp(2) + E_af(2) + E_af_ar(2) + E_af_n(2) + E_af_e(2) + E_dec(2) + E_ds(2) + E_build_sog(2)
    '<石炭(Mcal/y)>
    E_sog(3) = E_sio2(3) + E2_cp(3) + E_af(3) + E_af_ar(3) + E_af_n(3) + E_af_e(3) + E_dec(3) + E_ds(3) + E_build_sog(3)
    '<コスト(円/y)>
    C_sog = C_sio2 + C2_cp + C_af + C_p_af + C_af_ar + C_af_n + C_af_e + C_dec + C_p_dec + C_ds + C_p_ds + C_build_sog + C_oh_sog
    '<CO2(g-C/y)>
    CO2_sog = CO2_sio2 + CO2_cp + CO2_af + CO2_build_cement_sog + co2_elc * (E_af(1) + E_p_af_elc + E_af_ar(1) + E_af_n(1) + E_af_e(1) + E_dec(1) + E_p_dec_elc + E_ds(1) + E_p_ds_elc + E_build_sog(1) + E_oh_sog_elc) + co2_OIL * 10 ^ 6 * (E_af(2) + E_af_ar(2) + E_af_n(2) + E_af_e(2) + E_dec(2) + E_ds(2) + E_build_sog(2)) + co2_coal * 10 ^ 6 * (E_af(3) + E_af_ar(3) + E_af_n(3) + E_af_e(3) + E_dec(3) + E_ds(3) + E_build_sog(3))
'(4)SOG-Si製造原単位
    For n = 1 To 3
        ec_sog(n) = E_sog(n) / AW_sog
    Next n
    ec_sog(4) = C_sog / AW_sog
    uco2_sog = CO2_sog / AW_sog
    
'Debug.Print "<SOG-Si小計>"; "E_elc[kWh]="; E_sog(1); "E_oil[MJ]="; E_sog(2) * JCal; "E_coal[MJ]="; E_sog(3) * JCal; "Cost[Yen]="; C_sog; "CO2[t-C]="; CO2_sog / 10 ^ 6;

'*****B-5:poly-Siウェハ製造工程***********************************************

'B-5-1 インゴット鋳造工程---------------------------------------

'電磁鋳造炉最大処理能力
UCAP_em = 525.6 * s_ingot * r_em * sw_Si * CF_em_max
'電磁鋳造炉必要台数
N_em = Int(AW_sog / UCAP_em)
If N_em <> AW_sog / UCAP_em Then
    N_em = N_em + 1
End If
'電磁鋳造炉消費電力
s_ingot_0 = w_ingot_0 * l_ingot_0
UE_em = UE_em_0 * (s_ingot / s_ingot_0) ^ sf2
'電磁鋳造炉単位重量
UW_em = UW_em_0 * (UE_em / UE_em_0) ^ sf3
'電磁鋳造炉単価
UCAP_em_0 = 525.6 * s_ingot_0 * r_em_0 * sw_Si * CF_em_max_0
UC_em = UC_em_0 * (UCAP_em / UCAP_em_0) ^ sf4

'(1)装置製造関連
    '年間投入エネルギー
    For n = 1 To 3
        E_em(n) = UW_em * N_em * ec_ss(n) * f_w_wafer / L_equip
    Next n
    '年経費
    C_em = UC_em * N_em * f_c * r_c_equip

'(2)直接投入エネルギー
    '年間実稼働率
    CF_em = (CF_em_max / UCAP_em) * (AW_sog / N_em)
    E_p_em_elc = 8760 * CF_em * UE_em * N_em * f_e_wafer
    C_p_em = c_e * E_p_em_elc

'(3)人件費
    C_l_em = UL_em * N_em * f_l * c_l

'B-5-2 ウェハスライス工程---------------------------------------
    
'1バッチあたりウェハ生産枚数
UN_wafer = (h_ingot * n_ingot) / ((t_wafer + t_kerf) / 10 ^ 6)
'1台あたり年間最大ウェハ生産枚数
UCAP_mws = (8760 * CF_mws_max / T_mws_0) * UN_wafer
'マルチワイヤソー必要台数
N_mws = Int(AN_wafer / UCAP_mws)
If N_mws <> AW_wafer / UCAP_mws Then
    N_mws = N_mws + 1
End If
'消費電力
s_cell_0 = w_cell_0 * l_cell_0
UE_mws = UE_mws_0 * ((s_cell / s_cell_0) * (r_mws / r_mws_0) * (n_ingot / n_ingot_0)) ^ sf3
'単位重量
UW_mws = UW_mws_0 * (UE_mws / UE_mws_0) ^ sf2
'単価
UC_mws = UC_mws_0 * (UW_mws / UW_mws_0) ^ sf2

'(1)装置製造関連
    '年間投入エネルギー
    For n = 1 To 3
        E_mws(n) = UW_mws * N_mws * ec_ss(n) * f_w_wafer / L_equip
    Next n
    '年経費
    C_mws = UC_mws * N_mws * f_c * r_c_equip

'(2)直接投入エネルギー
    'MWS実働時間
    T_mws = T_mws_0 * (l_cell / l_cell_0) * (r_mws_0 / r_mws)
    E_p_mws_elc = UE_mws * T_mws * (8760 * CF_mws_max / T_mws_0) * N_mws * f_e_wafer
    C_p_mws = c_e * E_p_mws_elc

'(3)投入材料
    'SiC砥粒
    For n = 1 To 3
        E_mws_sic(n) = ec_sic(n) * w_sic * AW_ingot
    Next n
    C_mws_sic = ec_sic(4) * w_sic * AW_ingot
    'ピアノ線
    For n = 1 To 3
        E_mws_wire(n) = ec_wire(n) * w_wire * AW_ingot
    Next n
    C_mws_wire = ec_wire(4) * w_wire * AW_ingot
    'その他少量材料
    For n = 1 To 3
        E_mws_om(n) = 0.1 * (E_mws_sic(n) + E_mws_wire(n))
    Next n
    C_mws_om = 0.1 * (C_mws_sic + C_mws_wire)

'(4)人件費
    C_l_mws = UL_mws * N_mws * f_l * c_l

'B-5-3 ウェハ製造工程の建屋および空調照明-------------------------------------

'装置面積
AREA_wafer = uarea_wafer * N_mws * f_a
'建屋の年間投入エネルギー
For n = 1 To 3
    E_build_wafer(n) = ue_build(n) * AREA_wafer / L_build
Next n
'年経費
C_build_wafer = c_build * AREA_wafer * r_c_build
'セメント由来の年間CO2排出量 [g-C/y]
CO2_build_cement_wafer = c_cement_p * w_build_cement * AREA_wafer / L_build * 10 ^ 6
'照明・空調の年間エネルギー
E_oh_wafer_elc = e_oh * 8760 * AREA_wafer
'照明・空調の年経費
C_oh_wafer = c_e * E_oh_wafer_elc

'B-5-3 ウェハ製造工程の総年間エネルギー投入・CO2排出量および年経費--------------

'総年間エネルギー投入量
E_wafer(1) = E_em(1) + E_p_em_elc + E_mws(1) + E_p_mws_elc + E_mws_sic(1) + E_mws_wire(1) + E_mws_om(1) + E_build_wafer(1) + E_oh_wafer_elc
E_wafer(2) = E_em(2) + E_mws(2) + E_mws_sic(2) + E_mws_wire(2) + E_mws_om(2) + E_build_wafer(2)
E_wafer(3) = E_em(3) + E_mws(3) + E_mws_sic(3) + E_mws_wire(3) + E_mws_om(3) + E_build_wafer(3)

'総年経費
C_wafer = C_em + C_p_em + C_l_em + C_mws + C_p_mws + C_mws_sic + C_mws_wire + C_mws_om + C_l_mws + C_build_wafer + C_oh_wafer

'総年間CO2排出量
CO2_wafer = co2_elc * E_wafer(1) + co2_OIL * E_wafer(2) * 10 ^ 6 + co2_coal * E_wafer(3) * 10 ^ 6 + CO2_build_cement_wafer

'ウェハ1枚あたり原単位
    For n = 1 To 3
        ec_wafer(n) = (ec_sog(n) * AW_sog + E_wafer(n)) / AN_wafer
    Next n
    ec_wafer(4) = (ec_sog(4) * AW_sog + C_wafer) / AN_wafer
    uco2_wafer = (uco2_sog * AW_sog + CO2_wafer) / AN_wafer

'Debug.Print "<wafer小計>"; "E_elc[kWh]="; E_wafer(1); "E_oil[MJ]="; E_wafer(2) * JCal; "E_coal[MJ]="; E_wafer(3) * JCal; "Cost[Yen]="; C_wafer; "CO2[t-C]="; CO2_wafer / 10 ^ 6;
    
'*****B-6:セル製造工程***************************************************
    
'B-6-1 製造設備に関する諸量の計算--------------

For k = 1 To 8
    If k <> 4 Then
        
        '年間最大処理能力
        CAP_cell(k) = UCAP_cell_0(k) * r_cell_s * (T_cell_max / T_cell_max_0)
        
        '必要台数
        NN_cell(k) = Int(AN_cell / CAP_cell(k))
        If NN_cell(k) <> AN_cell / CAP_cell(k) Then
            NN_cell(k) = NN_cell(k) + 1
        End If
        
        '一台あたり重量
        UW_cell(k) = UW_cell_0(k) * (s_cell / s_cell_0) ^ sf4
        
        '一台あたり単価
        UC_cell(k) = UC_cell_0(k) * r_cell_s
                
    End If
Next k

'テフロン
AW_tef = uw_tef * 3 * AS_cell
For n = 1 To 3
    E_cell_tef(n) = ec_tef(n) * AW_tef / L_equip
Next n
C_cell_tef = ec_tef(4) * AW_tef * r_c_equip


'装置製造エネルギーおよび年経費
For n = 1 To 3
    E_cell_eq(n) = ec_ss(n) * f_w * (2 * NN_cell(1) * UW_cell(1) + NN_cell(2) * UW_cell(2) + NN_cell(3) * UW_cell(3) + 2 * NN_cell(5) * UW_cell(5) + 2 * NN_cell(6) * UW_cell(6) + 2 * NN_cell(7) * UW_cell(7) + NN_cell(8) * UW_cell(8)) / L_equip + E_cell_tef(n)
Next n
C_cell_eq = (2 * NN_cell(1) * UC_cell(1) + NN_cell(2) * UC_cell(2) + NN_cell(3) * UC_cell(3) + 2 * NN_cell(5) * UC_cell(5) + 2 * NN_cell(6) * UC_cell(6) + 2 * NN_cell(7) * UC_cell(7) + NN_cell(8) * UC_cell(8)) * f_c * r_c_equip + C_cell_tef


'B-6-2 直接投入エネルギー--------------
For k = 1 To 8
    If k <> 4 Then
        UE_cell_elc(k) = UE_cell_elc_0(k) * (s_cell / s_cell_0) ^ sf2
        T_cell(k) = T_cell_max * (AN_cell / (NN_cell(k) * CAP_cell(k)))
    End If
Next k

E_p_cell_elc = (UE_cell_elc(1) * T_cell(1) * NN_cell(1)) * 2 + UE_cell_elc(2) * T_cell(2) * NN_cell(2) + UE_cell_elc(3) * T_cell(3) * NN_cell(3) + 2 * (UE_cell_elc(5) * T_cell(5) * NN_cell(5) + UE_cell_elc(6) * T_cell(6) * NN_cell(6) + UE_cell_elc(7) * T_cell(7) * NN_cell(7)) + UE_cell_elc(8) * T_cell(8) * NN_cell(8)
C_p_cell_elc = c_e * E_p_cell_elc

'B-6-3 投入材料--------------

'KOH
uw_KOH = uw_KOH_0 * (s_cell / s_cell_0) ^ sf4
'IPA
uw_IPA = uw_IPA_0 * (s_cell / s_cell_0) ^ sf4
'蒸留水
uw_dw(1) = uw_dw_0(1) * (s_cell / s_cell_0) ^ sf4
uw_dw(2) = uw_dw_0(2) * (s_cell / s_cell_0) ^ sf4
'POCL3
uw_POCl3 = uw_POCl3_0 * (s_cell / s_cell_0) ^ sf4
'N2/O2
uw_N2O2(1) = uw_N2O2_0(1) * (s_cell / s_cell_0) ^ sf4
uw_N2O2(2) = uw_N2O2_0(2) * (s_cell / s_cell_0) ^ sf4
uw_N2O2(3) = uw_N2O2_0(3) * (s_cell / s_cell_0) ^ sf4
'石英管
uw_tube = uw_tube_0 * un_tube * (s_cell / s_cell_0) ^ sf4
'石英ホルダ
uw_folder = uw_folder_0 * un_folder * (s_cell / s_cell_0) ^ sf4
'TPT
uw_TPT = uw_TPT_0 * (s_cell / s_cell_0) ^ sf4
'Ag-Alペースト
uw_AgAl = uw_AgAl_0 * (s_cell / s_cell_0) ^ sf4
'Agペースト
uw_Ag = uw_Ag_0 * (s_cell / s_cell_0) ^ sf4
'スクリーン
uw_screen = uw_screen_0 * un_screen * (s_cell / s_cell_0) ^ sf4
'Xeランプ
un_Xe = un_Xe_0 * (s_cell / s_cell_0) ^ sf4
'ケース
un_case = un_case_0 * (s_cell / s_cell_0) ^ sf4


For n = 1 To 3
    E_cell_KOH(n) = 2 * uw_KOH * AN_cell * ec_KOH(n)
    E_cell_IPA(n) = 2 * uw_IPA * AN_cell * ec_IPA(n)
    E_cell_dw(n) = (2 * uw_dw(1) + uw_dw(2)) * AN_cell * ec_dw(n)
    E_cell_POCl3(n) = uw_POCl3 * AN_cell * ec_POCl3(n)
    E_cell_N2O2(n) = (uw_N2O2(1) + uw_N2O2(2) + uw_N2O2(3)) * AN_cell * ec_N2O2(n)
    E_cell_tube(n) = uw_tube * AN_cell * ec_tube(n)
    E_cell_folder(n) = uw_folder * AN_cell * ec_folder(n)
    E_cell_TPT(n) = uw_TPT * AN_cell * ec_TPT(n)
    E_cell_AgAl(n) = uw_AgAl * AN_cell * ec_AgAl(n)
    E_cell_Ag(n) = uw_Ag * AN_cell * ec_AgAl(n)
    E_cell_screen(n) = 2 * uw_screen * AN_cell * ec_screen(n)
    E_cell_Xe(n) = un_Xe * AN_cell * ec_Xe(n)
    E_cell_case(n) = un_case * AN_cell * ec_case(n)
Next n
C_cell_KOH = 2 * uw_KOH * AN_cell * ec_KOH(4)
C_cell_IPA = 2 * uw_IPA * AN_cell * ec_IPA(4)
C_cell_dw = (2 * uw_dw(1) + uw_dw(2)) * AN_cell * ec_dw(4)
C_cell_POCl3 = uw_POCl3 * AN_cell * ec_POCl3(4)
C_cell_N2O2 = (uw_N2O2(1) + uw_N2O2(2) + uw_N2O2(3)) * AN_cell * ec_N2O2(4)
C_cell_tube = un_tube * AN_cell * ec_tube(4)
C_cell_folder = un_folder * AN_cell * ec_folder(4)
C_cell_TPT = uw_TPT * AN_cell * ec_TPT(4)
C_cell_AgAl = uw_AgAl * AN_cell * ec_AgAl(4)
C_cell_Ag = uw_Ag * AN_cell * ec_AgAl(4)
C_cell_screen = 2 * un_screen * AN_cell * ec_screen(4)
C_cell_Xe = un_Xe * AN_cell * ec_Xe(4)
C_cell_case = un_case * AN_cell * ec_case(4)

'B-6-4 人件費--------------
C_l_cell = c_l * f_l * (UL_cell(1) * NN_cell(1) * 2 + UL_cell(2) * NN_cell(2) + UL_cell(3) * NN_cell(3) + UL_cell(4) * NN_cell(5) * 2 + UL_cell(5) * NN_cell(8))

'B-6-5 建屋建設および空調照明--------------

AREA_cell = f_a * (uarea_cell(1) * NN_cell(1) * 2 + uarea_cell(2) * NN_cell(2) + uarea_cell(3) * NN_cell(3) + uarea_cell(5) * NN_cell(5) * 2 + uarea_cell(6) * NN_cell(6) * 2 + uarea_cell(7) * NN_cell(7) * 2 + uarea_cell(8) * NN_cell(8))

'建屋の年間投入エネルギー
For n = 1 To 3
    E_build_cell(n) = ue_build(n) * AREA_cell / L_build
Next n
'年経費
C_build_cell = c_build * AREA_cell * r_c_build
'セメント由来の年間CO2排出量 [g-C/y]
CO2_build_cement_cell = c_cement_p * w_build_cement * AREA_cell / L_build * 10 ^ 6
'照明・空調の年間エネルギー
E_oh_cell_elc = e_oh * 8760 * AREA_cell
'照明・空調の年経費
C_oh_cell = c_e * E_oh_cell_elc

'B-6-6 セル製造工程の総年間エネルギー投入・CO2排出量および年経費--------------

E_cell(1) = E_cell_eq(1) + E_p_cell_elc + E_cell_KOH(1) + E_cell_IPA(1) + E_cell_dw(1) + E_cell_POCl3(1) + E_cell_N2O2(1) + E_cell_tube(1) + E_cell_folder(1) + E_cell_TPT(1) + E_cell_AgAl(1) + E_cell_Ag(1) + E_cell_screen(1) + E_cell_Xe(1) + E_cell_case(1) + E_build_cell(1) + E_oh_cell_elc
E_cell(2) = E_cell_eq(2) + E_cell_KOH(2) + E_cell_IPA(2) + E_cell_dw(2) + E_cell_POCl3(2) + E_cell_N2O2(2) + E_cell_tube(2) + E_cell_folder(2) + E_cell_TPT(2) + E_cell_AgAl(2) + E_cell_Ag(2) + E_cell_screen(2) + E_cell_Xe(2) + E_cell_case(2) + E_build_cell(2)
E_cell(3) = E_cell_eq(3) + E_cell_KOH(3) + E_cell_IPA(3) + E_cell_dw(3) + E_cell_POCl3(3) + E_cell_N2O2(3) + E_cell_tube(3) + E_cell_folder(3) + E_cell_TPT(3) + E_cell_AgAl(3) + E_cell_Ag(3) + E_cell_screen(3) + E_cell_Xe(3) + E_cell_case(3) + E_build_cell(3)
C_cell = C_cell_eq + C_p_cell_elc + C_cell_KOH + C_cell_IPA + C_cell_dw + C_cell_POCl3 + C_cell_N2O2 + C_cell_tube + C_cell_folder + C_cell_TPT + C_cell_AgAl + C_cell_Ag + C_cell_screen + C_cell_Xe + C_cell_case + C_l_cell + C_build_cell + C_oh_cell
CO2_cell = co2_elc * E_cell(1) + co2_OIL * E_cell(2) * 10 ^ 6 + co2_coal * E_cell(3) * 10 ^ 6 + CO2_build_cement_cell

'Debug.Print "<Cell小計>"; "E_elc[kWh]="; E_cell(1); "E_oil[MJ]="; E_cell(2) * JCal; "E_coal[MJ]="; E_cell(3) * JCal; "Cost[Yen]="; C_cell; "CO2[t-C]="; CO2_cell / 10 ^ 6;

'B-6-7 セル一枚の製造原単位--------------

    For n = 1 To 3
        ec_cell(n) = (ec_wafer(n) * AN_wafer + E_cell(n)) / AN_cell
    Next n
    ec_cell(4) = (ec_wafer(4) * AN_wafer + C_cell) / AN_cell
    uco2_cell = (uco2_wafer * AN_wafer + CO2_cell) / AN_cell


'*****B-7:モジュール組立工程**********************************************

'B-7-1 モジュール組立工程の設備製造--------------

'必要ライン数
LINE_module = Int(AS_module / (UCAP_module * CF_module_max * 365))
If LINE_module <> AS_module / (UCAP_module * CF_module_max * 365) Then
    LINE_module = LINE_module + 1
End If

'年間エネルギー投入量および年経費
For n = 1 To 3
    E_module_eq(n) = ec_ss(n) * f_w * LINE_module * UW_module / L_equip
Next n
C_module_eq = f_c * r_c_equip * LINE_module * UC_module

'B-7-2 直接投入エネルギー--------------

'実稼働率
CF_module = AS_module / (UCAP_module * LINE_module * 365)
'年間電力消費量およびその年経費
E_p_module_elc = UE_module * LINE_module * CF_module * 8760
C_p_module = c_e * E_p_module_elc

'B-7-3 投入材料--------------

'(1)保護ガラス
For n = 1 To 3
    E_glass(n) = ec_glass(n) * sw_glass * (t_glass * l_module * w_module) * AN_module
Next n
C_glass = ec_glass(4) * sw_glass * (t_glass * l_module * w_module) * AN_module
'(2)アルミフレーム
For n = 1 To 3
    E_Al(n) = ec_Al(n) * 2 * (w_Al_l * l_module + w_Al_w * w_module) * AN_module
Next n
C_Al = ec_Al(4) * 2 * (w_Al_l * l_module + w_Al_w * w_module) * AN_module
'(3)配線材料
For n = 1 To 3
    E_wire(n) = ec_Cu(n) * w_Cu * AS_module
Next n
C_wire = ec_Cu(4) * w_Cu * AS_module
'(4)充填材
For n = 1 To 3
    E_EVA(n) = ec_EVA(n) * w_EVA * AS_module
Next n
C_EVA = ec_EVA(4) * w_EVA * AS_module
'(5)裏面シート
For n = 1 To 3
    E_BS(n) = (ec_Tedlar(n) * w_Tedlar + ec_Al(n) * w_AS) * AS_module
Next n
C_BS = ec_BS(4) * (w_Tedlar + w_AS) * AS_module
'(6)周辺シール
For n = 1 To 3
    E_rubber(n) = (ec_rubber(n) * w_rubber) * AS_module
Next n
C_rubber = ec_rubber(4) * w_rubber * AS_module

'B-7-4 人件費--------------
C_l_module = c_l * f_l * UL_module * LINE_module

'B-7-5 建屋建設および空調照明--------------

AREA_module = f_a * uarea_module * LINE_module

'建屋の年間投入エネルギー
For n = 1 To 3
    E_build_module(n) = ue_build(n) * AREA_module / L_build
Next n
'年経費
C_build_module = c_build * AREA_module * r_c_build

'セメント由来の年間CO2排出量 [g-C/y]
CO2_build_cement_module = c_cement_p * w_build_cement * AREA_module / L_build * 10 ^ 6
'照明・空調の年間エネルギー
E_oh_module_elc = e_oh * 8760 * AREA_module
'照明・空調の年経費
C_oh_module = c_e * E_oh_module_elc

'B-7-6 モジュール組立工程の総年間エネルギー投入・CO2排出量および年経費--------------

E_module(1) = E_module_eq(1) + E_p_module_elc + (E_glass(1) + E_Al(1) + E_wire(1) + E_EVA(1) + E_BS(1) + E_rubber(1)) + E_build_module(1) + E_oh_module_elc
E_module(2) = E_module_eq(2) + (E_glass(2) + E_Al(2) + E_wire(2) + E_EVA(2) + E_BS(2) + E_rubber(2)) + E_build_module(2)
E_module(3) = E_module_eq(3) + (E_glass(3) + E_Al(3) + E_wire(3) + E_EVA(3) + E_BS(3) + E_rubber(3)) + E_build_module(3)
C_module = C_module_eq + C_p_module + (C_glass + C_Al + C_wire + C_EVA + C_BS + C_rubber) + C_l_module + C_build_module + C_oh_module
CO2_module = co2_elc * E_module(1) + co2_OIL * E_module(2) * 10 ^ 6 + co2_coal * E_module(3) * 10 ^ 6 + CO2_build_cement_module

'Debug.Print "<module小計>"; "E_elc[kWh]="; E_module(1); "E_oil[MJ]="; E_module(2) * JCal; "E_coal[MJ]="; E_module(3) * JCal; "Cost[Yen]="; C_module; "CO2[t-C]="; CO2_module / 10 ^ 6;

'B-7-7 モジュールの製造原単位--------------

    'モジュールの総投入量
    For n = 1 To 3  '[kWh/y],[Mcal/y]
        Total_ec_module(n) = ec_cell(n) * AN_cell + E_module(n)
    Next n
    Total_ec_module(4) = ec_cell(4) * AN_cell + C_module    '[Yen/y]
    Total_co2_module = uco2_cell * AN_cell + CO2_module     '[g-C/y]
    
    'モジュール1枚あたり
    For n = 1 To 3  '[kWh/module],[Mcal/module]
        ec_module(n) = Total_ec_module(n) / AN_module
    Next n
    ec_module(4) = Total_ec_module(4) / AN_module   '[Yen/module]
    uco2_module = Total_co2_module / AN_module      '[g-C/module]
    
    'モジュールWあたり
    For n = 1 To 3  '[kWh/W_module],[Mcal/W_module]
        ec_perW_module(n) = Total_ec_module(n) / AC_module
    Next n
    ec_perW_module(4) = Total_ec_module(4) / AC_module  '[Yen/W_module]
    uco2_perW_module = Total_co2_module / AC_module     '[g-C/W_module]
'    Debug.Print "<module total>"; "elc[kWh/W]"; ec_perW_module(1); "oil[MJ/W]"; ec_perW_module(2) * JCal; "coal[MJ/W]"; ec_perW_module(3) * JCal; "cost[Yen/W]"; ec_perW_module(4); "co2[g-C/W]"; uco2_perW_module

'Debug.Print AN_module
'Debug.Print "<module total>"; "elc[kWh/W]"; ec_perW_module(1); "oil[MJ/W]"; ec_perW_module(2) * JCal; "coal[MJ/W]"; ec_perW_module(3) * JCal; "cost[Yen/W]"; ec_perW_module(4); "co2[g-C/W]"; uco2_perW_module


End Sub '==========多結晶SiモジュールのLCA計算終了==========



'===BOSのLCA計算====
Sub system_LCA()

'技術前提条件ファイルの読み込み
filename3 = "RooftopBOS.csv"
Open filename3 For Input As #3
n = 0
Do While Not EOF(3)
    n = n + 1
    For m = 1 To 9
        Input #3, dbs3(n, m)
        'Debug.Print dbs(n, m)
    Next m
    'Debug.Print n, dbs1(n, 7)
Loop
Close #3

'*****データ配列から変数名への格納*****
'架台
AS_steal = dbs3(2, 3)
AS_Al = dbs3(2, 4)
AS_Cu = dbs3(2, 5)
AS_in = dbs3(2, 6)
AS_elc = dbs3(2, 7)
AS_OIL = dbs3(2, 8)
AS_cost = dbs3(2, 9)
'パワーコンディショナ
PC_steal = dbs3(3, 3)
PC_Al = dbs3(3, 4)
PC_Cu = dbs3(3, 5)
PC_in = dbs3(3, 6)
PC_elc = dbs3(3, 7)
PC_OIL = dbs3(3, 8)
PC_cost = dbs3(3, 9)
'配線材料
WM_steal = dbs3(4, 3)
WM_Al = dbs3(4, 4)
WM_Cu = dbs3(4, 5)
WM_in = dbs3(4, 6)
WM_elc = dbs3(4, 7)
WM_OIL = dbs3(4, 8)
WM_cost = dbs3(4, 9)
'設置工事
LC_steal = dbs3(5, 3)
LC_Al = dbs3(5, 4)
LC_Cu = dbs3(5, 5)
LC_in = dbs3(5, 6)
LC_elc = dbs3(5, 7)
LC_OIL = dbs3(5, 8)
LC_cost = dbs3(5, 9)

capacity = PV_LCA8.Label29.Caption  '[kW/system]
Life_module = PV_LCA8.Label31.Caption '[年]
Life_pc = PV_LCA8.Label33.Caption '[年]

'モジュールの投入量 [kWh] [Mcal] [Yen] [g-C]
    For n = 1 To 4
        Req_sys_module(n) = ec_perW_module(n) * (capacity * 10 ^ 3)
    Next n
    Req_sys_module(5) = uco2_perW_module * (capacity * 10 ^ 3)

'パワーコンディショナの投入量
      
    N_pc = Int(Life_module / Life_pc)   'システム耐用期間中の必要台数
    If N_pc <> Life_module / Life_pc Then
        N_pc = N_pc + 1
    End If

    '(1)鉄
    Steal_pc = capacity * PC_steal * N_pc            '鉄必要量[t]
    EC_pc_steal(1) = Steal_pc * ec_steal2(1)    '鉄製造の電力[kWh]
    EC_pc_steal(2) = Steal_pc * ec_steal2(2)    '鉄製造の石油[Mcal]
    EC_pc_steal(3) = Steal_pc * ec_steal2(3)    '鉄製造の石炭[Mcal]
    '(2)アルミ
    Al_pc = capacity * PC_Al * N_pc       'アルミ必要量[t]
    EC_pc_Al(1) = Al_pc * ec_Al(1)  'アルミ製造の電力[kWh]
    EC_pc_Al(2) = Al_pc * ec_Al(2)  'アルミ製造の石油[Mcal]
    EC_pc_Al(3) = Al_pc * ec_Al(3)  'アルミ製造の石炭[Mcal]
    '(3)銅
    Cu_pc = capacity * PC_Cu * N_pc       '銅必要量[t]
    EC_pc_Cu(1) = Cu_pc * ec_Cu(1)  '銅製造の電力[kWh]
    EC_pc_Cu(2) = Cu_pc * ec_Cu(2)  '銅製造の石油[Mcal]
    EC_pc_Cu(3) = Cu_pc * ec_Cu(3)  '銅製造の石炭[Mcal]
    '(4)絶縁材料
    In_pc = capacity * PC_in * N_pc       '絶縁材料必要量[t]
    EC_pc_in(1) = In_pc * ec_in(1)  '絶縁材料製造の電力[kWh]
    EC_pc_in(2) = In_pc * ec_in(2)  '絶縁材料製造の石油[Mcal]
    EC_pc_in(3) = In_pc * ec_in(3)  '絶縁材料製造の石炭[Mcal]
    
    '総投入電力 [kWh]
    Req_sys_pc(1) = (EC_pc_steal(1) + EC_pc_Al(1) + EC_pc_Cu(1) + EC_pc_in(1)) + PC_elc * capacity
    '総投入石油 [Mcal]
    Req_sys_pc(2) = (EC_pc_steal(2) + EC_pc_Al(2) + EC_pc_Cu(2) + EC_pc_in(2)) + PC_OIL * capacity
    '総投入石炭 [Mcal]
    Req_sys_pc(3) = (EC_pc_steal(3) + EC_pc_Al(3) + EC_pc_Cu(3) + EC_pc_in(3))
    '総コスト [Yen]　※コストは初期コスト(1台目のみ)
    Req_sys_pc(4) = PC_cost * capacity
    '総CO2排出量 [g-C]
    Req_sys_pc(5) = co2_elc * Req_sys_pc(1) + co2_OIL * Req_sys_pc(2) + co2_coal * Req_sys_pc(3)

'架台の投入量
    eata_module = (eata_cell * eata_r) * (n_cell * s_cell) / (l_module * w_module) 'モジュール効率
    AREA_as = (capacity * 10 ^ 3) / (eata_module * 10 ^ 3) '架台所要面積
    '(1)鉄
    Steal_as = AREA_as * AS_steal              '鉄必要量[t]
    EC_as_steal(1) = Steal_as * ec_steal1(1)    '鉄製造の電力[kWh]
    EC_as_steal(2) = Steal_as * ec_steal1(2)    '鉄製造の石油[Mcal]
    EC_as_steal(3) = Steal_as * ec_steal1(3)    '鉄製造の石炭[Mcal]

    Req_sys_as(1) = EC_as_steal(1)  '総投入電力 [kWh]
    Req_sys_as(2) = EC_as_steal(2)  '総投入石油 [Mcal]
    Req_sys_as(3) = EC_as_steal(3)  '総投入石炭 [Mcal]
    Req_sys_as(4) = AREA_as * AS_cost '総コスト [Yen]
    Req_sys_as(5) = co2_elc * Req_sys_as(1) + co2_OIL * Req_sys_as(2) + co2_coal * Req_sys_as(3)

'配線材料の投入量
    EC_wm_cu(1) = WM_Cu * ec_Cu(1)    '銅製造の電力[kWh]
    EC_wm_cu(2) = WM_Cu * ec_Cu(2)    '銅製造の石油[Mcal]
    EC_wm_cu(3) = WM_Cu * ec_Cu(3)    '銅製造の石炭[Mcal]
    
    Req_sys_wm(1) = EC_wm_cu(1)     '総投入電力 [kWh]
    Req_sys_wm(2) = EC_wm_cu(2)     '総投入石油 [Mcal]
    Req_sys_wm(3) = EC_wm_cu(3)     '総投入石炭 [Mcal]
    Req_sys_wm(4) = WM_cost         '総コスト [Yen]
    Req_sys_wm(5) = co2_elc * Req_sys_wm(1) + co2_OIL * Req_sys_wm(2) + co2_coal * Req_sys_wm(3)

'設置工事費用
    Req_sys_LC = LC_cost
    
End Sub

'===計算の総括====

Sub summary()

Cell_type = PV_LCA8.Label3.Caption
Cell_efficiency = PV_LCA8.Label18.Caption
Prod_scale = PV_LCA8.Label5.Caption
system_type = PV_LCA8.Label26.Caption
grid = PV_LCA8.Label27.Caption
Irradiation = PV_LCA8.Label8.Caption
PR = PV_LCA8.Label23.Caption
r = PV_LCA8.Label21.Caption
om = PV_LCA8.Label22.Caption
n_ds = PV_LCA8.Label35.Caption

Open "LCA_summary.csv" For Output As #4
Write #4, "Selected Assumption"
Write #4, "", "太陽電池種類", Cell_type
Write #4, "", "セル効率 [%]", Cell_efficiency
Write #4, "", "年間セル生産規模 [MW/年]", Prod_scale
Write #4, "", "製造技術オプション", Tech_type
Write #4, "", "方向性凝固回数", n_ds
Write #4, "", "評価対象とするシステム形態", system_type
Write #4, "", "系統連系の有無", grid
Write #4, "", "システム容量 [kW]", capacity
Write #4, "", "モジュール耐用年数 [年]", Life_module
Write #4, "", "パワーコンディショナ耐用年数 [年]", Life_pc
Write #4, "", "年間日射量 [kWh/m2/年]", Irradiation
Write #4, "", "システム出力係数", PR
Write #4, "", "系統電力の熱効率", eata_elc
Write #4, "", "系統電力のCO2排出原単位 [g-C/kWh]", co2_elc
Write #4, "", "年金利 [%/年]", r
Write #4, "", "年間修繕・保守費率 [%]", om
Write #4, "", "出荷価格係数", alpha

Power = capacity * Irradiation * PR
Write #4, "Annual Power Generation", "", Int(Power + 0.5), "[kWh/年]"

Power_MJ = Power / eata_elc * kWhMJ
Write #4, "", "(一次エネルギー換算)", Int(Power_MJ / 10 ^ 3 + 0.5), "[GJ/年]", "平均電源構成基準"

CO2_recover = Power * co2_elc   '[g-C/年]
Write #4, "", "(年間CO2排出回避量)", Int(CO2_recover / 10 ^ 3 + 0.5), "[kg-C/年]", "平均電源構成基準"

Write #4, "LCA Result", "Process", "ELC Input [kWh/system]", "OIL Input [MJ/system]", "COAL Input [MJ/system]", " Total Primary Energy Input [GJ/system]", "CO2 emission [kg-C/system]", "Price [千円/system]"

'Debug.Print E_sog(1), AC_module, capacity

'システムあたりのSOG-Si製造工程における投入エネルギーおよびCO2、初期価格
ELC_sys_sog = E_sog(1) / AC_module * capacity * 10 ^ 3   '[kWh/system]
OIL_sys_sog = E_sog(2) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
COAL_sys_sog = E_sog(3) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
PER_sys_sog = ELC_sys_sog / eata_elc * kWhMJ + OIL_sys_sog + COAL_sys_sog '[MJ/system]
CO2_sys_sog = CO2_sog / AC_module * capacity * 10 ^ 3  '[g-C/system]
COST_sys_sog = alpha * (C_sog / AC_module * capacity * 10 ^ 3) '[Yen/system]
Write #4, "", "SOG-Si Production", Int(ELC_sys_sog + 0.5), Int(OIL_sys_sog + 0.5), Int(COAL_sys_sog + 0.5), Int(PER_sys_sog / 10 ^ 3 + 0.5), Int(CO2_sys_sog / 10 ^ 3 + 0.5), Int(COST_sys_sog / 10 ^ 3 + 0.5)

'システムあたりのウェハ製造工程における投入エネルギーおよびCO2、初期価格
ELC_sys_wafer = E_wafer(1) / AC_module * capacity * 10 ^ 3   '[kWh/system]
OIL_sys_wafer = E_wafer(2) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
COAL_sys_wafer = E_wafer(3) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
PER_sys_wafer = ELC_sys_wafer / eata_elc * kWhMJ + OIL_sys_wafer + COAL_sys_wafer '[MJ/system]
CO2_sys_wafer = CO2_wafer / AC_module * capacity * 10 ^ 3  '[g-C/system]
COST_sys_wafer = alpha * (C_wafer / AC_module * capacity * 10 ^ 3) '[Yen/system]
Write #4, "", "Wafer Production", Int(ELC_sys_wafer + 0.5), Int(OIL_sys_wafer + 0.5), Int(COAL_sys_wafer + 0.5), Int(PER_sys_wafer / 10 ^ 3 + 0.5), Int(CO2_sys_wafer / 10 ^ 3 + 0.5), Int(COST_sys_wafer / 10 ^ 3 + 0.5)

'システムあたりのセル製造工程における投入エネルギーおよびCO2、初期価格
ELC_sys_cell = EC_as_steal(1)
OIL_sys_cell = E_cell(2) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
COAL_sys_cell = E_cell(3) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
PER_sys_cell = ELC_sys_cell / eata_elc * kWhMJ + OIL_sys_cell + COAL_sys_cell '[MJ/system]
CO2_sys_cell = CO2_cell / AC_module * capacity * 10 ^ 3  '[g-C/system]
COST_sys_cell = alpha * (C_cell / AC_module * capacity * 10 ^ 3) '[Yen/system]
Write #4, "", "Cell Processing", Int(ELC_sys_cell + 0.5), Int(OIL_sys_cell + 0.5), Int(COAL_sys_cell + 0.5), Int(PER_sys_cell / 10 ^ 3 + 0.5), Int(CO2_sys_cell / 10 ^ 3 + 0.5), Int(COST_sys_cell / 10 ^ 3 + 0.5)

'システムあたりのモジュール工程における投入エネルギーおよびCO2、初期価格
ELC_sys_module = E_module(1) / AC_module * capacity * 10 ^ 3   '[kWh/system]
OIL_sys_module = E_module(2) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
COAL_sys_module = E_module(3) * JCal / AC_module * capacity * 10 ^ 3  '[MJ/system]
PER_sys_module = ELC_sys_module / eata_elc * kWhMJ + OIL_sys_module + COAL_sys_module '[MJ/system]
CO2_sys_module = CO2_module / AC_module * capacity * 10 ^ 3 '[g-C/system]
COST_sys_module = alpha * (C_module / AC_module * capacity * 10 ^ 3) '[Yen/system]
Write #4, "", "Module Fabrication", Int(ELC_sys_module + 0.5), Int(OIL_sys_module + 0.5), Int(COAL_sys_module + 0.5), Int(PER_sys_module / 10 ^ 3 + 0.5), Int(CO2_sys_module / 10 ^ 3 + 0.5), Int(COST_sys_module / 10 ^ 3 + 0.5)

'モジュール分小計
ELC_sys_1 = ELC_sys_sog + ELC_sys_wafer + ELC_sys_cell + ELC_sys_module '[kWh/system]
OIL_sys_1 = OIL_sys_sog + OIL_sys_wafer + OIL_sys_cell + OIL_sys_module '[MJ/system]
COAL_sys_1 = COAL_sys_sog + COAL_sys_wafer + COAL_sys_cell + COAL_sys_module '[MJ/system]
PER_sys_1 = PER_sys_sog + PER_sys_wafer + PER_sys_cell + PER_sys_module '[MJ/system]
CO2_sys_1 = CO2_sys_sog + CO2_sys_wafer + CO2_sys_cell + CO2_sys_module '[g-C/system]
COST_sys_1 = COST_sys_sog + COST_sys_wafer + COST_sys_cell + COST_sys_module '[Yen/system]
Write #4, "", "Module subtotal", Int(ELC_sys_1 + 0.5), Int(OIL_sys_1 + 0.5), Int(COAL_sys_1 + 0.5), Int(PER_sys_1 / 10 ^ 3 + 0.5), Int(CO2_sys_1 / 10 ^ 3 + 0.5), Int(COST_sys_1 / 10 ^ 3 + 0.5)

'システムあたりのBOS投入エネルギーおよびCO2、初期価格
ELC_sys_BOS = Req_sys_pc(1) + Req_sys_as(1) + Req_sys_wm(1) '[kWh/system]
OIL_sys_BOS = (Req_sys_pc(2) + Req_sys_as(2) + Req_sys_wm(2)) * JCal '[MJ/system]
COAL_sys_BOS = (Req_sys_pc(3) + Req_sys_as(3) + Req_sys_wm(3)) * JCal '[MJ/system]
PER_sys_BOS = ELC_sys_BOS / eata_elc * kWhMJ + OIL_sys_BOS + COAL_sys_BOS '[MJ/system]
CO2_sys_BOS = Req_sys_pc(5) + Req_sys_as(5) + Req_sys_wm(5)  '[g-C/system]
COST_sys_BOS = alpha * (Req_sys_pc(4) + Req_sys_as(4) + Req_sys_wm(4)) + Req_sys_LC '[Yen/system]
    '※初期コストに含まれるのは、架台、配線材料、PC一台分、設置工事
Write #4, "", "Balance of System", Int(ELC_sys_BOS + 0.5), Int(OIL_sys_BOS + 0.5), Int(COAL_sys_BOS + 0.5), Int(PER_sys_BOS / 10 ^ 3 + 0.5), Int(CO2_sys_BOS / 10 ^ 3 + 0.5), Int(COST_sys_BOS / 10 ^ 3 + 0.5)

Write #4, "", "System Total", Int(ELC_sys_1 + ELC_sys_BOS + 0.5), Int(OIL_sys_1 + OIL_sys_BOS + 0.5), Int(COAL_sys_1 + COAL_sys_BOS + 0.5), Int((PER_sys_1 + PER_sys_BOS) / 10 ^ 3 + 0.5), Int((CO2_sys_1 + CO2_sys_BOS) / 10 ^ 3 + 0.5), Int((COST_sys_1 + COST_sys_BOS) / 10 ^ 3 + 0.5)

Write #4,   '空白行挿入

EPT = (PER_sys_1 + PER_sys_BOS) / Power_MJ '[年]
CO2_emission_rate = (CO2_sys_1 + CO2_sys_BOS) / (Power * Life_module) '[g-C/kWh]
CO2PT = (CO2_sys_1 + CO2_sys_BOS) / CO2_recover '[年]
Write #4, "Energy Payback Time [年]", EPT, "平均電源構成基準"
Write #4, "CO2 Payback Time [年]", CO2PT, "平均電源構成基準"
Write #4, "CO2 emission rate [g-C/kWh]", CO2_emission_rate


'発電コストの計算

r = Val(r) / 100
D_module = r * (1 + r) ^ Life_module / ((1 + r) ^ Life_module - 1)
D_pc = r * (1 + r) ^ Life_pc / ((1 + r) ^ Life_pc - 1)

 '(1)太陽電池モジュール
    ANNCOST_module = D_module * (COST_sys_sog + COST_sys_wafer + COST_sys_cell + COST_sys_module)
 '(2)架台
    ANNCOST_as = D_module * alpha * Req_sys_as(4)
 '(3)パワーコンディショナ
    If N_pc = 1 Then
        ANNCOST_pc = D_pc * alpha * Req_sys_pc(4)
    Else
        ANNCOST_pc = D_pc * alpha * Req_sys_pc(4) * (Life_pc + (Life_module - Life_pc) ^ 2 / ((1 + r) ^ Life_pc * Life_pc)) / Life_module
    End If
 '(4)配線材料
    ANNCOST_wm = D_module * alpha * Req_sys_wm(4)
 '(5)設置工事
    ANNCOST_lc = D_module * Req_sys_LC
 '(6)修繕保守
    ANNCOST_om = (COST_sys_1 + COST_sys_BOS) * Val(om) / 100


ANNCOST = ANNCOST_module + ANNCOST_as + ANNCOST_pc + ANNCOST_wm + ANNCOST_lc + ANNCOST_om
gc = ANNCOST / Power

Write #4,   '空白行挿入

Write #4, "Annual Cost [円/年]", "PV module", Int(ANNCOST_module + 0.5)
Write #4, "", "Array support & wire", Int(ANNCOST_as + ANNCOST_wm + 0.5)
Write #4, "", "Power conditioner", Int(ANNCOST_pc + 0.5)
Write #4, "", "Installation", Int(ANNCOST_lc + 0.5)
Write #4, "", "Operation & Maintenance", Int(ANNCOST_om + 0.5)
Write #4, "", "Total Annual Cost", Int(ANNCOST + 0.5)
Write #4, "Generation Cost [円/kWh]", gc

'Debug.Print gc

End Sub