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. OK, I Understand
 
Top