LAMP
Lampiran 1. Diaram AIir Model Pertumbuha~~ dan Perltembal~gz~n Tanalnan Ke11te11g
s > 0.8
pB = 0.161 i-(-0.1286 s Fase) pA = 0.1028 t (-0.095 s I7ase) pD = 0.161 + (-0.3278 x Fase) pU = I
sla = sla + Rsla s s p dlLD = sln s d W D ILD=dll.D + ILD IfII,D < O T h e n I L D = O
1 Gdmp = loc s Sin1 s 1 0 4
- (pB
+
pR +I'D)
Lonlpirnll 2. Di;lgrnm Alir Model Prediksi I < e ~ a d ~ aPer~ynli~l n llarrd! I);IIIII I < ~ e t ; ~ o g
0 Start
WlJi ILDi
ILD = ILD + dlLD
0 Sclesai
= WU = ILD
(11-I)
- I'L13
21
~
Lampiran 5. Modul Public Pada Dua Model Y m g Dibuat Public rad(1000) As Single 'radiasi Public hst(1000) A s Single 'hari setelah tanam Public ch(1000) As Single Public Tmin(1000) As Single Public Tmax(l000) As Single Public tr(1000) As Single 'h.ata2 Public Td(1000) As Single 'titik e n ~ b u n Public RH(1000) As Single Public tac(l000) A s Single 'suhu kolom a dikurangi kolom c ) Single ' fase Public ~ ( 1 0 0 0 As Public Dlai(1000) As Single 'Peruhahan LA1 Public Lail(1000) As Single Public Lai(1000) As Single 'Leaf Area Indeks Public Laip(1000) As Single ' lai setelah terkcna penyakit Public sla(1000) As Single'spesific leaf area Public Dflai(lOO0) A s Single ' pcrubaban LA1 akibat penyakit Public Hin(1000) As Single ' hari dcngan kemungkinan infeksi terjadi paling tinggi Public Hinl(l000) As Single Public wt(10000) As Single 'berat total potensial Public wta(1000) As Single 'berat total aktual Public stom(l0000) As Single 'stomata Public Bp(100000) As Single 'radiasi terinterseysi Public GDMp(10000) As Single 'potential gross dry mutter production Public GDMa(10000) As Single 'actual gross dry matter production (kglhaid) Public QlO(10000) A s Single Public Wdaun(1000) As Single 'berat daun Public Wbat(1000) As Single ' berat batang Public Wak(1000) As Single ' berat akar Public Wurnbi(1090) As Single 'berat ulnbi Public dl(1000) As Single Public ds(1000) As Single Public dr(1000) As Single Public dt(1000) As Single 'perubahan niasiiig niasing organ, daiin,batang;aliar,umbi Public rrnd(10000) As Single Poblic rmb(10000) As Single Public rma(10000) As Single Public r1nt(10000) As Single Public HU(1000) As Single 'Neat Unit Public pddcol 'untuk pengolalian gratik Public pddrow 'ontuk pengolalian gratik Public tbk(500) 'Suiiu bola kering (input) Public radiasi(500) 'radiasi (input) Public pdaun(1000) As Singlc ' proporsi bioni:tssa itntuk daun Public pbatang(1000) As Single ' proporsi biomassa untuk batang Public pakar(1000) As Single 'propoprsi bioniassa untuk akar Public pumbi(1000) As Single 'proporsi biomnssa untuk umbi I'oblic plb(1000) As Single 'pcrtalnbalian ILI:IS bercak
Public dplb(l000) A s Single ' perububahan lai karena luas bercak Public dpdaun(1000) As Single Public wdaunp(1000) As Single ' berat organ setelah terkcna setwngan penyakit Public wbatap(1000) As Single Public wakap(1000) As Single Public wumbip(1000) As Single Public tanggal 'tanggal tanam Public tanam 'tanggal mulai menanarn Public panen ' panen Public tunas Public bentukun?bi Public isiumbi Public matangumbi Public panonumbi Public syarattotal(1000) As Single Public HUl(1000) A s Single Public dflail(1000) As Single Public Ihari As lnteger Public HTanam As Integer Public Sub Inisialisasi() For g = 1 To 1000 tr(g) = 0 rad(g) = 0 tac(g) = 0 HU(g) = 0 syarattotal(g) = 0 =0 Dlai(g) = 0 Lai(g) = 0 wta(g) = 0 Wdaun(g) = 0 Wbat(g) = 0 Wak(g) = 0 Wumbi(g) = 0 sla(g) = 0 Laip(g) = 0 Next g End Sub
Lampiran 6 . Model Silnulasi P c r t ~ ~ m b u Dnn l ~ a Perke~nballgnr~ ~~ Tntlnrnae I<e~ttnsg Open txtfileinput.Text For Input As #I Open txtfileOutput.Text For Output As #2 a=O While Not EOF(1) a=a+ l Input #I, &(a), radiasi(a) Wend Close #I 'menghitung Parameter For b = 1bltanam.Caption To 500 'Labelmatang.Caption ' 1 To a rad(b) = radiasi(b) '* 0.042 tac(ij = tr(b) HUl(b) = HUl(b - 1) + (tac(a + I ) + tac(b)) Next b
iO
s(c) = s(c - I ) + (tac(c) 1 1000) ' s adalah tingkat perkembangan fenologi tanaman HU(c) = s(c) * HUl(c) Next c Ford= 1 T o a If 0.16 >= s(d) Then pbatang(d) = 0.121 * Exp(3.6797 * s(d)) pakar(d) = 0.2937 * Exp(-3.2141 * s(d)) pdaun(d) = l - (pakar(d) + pbatang(d)) pnmbi(d) = 0 End If If 0.16 < s(d) And s(d) <= 0.33 Then pbatang(d) = 0. IS43 * Exp(-1.1874 * s(d)) pakar(d) = 0.2187 * Exp(-2.6511 's(d)) pdaun(d) = 0.7159 * Expi-2.662 * s(d)) pumbi(d) = l - (pbatang(d) + pakar(d) + pdaun(d)) Elld If If 0.33 < s(d) And s(d) <= 0.44 Then pbatang(d) = 0.0053 + (0.2808 * s(d)) pakar(d) = 0.0518 + (0.0572 * s(d)) pdaun(d) = 0.5392 + (-0.7007 * s(d)) pumbi(d) = l - (pbatang(d) + pakar(d) + pdaun(d)) End If
If 0.44 < ~ ( d And ) s(d) <= 0.8 Then pbatang(d) = 0.2572 + (-0.2441 * s(d)) pakar(d) = 0.149 + (-0.1503 * s(d)) pdaon(d) = 0.3404 + (-0.2689 * s(d)) pumbi(d) = I - (pbatang(d) + pakar(d) + pdauti(d)) End If
purnbi(d) = 1 - (pbatang(d) End If Next d
+ pakar(d) + pdaun(d))
Fore= 1Toa Bp(e) = (1 - Exp(-0.35 * 0.0323)) * rad(e) 'produksi biornassa potensial (13p) rr=140* 12*(3600/10A6) If Bp(e) > rr Then Bp(e) = rr GDMp(e) = 2.76 'Bp(e) * (10 A 4) GDMa(e) = (I - 0.13) * GDMp(e) QlO(e) = 2 " (itr(e) - 20) 1 10)
rmdie) = kni * QlO(e) * Wdaun(e) rnlbie) = Imi * QlO(e) * Wbat(e) rma(e) = km * QIO(e) * Wak(e) rmt(e) = km * QlO(e) * Wurnbi(e)
sla(e) = -0.047 * s(e) + 0.0649 Dlai(e) = (sla(e) * dl(e)) Lai l(e) = (Lail(e - I) + Dlai(e - I)) Lai(e) = Lai i(e) 1 1000 Next e
Lainpiran 7. Model Prediloi I
0 And Td(b) > Tmin(d) Then syaratl(b) = 1 Else syaratl(b)= 0 'syarat suhu minimum If Tmin(b) >= 10 Then syarat2(b) = 0.5 Else syaratZ(b) = 0 'syarat suhu maks If Tmax(b) >= 24 And Tmax(b) <= 26 Then syarat3(b) = 1 Else syarat3(b) = 0 'syarat kelembaban If RH(b) > 90 Then syarat4(b) = 1.5 Else syarat4(b) = 0 'syarat serang If hst(b) > s(b) And s(b) <= 0.16 Then syaretS(b) = 1 Else syarat5(b) = 0 'syarat total syarattotal(b) = ((syaratl(b) Next b
+ syaratZ(b) + syarat3(b) + syarat4(b) + syarat5(b)) 15) * 100
s(l) = tac(l) 1 1000
Farc=2Toa s(c) = s(c - I ) + (tac(c) 1 1000) ' s adalah tingkat perkcmbangan fenolog~1:inatnall HU(C) S(C) * I-IUl(c) Next c
-
F w d = 1Toa If s(d) < 0.16 Then syarattotal(d) = 0 If 0.16 r~ ( d Then )
pbatang(d) = 0.121 * Exp(3.6797 * s(d)) pakar(d) = 0.2937 * Exp(-3.2141 * s(d)) pdaun(d) = l (pakar(d) + pbatang(d)) pumbi(d) = 0 End If
-
If 0.16 < s(d) And s(d) <= 0.33 Then pbatang(d) = 0.1843 * Exp(-1.1874 * s(d)) pakar(d) = 0.2187 * Exp(-2.6511 * s(d)) pdaun(d) = 0.7159 * Exp(-2.662 * s(d)) pumbi(d) = 1 - (pbatang(d) + pakar(d) + pdaun(d)) End If
If 0.33 < s(d) And s(d) <= 0.44 Then pbatang(d) = 0.0053 + (0.2808 * s(d)) pakar(d) = 0.05 18 + (0.0572 * s(d)) pdaun(d) = 0.5392 + (-0.7007 * s(d)) pumbi(d) = 1 (pbatang(d) + pakar(d) + pdaun(d)) End If
-
If 0.44 < s(d) And s(d) <= 0.8 Then pbatang(d) = 0.2572 + (-0.2441 * s(d)) pakar(d) = 0.149 + (-0.1503 * s(d)) pdaun(d) = 0.3404 + (-0.2689 * s(d)) pumbi(d) = l - (pbatang(d) + pakar(d) i- pdaun(d)) End If If s(d) > 0.8 Then pbatang(d) = 0.161 + (-0.1286 * s(d)) pakar(d) = 0.1028 + (-0.095 * s(d)) pdaun(d) = 0.161 + (-0.3278 * s(d)) pumbi(d) = 1 - (pbatang(d) + pakar(d) + pdaun(d)) End If Next d Fore= 1 T o a Bp(e) = (I Exp(-0.35 * 0.0323)) * rad(e) 'produksi biomassa potcnsial (Bp) rr=140*12*(36O011OA6) If Bp(e) > rr Then Bp(e) = rr GDMp(e) = 2.76 * Bp(e) * (10 " 4) GDMa(e) = (I - 0.13) * GDMp(e) QI 0(e) = 2 " ((tr(e) - 20) 1 10)
-
If syarattotal(e) >= SO Then X = X + I
If syorattotal(e) >= SO And X = 1 Then Hin(e) = e If syarattotal(e) >= SO And X = 1 Then Txtinf.Text = Hin(e)
plb(o) = (2 * phi) * (e * (4 * I0 " -3)) " 2 If syarattotal(e) >= 60 Then plbl(e) = plbl(e - 1) + plb(e - 1) Else pibl(e) = plbl(e - 1) If Lai(e) > 0 Then dplb(e) = (plb(e) / Lai(e)) * 100 If syarattotal(e) >= 80 Then Laip(e) = Lai(e) - plbl(e) Else Laip(e) = Laip(e - I) + (Dlai(e - I) 1 1000) If Laip(e) <= 0 Then Laip(e) = 0 If Laip(e I) <= 0 And s(e) > 0.44 Then Laip(e) = 0
-
wdaunp(e) = Wdaun(e) wnmbip(e) = Wumbi(e) Ifsyarattotal(e) >= 80 Then wdaunp(e) = Wdaun(e) - ylbl(e) Else wdaun(e) (dl(e - 1) 1 10000) If wdaunp(e) <= 0 Then wdaunp(e) = 0 If Laip(e) = 0 And s(e) >= 0.44 Then wumbip(e) = wumbip(e - 1) Dflai(e) = 0.012 * (Val(TxtinEText) " l.lG22) Next e
= wdaunp(c -
1) +