8.5. Questions and Tasks
161
REGEN_Si3 = CONVEYOR OUTFLOW
SILICATE(t)
SILICATE (t-dt) + (REGEN_Si3-Si_LOSS) * dt
INIT SILICATE = 6e-3
REGEN_Si3 = CONVEYOR OUTFLOW
Si_LOSS = CHECK*(TOTAL_DIATOMS)
DIA_UPTAKE_Si = Vrnax_Si*(SILICATE /(SILICATE+Ks_Si))
CHECK = if DIA_UPTAKE_Si /Vrnax_Si > = . 9 9 then Vrnax_Si
else DIA_UPTAKE_Si
Ks_Si = 0
Vrnax_Si = DIA_UPTAKE_N*1 .0
Z = if Length_of_Rotation_Period = 0 then
Mixed_Layer_Depth else
(Mixed_Layer_Depth /2*(1+SIN(2*PI*TIME /Length_of_Rotati°
n_Period) ) )
Julian_Day = IF TIME< = 365 THEN TIME ELSE (TIME365*(INT(TIME /365)))
Length_of_Rotation_Period = GRAPH(Julian_Day)
(0.00, 0.00),
(146 , 0.00),
(292, 0.00),
(36 .5, 0.00), (73.0, 0.00), (110, 0.00),
(183, 0.00),
(329, 0 .00) ,
(219, 0 .00),
(365 , 0 .00)
(256, 0.00) ,
Mixed_Layer_Depth = GRAPH(Julian_Day)
(0.00, 20.0),
(146, 20 .0),
(292, 20.0),
(36.5, 20.0), (73 .0, 20.0), (110, 20 .0),
(183, 20 .0),
(329, 20.0),
(219, 20.0),
(365, 20.0)
(256, 20.0),
DIATOMS(t) = DIATOMS(t-dt) + (DIATOM_DB-DIA_HOLD1) * dt
INIT DIATOMS = 50
DIATOM_DB = DIA_GROWTH*(DIATOMS+DIA_HOLD2)
DIA_HOLD1 = DIATOMS-50
DIA_HOLD2(t) = DIA_HOLD2(t-dt) + (DIA_HOLD1-DIA_HOLD3)
* dt
INIT DIA_HOLD2 = 0
TRANSIT TIME = 4
INFLOW LIMIT = INF
CAPACITY = INF
DIA_HOLD1 = DIATOMS-50
DIA_HOLD3 = CONVEYOR OUTFLOW
DIA_HOLD4(t) = DIA_HOLD4(t-dt) + (DIA_HOLD3-DIA_HOLD5)
* dt