新安江VB代码.doc

Option Explicit Private Sub Command1_Click() Dim P(366) As Single, E(366) As Single, PE(366) As Single Dim P1(366) As Single, P2(366) As Single, P3(366) As Single, P4(366) As Single Dim E0(366) As Single, Ep(366) As Single Dim EU(366) As Single, EL(366) As Single, ED(366) As Single Dim WU(366) As Single, WL(366) As Single, WD(366) As Single, W(366) As Single Dim Q(366) As Single, R(366) As Single, sumR As Single, sumQ As Single Dim i As Integer, Kc As Single, WMM As Single, a(366) As Single Const WUM = 20, WLM = 60, WDM = 60 Const C = 0.16 Const b = 0.3 Const WM = 140 Const FC = 24 Kc = Val(Text1.Text) If Option1.Value = True Then Open 89data.txt For Input As #1 打开文件 For i = 1 To 365 Input #1, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i) 引用89数据 P(i) = 0.33 * P1(i) + 0.14 * P2(i) + 0.33 * P3(i) + 0.2 * P4(i) sumQ = sumQ + Q(i) * 24 * 3.6 / 553 Text2.Text = sumQ Next i Close #1 End If If Option2.Value = True Then Open 90data.txt For Input As #2 打开文件 For i = 1 To 365 Input #2, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i) 引用90数据 P(i) = 0.33 * P1(i) + 0.14 * P2(i) + 0.33 * P3(i) + 0.2 * P4(i) sumQ = sumQ + Q(i) * 24 * 3.6 / 553 Text2.Text = sumQ Next i Close #2 End If WMM = WM * (1 + b) WU(1) = 10: WL(1) = 40: WD(1) = 60 W(1) = WU(1) + WL(1) + WD(1) a(1) = WMM * (1 - (1 - W(1) / WM) ^ (1 / (1 + b))) 三层蒸散发计算 For i = 1 To 365 Ep(i) = E0(i) * Kc Next i For i = 1 To 365 If WU(i) + P(i) = Ep(i) Then EU(i) = Ep(i) EL(i) = 0 ED(i) = 0 Else If WL(i) = C * WLM Then EU(i) = WU(i) + P(i) EL(i) = (Ep(i) - EU(i)) * WL(i) / WLM ED(i) = 0 Else If WL(i) = C * (Ep(i) - EU(i)) Then EU(i) = WU(i) + P(i) EL(i) = C * (Ep(i) - EU(i)) ED(i) = 0 Else EU(i) = WU(i) + P(i) EL(i) = WL(i) ED(i) = C * (Ep(i) - EU(i)) - EL(i) End If End If End If E(i) = EU(i) + EL(i) + ED(i) PE(i) = P(i) - E(i) 产流量计算 If PE(i) 0 Then 产流 If a(i) + PE(i) WMM Then 蓄满前 R(i) = PE(i) + W(i) - WM + WM * (1 - (PE(i) + a(i)) / WMM) ^ (b + 1) a(i + 1) = PE(i) + a(i) W(i + 1) = W(i) + PE(i) - R(i) Else 蓄满 R(i) = PE(i) + W(i) - WM a(i +