bio_gsj.nml 5.68 KB
#$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]
!  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.0
  mu4     = 0.02
  mu5     = 0.00
  w_d     = -5.0
  kc      = 0.03
 /