Public Class CZoneOD

    Public Event PBarSet(ByVal min As Integer, ByVal max As Integer) '//프로그레시브바 초기화
    Public Event Message(ByVal Msg As String)    '//메세지표시
    Public Event PBarVal(ByVal Value As Integer) '//프로그레시브 값설정

    ''' <summary>
    ''' (main)에너지요구량 계산
    ''' </summary>
    ''' <param name="Calc_Use">에너지소요량을 추가로 계산하려면 True 를 입력</param>
    ''' <param name="Log">로그가 입력될 스트링빌더</param>
    ''' <remarks></remarks>
    Public Overridable Function Calc_LoadOfZones() As Boolean
        Dim Title As String = "에너지요구량을 계산합니다"
        RaiseEvent Message(Title & vbCrLf & "잠시만 기다려주세요")
        RaiseEvent PBarSet(0, 100)

        RaiseEvent Message(Title & vbCrLf & "데이터저장소를 확보하는중.")
        'i_count_OpWe As Integer, d_we_month As Integer, Q_h_b_we as decimal,  ==================================2010.2.21 삭제

        Dim t_24 As Integer, f_wind As Integer, t_NA As Integer, t_h_op_d As Integer, i_count_Is As Integer
        Dim Delta_theta_er As Decimal, Delta_theta_i_NA As Decimal, Delta_theta_i_WE As Decimal, theta_i As Decimal, theta_i_h As Decimal, theta_i_c As Decimal, F_f As Decimal, F__F As Decimal, F_W As Decimal, F_S As Decimal, F_V As Decimal, h_r As Decimal, c_p_arhoa As Decimal, e_wind As Decimal, n_50 As Decimal, n_ue As Decimal, n_mech As Decimal, n_inf As Decimal, A_B As Decimal, A_u As Decimal, h_u As Decimal, V As Decimal, V_ue As Decimal
        Dim H_V_inf As Decimal, H_V_mech As Decimal, H_V_ue As Decimal, H_V_z As Decimal, theta_z As Decimal, theta_u As Decimal, H_T_D As Decimal, H_T_iu As Decimal, H_T_iz As Decimal, H_T_s As Decimal, Delta_U_WB As Decimal, H_V_mech_theta As Decimal, theta_i_h_soll As Decimal, theta_i_c_soll As Decimal, theta_V_mech As Decimal ', theta_e_min as decimal, theta_e_max as decimal
        Dim C_wirk As Decimal, tau As Decimal, a As Decimal, a_tb As Decimal, A_beheizt As Decimal, f_we As Decimal, f_NA As Decimal, f_tb As Decimal, Q_h_max As Decimal, theta_i_NA As Decimal, eta_V_mech As Decimal, theta_V_mech_RLT As Decimal
        Dim Q_V_inf_sink As Decimal, Q_V_mech_sink As Decimal, Q_V_z_sink As Decimal, Q_V_sink As Decimal, Q_T_e_sink As Decimal, Q_T_u_sink As Decimal, Q_T_z_sink As Decimal, Q_T_s_sink As Decimal, Q_T_sink As Decimal, Q_S_sink As Decimal, R_se As Decimal, Q_source As Decimal, Q_sink As Decimal, Q_I_source As Decimal
        Dim Q_V_inf_source As Decimal, Q_V_mech_source As Decimal, Q_V_z_source As Decimal, Q_V_source As Decimal, Q_T_e_source As Decimal, Q_T_u_source As Decimal, Q_T_z_source As Decimal, Q_T_s_source As Decimal, Q_T_source As Decimal, Q_S_source As Decimal
        Dim Q_I_source_p As Decimal, Q_I_source_fac As Decimal, Q_I_L As Decimal, Q_l_b As Decimal, Q_w_b As Decimal, gamma As Decimal, eta As Decimal, Q_h_b As Decimal, Q_h_b_op As Decimal, g_TI As Decimal
        Dim V_mech As Decimal
        'Dim d_op_month As Integer, theta_e_month as decimal
        'Dim A_H as decimal, A_C as decimal, A_L as decimal, A_W as decimal, A_V as decimal
        Dim Q_vh_b_op As Decimal
        Dim g_blinds As Decimal, tau_e_B As Decimal, rho_e_B As Decimal, U_blinds As Decimal, tau_e_B_corr As Decimal, rho_e_B_corr As Decimal 'EN 13363-1 적용
        Dim g As Decimal, alpa As Decimal
        'Dim i_sub_wall As Integer
        Dim V_mech_outdoor As Decimal
        Dim T As String
        T = "건축"
        '//2010.02.04추가
        Dim i_count_AHU As Integer
        Dim Phi As Decimal
        '//2010.02.21추가
        Dim d_op_mth As Integer

        '//2017.06.19 추가
        Dim Q_V_win_sink As Decimal, Q_V_win_source As Decimal, H_V_win As Decimal, n_win As Decimal

        Dim d_full_mth As Integer

        RaiseEvent Message(Title & vbCrLf & "저장소 기본값을 설정하는중.")
        t_24 = 24
        F_f = 0.5 '구조체와 하늘과의 형태계수
        F__F = 0.7 '창 프레임 감소계수 __필요
        F_W = 0.9 '수직입사가 아닌경우에 대한 감소계수
        F_S = 1 '주변환경에 의한 그늘짐
        F_V = 1 '0.9 '우선 먼지팩터는 없는 것으로 가정
        h_r = 5 * 0.9 '= 5*e(모를경우 e=0.9)
        Delta_theta_er = 10
        Delta_theta_i_NA = 4 '(용도별조건)
        'theta_i = 21
        'theta_i_h = 21
        theta_i_h_soll = 20
        theta_i_c = 24
        theta_i_c_soll = 26
        theta_V_mech_RLT = 27
        'theta_e_min = -11.3
        'theta_e_max = 31.5
        c_p_arhoa = 0.34    '//먼지모름
        'n_50 = 1.5
        n_ue = 0.6
        e_wind = 0.07
        f_wind = 15
        R_se = 0.043 '외기에 직접 열전달저항 확인요망
        g_TI = 0.35 '(2-113)
        alpa = 0.4 '우선 0.4으로 설정 표면흡수율 도표 2-6 참조
        n_mech = 1 '임시
        A_u = 1
        h_u = 1 '임시

        RaiseEvent Message(Title & vbCrLf & "결과저장소 초기화")

        Work_msg("요구량계산중...")
        Frm_Work.ProgressBar1.Maximum = 26
        Frm_Work.ProgressBar1.Minimum = 0
        Frm_Work.ProgressBar1.Value = 0
        Frm_Work.TopMost = True
        Frm_Work.Show()




       











        '/////////////////////////////////////////////////////
        For i_count_HeatCool As Integer = 1 To 2 '(1)   루프회전 난방/냉방 구분해서 회전한다. 둘의 산식이 거의 비슷하므로 루프를 이용해서 계산

         

            For i_count_month As Short = 0 To 12 '(2)    --월전체의 루프?

                'If i_count_HeatCool = 1 AndAlso i_count_month = 1 Then
                '    Console.WriteLine("col")
                'End If

                Work_msg("요구량계산중..." & CStr(IIf(i_count_HeatCool = 1, "난방", "냉방")) & "(" & i_count_month & "/12)") : Work_plus()
                Dim Mon As String = "M" & Format(i_count_month, "00")   '//저장소의 필드명이 M00 식으로되어있다.

                If i_count_month = 0 Then '//난방이라면? '//최소외기온도 난방기(theta_e_min)
                    CalcO.Pub.theta_e(i_count_month) = TOSG(T, IIf(i_count_HeatCool = 1, CalcO.Pub.theta_e_min, CalcO.Pub.theta_e_max))
                End If



                '//존 루프
                For Each DR존 As DS.tbl_zoneRow In DSET1.tbl_zone.Select("code = '0001'", "code") ' DSET1.tbl_zone.Rows.Count '(3)
                    '//방식이나 공조는 ZZ는 사용안함이다 V=환기  외기부하처리여부의경우 0001 은 예이다.
                    If (DR존.냉난방방식 <> "기능없음" Or (DR존.냉난방공조 <> "기능없음" And (DR존.냉난방공조 <> "환기" Or DR존.외기부하처리여부 = "예"))) And DR존.면적 <> "0" Then
                        Q_vh_b_op = 0

                        '//i_count_OpWe 1은 평일 2는 주말이다ㅡㅡ/
                        '======================================2010.2.21 삭제
                        'For i_count_OpWe = 1 To 2
                        '======================================2010.2.21 삭제
                        Q_V_inf_sink = 0
                        Q_V_inf_source = 0
                        Q_V_mech_sink = 0
                        Q_V_mech_source = 0
                        Q_V_z_sink = 0
                        Q_V_z_source = 0
                        Q_V_win_sink = 0 '//2011.04.04 추가
                        Q_V_win_source = 0 '//2011.04.04 추가

                        Q_T_e_sink = 0
                        Q_T_e_source = 0
                        Q_T_u_sink = 0
                        Q_T_u_source = 0
                        Q_T_z_sink = 0
                        Q_T_z_source = 0
                        Q_T_s_sink = 0
                        Q_T_s_source = 0

                        A_B = TOSG(T, DR존.면적) ' Worksheets("입력존").Range("전용면적1").Offset(0, i_count_zone - 1)
                        V = TOSG(T, DR존.순실체적) ' Worksheets("입력존").Range("체적1").Offset(0, i_count_zone - 1)
                        V_ue = A_u * h_u

                        'Log.AppendLine("               >>  V_ue = A_u * h_u =  // au와 hu 는 선언부에서 1값으로 고정되어있음 v_ue는 항상 1")

                        '//사용한 입력면 (존분류)

                        '//사용한 입력면의 용도프로필 확인

                        '//존의 프로필설정에대한 실제 프로필의 데이터를 가져옴
                        Dim DR프로필() As DS.tbl_profile_odRow = CType(DSET1.tbl_profile_od.Select("code='0'"), DS.tbl_profile_odRow())
                        If DR프로필.GetUpperBound(0) = -1 Then
                            MsgBox("존=" & DR존.code & " 의 프로필=" & DR존.프로필 & " 의 데이터가 존재하지않습니다", MsgBoxStyle.Critical, "확인")
                        Else
                            'Log.AppendLine("               >>  존의 프로필코드(" & DR존.프로필 & ") 에 해당하는 프로필정보를 가져옵니다")
                            t_h_op_d = TOIT(DR프로필(0).운전시간난방) '  Worksheets("설정조건").Range("사용시작시간1").Offset(5, i_count_profile)
                            '===================================================================================2010.2.21 추가
                            If i_count_month = 0 Then
                                d_op_mth = 0
                            Else
                                d_op_mth = CInt(TOSG(T, DR프로필(0)("day" & Format(i_count_month, "00"))))
                            End If
                            '===================================================================================2010.2.21 추가
                            theta_i_h_soll = TOSG(T, DR프로필(0).난방설정온도) ' Worksheets("설정조건").Range("난방설정온도1").Offset(0, i_count_profile)
                            theta_i_c_soll = TOSG(T, DR프로필(0).냉방설정온도) ' Worksheets("설정조건").Range("냉방설정온도1").Offset(0, i_count_profile)
                            Delta_theta_i_NA = TOSG(T, DR프로필(0).야간최저허용온도) ' Worksheets("설정조건").Range("야간최저허용온도1").Offset(0, i_count_profile)
                            Delta_theta_i_WE = TOSG(T, DR프로필(0).주말최저허용온도) ' Worksheets("설정조건").Range("주말최저허용온도1").Offset(0, i_count_profile)
                            '  If i_count_OpWe = 1 Then=============================================================2010.2.21 삭제
                            'Log.AppendLine("               >>  i_count_OpWe = 1 >> Q_I_source_p/Q_I_source_fac/V_mech 계산됨")
                            Q_I_source_p = TOSG(T, DR프로필(0).사람) * A_B '  Worksheets("설정조건").Range("사람1").Offset(0, i_count_profile) * A_B '(2-118)
                            Q_I_source_fac = TOSG(T, DR프로필(0).작업보조기기) * A_B ' Worksheets("설정조건").Range("작업보조기기1").Offset(0, i_count_profile) * A_B '(2-119)
                            V_mech = TOSG(T, DR프로필(0).최소도입외기량) * A_B '  Worksheets("설정조건").Range("최소도입외기량1").Offset(0, i_count_profile) * A_B
                            '    Else   '=============================================================2010.2.21 삭제
                            '    Log.AppendLine("               >>  i_count_OpWe <> 1 >> Q_I_source_p/Q_I_source_fac 는 0입니다")
                            '    Q_I_source_p = 0 '주말 사람없음
                            '    Q_I_source_fac = 0 '주말 기기발열 없음
                            'End If
                            'Q_I_L = 0 'mu_L*Q_I_L_elektr (2-123)추후 추가
                        End If

                        V_mech_outdoor = 1 '외기처리 관련 (ZZ=사용안함 Y=예)
                        If DR존.냉난방방식 <> "기능없음" And DR존.냉난방공조 <> "기능없음" And DR존.외기부하처리여부 = "예" Then
                            If i_count_HeatCool = 1 Then    '//난방이라면
                                If DR존.냉난방공조 <> "난방" AndAlso DR존.냉난방공조 <> "냉난방" Then '난방/냉난방 이 아니면!
                                    If theta_i_h_soll > CalcO.Pub.theta_e(i_count_month) Then    '//난방설정온도 > 외기온도(난방기)
                                        Q_vh_b_op = V_mech * (theta_i_h_soll - CalcO.Pub.theta_e(i_count_month)) * t_h_op_d * c_p_arhoa * d_op_mth '===============d_op(i_count_month) 2010.2.21 수정
                                        V_mech_outdoor = 0 ' 외기처리된 상태인 경우
                                        'Log.AppendLine("               >>  외기처리됨")
                                    End If
                                End If
                            Else
                                If DR존.냉난방공조 <> "냉방" AndAlso DR존.냉난방공조 <> "냉난방" Then '냉방/냉난방이 아닐경우
                                    If theta_i_c_soll < CalcO.Pub.theta_e(i_count_month) Then
                                        Q_vh_b_op = V_mech * (CalcO.Pub.theta_e(i_count_month) - theta_i_c_soll) * t_h_op_d * c_p_arhoa * d_op_mth '===============d_op(i_count_month) 2010.2.21 수정
                                        V_mech_outdoor = 0 ' 외기처리된 상태인 경우
                                        'Log.AppendLine("               >>  외기처리됨")
                                    End If
                                End If
                            End If
                        End If



                        n_50 = TOSG(T, DR존.침기율)    '//tosg 는 빈값일경우 0을 반환
                        n_inf = n_50 * e_wind '(2-59) 환기기기 없는 경우
                        'n_inf = n_50 * e_wind * (1 + f_V_mech * t_V_mech / 24) '(2-60) 급기 및 배기에 대한 환기회수 필요


                        n_win = 0.1 ' 최소외기 도입량은 환기기기 유무와 상관없이 적용되는 것으로 평가하기 때문에 자연환기 기본값만 적용

                        H_T_D = 0
                        H_T_iu = 0
                        H_T_iz = 0


                        Dim Drow프로필 As DS.tbl_profile_odRow = CType(DSET1.tbl_profile_od.Select("code='0'")(0), DS.tbl_profile_odRow)
                        'If DR존.조명부하산출방법 = "계산치" Then
                        'Q_I_L = TOSG(T,T,DR존.조명에너지부하율계산치) * TOSG(T,T,Drow프로필.일일사용시간) * A_B
                        'Else
                        'Q_I_L = TOSG(T,T,DR존.조명에너지부하율입력치) * TOSG(T,T,Drow프로필.일일사용시간) * A_B
                        'End If

                        Q_I_L = TOSG(T, DR존.조명에너지부하율입력치) * TOSG(T, Drow프로필.일일사용시간) * A_B
                        'Q_I_L = TOSG(T,T,IIf(DR존.조명부하산출방법 = "계산치", TOSG(T,T,DR존.조명에너지부하율계산치), TOSG(T,T,DR존.조명에너지부하율입력치))) * TOSG(T,T,Drow프로필.일일사용시간) * A_B
                        Q_w_b = TOSG(T, Drow프로필.일일급탕요구량) * A_B * d_op_mth
                        '===================================2010.2.21 삭제
                        'If DR존.조명부하산출방법 = "CALC" Then    '//OpWe의 2번은 주말인거 같다.
                        '    '//부하율계산치?입력치구분
                        '    Q_w_b = TOSG(T,T,DR존.급탕요구량) * A_B * Calco.Pub.d_op_mth(i_count_month)
                        'Else
                        '    Q_I_L = 0 '주말 소등으로 설정
                        '    Q_w_b = 0 '주말은 급탕사용안함
                        'End If
                        '===================================2010.2.21 삭제

                        'Log.AppendLine("               >>  현재존을 사용하는 각 입력면(" & DSET1.tbl_myoun.Select("존분류='" & DR존.code & "'").GetUpperBound(0) & ")으로부터 데이터를 누적")
                        For Each Dr입력면 As DS.tbl_myounRow In DSET1.tbl_myoun.Select("code <> '0'")


                            Dim CodeYk As String = Dr입력면.열관류율2
                            Dim DrYKRows() As DS.tbl_ykRow = DSET1.tbl_yk.Select("code='" + CodeYk + "'")
                            If DrYKRows.Length = 1 Then

                                '//170601
                                Dim DrYKRow As DS.tbl_ykRow = DrYKRows(0)
                                If DrYKRow.열교방지구조 Then
                                    Delta_U_WB = 0.1
                                Else
                                    Delta_U_WB = 0.15
                                End If

                                'If DR존.열교가산치 = "내단열" Then
                                '    Delta_U_WB = 0.15                            '내단열(=0001) Worksheets("입력존").Range("열교1").Offset(0, i_count_zone - 1) = "내단열" Then
                                'Else
                                '    Delta_U_WB = 0.1                               '내단열(=0001) Worksheets("입력존").Range("열교1").Offset(0, i_count_zone - 1) = "내단열" Then
                                'End If
                            Else
                                Delta_U_WB = 0 '//입력면에 연결된 열관류율 데이터가 없는 경우
                            End If



                            'If DR존.열교가산치 = "내단열" Then
                            '    Delta_U_WB = 0.15                            '내단열(=0001) Worksheets("입력존").Range("열교1").Offset(0, i_count_zone - 1) = "내단열" Then
                            'Else
                            '    Delta_U_WB = 0.1                               '내단열(=0001) Worksheets("입력존").Range("열교1").Offset(0, i_count_zone - 1) = "내단열" Then
                            'End If



                            Select Case Dr입력면.건축부위방식.Trim ' Case Worksheets("입력면").Range("입력부위" & i_sub_wall & "").Offset(0, i_count_wall)
                                Case "외벽" '"외벽"
                                    'Dr입력면.열관류율 >> Calco.Get_열관류율(Dr입력면)
                                    H_T_D = H_T_D + CalcO.Get_열관류율(Dr입력면) * TOSG(T, Dr입력면.건축부위면적) + Delta_U_WB * TOSG(T, Dr입력면.건축부위면적) ' Worksheets("입력면").Range("부위면적" & i_sub_wall & "").Offset(0, i_count_wall) '(2-44) Delta_U_WB는 외단열 0.1, 내단열 0.15
                                Case "외부창", "내부창" '=================2010.03.20 수정
                                    If Dr입력면.블라인드유무 = "유" Then ' '0001(=유) Range("에너지투과율" & i_sub_wall & "").Offset(3, i_count_wall) = "유" Then   
                                        If CalcO.Get_열관류율(Dr입력면) = 0 Then
                                            MsgBox("블라인드 사용이지만 열관류율 수치가 입력되지 않았습니다", MsgBoxStyle.Critical, "확인")
                                        Else
                                            Select Case Dr입력면.블라인드위치.Trim ' Range("에너지투과율" & i_sub_wall & "").Offset(4, i_count_wall)
                                                Case "내부" '"내부"
                                                    U_blinds = TOSG(T, 1 / (1 / CalcO.Get_열관류율(Dr입력면) + 1 / 30))
                                                Case "중간" '중간"
                                                    U_blinds = TOSG(T, 1 / (1 / CalcO.Get_열관류율(Dr입력면) + 1 / 3))
                                                Case "외부" ' '외부"
                                                    U_blinds = TOSG(T, 1 / (1 / CalcO.Get_열관류율(Dr입력면) + 1 / 5 + 1 / 10))
                                            End Select
                                        End If

                                    Else
                                        U_blinds = CalcO.Get_열관류율(Dr입력면) 'Worksheets("입력면").Range("열관류율" & i_sub_wall & "").Offset(0, i_count_wall)
                                    End If

                                    '=========================================================================================================2010.03.20 수정
                                    If Dr입력면.건축부위방식.Trim = "외부창" Then
                                        H_T_D = H_T_D + U_blinds * TOSG(T, Dr입력면.건축부위면적) + Delta_U_WB * TOSG(T, Dr입력면.건축부위면적) '(2-44) Delta_U_WB는 외단열 0.1, 내단열 0.15
                                    Else
                                        H_T_iu = H_T_iu + U_blinds * TOSG(T, Dr입력면.건축부위면적)
                                    End If
                                    '=========================================================================================================2010.03.20 수정
                                Case "내벽" ' "내벽"
                                    H_T_iu = H_T_iu + CalcO.Get_열관류율(Dr입력면) * TOSG(T, Dr입력면.건축부위면적) '(2-48)
                                Case "간벽" '"간벽"
                                    H_T_iz = H_T_iz + CalcO.Get_열관류율(Dr입력면) * TOSG(T, Dr입력면.건축부위면적) '(2-52)
                                Case Else
                                    MsgBox("입력면 건축부위방식이 없습니다", MsgBoxStyle.Critical, "확인")
                                    Return False
                            End Select
                        Next


                        If H_T_D = 0 Then
                            H_V_inf = 0
                        Else
                            H_V_inf = n_inf * V * c_p_arhoa '(2-58) c_p_arhoa는 0.34 V는 존의 체적
                        End If
                        H_V_win = n_win * V * c_p_arhoa '//2011.04.04 추가
                        H_V_mech = V_mech * c_p_arhoa 'n_mech * V * c_p_arhoa '(2-83)
                        H_V_ue = n_ue * V_ue * c_p_arhoa '(2-94)
                        'H_V_z = n_z_d * c_p_arhoa '(2-99)

                        If DR존.열저장능력 = "" Then
                            MsgBox(DR존.설명 & "의 열저장능력이 입력되어 있지않습니다", MsgBoxStyle.Critical, "확인")
                            C_wirk = 0
                        Else
                            C_wirk = TOSG(T, DR존.열저장능력) * A_B ' .Worksheets("입력존").Range("열저장능력1").Offset(0, i_count_zone - 1) * A_B
                        End If

                        If i_count_HeatCool = 1 Then
                            H_V_mech_theta = 0 'H_V_mech * (theta_i_h_soll - theta_V_mech) / 6 '(2-130) 우선 고려안함
                        Else
                            H_V_mech_theta = 0 'H_V_mech * (theta_i_c_soll - theta_V_mech) / 6 '(2-130)
                        End If
                        tau = C_wirk / ((H_T_D + H_T_iu + H_T_iz) + (H_V_inf + H_V_mech + H_V_ue) + H_V_mech_theta) '(2-22, 129) 추후 H_V_z 추가
                        If tau < 48 Then tau = 48

                        a = 1 + tau / 16 '(2-25)

                        a_tb = 1 ' =A_beheizt / A_B 'A_beheizt는 난방면적, A_B는 전용면적

                        t_NA = 24 - t_h_op_d
                        f_tb = TOSG(T, 0.8 * (1 - Math.Exp(-Q_h_max / (A_B * 35))) * a_tb ^ 2) '(2-34)

                        '================================================2010.2.21 삭제
                        'If i_count_OpWe = 1 Then
                        '================================================2010.2.21 삭제
                        Select Case DR존.야간운전방식.Trim '공통1002 Worksheets("입력존").Range("야간가동1").Offset(0, i_count_zone - 1)
                            Case "정상가동" '정상가동"
                                f_NA = 0
                            Case "감소가동" '감소가동"
                                f_NA = TOSG(T, 0.13 * t_NA * Math.Exp(-tau / 250) / 24) '(2-28) 감소가동
                            Case "가동정지" '가동정지"
                                f_NA = TOSG(T, 0.26 * t_NA * Math.Exp(-tau / 250) / 24) '(2-29) 가동정지
                        End Select
                        theta_i_h = theta_i_h_soll - f_NA * (theta_i_h_soll - CalcO.Pub.theta_e(i_count_month))
                        If theta_i_h < theta_i_h_soll - Delta_theta_i_NA * t_NA / 24 Then theta_i_h = theta_i_h_soll - Delta_theta_i_NA * t_NA / 24 '(2-27) 야간감소
                        '================================================2010.2.21 삭제
                        'Else
                        '================================================2010.2.21 삭제
                        Select Case DR존.주말운전방식.Trim ' Worksheets("입력존").Range("주말가동1").Offset(0, i_count_zone - 1)
                            Case "정상가동" '정상가동"
                                f_we = 0
                            Case "감소가동" '감소가동"
                                f_we = TOSG(T, 0.2 * (1 - 0.4 * tau / 250)) '(2-31) 난방감소
                            Case "가동정지" '가동정지"
                                f_we = TOSG(T, 0.3 * (1 - 0.2 * tau / 250)) '(2-32) 난방정지식
                        End Select
                        theta_i_h = theta_i_h_soll - f_we * (theta_i_h_soll - CalcO.Pub.theta_e(i_count_month))
                        If theta_i_h < theta_i_h_soll - Delta_theta_i_WE Then theta_i_h = theta_i_h_soll - Delta_theta_i_WE '(2-30) 주말감소식 변형
                        '================================================2010.2.21 삭제
                        'End If
                        '================================================2010.2.21 삭제

                        If i_count_HeatCool = 1 Then
                            theta_i = theta_i_h
                            theta_z = theta_i_h
                            theta_u = TOSG(T, theta_i_h_soll - 0.7 * (theta_i_h_soll - CalcO.Pub.theta_e(i_count_month)))
                        Else
                            theta_i = theta_i_c_soll '24'에너지 사용량 계산할때 24 사용
                            theta_z = theta_i_c_soll
                            theta_u = TOSG(T, theta_i_c_soll - 0.7 * (theta_i_c_soll - CalcO.Pub.theta_e(i_count_month)))
                        End If

                        'theta_i_h = theta_i_h_soll - f_tb * (theta_i_h_soll - theta_e(i_count_month)) '(2-33) 공간적제한
                        'theta_i_h = theta_i_NA - f_tb * (theta_i_NA - theta_e(i_count_month)) '(2-35) 공간적 시간적 제한

                        'theta_V_mech = theta_e(i_count_month) '(2-90) 공조처리 없는 경우
                        'theta_V_mech = theta_e(i_count_month) + eta_V_mech * (theta_i - theta_e(i_count_month)) '(2-91) 비제어적 열교환 eta_V_mech은 사용자입력
                        theta_V_mech_RLT = TOSG(T, IIf(i_count_HeatCool = 1, 27, 16))

                        theta_V_mech = CalcO.Pub.theta_e(i_count_month) 'theta_V_mech_RLT '(2-92) 공조처리된 급기온도(도표 3-3, 3-4, 7-5) 난방 27, 냉방 16으로 설정

                        '=======================================================================2010.02.04 추가부분 시작
                        Dim 공조기기갯수 As Integer = DSET1.tbl_kongjo.Select("code <> '0'").Length
                        If 공조기기갯수 > 0 Then   '//공조기기가 있는 경우 처리  171225
                            'If DR존.냉난방공조처리시스템 <> "0" Or DR존.냉난방공조 = "환기" Then   '//냉난방공조처리시스템이 선택되어있고 공조가 환기라면
                            '//공조처리스템과 일치하는 공조기기(tbl_kongjo) 를 찾아서 계산한다.

                            Dim 풍량합_총 As Double = 0
                            Dim 가중풍량 As Double = 0
                            Dim 가중효율 As Double = 0
                            For Each drKJ As DS.tbl_kongjoRow In DSET1.tbl_kongjo.Select("code <> '0'")
                                Dim 풍량합 As Double = (TOSG(drKJ.급기풍량) + TOSG(drKJ.배기풍량)) * TOSG(drKJ.대수)  '==========2017.12.28 공조기 대수 추가
                                풍량합_총 += 풍량합
                                Dim 효율 As Double = 0
                                If drKJ.열교환기유형 = "전열교환" Or drKJ.열교환기유형 = "현열교환" Then
                                    If theta_i > CalcO.Pub.theta_e(i_count_month) Then
                                        '//난방
                                        효율 = TOSG(drKJ.열회수율) / 100
                                    Else
                                        '//냉방
                                        효율 = TOSG(drKJ.열회수율_냉방) / 100
                                    End If
                                Else
                                    효율 = 0
                                End If
                                가중풍량 += (풍량합 * 효율)
                            Next
                            가중효율 = 가중풍량 / 풍량합_총


                            ' Dim DR() As DS.tbl_kongjoRow = CType(DSET1.tbl_kongjo.Select("code='" & DR존.냉난방공조처리시스템 & "' and code <> '0'"), DS.tbl_kongjoRow())
                            Phi = 가중효율 ' 0.8
                            theta_V_mech = CalcO.Pub.theta_e(i_count_month) + Phi * (theta_i - CalcO.Pub.theta_e(i_count_month))
                            'End If
                        Else
                            Phi = 0
                        End If

                        '=======================================================================2010.02.04 추가부분 끝

                        If theta_i > CalcO.Pub.theta_e(i_count_month) Then
                            Q_V_inf_sink = H_V_inf * (theta_i - CalcO.Pub.theta_e(i_count_month)) * t_24 '(2-56)
                            Q_V_win_sink = H_V_win * (theta_i - CalcO.Pub.theta_e(i_count_month)) * t_24 '//2011.04.04 추가
                        Else
                            Q_V_inf_source = H_V_inf * (CalcO.Pub.theta_e(i_count_month) - theta_i) * t_24 '(2-
                            Q_V_win_source = H_V_win * (CalcO.Pub.theta_e(i_count_month) - theta_i) * t_24 '//2011.04.04 추가
                        End If



                        '================================================2010.2.21 삭제
                        'If i_count_OpWe = 1 Then
                        '================================================2010.2.21 삭제
                        If theta_i > theta_V_mech Then
                            Q_V_mech_sink = H_V_mech * V_mech_outdoor * (theta_i - theta_V_mech) * t_h_op_d '(2-81)
                        Else
                            Q_V_mech_source = H_V_mech * V_mech_outdoor * (theta_V_mech - theta_i) * t_h_op_d '(2-
                        End If
                        '================================================2010.2.21 삭제
                        'Else
                        'If theta_i > theta_V_mech Then
                        '    Q_V_mech_sink = 0
                        'Else
                        '    Q_V_mech_source = 0
                        'End If
                        'End If
                        '================================================2010.2.21 삭제

                        If theta_i > theta_z Then
                            Q_V_z_sink = H_V_z * (theta_i - theta_z) * t_24 '(2-97)
                        Else
                            Q_V_z_source = H_V_z * (theta_z - theta_i) * t_24 '(2-
                        End If

                        Q_V_sink = Q_V_inf_sink + Q_V_mech_sink + Q_V_z_sink + Q_V_win_sink '(2-13)

                        '//Console.WriteLine(String.Format("mon={0},QV_sink={1},{2},{3},{4}", i_count_month, Q_V_sink, Q_V_inf_sink, Q_V_mech_sink, Q_V_z_sink))

                        Q_V_source = Q_V_inf_source + Q_V_mech_source + Q_V_z_source + Q_V_win_source '(2-

                        If theta_i > CalcO.Pub.theta_e(i_count_month) Then
                            Q_T_e_sink = H_T_D * (theta_i - CalcO.Pub.theta_e(i_count_month)) * t_24 '(2-42)
                        Else
                            Q_T_e_source = H_T_D * (CalcO.Pub.theta_e(i_count_month) - theta_i) * t_24 '(2-
                        End If

                        If theta_i > theta_u Then
                            Q_T_u_sink = H_T_iu * (theta_i - theta_u) * t_24 '(2-46)
                        Else
                            Q_T_u_source = H_T_iu * (theta_u - theta_i) * t_24 '(2-
                        End If

                        If theta_i > theta_z Then
                            Q_T_z_sink = H_T_iz * (theta_i - theta_z) * t_24 '(2-50)
                        Else
                            Q_T_z_source = H_T_iz * (theta_z - theta_i) * t_24 '(2-
                        End If

                        If theta_i > CalcO.Pub.theta_e(i_count_month) Then
                            Q_T_s_sink = H_T_s * (theta_i - CalcO.Pub.theta_e(i_count_month)) * t_24 '(2-53) 추후수정
                        Else
                            Q_T_s_source = H_T_s * (CalcO.Pub.theta_e(i_count_month) - theta_i) * t_24 '(2-53) 추후수정
                        End If

                        Q_T_sink = Q_T_e_sink + Q_T_u_sink + Q_T_z_sink + Q_T_s_sink '(2-12)
                        Q_T_source = Q_T_e_source + Q_T_u_source + Q_T_z_source + Q_T_s_source '(2-12)

                        Q_S_sink = 0
                        Q_S_source = 0

                        For Each DRF입력면 As DS.tbl_myounRow In DSET1.tbl_myoun.Select("code <> '0'")    '//해당존을 사용하는 모든 입력면을 가져옴
                            Select Case DRF입력면.방위.Trim ' Worksheets("입력면").Range("방위" & i_sub_wall & "").Offset(0, i_count_wall)
                                Case "수평" '"수평"
                                    i_count_Is = 0
                                Case "남" '"남"
                                    i_count_Is = 1
                                Case "남동" '"남동"
                                    i_count_Is = 2
                                Case "남서" '"남서"
                                    i_count_Is = 3
                                Case "동" '"동"
                                    i_count_Is = 4
                                Case "서" '"서"
                                    i_count_Is = 5
                                Case "북동" '"북동"
                                    i_count_Is = 6
                                Case "북서" '"북서"
                                    i_count_Is = 7
                                Case "북" '"북"
                                    i_count_Is = 8
                                    '==========================================================2010.03.20 추가
                                Case Else
                                    i_count_Is = -1
                                    '==========================================================2010.03.20 추가
                            End Select

                            If i_count_Is <> -1 Then '=====================================2010.03.20 추가


                                Select Case DRF입력면.건축부위방식.Trim ' Range("입력부위" & i_sub_wall & "").Offset(0, i_count_wall)
                                    Case "외벽" '외벽
                                        If alpa * CalcO.Pub.I_s(i_count_Is, i_count_month) < F_f * h_r * Delta_theta_er Then
                                            Q_S_sink = Q_S_sink + R_se * TOSG(T, DRF입력면.열관류율) * TOSG(T, DRF입력면.건축부위면적) * (F_f * h_r * Delta_theta_er - alpa * CalcO.Pub.I_s(i_count_Is, i_count_month)) * t_24 '(2-111)
                                        Else
                                            'Log.AppendLine("외벽 if else = " & Q_S_source)
                                            Q_S_source = Q_S_source + R_se * TOSG(T, DRF입력면.열관류율) * TOSG(T, DRF입력면.건축부위면적) * (alpa * CalcO.Pub.I_s(i_count_Is, i_count_month) - F_f * h_r * Delta_theta_er) * t_24 '(2-
                                            'Log.AppendLine(">>" & Q_S_source)
                                        End If

                                    Case "외부창", "내부창" '=====================================2010.03.20 수정
                                        If TOSG(T, DRF입력면.수평차양각) < 67.5 Then 'Range("에너지투과율" & i_sub_wall & "").Offset(1, i_count_wall) < 67.5 Then

                                            If DRF입력면.블라인드유무 = "유" Then ' Range("에너지투과율" & i_sub_wall & "").Offset(3, i_count_wall) = "유" Then
                                                Select Case DRF입력면.블라인드빛종류.Trim ' Range("에너지투과율" & i_sub_wall & "").Offset(6, i_count_wall)
                                                    Case "불투과(t=0.0)" '불투과(t=0.0)"
                                                        tau_e_B = 0
                                                        Select Case DRF입력면.블라인드색상.Trim ' Range("에너지투과율" & i_sub_wall & "").Offset(7, i_count_wall)
                                                            Case "흰색" '흰색"
                                                                rho_e_B = 0.7
                                                            Case "밝은색" '밝은색"
                                                                rho_e_B = 0.5
                                                            Case "어두운색" '어두운색"
                                                                rho_e_B = 0.3
                                                            Case "검은색" '검은색"
                                                                rho_e_B = 0.1
                                                            Case Else
                                                                MsgBox("불투과 블라인드색상이 없습니다", MsgBoxStyle.Critical, "확인")
                                                        End Select
                                                    Case "약투과(t=0.2)" '약투과(t=0.2)"
                                                        tau_e_B = 0.2
                                                        Select Case DRF입력면.블라인드색상.Trim 'Range("에너지투과율" & i_sub_wall & "").Offset(7, i_count_wall)
                                                            Case "흰색"
                                                                rho_e_B = 0.6
                                                            Case "밝은색"
                                                                rho_e_B = 0.4
                                                            Case "어두운색"
                                                                rho_e_B = 0.2
                                                            Case "검은색"
                                                                rho_e_B = 0.1
                                                            Case Else
                                                                MsgBox("약투과 블라인드색상이 없습니다", MsgBoxStyle.Critical, "확인")
                                                        End Select
                                                    Case "반투과(t=0.4)" '반투과(t=0.4)"
                                                        tau_e_B = 0.4
                                                        Select Case DRF입력면.블라인드색상.Trim 'Range("에너지투과율" & i_sub_wall & "").Offset(7, i_count_wall)
                                                            Case "흰색"
                                                                rho_e_B = 0.4
                                                            Case "밝은색"
                                                                rho_e_B = 0.3
                                                            Case "어두운색"
                                                                rho_e_B = 0.2
                                                            Case "검은색"
                                                                rho_e_B = 0.1
                                                            Case Else
                                                                MsgBox("반투과 블라인드색상이 없습니다", MsgBoxStyle.Critical, "확인")
                                                        End Select
                                                End Select

                                                If DRF입력면.블라인드각도 = "45도" Then ' Range("에너지투과율" & i_sub_wall & "").Offset(5, i_count_wall) = "45도" Then
                                                    tau_e_B_corr = TOSG(T, 0.65 * tau_e_B + 0.15 * rho_e_B)
                                                    rho_e_B_corr = TOSG(T, rho_e_B * (0.75 + 0.7 * tau_e_B))
                                                    tau_e_B = tau_e_B_corr
                                                    rho_e_B = rho_e_B_corr
                                                End If

                                                g = TOSG(T, DRF입력면.투과율) ' Range("에너지투과율" & i_sub_wall & "").Offset(0, i_count_wall)
                                                If DRF입력면.블라인드위치.Trim <> "" AndAlso TOSG(T, DRF입력면.열관류율) = 0 Then
                                                    MsgBox(DRF입력면.설명 & "의 열관류율이 0입니다", MsgBoxStyle.Critical, "확인")
                                                Else
                                                    Select Case DRF입력면.블라인드위치.Trim 'Range("에너지투과율" & i_sub_wall & "").Offset(4, i_count_wall)
                                                        Case "내부" '내부"
                                                            U_blinds = TOSG(T, 1 / (1 / TOSG(T, DRF입력면.열관류율) + 1 / 30))
                                                            g_blinds = g * (1 - g * rho_e_B - (1 - tau_e_B - rho_e_B) * U_blinds / 30)
                                                        Case "중간" '중간"
                                                            U_blinds = TOSG(T, 1 / (1 / TOSG(T, DRF입력면.열관류율) + 1 / 3))
                                                            g_blinds = g * tau_e_B + g * ((1 - tau_e_B - rho_e_B) + (1 - g) * rho_e_B) * U_blinds / 3
                                                        Case "외부" '외부"
                                                            U_blinds = TOSG(T, 1 / (1 / TOSG(T, DRF입력면.열관류율) + 1 / 5 + 1 / 10))
                                                            g_blinds = tau_e_B * g + (1 - tau_e_B - rho_e_B) * U_blinds / 10 + tau_e_B * (1 - g) * g / 5
                                                        Case Else
                                                            MsgBox("블라인드위치가 없습니다", MsgBoxStyle.Critical, "확인")
                                                    End Select
                                                End If
                                                '===========================================================================2010.02.27 추가
                                                If i_count_month = 0 Then
                                                    If i_count_HeatCool = 1 Then
                                                        g_blinds = (1 - CalcO.Pub.F_sh_gl(i_count_Is - 1, 1 - 1)) * g + CalcO.Pub.F_sh_gl(i_count_Is - 1, 1 - 1) * g_blinds
                                                    Else
                                                        g_blinds = (1 - CalcO.Pub.F_sh_gl(i_count_Is - 1, 7 - 1)) * g + CalcO.Pub.F_sh_gl(i_count_Is - 1, 7 - 1) * g_blinds
                                                    End If
                                                Else
                                                    g_blinds = (1 - CalcO.Pub.F_sh_gl(i_count_Is - 1, i_count_month - 1)) * g + CalcO.Pub.F_sh_gl(i_count_Is - 1, i_count_month - 1) * g_blinds
                                                End If
                                                '===========================================================================2010.02.27 추가

                                            Else
                                                g_blinds = TOSG(T, DRF입력면.투과율) ' Range("에너지투과율" & i_sub_wall & "").Offset(0, i_count_wall)

                                            End If
                                            Q_S_source = TOSG(T, Q_S_source + F__F * TOSG(T, DRF입력면.건축부위면적) * F_S * F_W * F_V * g_blinds * Math.Cos(1.33 * TOSG(T, DRF입력면.수평차양각) * 3.141592 / 180) * (1 - TOSG(T, DRF입력면.수직차양각) / 300) * CalcO.Pub.I_s(i_count_Is, i_count_month) * t_24) '(2-113)
                                        End If
                                End Select
                            End If '==========================================================2010.03.20 추가
                        Next

                        Q_I_source = Q_I_source_p + Q_I_source_fac + Q_I_L '(2-20)
                        Q_source = Q_S_source + Q_T_source + Q_V_source + Q_I_source '(2-16)
                        If i_count_month = 0 AndAlso i_count_HeatCool = 1 Then Q_source = 0 ' 난방 최대부하 계산
                        Q_sink = Q_T_sink + Q_V_sink + Q_S_sink '(2-11) Q_I_sink는 현재 고려안함

                        If Q_sink = 0 Then
                            eta = 0
                        Else
                            gamma = Q_source / Q_sink '(2-21)
                            'If gamma > 1.1 And a > 100 Then
                            If gamma ^ a > 1000000 Then
                                eta = 1 / gamma
                            Else
                                If gamma = 1 Then
                                    eta = a / (a + 1) '(2-24)
                                Else
                                    eta = TOSG(T, (1 - gamma ^ a) / (1 - gamma ^ (a + 1))) '(2-23)
                                End If
                            End If

                            If 1 - (eta * gamma) < 0.01 Then eta = 1 / gamma '(2-137)
                            If (1 - eta) * gamma < 0.05 Then eta = 1 '(2-138)   '//211123
                            'If V_mech >= Q_C_max / (c_p_arhoa * (theta_i - theta_mech)) Then eta = 1 '(2-139) 추후 고민
                        End If

                        If i_count_HeatCool = 1 Then
                            Q_h_b = Q_sink - eta * Q_source '(2-1)
                        Else
                            Q_h_b = (1.0 - eta) * Q_source '(2-7)
                        End If

                        If i_count_month = 3 Then
                            Console.WriteLine(String.Format(":::::: Q_h_b={3},Q_sink={0},eta={1},Q_source={2},i_count_HeatCool={4}", _
                                              Q_sink, eta, Q_source, Q_h_b, i_count_HeatCool))
                        End If

                        If Q_h_b < 0.001 Then Q_h_b = 0 '1W 미만 제거

                        '========================================================2010.2.21 수정 And i_count_OpWe = 1  삭제
                        If i_count_month = 0 Then
                            If t_h_op_d <> 0 Then
                                Q_h_b_op = Q_h_b / t_h_op_d
                            Else
                                Q_h_b_op = Q_h_b
                            End If
                        Else
                            '=================================================2010.2.21 삭제
                            'If i_count_OpWe = 1 Then
                            '=================================================2010.2.21 삭제
                            Q_h_b_op = Q_h_b * d_op_mth

                            If i_count_month = 3 Then
                                Console.WriteLine(String.Format("======== Q_h_b_op={0},Q_h_b={1},d_op_mth={2}", _
                                                  Q_h_b_op, Q_h_b, d_op_mth))
                            End If
                            ' MsgBox(i_count_month & "    " & Q_h_b & "     " & d_op(i_count_month) & "     " & Q_h_b_op)
                            '=================================================2010.2.21 삭제
                            'Else
                            'Q_h_b_we = Q_h_b * Calco.Pub.d_we(i_count_month)
                            '=================================================2010.2.21 삭제
                            Dim DRProf As DS.tbl_profile_odRow = CType(DSET1.tbl_profile_od.Select("code='0'")(0), DS.tbl_profile_odRow)
                            'If DR존.조명부하산출방법 = "계산치" Then 'Worksheets("조명").Range("조명부하산출방법1").Offset(0, i_count_zone - 1) = "계산치" Then
                            'Q_l_b = TOSG(T,T,DR존.조명에너지부하율계산치) * TOSG(T,T,DRProf.일일사용시간) * A_B * d_op_mth
                            'Else
                            'Q_l_b = TOSG(T,T,DR존.조명에너지부하율입력치) * TOSG(T,T,DRProf.일일사용시간) * A_B * d_op_mth
                            'End If
                            Q_l_b = TOSG(T, DR존.조명에너지부하율입력치) * TOSG(T, DRProf.일일사용시간) * A_B * d_op_mth
                            Q_w_b = TOSG(T, DRProf.일일급탕요구량) * A_B * d_op_mth 'Worksheets("입력존").Range("일일급탕요구량1").Offset(0, i_count_zone - 1) * A_B * d_op(i_count_month)
                            '=================================================2010.2.21 삭제
                            'End If
                            '=================================================2010.2.21 삭제
                        End If
                        'Q_h_b_mth = d_op(i_count_month) * (Q_sink_op - eta_op * Q_source_op) + d_we(i_count_month) * (Q_sink_we - eta_we * Q_source_we) '(2-6)
                        '======================================2010.2.21 삭제
                        'Next
                        '======================================2010.2.21 삭제
                    Else

                        A_B = TOSG(T, DR존.면적) ' Worksheets("입력존").Range("전용면적1").Offset(0, i_count_zone - 1)

                        Q_vh_b_op = 0
                        Q_h_b_op = 0
                        'Q_h_b_we = 0

                        '//조명부하산출방법이 계산치(0001)일경우 계산치로 그렇지않을경우 입력치로 계산함
                        Dim DRProf As DS.tbl_profile_odRow = CType(DSET1.tbl_profile_od.Select("code='0'")(0), DS.tbl_profile_odRow)

                        'If DR존.조명부하산출방법 = "계산치" Then
                        'Q_l_b = TOSG(T,T,DR존.조명에너지부하율계산치) * TOSG(T,T,DRProf.일일사용시간) * A_B * d_op_mth
                        'Else
                        'Q_l_b = TOSG(T,T,DR존.조명에너지부하율입력치) * TOSG(T,T,DRProf.일일사용시간) * A_B * d_op_mth
                        'End If

                        Q_l_b = TOSG(T, DR존.조명에너지부하율입력치) * TOSG(T, DRProf.일일사용시간) * A_B * d_op_mth
                        Q_w_b = TOSG(T, DRProf.일일급탕요구량) * A_B * d_op_mth '===============d_op(i_count_month) 2010.2.21 수정
                    End If

                    Console.WriteLine(String.Format("mon={0},Q_l_b ={1},Q_w_b ={2}", i_count_month, Q_l_b, Q_w_b))   '//191120


                    Dim Dr최대부하 As New C최대냉난방부하OD.Citem(DR존.code) ' DSR.최대냉난방부하Row = CType(DSETR1.최대냉난방부하.Select("zone='" & DR존.code & "'")(0), DSR.최대냉난방부하Row)

                    If i_count_HeatCool = 1 Then
                        If i_count_month = 0 Then   '//0번시트는최대 냉난방부하 시트이다.
                            Dr최대부하.최대난방기기부하 = TOSG(T, (Q_h_b_op) / 1000) ''=================== Q_h_b_we 2010.2.21 삭제
                            'Worksheets("" & i_count_month & "").Range("최대난방부하1").Offset(0, i_count_zone - 1) = (Q_h_b_op + Q_h_b_we) / 1000
                            If DR존.야간운전방식 = "가동정지" Then ' Worksheets("입력존").Range("야간가동1").Offset(0, i_count_zone - 1) = "가동정지" Then
                                Dr최대부하.일일난방최대가동시간 = t_h_op_d
                                'Dr최대냉난방부하.난방최대가동시간 = t_h_op_d '  Worksheets("" & i_count_month & "").Range("최대난방부하1").Offset(2, i_count_zone - 1) = t_h_op_d
                            Else
                                Dr최대부하.일일난방최대가동시간 = 24
                                'Dr최대냉난방부하.난방최대가동시간 = 24 ' Worksheets("" & i_count_month & "").Range("최대난방부하1").Offset(2, i_count_zone - 1) = 24
                            End If
                        Else
                            Dim Dr월에너지요구량 As New C에너지분석OD(i_count_month, DR존.code)

                            'Dim Dr월에너지요구량 As DSR.월별에너지분석Row = CType(DSETR1.월별에너지분석.Select("월='" & Mon & "' and 존='" & DR존.code & "'")(0), DSR.월별에너지분석Row)
                            'Dr월에너지요구량 = DSETR1.에너지분석및요구량.Rows(i_count_month - 1) '//각 월별 시트에 값을 셋트(에너지분석및요구량 테이블에 들어있다)
                            If DR존.냉난방방식 = "냉난방" Or DR존.냉난방방식 = "난방" Then    '//냉난방,난방
                                Dr월에너지요구량.난방에너지요구량열 = TOSG(T, (Q_h_b_op) / 1000) '=================== Q_h_b_we 2010.2.21 삭제
                                Dr월에너지요구량.난방에너지요구량공조환기 = TOSG(T, Q_vh_b_op / 1000)
                                '  If i_count_month = 1 Then MsgBox("1")
                                'Log.AppendLine("공조값설정=" & CSng(Q_vh_b_op / 1000))
                            ElseIf DR존.냉난방공조 = "냉난방" Or DR존.냉난방공조 = "난방" Then    '//냉난방,난방
                                Dr월에너지요구량.난방에너지요구량열 = 0
                                Dr월에너지요구량.난방에너지요구량공조환기 = TOSG(T, (Q_h_b_op) / 1000) '=================== Q_h_b_we 2010.2.21 삭제

                                '//오류확인용 존1만 확인한다.
                                'If DR존.code = "0001" Then MsgBox(DR존.냉난방방식 & vbCrLf & DR존.냉난방공조 & vbCrLf & CSng((Q_h_b_op + Q_h_b_we) / 1000), MsgBoxStyle.Information, DR존.code)

                            Else
                                If (DR존.냉난방공조 = "냉방" Or DR존.냉난방공조 = "환기") And DR존.외기부하처리여부 = "예" Then '냉방/환기/"예" Then
                                    Dr월에너지요구량.난방에너지요구량열 = 0
                                    ' If i_count_month = 1 Then MsgBox("3")
                                    Dr월에너지요구량.난방에너지요구량공조환기 = TOSG(T, Q_vh_b_op / 1000)
                                Else
                                    ' If i_count_month = 1 Then MsgBox("4")
                                    Dr월에너지요구량.난방에너지요구량열 = 0
                                    Dr월에너지요구량.난방에너지요구량공조환기 = 0
                                End If
                            End If

                            Dr월에너지요구량.조명에너지요구량 = Q_l_b / 1000

                            ' Log.AppendLine("i_count_HeatCool=" & i_count_HeatCool & "i_count_month=" & i_count_month & ",급탕에너지요구량=" & CSng(Q_w_b / 1000))
                            Dr월에너지요구량.급탕에너지요구량 = Q_w_b / 1000
                            ' MsgBox("난방요구량" & Dr월에너지요구량.난방요구량공조, MsgBoxStyle.Information, i_count_month & "CZONE")
                        End If

                    Else
                        If i_count_month = 0 Then
                            Dr최대부하.최대냉방기기부하 = TOSG(T, (Q_h_b_op) / 1000) '=================== Q_h_b_we 2010.2.21 삭제
                            'Worksheets("" & i_count_month & "").Range("최대냉방부하1").Offset(0, i_count_zone - 1) = (Q_h_b_op + Q_h_b_we) / 1000
                            If DR존.야간운전방식 = "가동정지" Then  '//가동정지
                                Dr최대부하.일일냉방최대가동시간 = t_h_op_d ' Worksheets("" & i_count_month & "").Range("최대냉방부하1").Offset(2, i_count_zone - 1) = t_h_op_d
                            Else
                                Dr최대부하.일일냉방최대가동시간 = 24 'Worksheets("" & i_count_month & "").Range("최대냉방부하1").Offset(2, i_count_zone - 1) = 24
                            End If

                        Else
                            Dim Dr월에너지요구량 As DSR.월별에너지분석Row = CType(DSETR1.월별에너지분석.Select("월='" & Mon & "' and 존='" & DR존.code & "'")(0), DSR.월별에너지분석Row)

                            If DR존.프로필 = "5000" Then '주거부문일 경우 냉방 가동시간 3.8시간 250625
                                Q_h_b_op = TOSG(Q_h_b_op) * 3.8 / 24
                                Q_vh_b_op = TOSG(Q_vh_b_op) * 3.8 / 24
                            End If

                            If DR존.냉난방방식 = "냉난방" Or DR존.냉난방방식 = "냉방" Then    '//냉난방/냉방
                                Dr월에너지요구량.냉방요구량열 = TOSG(T, (Q_h_b_op) / 1000) '=================== Q_h_b_we 2010.2.21 삭제
                                Dr월에너지요구량.냉방요구량공조 = TOSG(T, Q_vh_b_op / 1000)
                            Else
                                If DR존.냉난방공조 = "냉난방" Or DR존.냉난방공조 = "냉방" Then    '//냉난방/냉방
                                    Dr월에너지요구량.냉방요구량열 = 0
                                    Dr월에너지요구량.냉방요구량공조 = TOSG(T, (Q_h_b_op) / 1000) '=================== Q_h_b_we 2010.2.21 삭제
                                    Console.WriteLine(String.Format("+++++++ 냉방요구량공조={0},Q_h_b_op={1}", Dr월에너지요구량.냉방요구량공조, Q_h_b_op))
                                Else
                                    If (DR존.냉난방공조 = "난방" Or DR존.냉난방공조 = "환기") And DR존.냉난방공조 = "예" Then '//난방/환기/예
                                        Dr월에너지요구량.냉방요구량열 = 0
                                        Dr월에너지요구량.냉방요구량공조 = TOSG(T, Q_vh_b_op / 1000)
                                    Else
                                        Dr월에너지요구량.냉방요구량열 = 0
                                        Dr월에너지요구량.냉방요구량공조 = 0
                                    End If

                                End If
                            End If

                            'MsgBox("냉방요구량" & Dr월에너지요구량.냉방요구량공조, MsgBoxStyle.Information, i_count_month & "CZONE")
                    End If
                    End If
                Next '(3)

                Dim CAir As New CAirHandlingOD
                If i_count_month = 0 Then
                    CAir.Calc_Max_Power(i_count_HeatCool) '//AIR HANDLING 모듈 에 존재함..역시 내용이 크다 ㅠ_ㅠ
                End If

            Next '(2)
        Next '(1)

        Work_msg("에너지분석데이터 집계중...")
        Frm_Work.ProgressBar1.Maximum = 13
        Frm_Work.ProgressBar1.Minimum = 0
        Frm_Work.ProgressBar1.Value = 0
        Frm_Work.TopMost = True
        Frm_Work.Show()
        '//요구량 계산이 완료되었으므로 Summary 계산한다. (이것이 없어서 소요량이 오류났음)
        For Each Dr As C에너지분석OD In Result1o.E분석 '//12달치의 각존들의 데이터를 가지고 합계내역을 추린다.
            Work_plus()
            Dr.Calc()
        Next
        DSETR1.AcceptChanges()
        Return True
    End Function


    ''' <summary>
    ''' (main2)소요량 계산
    ''' </summary>
    ''' <remarks></remarks>
    Public Function Calc_Useage() As Boolean
        'Dim d_op_month As Integer
        Dim theta_e_month As Decimal
        Dim d_full_mth As Integer

        Work_msg("소요량계산중...")
        Frm_Work.ProgressBar1.Maximum = 12
        Frm_Work.ProgressBar1.Minimum = 0
        Frm_Work.ProgressBar1.Value = 0
        Frm_Work.TopMost = True
        Frm_Work.Show()

        For i_count_month As Short = 1 To 12 '//월별로 난/냉방 소요량을 계산
            Work_msg("소요량계산중..." & "(" & i_count_month & "/12)") : Work_plus()
            '===============================================2010.2.21 수정
            'd_op_month = d_op(i_count_month)
            d_full_mth = CInt(Calco.Pub.d_full(i_count_month))
            '===============================================2010.2.21 수정
            theta_e_month = Calco.Pub.theta_e(i_count_month)

            '///난방
            Calco.Heating.Calc_WaterheatingSystem(i_count_month, d_full_mth, theta_e_month)  '//급탕소요량 '========2010.2.21 수정
            Calco.AirHandling.Calc_AirHandling(i_count_month, d_full_mth) '//환기? '========2010.2.21 수정
            If Result1o.E분석(i_count_month).난방에너지요구량 <> 0 Then Calco.Heating.Calc_HeatingSystem(i_count_month, theta_e_month) '//난방소요량

            '///냉방
            If Result1o.E분석(i_count_month).냉방에너지요구량 <> 0 Then Calco.Cooling.Calc_CoolingSystem(i_count_month, theta_e_month)
        Next
        DSETR1.AcceptChanges()
        Return True
    End Function



End Class