MINI-MODELOS DE ECOSSISTEMAS 
SIMULAÇÕES EM COMPUTADOR

Parte 1. Modelos básicos (quinto conjunto)

(j)  Banhado ou marisma (zona com pradaria resistente a água salgada e formação de turfa)
(k) Alagadiço (terras alagada com água doce, com pastos e formação de turfa)
(l)  Lagoa natural de criação de camarão com rios, mangues e marés

Principais tipos de Fontes e formas de Captura e Consumo de Energia
 
20. Banhado ou Marisma (pradaria resistente a água salgada, os nutrientes geram turfa)


Programa saltmrsh

Exemplo de programação em Basic: saltmarsh.bas
3 REM IBM
4 REM SALTMRSH, Salt marsh & input rhythm
5 SCREEN 1, 0:  COLOR 0, 0
12 LINE (0, 0)-(319, 180), 3, B
15 LINE (0, 50)-(319, 50)
20 LINE (0, 30)-(319, 30)
23 LINE (0, 70)-(319, 70)
24 LINE (0, 90)-(319, 90)
25 I0 = 200
30 M = 7
35 N = 6.5
37 E0 = .2
38 J0 = .002
40 Q1 = .2
50 Q2 = 10
55 Q3 = 500
60 J = .02
70 E = 1
75 P1 = .05
77 P2 = 30
80 P3 = 30
100 K0 = 20
110 K1 = .02
120 K2 = .0004
130 K3 = .0000004
 
140 K4 = .05
150 K5 = .00002
160 K6 = .00004
161 K7 = .0036
162 K8 = 6.000001E-04
163 K9 = 3.2E-09

165 I = 2000 + (1 + SIN(T / 58)) * 2000
166 J = .02 * (1 + COS(T / 58))
167 E = 2 + SIN(T / 27)
170 R = I / (1 + K0 * Q1)
180 D1 = J - K5 * Q1 * R + K3 * Q2 + K9 * Q3 - K4 * Q1 * E
190 D2 = K1 * Q1 * R - K7 * Q2 - K2 * Q2
195 D3 = K7 * Q2 - K8 * Q3 - K6 * Q3 * E
200 Q1 = Q1 + M * D1
205 IF Q1 < .01 THEN Q1 = .01
210 Q2 = Q2 + M * D2
215 Q3 = Q3 + M * D3
230 PSET (T / N, 90 - Q1 / P1), 3
235 PSET (T / N, 70 - J / J0), 3
250 PSET (T / N, 180 - (Q2 / P2)), 1
253 PSET (T / N, 30 - I / I0), 3
257 PSET (T / N, 180 - (Q3 / P3)), 2
259 PSET (T / N, 50 - E / E0), 3
260 T = T + M
270 IF (T / N) < 319 GOTO 165


 
21. Alagadiço (zonas de terra alagadas com crescimento de pastos e produção de turfa)


Programa wetland.

Exemplo de programação em Basic: wetland.bas
2 REM WETLAND.bas (Wet grassland)
3 SCREEN 1, 0: COLOR 7, 0
4 LINE (0, 0)-(319, 180), 3, B
5 LINE (0, 50)-(319, 50), 3
6 LINE (0, 120)-(319, 120), 3
10  DIM A(12)

11  DATA  2300,2000,2800,3400,4000,4200,4400,3600,3700,3300,3000,3200
12  FOR M = 0 TO 11
13  READ A(M)
14  NEXT

20 J = .005
30 Q1 = .05
40 Q2 = 100
50 Q3 = 200
60 K0 = 113
65 K1 = .268
70 K2 = .005
75 K3 = .0000016
80 K4 = .03
85 K5 = .000268
90 K6 = .0016
95 K7 = .005
100 K8 = 1.34E-06
 

110 N0 = 20
120 Q0 = .02
125 P0 = .02
130 I0 = .01
135 TX = 3560
140 T0 = 320 / TX
145 DT = 4

