• API
• FAQ
• Tools
• Archive
SHARE
TWEET

Neraca air MSP GFM

a guest Dec 13th, 2018 154 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. Public Sub neraca_air()
2.         'Intersepsi Tajuk Tanaman :
3.         If lai(i) < 3 Then
4.             ic(i) = 0.4233 * lai(i)
5.         Else
6.             ic(i) = 1.27
7.         End If
8.
9.         If ch(i) < ic(i) Then ic(i) = ch(i)
10.
11.         'Curah Hujan Effektif/netto
12.         ch_ef(i) = ch(i) - ic(i)
13.
14.         'Run off fungsi kemiringan dan CHeff:
15.         ro(i) = Math.Sin(3.14 * miring / 180) * ch_ef(i)
16.
17.         'Evapotranspirasi (Penman):
18.         f1 = 0.64 * (1.054 * angin(i))
19.         del = 208.84
20.         ETP(i) = (del * 0.5 * rad(i) + f1 * (100 - rh(i)) / 100 * 2000) / (del * 66.1)
21.
22.         'Evaporasi (Ea) dan Transpirasi (Ta):
23.         Em(i) = ETP(i) * Math.Exp(-0.5 * lai(i))
24.         Tm(i) = ETP(i) - Em(i)
25.
26.         If KAT(i - 1) < KL Then
27.             If KAT(i - 1) > TLP Then Ta(i) = Tm(i) * (KAT(i - 1) - TLP) / (KL - TLP) Else Ta(i) = 0
28.         Else
29.             Ta(i) = Tm(i)
30.         End If
31.         If KAT(i - 1) < (TLP / 2) Then Ea(i) = 0 Else Ea(i) = Em(i)
32.
33.         'Kadar Air Tanah (mm) :
34.         KAT(i) = KAT(i - 1) + ch_ef(i) - ro(i) - Ea(i) - Ta(i)
35.
36.         'Perkolasi :
37.         If KAT(i) > KL Then
38.             Pc(i) = KAT(i) - KL
39.             KAT(i) = KL
40.         Else
41.             Pc(i) = 0
42.         End If
43.         If KAT(i) < 0 Then KAT(i) = 0
44.
45.         'water defisit factor
46.         wdf(i) = Ta(i) / Tm(i)
47.
48.     End Sub
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top