blob: fe05881eb9025bdc0f0a8d68d777b68a6f8319bb (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
n t=15
note*** wasserökosystem nach abel und reich
note*** in: microextra 4/83 seite 34 ff
note************************************************************************
note* hilfsgleichung fuer temperatur t
note* die zeit time in wochen
a t.k=15+4*sin((time.k-10)*2*pi/52) temperatur t; time in wochen
c pi=3.1415
note gleichung fuer phytoplankton p
l p.k=p.j+dt*(p.j*(p1*n.j*t.j-p2*z.j)(100-p.j)/100) phytoplankton p
note gleichung fuer zooplankton z
l z.k=z.j+dt*(z.j*(p3*t.j*p.j+p4*n.j-(p5*f.j+p6*b.j)-1/p.j)(30-z.j)/30)
note gleichung fuer fische f
l f.k=f.j+dt*(f.j*(p7*z.j-p8*b.j-p9/(z.j+p.j))(10-f.j)/10)
note gleichung fuer raubfisch barsch b
l b.k=b.j+dt*(b.j*(p10*f.j+p11*z.j-1/(p12*f.j))(0.1-b.j)/0.1)
note **** gleichung fuer naehrstoffmenge n
l n.k=n.j+dt*(p13-n.j*(p14*p.j-p15*z.j))
note **** anfangswerte ****************************************************
n p=p0
n z=z0
n f=f0
n b=b0
n n=n0
c p0=10
c z0=3
c f0=1
c b0=0.01
c n0=30 in kg/volumeneinheit bzw. Stück/volumeneinhe�[
note ***** konstanten ********************************************************
c p1=0.006
c p2=1
c p3=0.006
c p4=0.03
c p5=1
c p6=100
c p7=0.33
c p8=100
c p9=1E-4
c p10=1
c p11=1
c p12=0.25
c p13=10
c p14=0.1
c p15=0.2
note **** simulationskonstanten *********************************************
c dt=0.5
c length=60
c pltper=1
note***** outputvariablen****************************************************
a lp.k=ln(p.k/p0)
a lz.k=ln(z.k/z0)
a lf.k=ln(f.k/f0)
a lb.k=ln(b.k/b0)
a logn.k=ln(n.k/n0)
plot lp=p,lz=z,lf=f,lb=b,logn=n(-4,4)
|
