【正文】 (n) Eu) Then Eu = wu + p(n) El = wl Ed = C * (Ep(n) Eu) El End If et = Eu + El + Ed 39。三層蒸發(fā)模式計算蒸發(fā) pe = p(n) et If pe = 0 Then QIM(n) = 0 Else QIM(n) = pe * IM * 553 / ( * 3) End If If pe = 0 Then r(n) = 0 ElseIf pe 0 And (pe + a) = WMM Then r(n) = pe wm * (1 a / WMM) ^ (B + 1) + wm * (1 (pe + a) / WMM) ^ (B + 1) ElseIf pe 0 And (pe + a) WMM Then r(n) = pe wm + W End If If r(n) 0 Then If pe = fc Then rg(n) = fc * r(n) / pe rd(n) = (pe fc) * r(n) / pe ElseIf pe = 0 And pe fc Then rg(n) = pe * r(n) / pe rd(n) = 0 End If Else rg(n) = 0 rd(n) = 0 End If 39。產(chǎn)流計算以及水源劃分 If p(n) 0 And wu = 0 Then wu = 0 ElseIf p(n) = 0 And wu = 0 Then wu = 0 ElseIf p(n) = 0 And wu 0 Then wu = wu Eu ElseIf p(n) 0 And p(n) Eu And wu 0 Then wu = wu (Eu p(n)) ElseIf p(n) 0 And p(n) Eu And wu 0 Then wu = wu End If If El 0 And wl = 0 Then wl = 0 ElseIf wl 0 Then wl = wl El End If If Ed 0 And wd = 0 Then wd = 0 ElseIf wd 0 Then wd = wd Ed End If If pe 0 And pe r(n) = wum wu Then wu = pe r(n) + wu wl = wl wd = wd ElseIf pe 0 And pe r(n) wum wu And (pe r(n) (wum wu) = (wlm wl)) Then wl = wl + pe r(n) (wum wu) wu = wum wd = wd ElseIf pe 0 And pe r(n) wum wu And (pe r(n) (wum wu) (wlm wl)) And ((pe r(n) (wum wu) (wlm wl)) = (wdm wd)) Then wd = wd + pe r(n) (wum wu) (wlm wl) wu = wum wl = wlm ElseIf pe 0 And pe r(n) wum wu And (pe r(n) (wum wu) (wlm wl)) And ((pe r(n) (wum wu) (wlm wl)) (wdm wd)) Then wu = wum wl = wlm wd = wdm End If 39。下一次初始包氣帶含水量變化計算 n = n + 1 Loop For i = 1 To 28 QD(i) = 0 For n = 1 To i QD(i) = QD(i) + rd(n) * UH(i + 1 n) / 10 39。直接徑流匯流計算 Next n Next i For i = 1 To 28 QD(i) = QD(i) * (1 IM) + QIM(i) Next i QG(1) = CG * QG0 + (1 CG) * rg(1) * 553 / ( * 3) For i = 2 To 28 QG(i) = CG * QG(i 1) + (1 CG) * rg(i) * 553 / ( * 3) 39。地下徑流匯流計算 Next i For i = 1 To 28 Q(i) = QG(i) + QD(i) 39?？倧搅鲄R流計算 Next i For i = 1 To 28 Call fact(rd(i)) Call fact(rg(i)) Call fact(r(i)) Call fact(QD(i)) Call fact(QG(i)) Call fact(Q(i)) 第 amp。 i amp。 時段 amp。 P= amp。 p(i) amp。 RD= amp。 rd(i) amp。 RG= amp。 rg(i) amp。 R= amp。 r(i) amp。 QD= amp。 QD(i) amp。 QG= amp。 QG(i) amp。 Q= amp。 Q(i) Next iWith (0) For i = 1 To 28 .Add Q(i), i, vbBlue Next iEnd WithWith (1) For i = 1 To 28 .Add QD(i), i, vbRed Next iEnd WithWith (2) For i = 1 To 28 .Add QG(i), i, vbBlue Next iEnd With End SubPrivate Sub Command2_Click()End SubPrivate Sub Command3_Click() End SubPrivate Function fact(a As Single)Dim d As Single d = CInt(a * 10) / 10 a = dEnd Function八、心得體會水文預(yù)報在防洪、抗旱、水資源開發(fā)利用、國民經(jīng)濟建設(shè)和國防等領(lǐng)域都有廣泛的應(yīng)用，經(jīng)濟效益巨大，應(yīng)用單位眾多。水文預(yù)報是與國家社會經(jīng)濟發(fā)展與安全生活生產(chǎn)等密切相關(guān)的基礎(chǔ)業(yè)務(wù)，它看似單一但卻與我么的生產(chǎn)生活息息相關(guān)，十分重要。通過本次課程設(shè)計，有效的溫習(xí)了產(chǎn)、匯流和蓄滿產(chǎn)流模型（三層模型）相關(guān)的知識，經(jīng)過相應(yīng)的程序的編寫和計算使自己對產(chǎn)、匯流的過程和蓄滿產(chǎn)流模型的機理更加的熟練。通過對程序?qū)c的率定和檢驗，結(jié)果顯示各年的計算徑流和實際徑流之間的相對誤差都較大。這與率定參數(shù)時資料序列較短都有關(guān)，另一方面也顯示出了蓄滿產(chǎn)流模型本身的不確定，以及穩(wěn)定性和適應(yīng)性差得缺點。當然流域氣候和地質(zhì)的時空不均勻性，由于流域本身的地形、土壤、植被、氣候等條件各不相同，且同一流域不同地區(qū)、不同時間各要素也不并不相同，都處于一個時變的狀態(tài)，這是導(dǎo)致水文模型的計算結(jié)果誤差較大的本質(zhì)原因。雖然當前的水文模型仍然存在很多的弊端和不足，但隨著計算機科學(xué)和信息技術(shù)的發(fā)展和在水文與水資源中的應(yīng)用不斷地推廣，使得流域水文模型的研究在近年取得了迅速的發(fā)展，尤其是3s技術(shù)在水文模型中的應(yīng)用，分布式水文模型的開發(fā)等，相信今后的水文模型將會取得突破性的進展，將在解決生產(chǎn)實際問題中發(fā)揮著更加重要的作用。 最后要感謝這次課設(shè)當中給予我?guī)椭睦蠋熀屯瑢W(xué)表示感謝！24