Commit dd33917e authored by Philippe Klotz's avatar Philippe Klotz

Merge branch 'master' of gitlasso.uqar.ca:dumoda01/gotm_ismer

parents 54921826 b9e7fd9a
This diff is collapsed.
2007-06-01 00:00:00 4 2
-1.0 0.5 0.0
-20.0 0.5 0.0
-30.0 -0.5 0.0
-40.0 -0.5 0.0
-50.0 0.0 0.0
2007-06-15 00:00:00 4 2
-1.0 0.5 0.0
-20.0 0.5 0.0
-30.0 -0.5 0.0
-40.0 -0.5 0.0
-50.0 0.0 0.0
This diff is collapsed.
This diff is collapsed.
......@@ -575,6 +575,7 @@
case (5)
call surface_fluxes_mab(nlev,t(nlev),s(nlev))
case (6)
case (7)
end select
if (bio_eulerian) then
......@@ -794,6 +795,7 @@
if (allocated(lumlim2)) deallocate(lumlim2)
if (allocated(nitlim2)) deallocate(nitlim2)
if (allocated(ammlim2)) deallocate(ammlim2)
if (allocated(ppnet)) deallocate(ppnet)
! The external provide arrays
if (allocated(h)) deallocate(h)
......@@ -898,6 +900,9 @@
allocate(ammlim2(0:nlev),stat=rc)
if (rc /= 0) stop 'init_bio(): Error allocating (ammlim2)'
allocate(ppnet(0:nlev),stat=rc)
if (rc /= 0) stop 'init_bio(): Error allocating (ppnet)'
! The external provide arrays
allocate(h(0:nlev),stat=rc)
if (rc /= 0) stop 'init_bio(): Error allocating (h)'
......
......@@ -595,6 +595,8 @@
dd(a,b,ci)=vb*min67/(k4+min67+cc(l,ci))*(cc(b,ci)+b0)
dd(l,b,ci)=vb*cc(l,ci)/(k4+min67+cc(l,ci))*(cc(b,ci)+b0)
ppnet(ci) =(dd(a,p,ci)+dd(n,p,ci)-dd(p,d,ci)-dd(p,z,ci)-dd(p,l,ci))*secs_pr_day
do i=1,numc
do j=1,numc
pp(i,j,ci)=dd(j,i,ci)
......
......@@ -707,6 +707,9 @@
ws(p1,ci) = w_p1
ws(p2,ci) = w_p2
ppnet(ci) =(dd(a,p1,ci)+dd(n,p1,ci)-dd(p1,d,ci)-dd(p1,z1,ci)-dd(p1,z2,ci)-dd(p1,l,ci) &
+dd(a,p2,ci)+dd(n,p2,ci)-dd(p2,d,ci)-dd(p2,z1,ci)-dd(p2,z2,ci)-dd(p2,l,ci))*secs_pr_day
do i=1,numc
do j=1,numc
pp(i,j,ci)=dd(j,i,ci)
......
......@@ -107,6 +107,8 @@
iret = set_attributes(ncid,nitlim2_id,units='',long_name='nitrate limitation for microphyto')
iret = new_nc_variable(ncid,'ammlim2',NF_REAL,4,dims,ammlim2_id)
iret = set_attributes(ncid,ammlim2_id,units='',long_name='ammonium limitation for microphyto')
iret = new_nc_variable(ncid,'ppnet',NF_REAL,4,dims,ppnet_id)
iret = set_attributes(ncid,ppnet_id,units='1/day',long_name='net primary production rate')
!DD Diagnostic de npar (nb de particules lagrangiennes) pour bebogage
!iret = new_nc_variable(ncid,'npar',NF_REAL,4,dims,npar_id)
......@@ -138,6 +140,7 @@
iret = store_data(ncid,lumlim2_id,XYZT_SHAPE,nlev,array=lumlim2(:))
iret = store_data(ncid,nitlim2_id,XYZT_SHAPE,nlev,array=nitlim2(:))
iret = store_data(ncid,ammlim2_id,XYZT_SHAPE,nlev,array=ammlim2(:))
iret = store_data(ncid,ppnet_id,XYZT_SHAPE,nlev,array=ppnet(:))
!DD
!iret = store_data(ncid,npar_id,XYZT_SHAPE,nlev,array=npar(:))
......
......@@ -22,6 +22,7 @@
REALTYPE, dimension(:), allocatable :: zlev
REALTYPE, dimension(:), allocatable :: par
REALTYPE, dimension(:), allocatable :: lumlim1,nitlim1,ammlim1,lumlim2,nitlim2,ammlim2 !DD
REALTYPE, dimension(:), allocatable :: ppnet
REALTYPE, dimension(:,:), allocatable :: cc,ws
integer :: surface_flux_method=-1
integer :: n_surface_fluxes=-1
......@@ -37,6 +38,7 @@
integer :: wp_id !CHG2
integer :: par_id !CHG4
integer :: lumlim1_id,nitlim1_id,ammlim1_id,lumlim2_id,nitlim2_id,ammlim2_id !DD
integer :: ppnet_id
character(len=64), dimension(:), allocatable :: var_names
character(len=64), dimension(:), allocatable :: var_units
character(len=64), dimension(:), allocatable :: var_long
......
!$Id: airsea.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
!-------------------------------------------------------------------------------
!
!-------------------------------------------------------------------------------
! air-sea interaction (heat, momentum and freshwater fluxes and solar radiation)
!
! calc_fluxes -> surface fluxes calculated by means of bulk formulae
! (.true./.false.). Solar radiation is calculated from
! time, latitude, longitude and clouds. If (.true.),
! meteo_file must be given and wet_mode must be specified.
! If (.false.), surface fluxes and solar radiation are
! prescribed.
!
! meteo_file -> file with meteo data (for calc_fluxes=.true.) with
! - date (yyyy-mm-dd hh:mm:ss)
! - x-comp. of wind (10 m) in m/s
! - y-comp. of wind (10 m) in m/s
! - air pressure ( 2 m) in hectopascal
! - dry air temp. ( 2 m) in Celsius
! - relative humidity in % or wet bulb temperature in C
! or dew point temperature in C (depending on wet_mode)
! - cloud cover in 1/10
!
! wet_mode -> decides what is given in 7. column in meteo_file
! 1: relative humidity
! 2: wet bulb temperature
! 3: dew point temperature
!
! heat_method -> method to provide short wave radiation (swr) and
! surface heat flux (qh)
! (only for calc_fluxes=.false.)
! 0: heat flux not prescribed
! 1: constant "const_swr" and "const_qh" given (see below)
! 2: swr and qh are read from heatflux_file
! const_swr -> constant value of incoming short wave radiation in W/m^2
! (always positive)
!
! const_qh -> constant value of surface heat flux in W/m^2
! (negative for heat loss)
!
! heatflux_file -> file with qin and qout given in W/m^2
! (negative for net outgoing)
!
! momentum_method -> method how momentum fluxes are given
! (only for calc_fluxes=.false.)
! 0: momentum flux not prescribed
! 1: constant surface momentum fluxes given
! 2: surface momentum fluxes given from file momentumflux_file
!
! const_tx -> x-component of surface momentum flux in N/m^2
! const_ty -> y-component of surface momentum flux in N/m^2
!
! momentumflux_file-> file with tx and ty given in N/m^2
!
! p_e_method -> method how fresh water fluxes (P-E) are given
! 0: P-E not used
! 1: constant value for P-E (in m/s) used
! (P-E = precipitation-evaporation)
! 2: values for P-E read from file
!
! const_p_e -> constant value for P-E in m/s (positive for P>E)
!
! p_e_flux_file -> file with value for P-E (positive for P>E)
! used if p_e_method=2
!
! sst_method -> method how sea surface temperature (SST) is given
! 0: no independent SST observation is read from file
! 2: independent SST observation is read from file,
! only for output
!
! sst_file -> file with independent SST observation
!
! sss_method -> method how sea surface salinity (SSS) is given
! 0: no independent SSS observation is read from file
! 2: independent SSS observation is read from file,
! only for output
!-------------------------------------------------------------------------------
&airsea
calc_fluxes= .false.
meteo_file= 'meteo.dat'
wet_mode= 1
heat_method= 0
const_swr= 200.0
const_heat= 0.0
heatflux_file= 'heatflux.dat'
momentum_method= 0
const_tx= 0.1
const_ty= 0.0
momentumflux_file='momentumflux.dat'
p_e_method= 0
const_p_e= 0.
p_e_flux_file= 'p_e.dat'
sst_method= 0
sst_file= 'sst.dat'
sss_method= 0
sss_file= 'sss.dat'
/
!$Id$
!-------------------------------------------------------------------------------
! Basic settings for biogeochemical model
!
! bio_calc -> calculation of the bio model 'bio_model' (.true./.false.)
! bio_model -> choice of the bio model:
! 1: NPZD (4 variables)
! 2: IOW-ERGOM (9 variables)
! 3: Suspended matter only (1 variable)
! 4: Fasham et al. 1990 (7 variables)
! 5: IOW-ERGOM MaBenE version (9 variables)
! 6: ISMER model (9 variables)
!
! bio_eulerian -> state variables are Eulerian (.true./.false.)
!
! cnpar -> Cranck-Nicolson parameter for vertical diffusion
!
! w_adv_discr -> advection scheme for vertical motion
! 1: first order upstream
! 2: not coded yet
! 3: third-order polynomial
! 4: TVD with Superbee limiter
! 5: TVD with MUSCL limiter
! 6: TVD with ULTIMATE QUICKEST
!
! ode_method -> ODE scheme for source and sink dynamics
! 1: first-order explicit (not positive)
! 2: second order explicit Runge-Kutta (not positive)
! 3: fourth-order explicit Runge-Kutta (not positive)
! 4: Patankar (first order, not conservative)
! 5: Patankar-RK (second order, not conservative)
! 6: Patankar-RK (does not work, not conservative)
! 7: Modified Patankar (1. order, conservat., posit.)
! 8: Modified Patankar-RK (2. order, conservat., posit.)
! 9: Modified Patankar-RK (does not work, conservat., posit.)
! 10: Ext. Modified Patankar (1. order, conservat., posit.)
! 11: Ext. Modified Patankar-RK (2. order, conservat., posit.)
!
! split_factor -> number of biogeochemical time steps per physical time step
!
! bioshade_feedback -> feedback of bio-turbidity to temp. eq. (.true./.false.)
!
! bio_lagrange_mean -> averaging Lagrangian conc. on output (.true./.false.)
!
! bio_npar -> total number of Lagrangian particles
!-------------------------------------------------------------------------------
&bio_nml
bio_calc= .true.
bio_model= 4
bio_eulerian= .true.
cnpar= 1.0
w_adv_discr= 4
ode_method= 11
split_factor= 1
bioshade_feedback= .true.
bio_lagrange_mean= .true.
bio_npar= 10000
/
#$Id$
!-------------------------------------------------------------------------------
! Fasham et al. biological model with modifications by Kuehn and Radach
!
! numc= number of compartments for geobiochemical model
!
! p_initial= initial phytoplankton concentration [mmol n/m3]
! z_initial= initial zooplakton concentration [mmol n/m3]
! b_initial= initial bacteria concentration [mmol n/m3]
! d_initial= initial detritus concentration [mmol n/m3]
! n_initial= initial nitrate concentration [mmol n/m3]
! a_initial= initial ammonium concentration [mmol n/m3]
! l_initial= 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]
! vp = maximum phytoplankton uptake rate [1/day]
! alpha = slope of the PI-curvea [m2/(W day)]
! k1 = half saturation constant nitrate uptake [mmol n/m3]
! k2 = half saturation constant ammonium uptake [mmol n/m3]
! mu1 = phytoplankton mortality rate [1/day]
! k5 = half saturation constant phytoplankton mortality [mmol n/m3]
! gamma = exudation fraction [-]
! w_p = phytoplankton settling velocity [m/day]
! gmax = maximum ingestion rate [1/day]
! k3 = half saturation constant ingestion [mmol n/m3]
! beta = grazing efficiency [-]
! mu2 = maximum zooplankton loss rate [1/day]
! k6 = half saturation zooplankton loss [mmol n/m3]
! delta = fractional zooplankton loss to LDON [-]
! epsi = fractional zooplankton loss to ammonium [-]
! r1 = grazing preference phytoplankton [-]
! r2 = grazing preference bacteria [-]
! r3 = grazing preference detritus [-]
! vb = maximum bacterial uptake rate [1/day]
! k4 = half saturation bacterial uptake [mmol n/m3]
! mu3 = bacteria excretion rate [1/day]
! eta = uptake ratio ammonium:LDON [-]
! mu4 = detritus breakdown rate [1/day]
! w_d = detritus settling velocity [m/day]
! kc = attenuation constant for the self shading effect [m**2/mmol N]
! I_min = minimum photosynthetically active radiation (PAR) [W/m**2]
! I_opt = optimal photosynthetically active radiation (PAR) [W/m**2] !CHG1
! inib = inhibition slope of the PI-curve (positive) [m2/(W day)] !CHG1
! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)] !CHG2
!-------------------------------------------------------------------------------
&bio_fasham_nml
numc= 7
p_initial= 0.0001
z_initial= 0.0001
b_initial= 0.0001
d_initial= 0.001
l_initial= 0.001
p0= 0.00001
z0= 0.00001
b0= 0.00001
vp= 1.5
alpha= 0.45
inib= 0.001
I_opt= 5.0
k1= 1.5
k2= 1.5
mu1= 0.05
k5= 0.3
gamma= 0.05
w_p= -0.38
theta= 0.0
w_pmin= -0.06
w_pmax= -0.38
gmax= 1.0
k3= 1.0
beta= 0.625
mu2= 0.3
k6= 0.2
delta= 0.1
epsi= 0.70
r1= 0.55
r2= 0.3
r3= 0.15
vb= 0.24
k4= 0.5
mu3= 0.03
eta= 0.0
mu4= 0.02
w_d= -5.0
kc= 0.03
/
#$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 [-]
! vb = maximum bacterial uptake rate [1/day]
! k4 = half saturation bacterial uptake [mmol n/m3]
! mu3 = bacteria excretion rate [1/day]
! eta = uptake ratio ammonium: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
vb = 0.24
k4 = 0.5
k10 = 0.5
w_h = 0.0
mu3 = 0.03
eta = 0.0
mu4 = 0.02
mu5 = 0.00
w_d = -5.0
kc = 0.03
/
#$Id$
!-------------------------------------------------------------------------------
! IOW biological model according to Neumann et al. 2002 and Burchard et al. 2005
!
! numc -> number of compartments for geobiochemical model
! p1_initial -> initial diatoms concentration [mmol n/m3]
! p2_initial -> initial flagellates concentration [mmol n/m3]
! p3_initial -> initial cyanobacterias concentration [mmol n/m3]
! zo_initial -> initial zooplankton concentration [mmol n/m3]
! de_initial -> initial detritus concentration [mmol n/m3]
! am_initial -> initial ammonium concentration [mmol n/m3]
! ni_initial -> initial nitrate concentration [mmol n/m3]
! po_initial -> initial phosphate concentration [mmol p/m3]
! o2_initial -> initial oxygen concentration [mmol n/m3]
! sfl_po -> constant surface phosphate flux [mmol n/m2/d]
! sfl_am -> constant surface ammonium flux [mmol n/m2/d]
! sfl_ni -> constant surface nitrate flux [mmol n/m2/d]
! surface_flux_method -> 0: constant; 2: from file 'bio_fluxes.dat'
! fluff -> .true.: caluclate fluff layer
! fl_initial -> initial fluff concentration [mmol n/m2]
! p10= -> minimum diatoms concentration [mmol n/m3]
! p20= -> minimum flagellates concentration [mmol n/m3]
! p30= -> minimum cyanobacterias concentration [mmol n/m3]
! zo0= -> minimum zooplankton concentration [mmol n/m3]
! w_p1= -> sinking speed diatom [m/d]
! w_p2= -> sinking speed flagellates [m/d]
! w_p3= -> sinking speed cyanobacteria [m/d]
! w_de= -> sinking speed detritus [m/d]
! kc= -> attenuation constant for the self shading effect [m2/mmol N]
! i_min= -> minimum photosynthetically active radiation (PAR) [W/m2]
! r1max= -> maximum growth rate diatoms [1/d]
! r2max= -> maximum growth rate flagellates [1/d]
! r3max= -> maximum growth rate cyanobacteria [1/d]
! alpha1= -> half saturation diatoms [mmol n/m3]
! alpha2= -> half saturation flagellates [mmol n/m3]
! alpha3= -> half saturation cyanobacteria [mmol n/m3]
! lpa= -> p --> a rate (p respiration) [1/d]
! lpd= -> p --> d rate (p mortality) [1/d]
! tf= -> temperature control flagellates [deg C]
! tbg= -> temperature control cyanos [deg C]
! beta_bg= -> temperature control cyanos [1/deg C]
! g1max= -> maximum grazing rate on diatoms [1/d]
! g2max= -> maximum grazing rate on flagellates [1/d]
! g3max= -> maximum grazing rate on cyanos [1/d]
! lza= -> z --> a rate (z exkudation) [m3/d/mmol n]
! lzd= -> z --> d rate (z mortality) [m3/d/mmol n]
! iv= -> Ivlev constant [1/(mmol n/m3)3]
! topt= -> optimum temperature [deg C]
!
! lan= -> nitrification constant [1/d]
! oan= -> oxygen parameter [ - ]
! beta_an= -> temperature control [1/deg C]
!
! lda= -> detritus mineralisation constant [1/d]
! tda= -> temperature control [deg C]
! lds= -> rate of detritus sinking into sediment [m/d]
! lsa= -> rate of sediment mineralisation [1/d]
! bsa= -> temperature control of sediment mineralisation [1/deg C]
! ph1= -> release of phosphate [ - ]
! ph2= -> release of phosphate [ - ]
! pvel= -> piston velocity [m/d]
!
! sr= -> redfield ratio(p/n) [ - ]
! s1= -> reduced nitrate/oxidized detritus [ - ]
! s2= -> oxygen produktion/recycled nitrogen [ - ]
! s3= -> oyxgen produktion related to n [ - ]
! s4= -> nitrification [ - ]
!
! a0= -> osat Parameter [mmol/m**3]
! a1= -> osat Parameter [ - ]
! a2= -> osat Parameter [1/ deg C]
!
! aa -> weighting factor for visible light fraction [-]
! g2 -> attenuation length for visible light [m]
!-------------------------------------------------------------------------------
&bio_iow_nml
numc= 9
p1_initial= 0.1
p2_initial= 0.1
p3_initial= 0.1
zo_initial= 0.1
de_initial= 5.0
am_initial= 0.1
ni_initial= 5.0
po_initial= 0.6
o2_initial= 85.0
sfl_po= -0.055
sfl_am= -0.36
sfl_ni= -0.45
surface_flux_method= 0
fluff= .true.
fl_initial= 0.0
p10= 0.0045
p20= 0.0045
p30= 0.0045
zo0= 0.0045
w_p1= -0.5
w_p2= 0.
w_p3= 0.1
w_de= -3.
kc= 0.03
i_min= 25.
r1max= 2.0
r2max= 0.7
r3max= 0.5
alpha1= 1.35
alpha2= 0.675
alpha3= 0.5
lpa= 0.01
lpd= 0.02
tf= 10.
tbg= 14.
beta_bg= 1.
g1max= 0.5
g2max= 0.5
g3max= 0.25
lza= 0.0666666666
lzd= 0.1333333333
iv= 0.24444444
topt= 20.
lan= 0.1
oan= 0.01
beta_an= 0.11
lda= 0.003
tda= 13.
beta_da= 20.
lds= 3.5
lsa= 0.001
bsa= 0.15
ph1= 0.15
ph2= 0.1
pvel= 5.
sr= 0.0625
s1= 5.3
s2= 6.625
s3= 8.125
s4= 0.6666666666666
a0= 31.25
a1= 14.603
a2= 0.4025
aa= 0.78
g2= 7.9
/
#$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]
! n_init = *** see obs.nml *** [mmol n/m3]
! a_init = *** see obs.nml *** [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]
! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)] !CHG2
! 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)] !CHG1
! inib2 = inhibition slope of the PI-curve (pos.) [m2/(W day)] !CHG1
! kn1 = half sat. constant nitrate uptake by fla [mmol n/m3]
! ka1 = half sat. constant ammonium uptake by fla [mmol n/m3]
! kn2 = half sat. constant nitrate uptake by diatoms [mmol n/m3]
! ka2 = half sat. constant ammonium uptake by diatoms [mmol n/m3]