1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
2 | !! PISCES (key_pisces) reference namelist (see below for key_pisces_reduced) |
---|
3 | !! 1 - air-sea exchange (nampisext) |
---|
4 | !! 2 - biological parameters (nampisbio) |
---|
5 | !! 3 - parameters for nutrient limitations (nampislim) |
---|
6 | !! 4 - parameters for phytoplankton (nampisprod,nampismort) |
---|
7 | !! 5 - parameters for zooplankton (nampismes,nampiszoo) |
---|
8 | !! 6 - parameters for remineralization (nampisrem) |
---|
9 | !! 7 - parameters for calcite chemistry (nampiscal) |
---|
10 | !! 8 - parameters for inputs deposition (nampissed) |
---|
11 | !! 9 - parameters for Kriest parameterization (nampiskrp, nampiskrs) |
---|
12 | !! 10 - additional 2D/3D diagnostics (nampisdia) |
---|
13 | !! 11 - Damping (nampisdmp) |
---|
14 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
15 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
16 | &nampisext ! air-sea exchange |
---|
17 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
18 | ln_co2int = .false. ! read atm pco2 from a file (T) or constant (F) |
---|
19 | atcco2 = 280. ! Constant value atmospheric pCO2 - ln_co2int = F |
---|
20 | clname = 'atcco2.txt' ! Name of atm pCO2 file - ln_co2int = T |
---|
21 | nn_offset = 0 ! Offset model-data start year - ln_co2int = T |
---|
22 | ! ! If your model year is iyy, nn_offset=(years(1)-iyy) |
---|
23 | ! ! then the first atmospheric CO2 record read is at years(1) |
---|
24 | / |
---|
25 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
26 | &nampisatm ! Atmospheric prrssure |
---|
27 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
28 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
29 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
30 | sn_patm = 'presatm' , -1 , 'patm' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
31 | sn_atmco2 = 'presatmco2_reg' , -1 , 'xco2' , .true. , .true. , 'yearly' , 'weights_AIRS_bilinear.nc' , '' , '' |
---|
32 | cn_dir = './' ! root directory for the location of the dynamical files |
---|
33 | ! |
---|
34 | ln_presatm = .false. ! constant atmopsheric pressure (F) or from a file (T) |
---|
35 | ln_presatmco2= .false. ! Read spatialized atm co2 files [ppm] if TRUE |
---|
36 | / |
---|
37 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
38 | &nampisbio ! biological parameters |
---|
39 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
40 | nrdttrc = 2 ! time step frequency for biology |
---|
41 | wsbio = 2. ! POC sinking speed |
---|
42 | xkmort = 1.E-8 ! half saturation constant for mortality |
---|
43 | ferat3 = 10.E-6 ! Fe/C in zooplankton |
---|
44 | no3rat3 = 0.182 ! N/C ratio in zooplankton |
---|
45 | po4rat3 = 0.0094 ! P/C ratio in zooplankton |
---|
46 | wsbio2 = 50. ! Big particles sinking speed |
---|
47 | wsbio2max = 50. ! Big particles maximum sinking speed |
---|
48 | wsbio2scale= 5000. ! Big particles length scale of sinking |
---|
49 | niter1max = 1 ! Maximum number of iterations for POC |
---|
50 | niter2max = 1 ! Maximum number of iterations for GOC |
---|
51 | ! wfep = 0.2 ! FeP sinking speed |
---|
52 | ! ldocp = 1.E-5 ! Phyto ligand production per unit doc |
---|
53 | ! ldocz = 1.E-5 ! Zoo ligand production per unit doc |
---|
54 | ! lthet = 0.5 ! Proportional loss of ligands due to Fe uptake |
---|
55 | / |
---|
56 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
57 | &nampislim ! parameters for nutrient limitations |
---|
58 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
59 | concnno3 = 4e-6 ! Nitrate half saturation of nanophytoplankton |
---|
60 | concpno3 = 1e-6 |
---|
61 | concdno3 = 4E-6 ! Phosphate half saturation for diatoms |
---|
62 | concnnh4 = 1.2E-6 ! NH4 half saturation for phyto |
---|
63 | concpnh4 = 3E-7 |
---|
64 | concdnh4 = 1.2E-6 ! NH4 half saturation for diatoms |
---|
65 | concnpo4 = 4E-6 ! PO4 half saturation for phyto |
---|
66 | concppo4 = 1E-6 |
---|
67 | concdpo4 = 4E-6 ! PO4 half saturation for diatoms |
---|
68 | concnfer = 4E-9 ! Iron half saturation for phyto |
---|
69 | concpfer = 1E-9 |
---|
70 | concdfer = 4E-9 ! Iron half saturation for diatoms |
---|
71 | concbfe = 1.E-11 ! Half-saturation for Fe limitation of Bacteria |
---|
72 | concbnh4 = 1.E-7 ! NH4 half saturation for phyto |
---|
73 | concbno3 = 5.E-7 ! Phosphate half saturation for diatoms |
---|
74 | concbpo4 = 1E-7 ! Phosphate half saturation for bacteria |
---|
75 | xsizedia = 1.E-6 ! Minimum size criteria for diatoms |
---|
76 | xsizephy = 1.E-6 ! Minimum size criteria for phyto |
---|
77 | xsizepic = 1.E-6 |
---|
78 | xsizern = 3.0 ! Size ratio for nanophytoplankton |
---|
79 | xsizerp = 1.0 |
---|
80 | xsizerd = 3.0 ! Size ratio for diatoms |
---|
81 | xksi1 = 8.E-6 ! half saturation constant for Si uptake |
---|
82 | xksi2 = 20E-6 ! half saturation constant for Si/C |
---|
83 | xkdoc = 417.E-6 ! half-saturation constant of DOC remineralization |
---|
84 | caco3r = 0.35 ! mean rain ratio |
---|
85 | / |
---|
86 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
87 | &nampisquota ! parameters for nutrient limitations |
---|
88 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
89 | qfnopt = 7.E-6 ! Optimal quota of phyto |
---|
90 | qfpopt = 7.E-6 |
---|
91 | qfdopt = 7.E-6 ! Optimal quota of diatoms |
---|
92 | qnnmin = 0.29 |
---|
93 | qnnmax = 1.39 |
---|
94 | qpnmin = 0.24 |
---|
95 | qpnmax = 1.06 |
---|
96 | qnpmin = 0.42 |
---|
97 | qnpmax = 1.39 |
---|
98 | qppmin = 0.19 |
---|
99 | qppmax = 0.7 |
---|
100 | qndmin = 0.25 |
---|
101 | qndmax = 1.39 |
---|
102 | qpdmin = 0.24 |
---|
103 | qpdmax = 1.32 |
---|
104 | qfnmax = 40E-6 |
---|
105 | qfpmax = 40E-6 |
---|
106 | qfdmax = 40E-6 |
---|
107 | / |
---|
108 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
109 | &nampisopt ! parameters for optics |
---|
110 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
111 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
112 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
113 | sn_par = 'par.orca' , 24 , 'fr_par' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
114 | cn_dir = './' ! root directory for the location of the dynamical files |
---|
115 | ln_varpar = .true. ! boolean for PAR variable |
---|
116 | parlux = 0.43 ! Fraction of shortwave as PAR |
---|
117 | / |
---|
118 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
119 | &nampisprod ! parameters for phytoplankton growth |
---|
120 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
121 | pislopen = 4. ! P-I slope |
---|
122 | pislopep = 4. ! P-I slope for picophytoplankton |
---|
123 | pisloped = 4. ! P-I slope for diatoms |
---|
124 | xadap = 0. ! Adaptation factor to low light |
---|
125 | excretn = 0.05 ! excretion ratio of phytoplankton |
---|
126 | excretp = 0.05 ! excretion ratio of picophytoplankton |
---|
127 | excretd = 0.05 ! excretion ratio of diatoms |
---|
128 | bresp = 0.02 ! Basal respiration rate |
---|
129 | thetannm = 0.3 ! Maximum Chl/N in nanophytoplankton |
---|
130 | thetanpm = 0.3 ! Maximum Chl/N in picophytoplankton |
---|
131 | thetandm = 0.35 ! Maximum Chl/N in diatoms |
---|
132 | chlcmin = 0.004 ! Minimum Chl/c in phytoplankton |
---|
133 | grosip = 0.39 ! mean Si/C ratio |
---|
134 | / |
---|
135 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
136 | &nampismort ! parameters for phytoplankton sinks |
---|
137 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
138 | wchl = 0.01 ! quadratic mortality of phytoplankton |
---|
139 | wchlp = 0.01 |
---|
140 | wchld = 0.01 ! maximum quadratic mortality of diatoms |
---|
141 | wchldm = 0.02 ! maximum quadratic mortality of diatoms |
---|
142 | mprat = 0.01 ! phytoplankton mortality rate |
---|
143 | mpratp = 0.01 |
---|
144 | mprat2 = 0.01 ! Diatoms mortality rate |
---|
145 | / |
---|
146 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
147 | &nampismes ! parameters for mesozooplankton |
---|
148 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
149 | part2 = 0.75 ! part of calcite not dissolved in mesozoo guts |
---|
150 | grazrat2 = 0.8 ! maximal mesozoo grazing rate |
---|
151 | bmetexc2 = .false. ! Metabolic use of excess carbon |
---|
152 | resrat2 = 0.005 ! exsudation rate of mesozooplankton |
---|
153 | mzrat2 = 0.02 ! mesozooplankton mortality rate |
---|
154 | xpref2d = 1. ! zoo preference for phyto |
---|
155 | xpref2p = 0.8 ! zoo preference for POC |
---|
156 | xpref2z = 1. ! zoo preference for zoo |
---|
157 | xpref2m = 0.2 ! meso preference for zoo |
---|
158 | xpref2c = 0.2 ! zoo preference for poc |
---|
159 | xthresh2zoo = 1E-8 ! zoo feeding threshold for mesozooplankton |
---|
160 | xthresh2dia = 1E-8 ! diatoms feeding threshold for mesozooplankton |
---|
161 | xthresh2phy = 1E-8 ! nanophyto feeding threshold for mesozooplankton |
---|
162 | xthresh2mes = 1E-8 ! meso feeding threshold for mesozooplankton |
---|
163 | xthresh2poc = 1E-8 ! poc feeding threshold for mesozooplankton |
---|
164 | xthresh2 = 1E-8 ! Food threshold for grazing |
---|
165 | xkgraz2 = 20.E-6 ! half sturation constant for meso grazing |
---|
166 | epsher2 = 0.5 ! Efficicency of Mesozoo growth |
---|
167 | ssigma2 = 0.5 ! Fraction excreted as semi-labile DOM |
---|
168 | srespir2 = 0.2 ! Active respiration |
---|
169 | unass2c = 0.3 ! non assimilated fraction of P by mesozoo |
---|
170 | unass2n = 0.3 ! non assimilated fraction of N by mesozoo |
---|
171 | unass2p = 0.3 ! non assimilated fraction of P by mesozoo |
---|
172 | grazflux = 3.e3 ! flux-feeding rate |
---|
173 | / |
---|
174 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
175 | &nampiszoo ! parameters for microzooplankton |
---|
176 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
177 | part = 0.5 ! part of calcite not dissolved in microzoo gutsa |
---|
178 | grazrat = 2.5 ! maximal zoo grazing rate |
---|
179 | bmetexc = .false. ! Metabolic use of excess carbon |
---|
180 | resrat = 0.03 ! exsudation rate of zooplankton |
---|
181 | mzrat = 0.005 ! zooplankton mortality rate |
---|
182 | xprefc = 0.1 ! Microzoo preference for POM |
---|
183 | xprefn = 1. ! Microzoo preference for Nanophyto |
---|
184 | xprefp = 1.4 |
---|
185 | xprefd = 0.8 ! Microzoo preference for Diatoms |
---|
186 | xprefz = 0.1 ! Microzoo preference for microzooplankton |
---|
187 | xthreshdia = 1.E-8 ! Diatoms feeding threshold for microzooplankton |
---|
188 | xthreshphy = 1.E-8 ! Nanophyto feeding threshold for microzooplankton |
---|
189 | xthreshpic = 1.E-8 |
---|
190 | xthreshzoo = 1.E-8 ! Nanophyto feeding threshold for microzooplankton |
---|
191 | xthreshpoc = 1.E-8 ! POC feeding threshold for microzooplankton |
---|
192 | xthresh = 1.E-8 ! Food threshold for feeding |
---|
193 | xkgraz = 20.E-6 ! half sturation constant for grazing |
---|
194 | epsher = 0.5 ! Efficiency of microzoo growth |
---|
195 | ssigma = 0.5 ! Fraction excreted as semi-labile DOM |
---|
196 | srespir = 0.2 ! Active respiration |
---|
197 | unassc = 0.3 ! non assimilated fraction of C by zoo |
---|
198 | unassn = 0.3 ! non assimilated fraction of C by zoo |
---|
199 | unassp = 0.3 ! non assimilated fraction of C by zoo |
---|
200 | / |
---|
201 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
202 | &nampisfer ! parameters for iron chemistry |
---|
203 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
204 | ln_fechem = .false. ! complex iron chemistry ( T/F ) |
---|
205 | ln_ligvar = .false. ! variable ligand concentration |
---|
206 | ln_fecolloid = .false. ! variable colloidal fraction |
---|
207 | xlam1 = 0.005 ! scavenging rate of Iron |
---|
208 | xlamdust = 150.0 ! Scavenging rate of dust |
---|
209 | ligand = 0.7E-9 ! Ligands concentration |
---|
210 | kfep = 0. ! Nanoparticle formation rate constant |
---|
211 | / |
---|
212 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
213 | &nampisrem ! parameters for remineralization |
---|
214 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
215 | xremikc = 0.25 ! remineralization rate of DOC |
---|
216 | xremikn = 0.35 ! remineralization rate of DON |
---|
217 | xremikp = 0.4 ! remineralization rate of DOP |
---|
218 | nitrif = 0.05 ! NH4 nitrification rate |
---|
219 | xsirem = 0.003 ! remineralization rate of Si |
---|
220 | xsiremlab = 0.03 ! fast remineralization rate of Si |
---|
221 | xsilab = 0.5 ! Fraction of labile biogenic silica |
---|
222 | oxymin = 1.E-6 ! Half-saturation constant for anoxia |
---|
223 | oxymin2 = 6.E-6 ! Minimum O2 concentration for oxic remin. |
---|
224 | feratb = 10E-6 ! Bacterial Fe/C ratio |
---|
225 | xkferb = 3E-10 ! Half-saturation constant for bact. Fe/C |
---|
226 | / |
---|
227 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
228 | &nampispoc ! parameters for organic particles |
---|
229 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
230 | xremipc = 0.03 ! remineralisation rate of POC |
---|
231 | xremipn = 0.035 ! remineralisation rate of PON |
---|
232 | xremipp = 0.04 ! remineralisation rate of POP |
---|
233 | jcpoc = 15 ! Number of lability classes |
---|
234 | rshape = 1.0 ! Shape of the gamma function |
---|
235 | / |
---|
236 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
237 | &nampiscal ! parameters for Calcite chemistry |
---|
238 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
239 | kdca = 6. ! calcite dissolution rate constant (1/time) |
---|
240 | nca = 1. ! order of dissolution reaction (dimensionless) |
---|
241 | / |
---|
242 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
243 | &nampissbc ! parameters for inputs deposition |
---|
244 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
245 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
246 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
247 | sn_dust = 'dust.orca' , -1 , 'dust' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
248 | sn_solub = 'solubility.orca' , -12 , 'solubility1' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
249 | sn_riverdic = 'river.orca' , 120 , 'riverdic' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
250 | sn_riverdoc = 'river.orca' , 120 , 'riverdoc' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
251 | sn_riverdin = 'river.orca' , 120 , 'riverdin' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
252 | sn_riverdon = 'river.orca' , 120 , 'riverdon' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
253 | sn_riverdip = 'river.orca' , 120 , 'riverdip' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
254 | sn_riverdop = 'river.orca' , 120 , 'riverdop' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
255 | sn_riverdsi = 'river.orca' , 120 , 'riverdsi' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
256 | sn_ndepo = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
257 | sn_ironsed = 'bathy.orca' , -12 , 'bathy' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
258 | sn_hydrofe = 'hydrofe.orca' , -12 , 'epsdb' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
259 | ! |
---|
260 | cn_dir = './' ! root directory for the location of the dynamical files |
---|
261 | ln_dust = .true. ! boolean for dust input from the atmosphere |
---|
262 | ln_solub = .true. ! boolean for variable solubility of atm. Iron |
---|
263 | ln_river = .true. ! boolean for river input of nutrients |
---|
264 | ln_ndepo = .true. ! boolean for atmospheric deposition of N |
---|
265 | ln_ironsed = .true. ! boolean for Fe input from sediments |
---|
266 | ln_ironice = .true. ! boolean for Fe input from sea ice |
---|
267 | ln_hydrofe = .false. ! boolean for from hydrothermal vents |
---|
268 | sedfeinput = 2.e-9 ! Coastal release of Iron |
---|
269 | dustsolub = 0.02 ! Solubility of the dusta |
---|
270 | mfrac = 0.035 ! Fe mineral fraction of dust |
---|
271 | wdust = 2.0 ! Dust sinking speed |
---|
272 | icefeinput = 15.e-9 ! Iron concentration in sea ice |
---|
273 | nitrfix = 1.e-7 ! Nitrogen fixation rate |
---|
274 | diazolight = 50. ! Diazotrophs sensitivity to light (W/m2) |
---|
275 | concfediaz = 1.e-10 ! Diazotrophs half-saturation Cste for Iron |
---|
276 | hratio = 1.e+7 ! Fe to 3He ratio assumed for vent iron supply |
---|
277 | ! fep_rats = 1. ! Fep/Fer ratio from sed sources |
---|
278 | ! fep_rath = 1. ! Fep/Fer ratio from sed hydro sources |
---|
279 | / |
---|
280 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
281 | &nampislig ! Namelist parameters for ligands, nampislig |
---|
282 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
283 | rfep = 0.001 ! Dissolution rate of FeP |
---|
284 | rlgw = 1. ! Lifetime (years) of weak ligands |
---|
285 | rlig = 1.E-4 ! Remin ligand production per unit C |
---|
286 | prlgw = 1.E-4 ! Photolysis of weak ligand |
---|
287 | rlgs = 1000. ! Lifetime (years) of strong ligands |
---|
288 | / |
---|
289 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
290 | &nampisice ! Prescribed sea ice tracers |
---|
291 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
292 | ! constant ocean tracer concentrations are defined in trcice_pisces.F90 (Global, Arctic, Antarctic and Baltic) |
---|
293 | ! trc_ice_ratio * betw 0 and 1: prescribed ice/ocean tracer concentration ratio |
---|
294 | ! * -1 => the ice-ocean tracer concentration ratio follows the |
---|
295 | ! ice-ocean salinity ratio |
---|
296 | ! * -2 => tracer concentration in sea ice is prescribed and |
---|
297 | ! trc_ice_prescr is used |
---|
298 | ! trc_ice_prescr * prescribed tracer concentration. used only if |
---|
299 | ! trc_ice_ratio = -2. equals -99 if not used. |
---|
300 | ! cn_trc_o * 'GL' use global ocean values making the Baltic distinction only |
---|
301 | ! 'AA' use specific Arctic/Antarctic/Baltic values |
---|
302 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
303 | ! sn_tri_ ! trc_ice_ratio ! trc_ice_prescr ! cn_trc_o |
---|
304 | sn_tri_dic = -1., -99., 'AA' |
---|
305 | sn_tri_doc = 0., -99., 'AA' |
---|
306 | sn_tri_tal = -1., -99., 'AA' |
---|
307 | sn_tri_oxy = -1., -99., 'AA' |
---|
308 | sn_tri_cal = 0., -99., 'AA' |
---|
309 | sn_tri_po4 = -1., -99., 'AA' |
---|
310 | sn_tri_poc = 0., -99., 'AA' |
---|
311 | sn_tri_goc = 0., -99., 'AA' |
---|
312 | sn_tri_bfe = 0., -99., 'AA' |
---|
313 | sn_tri_num = 0., -99., 'AA' |
---|
314 | sn_tri_sil = -1., -99., 'AA' |
---|
315 | sn_tri_dsi = 0., -99., 'AA' |
---|
316 | sn_tri_gsi = 0., -99., 'AA' |
---|
317 | sn_tri_phy = 0., -99., 'AA' |
---|
318 | sn_tri_dia = 0., -99., 'AA' |
---|
319 | sn_tri_zoo = 0., -99., 'AA' |
---|
320 | sn_tri_mes = 0., -99., 'AA' |
---|
321 | sn_tri_fer = -2., 15E-9, 'AA' |
---|
322 | sn_tri_sfe = 0., -99., 'AA' |
---|
323 | sn_tri_dfe = 0., -99., 'AA' |
---|
324 | sn_tri_nfe = 0., -99., 'AA' |
---|
325 | sn_tri_nch = 0., -99., 'AA' |
---|
326 | sn_tri_dch = 0., -99., 'AA' |
---|
327 | sn_tri_no3 = -1., -99., 'AA' |
---|
328 | sn_tri_nh4 = 1., -99., 'AA' |
---|
329 | / |
---|
330 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
331 | &nampiskrp ! Kriest parameterization : parameters "key_kriest" |
---|
332 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
333 | xkr_eta = 1.17 ! Sinking exponent |
---|
334 | xkr_zeta = 2.28 ! N content exponent |
---|
335 | xkr_ncontent = 5.7E-6 ! N content factor |
---|
336 | xkr_mass_min = 0.0002 ! Minimum mass for Aggregates |
---|
337 | xkr_mass_max = 1. ! Maximum mass for Aggregates |
---|
338 | / |
---|
339 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
340 | &nampiskrs ! Kriest parameterization : size classes "key_kriest" |
---|
341 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
342 | xkr_sfact = 942. ! Sinking factor |
---|
343 | xkr_stick = 0.5 ! Stickiness |
---|
344 | xkr_nnano = 2.337 ! Nbr of cell in nano size class |
---|
345 | xkr_ndiat = 3.718 ! Nbr of cell in diatoms size class |
---|
346 | xkr_nmeso = 7.147 ! Nbr of cell in mesozoo size class |
---|
347 | xkr_naggr = 9.877 ! Nbr of cell in aggregates size class |
---|
348 | / |
---|
349 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
350 | &nampisdia ! additional 2D/3D tracers diagnostics |
---|
351 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
352 | ! ! name ! title of the field ! units ! |
---|
353 | ! ! ! ! ! |
---|
354 | pisdia2d(1) = 'Cflx ' , 'DIC flux ', 'molC/m2/s ' |
---|
355 | pisdia2d(2) = 'Oflx ' , 'Oxygen flux ', 'molC/m2/s ' |
---|
356 | pisdia2d(3) = 'Kg ' , 'Gas transfer ', 'mol/m2/s/uatm' |
---|
357 | pisdia2d(4) = 'Delc ' , 'Delta CO2 ', 'uatm ' |
---|
358 | pisdia2d(5) = 'PMO ' , 'POC export ', 'molC/m2/s ' |
---|
359 | pisdia2d(6) = 'PMO2 ' , 'GOC export ', 'molC/m2/s ' |
---|
360 | pisdia2d(7) = 'ExpFe1 ' , 'Nano iron export ', 'molFe/m2/s ' |
---|
361 | pisdia2d(8) = 'ExpFe2 ' , 'Diatoms iron export ', 'molFe/m2/s ' |
---|
362 | pisdia2d(9) = 'ExpSi ' , 'Silicate export ', 'molSi/m2/s ' |
---|
363 | pisdia2d(10) = 'ExpCaCO3 ' , 'Calcite export ', 'molC/m2/s ' |
---|
364 | pisdia2d(11) = 'heup ' , 'euphotic layer depth ', 'm ' |
---|
365 | pisdia2d(12) = 'Fedep ' , 'Iron dep ', 'molFe/m2/s ' |
---|
366 | pisdia2d(13) = 'Nfix ' , 'Nitrogen Fixation ', 'molN/m2/s ' |
---|
367 | pisdia3d(1) = 'PH ' , 'PH ', '- ' |
---|
368 | pisdia3d(2) = 'CO3 ' , 'Bicarbonates ', 'mol/l ' |
---|
369 | pisdia3d(3) = 'CO3sat ' , 'CO3 saturation ', 'mol/l ' |
---|
370 | pisdia3d(4) = 'PAR ' , 'light penetration ', 'W/m2 ' |
---|
371 | pisdia3d(5) = 'PPPHY ' , 'Primary production of nanophyto ', 'molC/m3/s ' |
---|
372 | pisdia3d(6) = 'PPPHY2 ' , 'Primary production of diatoms ', 'molC/m3/s ' |
---|
373 | pisdia3d(7) = 'PPNEWN ' , 'New Primary production of nano ', 'molC/m3/s ' |
---|
374 | pisdia3d(8) = 'PPNEWD ' , 'New Primary production of diat ', 'molC/m3/s ' |
---|
375 | pisdia3d(9) = 'PBSi ' , 'Primary production of Si diatoms ', 'molSi/m3/s ' |
---|
376 | pisdia3d(10) = 'PFeN ' , 'Primary production of nano iron ', 'molFe/m3/s ' |
---|
377 | pisdia3d(11) = 'PFeD ' , 'Primary production of diatoms iron', 'molFe/m3/s ' |
---|
378 | / |
---|
379 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
380 | &nampisdmp ! Damping |
---|
381 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
382 | ln_pisdmp = .true. ! Relaxation fo some tracers to a mean value |
---|
383 | nn_pisdmp = 2190 ! Frequency of Relaxation |
---|
384 | / |
---|
385 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
386 | &nampismass ! Mass conservation |
---|
387 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
388 | ln_check_mass = .false. ! Check mass conservation |
---|
389 | / |
---|
390 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
391 | !! PISCES reduced (key_pisces_reduced, ex LOBSTER) : namelists |
---|
392 | !! 1 - biological parameters for phytoplankton (namlobphy) |
---|
393 | !! 2 - biological parameters for nutrients (namlobnut) |
---|
394 | !! 3 - biological parameters for zooplankton (namlobzoo) |
---|
395 | !! 4 - biological parameters for detritus (namlobdet) |
---|
396 | !! 5 - biological parameters for DOM (namlobdom) |
---|
397 | !! 6 - parameters from aphotic layers to sediment (namlobsed) |
---|
398 | !! 7 - general coefficients (namlobrat) |
---|
399 | !! 8 - optical parameters (namlobopt) |
---|
400 | |
---|
401 | !! 10 - biological diagnostics trends (namlobdbi) |
---|
402 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
403 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
404 | &namlobphy ! biological parameters for phytoplankton |
---|
405 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
406 | tmumax = 1.21e-5 ! maximal phytoplankton growth rate [s-1] |
---|
407 | rgamma = 0.05 ! phytoplankton exudation fraction [%] |
---|
408 | fphylab = 0.75 ! NH4 fraction of phytoplankton exsudation |
---|
409 | tmminp = 5.8e-7 ! minimal phytoplancton mortality rate [0.05/86400 s-1=20 days] |
---|
410 | aki = 33. ! light photosynthesis half saturation constant[W/m2] |
---|
411 | / |
---|
412 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
413 | &namlobnut ! biological parameters for nutrients |
---|
414 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
415 | akno3 = 0.7 ! nitrate limitation half-saturation value [mmol/m3] |
---|
416 | aknh4 = 0.001 ! ammonium limitation half-saturation value [mmol/m3] |
---|
417 | taunn = 5.80e-7 ! nitrification rate [s-1] |
---|
418 | psinut = 3. ! inhibition of nitrate uptake by ammonium |
---|
419 | / |
---|
420 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
421 | &namlobzoo ! biological parameters for zooplankton |
---|
422 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
423 | rppz = 0.8 ! zooplankton nominal preference for phytoplancton food [%] |
---|
424 | taus = 9.26E-6 ! specific zooplankton maximal grazing rate [s-1] |
---|
425 | ! ! 0.75/86400 s-1=8.680555E-6 1/86400 = 1.15e-5 |
---|
426 | aks = 1. ! half-saturation constant for total zooplankton grazing [mmolN.m-3] |
---|
427 | rpnaz = 0.3 ! non-assimilated phytoplankton by zooplancton [%] |
---|
428 | rdnaz = 0.3 ! non-assimilated detritus by zooplankton [%] |
---|
429 | tauzn = 8.1e-7 ! zooplancton specific excretion rate [0.1/86400 s-1=10 days] |
---|
430 | fzoolab = 0.5 ! NH4 fraction of zooplankton excretion |
---|
431 | fdbod = 0.5 ! zooplankton mortality fraction that goes to detritus |
---|
432 | tmminz = 2.31e-6 ! minimal zooplankton mortality rate [(mmolN/m3)-1 d-1] |
---|
433 | / |
---|
434 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
435 | &namlobdet ! biological parameters for detritus |
---|
436 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
437 | taudn = 5.80e-7 ! detritus breakdown rate [0.1/86400 s-1=10 days] |
---|
438 | fdetlab = 0. ! NH4 fraction of detritus dissolution |
---|
439 | / |
---|
440 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
441 | &namlobdom ! biological parameters for DOM |
---|
442 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
443 | taudomn = 6.43e-8 ! DOM breakdown rate [s-1] |
---|
444 | ! ! slow remineralization rate of semi-labile dom to nh4 (1 month) |
---|
445 | / |
---|
446 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
447 | &namlobsed ! parameters from aphotic layers to sediment |
---|
448 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
449 | sedlam = 3.86e-7 ! time coefficient of POC remineralization in sediments [s-1] |
---|
450 | sedlostpoc = 0. ! mass of POC lost in sediments |
---|
451 | vsed = 3.47e-5 ! detritus sedimentation speed [m/s] |
---|
452 | xhr = -0.858 ! coeff for martin''s remineralisation profile |
---|
453 | / |
---|
454 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
455 | &namlobrat ! general coefficients |
---|
456 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
457 | rcchl = 60. ! Carbone/Chlorophyl ratio [mgC.mgChla-1] |
---|
458 | redf = 6.56 ! redfield ratio (C:N) for phyto |
---|
459 | reddom = 6.56 ! redfield ratio (C:N) for DOM |
---|
460 | / |
---|
461 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
462 | &namlobopt ! optical parameters |
---|
463 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
464 | xkg0 = 0.0232 ! green absorption coefficient of water |
---|
465 | xkr0 = 0.225 ! red absorption coefficent of water |
---|
466 | xkgp = 0.074 ! green absorption coefficient of chl |
---|
467 | xkrp = 0.037 ! red absorption coefficient of chl |
---|
468 | xlg = 0.674 ! green chl exposant for absorption |
---|
469 | xlr = 0.629 ! red chl exposant for absorption |
---|
470 | rpig = 0.7 ! chla/chla+pheo ratio |
---|
471 | / |
---|
472 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
473 | &nampisdbi ! biological diagnostics trends |
---|
474 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
475 | ! ! 2D bio diagnostics units : mmole/m2/s ("key_trdmld_trc") |
---|
476 | ! ! name ! title of the field ! units ! |
---|
477 | pisdiabio(1) = 'NO3PHY' , 'Flux from NO3 to PHY ', 'mmole/m3/s' |
---|
478 | pisdiabio(2) = 'NH4PHY' , 'Flux from NH4 to PHY ', 'mmole/m3/s' |
---|
479 | pisdiabio(3) = 'PHYNH4' , 'Flux from PHY to NH4 ', 'mmole/m3/s' |
---|
480 | pisdiabio(4) = 'PHYDOM' , 'Flux from PHY to DOM ', 'mmole/m3/s' |
---|
481 | pisdiabio(5) = 'PHYZOO' , 'Flux from PHY to ZOO ', 'mmole/m3/s' |
---|
482 | pisdiabio(6) = 'PHYDET' , 'Flux from PHY to DET ', 'mmole/m3/s' |
---|
483 | pisdiabio(7) = 'DETZOO' , 'Flux from DET to ZOO ', 'mmole/m3/s' |
---|
484 | pisdiabio(8) = 'DETSED' , 'Flux from DET to SED ', 'mmole/m3/s' |
---|
485 | pisdiabio(9) = 'ZOODET' , 'Flux from ZOO to DET ', 'mmole/m3/s' |
---|
486 | pisdiabio(10) = 'ZOOBOD' , 'Zooplankton closure ', 'mmole/m3/s' |
---|
487 | pisdiabio(11) = 'ZOONH4' , 'Flux from ZOO to NH4 ', 'mmole/m3/s' |
---|
488 | pisdiabio(12) = 'ZOODOM' , 'Flux from ZOO to DOM ', 'mmole/m3/s' |
---|
489 | pisdiabio(13) = 'NH4NO3' , 'Flux from NH4 to NO3 ', 'mmole/m3/s' |
---|
490 | pisdiabio(14) = 'DOMNH4' , 'Flux from DOM to NH4 ', 'mmole/m3/s' |
---|
491 | pisdiabio(15) = 'DETNH4' , 'Flux from DET to NH4 ', 'mmole/m3/s' |
---|
492 | pisdiabio(16) = 'DETDOM' , 'Flux from DET to DOM ', 'mmole/m3/s' |
---|
493 | pisdiabio(17) = 'SEDNO3' , 'NO3 remineralization from SED', 'mmole/m3/s' |
---|
494 | / |
---|