bio_gsj_sansHCsurzoo.nml 5.75 KB
Newer Older
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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
#$Id$
!-------------------------------------------------------------------------------
! Fasham et al. biological model with modifications by Kuehn and Radach
! 
!  numc=        number of compartments for geobiochemical model
!
!  p1_init = initial flagellate concentration                    [mmol n/m3]
!  p2_init = initial diatom concentration                        [mmol n/m3]
!  z1_init = initial micro-zooplakton concentration              [mmol n/m3] 
!  z2_init = initial meso-zooplakton concentration               [mmol n/m3] 
!  b_init  = initial bacteria concentration                      [mmol n/m3]
!  d_init  = initial detritus concentration                      [mmol n/m3]
!  l_init  = initial LDON concentration                          [mmol n/m3]
!  p0      = minimum phytoplankton concentration                 [mmol n/m3]
!  z0      = minimum zooplakton concentration                    [mmol n/m3]
!  b0      = minimum bacteria concentration                      [mmol n/m3]
!  mte     = if .true. use temperature-dependent metabolic rates
!  ca1     = temp-dependence coeff for p1
!  ca2     = temp-dependence coeff for p2
!  ch1     = temp-dependence coeff for z1
!  ch2     = temp-dependence coeff for z2
!  amratio = Mass ratio between p2 and p1
!  hmratio = Mass ratio between z2 and z1
!  vp1     = maximum flagellate uptake rate by flagellates       [1/day] 
!  vp2     = maximum diatom uptake rate by diatoms               [1/day] 
!  alpha1  = slope of the flagellate PI-curve                    [m2/(W day)] 
!  alpha2  = slope of the diatom PI-curve                        [m2/(W day)] 
!  inib1   = inhibition slope of the flagellate PI-curve (pos.)  [m2/(W day)]
!  inib2   = inhibition slope of the PI-curve (pos.)             [m2/(W day)]
!  kn1     = half sat. constant nitrate uptake by p1             [mmol n/m3]
!  ka1     = half sat. constant ammonium uptake by p1            [mmol n/m3]
!  kn2     = half sat. constant nitrate uptake by p2             [mmol n/m3]
!  ka2     = half sat. constant ammonium uptake by p2            [mmol n/m3]
!  mu11    = mortality rate for p1                               [1/day] 
!  mu12    = mortality rate for p2                               [1/day] 
!  k5      = half sat. constant phy. mortality                   [mmol n/m3]
!  gamma   = exudation fraction                                  [-]
!  w_p1    = settling velocity for p1                            [m/day]
!  w_p2    = settling velocity for p2                            [m/day]
!  theta   = phytoplancton buoyancy parameter                    [m3 day/(mmol N)]
!  g1max   = maximum ingestion rate for z1                       [1/day]
!  g2max   = maximum ingestion rate for z2                       [1/day]
!  k3      = half saturation constant ingestion                  [mmol n/m3]
!  beta    = grazing efficiency                                  [-]
!  k6      = half saturation zooplankton loss (z1 & z2)          [mmol n/m3]
!  mu21    = maximum loss rate for z1                            [1/day]
!  mu22    = maximum loss rate for z2                            [1/day]
!  delta   = fractional zooplankton loss to LDON (z1 & z2)       [-]
!  epsi    = fractional zooplankton loss to ammonium (z1 & z2)   [-]
!  r11     = z1 grazing preference on p1                         [-]
!  r12     = z1 grazing preference on p2                         [-] 
!  r13     = z1 grazing preference on bacteria                   [-] 
!  r14     = z1 grazing preference on detritus                   [-] 
!  r21     = z2 grazing preference on p1                         [-]
!  r22     = z2 grazing preference on p2                         [-] 
!  r23     = z2 grazing preference on detritus                   [-] 
!  r24     = z2 grazing preference on z1                         [-] 
!  vb1     = maximum bac1 uptake rate                            [1/day]
!  vb2     = maximum bac2 uptake rate                            [1/day]
!  k4      = half saturation bacterial uptake                    [mmol n/m3]
!  mu3     = bacteria excretion rate                             [1/day]  
!  etaa    = uptake ratio ammonium:LDON                          [-]  
!  etah    = uptake ratio hydrocarbon:LDON                       [-]  
!  mu4     = detritus breakdown rate                             [1/day] 
!  mu5     = nitrification rate                                  [1/day] 
!  w_d     = detritus settling velocity                          [m/day]
!  w_h     = hydrocarbon settling velocity                       [m/day]
!  kc      = attenuation constant for the self shading effect    [m**2/mmol N]
!-------------------------------------------------------------------------------
 &bio_gsj_nml
  numc    = 11 
  p1_init = 0.012
  p2_init = 0.012
  z1_init = 0.012
  z2_init = 0.012
  b_init  = 0.001
  d_init  = 0.01
  l_init  = 0.1
  p0      = 0.0001
  z0      = 0.0001
  b0      = 0.0001
  mte     = .true.
  ca1     = 3.61
  ca2     = 14.58
  ch1     = 3.265
  ch2     = 24.923
  amratio = 200
  hmratio = 1000
  vp1     = 0.02
  vp2     = 0.8
  alpha1  = 0.02
  alpha2  = 0.04
  inib1   = 0.0
  inib2   = 0.006
  kn1     = 1.0
  ka1     = 0.8
  kn2     = 1.0
  ka2     = 0.8
  mu11    = 0.05
  mu12    = 0.05
  k5      = 0.2
  gamma   = 0.05
  w_p1    =-0.38
  w_p2    =-0.00
  theta   = 0.0
  w_p1min =-0.01
  w_p1max =-0.10
  w_p2min =-0.05
  w_p2max =-0.38
  g1max   = 1.0
  g2max   = 1.0
  k3      = 1.0
  beta    = 0.625
  mu21    = 0.3
  mu22    = 0.3
  k6      = 0.2
  delta   = 0.1
  epsi    = 0.70
  r11     = 0.55
  r12     = 0.30
  r13     = 0.05
  r14     = 0.10
  r21     = 0.50
  r22     = 0.30
  r23     = 0.05
  r24     = 0.15
  vb1     = 0.24
  vb2     = 0.24
  k4      = 0.5
  k10     = 0.5
  w_h     = 100.0
  mu3     = 0.03
  etaa    = 0.0
  etah    = 0.2
  mu4     = 0.02
  mu5     = 0.00
  w_d     = -5.0
  kc      = 0.03
 /