New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
namelist_medusa_ref in branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/CONFIG/SHARED – NEMO

source: branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/CONFIG/SHARED/namelist_medusa_ref @ 6057

Last change on this file since 6057 was 6057, checked in by jpalmier, 9 years ago

JPALM -- 15-12-2015 -- mini update

File size: 16.6 KB
Line 
1!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!
3! BIO namelist :  MEDUSA biological model option and parameter input
4! -------------
5!
6!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
7!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
8!       natbio   Shared parameters for dynamics/advection/thermo
9!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10!
11!! Primary production and chl related quantities
12!!       xxi         :  conversion factor from gC to mmolN
13!!       xaln        :  Chl-a specific initial slope of P-I curve for non-diatoms
14!!       xald        :  Chl-a specific initial slope of P-I curve for diatoms
15!!       jphy        :  phytoplankton T-dependent growth switch
16!!       xvpn        :  maximum growth rate for non-diatoms
17!!       xvpd        :  maximum growth rate for diatoms
18!!       xthetam     :  maximum Chl to C ratio for non-diatoms
19!!       xthetamd    :  maximum Chl to C ratio for diatoms
20!!
21!! Diatom silicon parameters
22!!       xsin0       :  minimum diatom Si:N ratio
23!!       xnsi0       :  minimum diatom N:Si ratio
24!!       xuif        :  hypothetical growth ratio at infinite Si:N ratio
25!!
26!! Nutrient limitation
27!!       jliebig     :  Liebig nutrient uptake switch
28!!       xnln        :  half-sat constant for DIN uptake by non-diatoms
29!!       xnld        :  half-sat constant for DIN uptake by diatoms
30!!       xsl         :  half-sat constant for Si uptake by diatoms
31!!       xfld        :  half-sat constant for Fe uptake by diatoms
32!!       xfln        :  half-sat constant for Fe uptake by non-datoms
33!!
34!! Grazing
35!!       xgmi        :  microzoo maximum growth rate
36!!       xgme        :  mesozoo maximum growth rate
37!!       xkmi        :  microzoo grazing half-sat parameter
38!!       xkme        :  mesozoo grazing half-sat parameter
39!!       xphi        :  micro/mesozoo grazing inefficiency
40!!       xbetan      :  micro/mesozoo N assimilation efficiency
41!!       xbetac      :  micro/mesozoo C assimilation efficiency
42!!       xkc         :  micro/mesozoo net C growth efficiency
43!!       xpmipn      :  grazing preference of microzoo for non-diatoms
44!!       xpmid       :  grazing preference of microzoo for diatoms
45!!       xpmepn      :  grazing preference of mesozoo for non-diatoms
46!!       xpmepd      :  grazing preference of mesozoo for diatoms
47!!       xpmezmi     :  grazing preference of mesozoo for microzoo
48!!       xpmed       :  grazing preference of mesozoo for detritus
49!!
50!! Metabolic losses
51!!       xmetapn     :  non-diatom metabolic loss rate
52!!       xmetapd     :  diatom     metabolic loss rate
53!!       xmetazmi    :  microzoo   metabolic loss rate
54!!       xmetazme    :  mesozoo    metabolic loss rate
55!!
56!! Metabolic losses
57!!       xmetapn     :  non-diatom metabolic loss rate
58!!       xmetapd     :  diatom     metabolic loss rate
59!!       xmetazmi    :  microzoo   metabolic loss rate
60!!       xmetazme    :  mesozoo    metabolic loss rate
61!!
62!! Mortality losses
63!!       jmpn        :  non-diatom mortality functional form
64!!       xmpn        :  non-diatom mortality rate
65!!       xkphn       :  non-diatom mortality half-sat constant
66!!       jmpd        :  diatom     mortality functional form
67!!       xmpd        :  diatom     mortality rate
68!!       xkphd       :  diatom     mortality half-sat constant
69!!       jmzmi       :  microzoo   mortality functional form
70!!       xmzmi       :  microzoo   mortality rate
71!!       xkzmi       :  microzoo   mortality half-sat constant
72!!       jmzme       :  mesozoo    mortality functional form
73!!       xmzme       :  mesozoo    mortality rate
74!!       xkzme       :  mesozoo    mortality half-sat constant
75!!
76!! Remineralisation
77!!       jmd         :  detritus T-dependent remineralisation switch
78!!       jsfd        :  accelerate seafloor detritus remin. switch
79!!       xmd         :  detrital nitrogen remineralisation rate
80!!       xmdc        :  detrital carbon   remineralisation rate
81!!
82!! Stochiometric ratios
83!!       xthetapn    :  non-diatom C:N ratio
84!!       xthetapd    :  diatom C:N ratio
85!!       xthetazmi   :  microzoo C:N ratio
86!!       xthetazme   :  mesozoo C:N ratio
87!!       xthetad     :  detritus C:N ratio (FOR IMPLICIT DTC ONLY)
88!!       xrfn        :  phytoplankton Fe:N ratio
89!!       xrsn        :  phytoplankton Si:N ratio (LEGACY; NOT USED)
90!!
91!! Iron parameters
92!!       jiron       :  iron scavenging submodel switch
93!!       xfe_mass    :  iron atomic mass
94!!       xfe_sol     :  aeolian iron solubility
95!!       xfe_sed     :  sediment iron input
96!!       xLgT        :  total ligand concentration (umol/m3)
97!!       xk_FeL      :  dissociation constant for (Fe + L)
98!!       xk_sc_Fe    :  scavenging rate of "free" iron
99!!
100!! Fast-sinking detritus parameters
101!!       jmartin     :  fast detritus remineralisation switch
102!!       jfdfate     :  fate of fast detritus at seafloor switch
103!!       jrratio     :  rain ratio switch
104!!       jocalccd    :  CCD switch
105!!       xridg_r0    :  Ridgwell rain ratio coefficient
106!!       xfdfrac1    :  fast-sinking fraction of diatom nat. mort. losses
107!!       xfdfrac2    :  fast-sinking fraction of mesozooplankton mort. losses
108!!       xfdfrac3    :  fast-sinking fraction of diatom silicon grazing losses
109!!       xcaco3a     :  polar (high latitude) CaCO3 fraction
110!!       xcaco3b     :  equatorial (low latitude) CaCO3 fraction
111!!       xmassc      :  organic C mass:mole ratio, C106 H175 O40 N16 P1
112!!       xmassca     :  calcium carbonate mass:mole ratio, CaCO3
113!!       xmasssi     :  biogenic silicon mass:mole ratio, (H2SiO3)n
114!!       xprotca     :  calcium carbonate protection ratio
115!!       xprotsi     :  biogenic silicon protection ratio
116!!       xfastc      :  organic C remineralisation length scale
117!!       xfastca     :  calcium carbonate dissolution length scale
118!!       xfastsi     :  biogenic silicon dissolution length scale
119!!
120!! Benthos parameters
121!!       jorgben     :  does   organic detritus go to the benthos?
122!!       jinorgben   :  does inorganic detritus go to the benthos?
123!!       xsedn       :  organic   nitrogen sediment remineralisation rate
124!!       xsedfe      :  organic   iron     sediment remineralisation rate
125!!       xsedsi      :  inorganic silicon  sediment dissolution      rate
126!!       xsedc       :  organic   carbon   sediment remineralisation rate
127!!       xsedca      :  inorganic carbon   sediment dissolution      rate
128!!       xburial     :  burial fraction of sediment material
129!!
130!! River parameters
131!!       jriver_n    :  riverine N          input?
132!!       jriver_si   :  riverine Si         input?
133!!       jriver_c    :  riverine C          input?
134!!       jriver_alk  :  riverine alkalinity input?
135!!       jriver_dep  :  depth of riverine   input?
136!!
137!! Miscellaneous
138!!       xsdiss      :  diatom frustule dissolution rate
139!!
140!! Gravitational sinking
141!!       vsed        :  detritus gravitational sinking rate
142!!       xhr         :  coeff for Martin's remineralisation profile
143!!
144!! Additional parameters
145!!       jpkb        :  vertical layer for diagnostic of the vertical flux
146!!
147!! UKESM1 - new diagnostics  !! Jpalm
148!!       jdms        :  include dms diagnostics
149!!
150!!
151&natbio
152!
153!! Primary production and chl related quantities
154      xxi       = 0.01257   !!
155      xaln      = 15.0      !! EKP (05/03/09)
156      xald      = 11.25     !! EKP (05/03/09)
157      jphy      = 1         !! 0 = T-independent; 1 = T-dependent
158      xvpn      = 1.060     !! EKP (05/03/09)
159      xvpd      = 1.000     !! EKP (05/03/09)
160      xthetam   = 0.05      !! EKP (05/03/09)
161      xthetamd  = 0.05      !! EKP (05/03/09)
162!!
163!! Diatom silicon parameters
164      xsin0     = 0.2       !! Mongin et al. (2003)
165      xnsi0     = 0.2       !! Mongin et al. (2003)
166      xuif      = 1.5       !! Mongin et al. (2003)
167!!
168!! Nutrient limitation
169      jliebig   = 0         !! 0 = multiplicative, 1 = Liebig
170      xnln      = 0.5       !! Fasham (1993)
171      xnld      = 0.75      !! EKP (05/03/09)
172      xsld      = 3.00      !! Fasham et al. (2006)-ish
173      xfln      = 0.00033   !! EKP (29/10/09); Parekh et al. (2005)
174      xfld      = 0.00067   !! EKP (29/10/09); Parekh et al. (2005); filtered through Moore et al. (2004)
175!!
176!! Grazing
177      xgmi      = 2.0       !! EKP (05/03/09)
178      xgme      = 0.5       !! AXY (21/12/09)
179      xkmi      = 0.8       !! EKP (26/02/09)
180      xkme      = 0.3       !! EKP (26/02/09)
181      xphi      = 0.20      !! Anderson & Pondaven (2003; 0.23); EKP (05/03/09; 0.20)
182      xbetan    = 0.77      !! Anderson & Pondaven (2003)
183      xbetac    = 0.64      !! Anderson & Pondaven (2003)
184!!    xbetan    = 0.69      !! compromise values inbetween those from ...
185!!    xbetac    = 0.69      !! Anderson & Pondaven (2003)for detrital carbon model
186      xkc       = 0.80      !! Anderson & Pondaven (2003)
187      xpmipn    = 0.75      !!
188      xpmid     = 0.25      !!
189      xpmepn    = 0.15      !! EKP (26/02/09)
190      xpmepd    = 0.35      !! EKP (26/02/09)
191      xpmezmi   = 0.35      !! EKP (26/02/09)
192      xpmed     = 0.15      !! EKP (26/02/09)
193!!
194!! Metabolic losses
195      xmetapn   = 0.02      !!
196      xmetapd   = 0.02      !!
197      xmetazmi  = 0.02      !!
198      xmetazme  = 0.02      !!
199!!
200!! Mortality losses
201      jmpn      = 3         !! 1=linear; 2=quadratic; 3=hyperbolic; 4=sigmoid
202      xmpn      = 0.1       !!
203      xkphn     = 0.5       !! EKP (26/02/09): new parameter
204      jmpd      = 3         !! 1=linear; 2=quadratic; 3=hyperbolic; 4=sigmoid
205      xmpd      = 0.1       !!
206      xkphd     = 0.5       !! EKP (26/02/09): new parameter
207      jmzmi     = 3         !! 1=linear; 2=quadratic; 3=hyperbolic; 4=sigmoid
208      xmzmi     = 0.1       !!
209      xkzmi     = 0.5       !!
210      jmzme     = 3         !! 1=linear; 2=quadratic; 3=hyperbolic; 4=sigmoid
211      xmzme     = 0.2       !! EKP (26/02/09)
212      xkzme     = 0.75      !! EKP (26/02/09)
213!!
214!! Remineralisation
215      jmd       = 1         !! 0 = T-independent; 1 = T-dependent
216      jsfd      = 0         !! 0 = unchanged seafloor remin.; 1 = accelerated seafloor remin.
217      xmd       = 0.0158    !! T-dependent; 0.05 at 18C (Yool et al., 2010)
218      xmdc      = 0.0127    !! T-dependent; 0.04 at 18C (Yool et al., 2010)
219!!
220!! Stochiometric ratios
221      xthetapn  = 6.625     !! Redfield
222      xthetapd  = 6.625     !! Redfield
223      xthetazmi = 5.625     !! Redfield
224      xthetazme = 5.625     !! Redfield
225!! AXY (26/11/08): implicit detrital carbon change
226      xthetad   = 6.625     !! Redfield
227      xrfn      = 30.0e-6   !! Parekh et al. (2005)
228      xrsn      = 1.0       !! not used here; retained for Lobster
229!!
230!! Iron parameters
231      jiron     = 1         !! iron scavenging submodel switch
232      xfe_mass  = 55.485    !! Fe atomic mass
233      xfe_sol   = 0.00532   !! scaled to get same input as Parekh et al. (2005)
234      xfe_sed   = 0.000228  !! AXY (10/07/12): Moore-ish
235      xLgT      = 1.0       !! from Diat-HadOCC; from Parekh et al. (2005)
236      xk_FeL    = 100.0     !! from Diat-HadOCC
237      xk_sc_Fe  = 1.e-3     !! from Dutkiewicz et al. (2005)
238!!
239!! Fast-sinking detritus parameters
240      jexport   = 1         !! 1 = ballast; 2 = Martin (1987); 3 = Henson (2011)
241      jfdfate   = 0         !! 0 = instant remin.; 1 = transmogrify into slow detritus
242      jrratio   = 2         !! 0 = Dunne (2005); 1 = Ridgwell (2007) surface; 2 = Ridgwell (2007) 3D
243      jocalccd  = 1         !! 0 = default, fixed; 1 = calculated, dynamic
244      xridg_r0  = 0.026     !! = 0.044, Ridgwell et al. (2007); scaled to 60% from preliminary sims
245      xfdfrac1  = 0.333     !! control parameter; no "real" value (AXY; 12/03/09)
246      xfdfrac2  = 1.00      !! control parameter; no "real" value (AXY; 12/03/09)
247      xfdfrac3  = 0.80      !! control parameter; no "real" value (AXY; 12/03/09)
248      xcaco3a   = 0.02      !! Dunne et al. (2005); not an exact value
249      xcaco3b   = 0.10      !! Dunne et al. (2005); not an exact value
250!! AXY (06/01/09): after contacting Dunne, new values are used for masses ...
251!!    xmassc    = 22.1188   !!
252!!    xmassca   = 100.086   !! original values; replaced by Dunne suggestion
253!!    xmasssi   = 78.0988   !!
254      xmassc    = 12.011    !! C atomic mass
255      xmassca   = 100.086   !! CaCO3; based on atomic masses
256      xmasssi   = 60.084    !! SiO2; based on atomic masses
257      xprotca   = 0.070     !! Dunne et al. (2005)
258      xprotsi   = 0.026     !! Dunne et al. (2005)
259      xfastc    = 188.0     !! Dunne et al. (2005)
260      xfastca   = 3500.0    !! Dunne et al. (2005)
261      xfastsi   = 2000.0    !! Dunne et al. (2005)
262!!
263!! Benthos parameters
264      jorgben   = 1         !! 0 = no; 1 = yes
265      jinorgben = 1         !! 0 = no; 1 = yes
266      xsedn     = 0.05      !! default ERSEM-ish value
267      xsedfe    = 0.05      !! default ERSEM-ish value
268      xsedsi    = 0.01      !! default ERSEM-ish value
269      xsedc     = 0.05      !! default ERSEM-ish value
270      xsedca    = 0.01      !! default ERSEM-ish value
271      xburial   = 0.0       !! set to zero
272!!
273!! River parameters
274      jriver_n  = 0         !! 0 = no; 1 = runoff-dependent fluxes; 2= fixed fluxes
275      jriver_si = 0         !! 0 = no; 1 = runoff-dependent fluxes; 2= fixed fluxes
276      jriver_c  = 0         !! 0 = no; 1 = runoff-dependent fluxes; 2= fixed fluxes
277      jriver_alk= 0         !! 0 = no; 1 = runoff-dependent fluxes; 2= fixed fluxes
278      jriver_dep= 5         !! depth to which river nutrients are added
279!!
280!! Miscellaneous
281      xsdiss    = 0.006     !! Mongin et al. (2003); 0.006 is an alternative
282!!
283!! Gravitational sinking
284      vsed      = 3.472e-5  !! EKP (05/03/09)
285!!
286!! Additional parameters
287      jpkb      = 16        !!
288!!
289/
290!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
291!       natroam
292!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
293!
294!!       xthetaphy :  oxygen evolution/consumption by phytoplankton
295!!                    = (32 / 16) + (6.625 * (119 / 106)) = (2 + 7.4375) mol O2 / mol N
296!!       xthetazoo :  oxygen consumption by zooplankton
297!!                    = (32 / 16) + (5.625 * (119 / 106)) = (2 + 6.3149) mol O2 / mol N
298!!       xthetanit :  oxygen consumption by nitrogen remineralisation
299!!                    = (32 / 16) = 2 mol O2 / mol N
300!!       xthetarem :  oxygen consumption by carbon remineralisation
301!!                    = (119 / 106) = 1.1226 mol O2 / mol C
302!!       xo2min    :  oxygen minimum concentration
303!
304&natroam
305   xthetaphy = 9.4375
306   xthetazoo = 8.3149
307   xthetanit = 2.0
308   xthetarem = 1.1226
309   xo2min    = 4.0
310/
311!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
312!       natopt
313!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
314!
315!   xkg0 : green water absorption coefficient   [m-1]
316!   xkr0 : red water absorption coefficient     [m-1]
317!  xkgp  : pigment green absorption coefficient [m-1]
318!  xkrp  : pigment red absorption coefficient   [m-1]
319!  xlg   : green chl exposant
320!  xlr   : red chl exposant
321!  rpig  : chla / (chla+phea) ratio
322!
323&natopt
324   xkg0 = 0.0232
325   xkr0 = 0.225
326   xkgp = 0.074
327   xkrp = 0.037
328   xlg  = 0.629
329   xlr  = 0.674
330   rpig = 0.7
331/
332!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
333!       natdbi used if key_trc_diabio is set
334!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
335!
336!  ctrbio    : array of biological trend numbers
337!  ctrbil    : character string describing each biological trend
338!  ctrbiu    : array of unit numbers used for trend output
339!  nwritebio : frequency of biological model output [timesteps]
340!
341&natdbi
342   ctrbio    = "BIO_01", "BIO_02", "BIO_03", "BIO_04", "BIO_05", "BIO_06", "BIO_07", "BIO_08", "BIO_09", "BIO_10", "BIO_11", "BIO_12", "BIO_13", "BIO_14", "BIO_15", "BIO_16", "BIO_17", "BIO_18", "BIO_19", "BIO_20", "BIO_21", "BIO_22", "BIO_23", "BIO_24", "BIO_25", "BIO_26", "BIO_27", "BIO_28", "BIO_29", "BIO_30", "BIO_31", "BIO_32", "BIO_33", "BIO_34"
343   ctrbil    = "BIO_01", "BIO_02", "BIO_03", "BIO_04", "BIO_05", "BIO_06", "BIO_07", "BIO_08", "BIO_09", "BIO_10", "BIO_11", "BIO_12", "BIO_13", "BIO_14", "BIO_15", "BIO_16", "BIO_17", "BIO_18", "BIO_19", "BIO_20", "BIO_21", "BIO_22", "BIO_23", "BIO_24", "BIO_25", "BIO_26", "BIO_27", "BIO_28", "BIO_29", "BIO_30", "BIO_31", "BIO_32", "BIO_33", "BIO_34"
344   ctrbiu    = "m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m","m"
345   nwritebio = 10
346!/
347!&natopt
348/
349!=================================================
350! nammeddia    Passive tracers additional diagnostics
351!=================================================
352!
353!      ctrc2d  : 2d output field name
354!      ctrc2u  : 2d output field unit
355!      ctrc2l  : 2d output field long name
356!
357!      nwriteadd: frequency of additional arrays outputs(namelist)
358&nammeddia
359/
Note: See TracBrowser for help on using the repository browser.