150 I = A(M)
200 R = I / (1 + K0 * Q1)
205 D1 = J + K8 * Q2 - K4 * Q1 + K3 * Q3 - K5 * Q1 * R
210 D2 = K1 * Q1 * R - K2 * Q2 - K7 * Q2
220 D3 = K2 * Q2 - K6 * Q3
230 Q1 = Q1 + D1 * DT
235  IF Q1 < .0000001 THEN Q1 = .0000001
240 Q2 = Q2 + D2 * DT
300 Q3 = Q3 + D3 * DT
310 T = T + DT
315 Y = INT(T / 365)
318 W = 12 * Y
320 M = INT(T / 30.4) - W
400  PSET (T0 * T, 180 - Q1 * N0), 1
410  PSET (T0 * T, 120 - Q2 * Q0), 2
420  PSET (T0 * T, 120 - Q3 * P0), 1
430  PSET (T0 * T, 50 - I * I0), 1
440  IF T * T0 < 319 GOTO 150               


 
22. Lagoa beira mar (área natural com mangue, rios e mares, que é um criadouro natural de camarão)


Programa lagoon

Exemplo de programação em Basic: lagoon.bas
 
3 REM LAGOON.BAS (Tropical marine lagoon)
5 SCREEN 12, 0: REMCOLOR 0, 0
6 LINE (0, 0)-(320, 60), , B
7 LINE (0, 70)-(320, 130), , B
8 LINE (0, 140)-(320, 200), , B
9 LINE (0, 210)-(320, 270), , B
20 REM Scaling Factors
21 DT = .002
23 T0 = 35
27 C0 = .4
30 P0 = .02
35 D0 = .02
37 M0 = .005
40 Y0 = 4
43 N0 = 10
44 N = .1
47 REM Sources and Starting Conditions
52 L = .1
55 NG = 1
60 NR = 1
65 J = 1
66 A = .1
67 M = 10000
68 D = 100
69 C = 2
70 S = 2
80 TD = 1
100REMCoefficients
110 K0 = 45
120 K1 = 90
130 K2 = 200000!
140 K3 = .167
150 K4 = .01
155 K5 = .067
157 K6 = .0167
160 K7 = 10
165 K8 = 1.5
168 K9 = .33
170 L1 = 50
175 L2 = 5000!
178 L3 = .01
180 L4 = .3
190 L5 = 240
195 L6 = .05
197 L7 = 3
198 L8 = .02
199 L9 = .002

200 REM Plotting
220 REM PSET (T*T0, 270 - M*M0)
230 PSET (T * T0, 270 - P * P0), 1
250 PSET (T * T0, 200 - D * D0), 2
285 PSET (T * T0, 60 - Y * Y0), 1
290 PSET (T * T0, 130 - C * C0), 2
300 REM EQUATIONS
305 R1 = J / (1 + K0 * A * N0)
310 R2 = J / (1 + K1 * N)
317 N0 = (NR + NG) / (1 + L1 * A * R1)
320 N = (N0 + L6 * C) / (1 + L5 * R2)
325 IF N < .00001 THEN N = .00001
330 P = K2 * R2 * N
340 Y = K9 * S * TD
350 REM NEW VALUES OF STORAGES
352 DD = K2 * N * R2 - K3 * D - K5 * C * D + K4 * M - L8 * D
355 DC = K6 * D * C - K7 * C - K8 * C * S + L9 * D
360 DS = L4 * C * S - K9 * S - L7 * S + L * TD
365 DM = L2 * A * N0 * R1 - L3 * M
370 D = D + DD * DT
375 C = C + DC * DT
377 IF C < .0001 THEN C = .0001
380 S = S + DS * DT
390 M = M + DM * DT
392 T = T + DT
400 IF T * T0 < 320 GOTO 200



Index Java  | Index Livro | Odum | LEIA | Top