3*************************************************** hernan g. arango ***
4** copyright(c) 2002-2025 the roms group **
5** licensed under a mit/x style license **
6** see license_roms.md **
7************************************************************************
9** defines ecosim bio-optical model input parameters in
output **
10** netcdf files. it is included in routine
"def_info.F". **
12************************************************************************
19 vinfo( 2)=
'maximum number of iterations to achieve convergence'
20 status=def_var(ng, model, piofile, piovar, pio_int, &
21 & 1, (/0/), aval, vinfo, ncname, &
22 & setparaccess = .false.)
23 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
25 vinfo( 1)=
'RtUVR_flag'
26 vinfo( 2)=
'switch to calculate CDOC UV photolysis.'
29 status=def_var(ng, model, piofile, piovar, pio_int, &
30 & 1, (/0/), aval, vinfo, ncname, &
31 & setparaccess = .false.)
32 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
35 vinfo( 2)=
'switch to calculate temperature based N fixation'
38 status=def_var(ng, model, piofile, piovar, pio_int, &
39 & 1, (/0/), aval, vinfo, ncname, &
40 & setparaccess = .false.)
41 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
43 vinfo( 1)=
'Regen_flag'
44 vinfo( 2)=
'switch to calculate fecal matter regeneration'
47 status=def_var(ng, model, piofile, piovar, pio_int, &
48 & 1, (/0/), aval, vinfo, ncname, &
49 & setparaccess = .false.)
50 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
53 vinfo( 2)=
'center wavelength of spectral irradiance bands'
55 status=def_var(ng, model, piofile, piovar, pio_type, &
56 & 1, (/lightdim/), aval, vinfo, ncname, &
57 & setparaccess = .false.)
58 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
61 vinfo( 2)=
'half-saturation for phytoplankton NO3 uptake'
62 vinfo( 3)=
'micromole_NO3 liter-1'
63 status=def_var(ng, model, piofile, piovar, pio_type, &
64 & 1, (/phydim/), aval, vinfo, ncname, &
65 & setparaccess = .false.)
66 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
69 vinfo( 2)=
'half-saturation for phytoplankton NH4 uptake'
70 vinfo( 3)=
'micromole_NH4 liter-1'
71 status=def_var(ng, model, piofile, piovar, pio_type, &
72 & 1, (/phydim/), aval, vinfo, ncname, &
73 & setparaccess = .false.)
74 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
77 vinfo( 2)=
'half-saturation for phytoplankton SiO uptake'
78 vinfo( 3)=
'micromole_SiO liter-1'
79 status=def_var(ng, model, piofile, piovar, pio_type, &
80 & 1, (/phydim/), aval, vinfo, ncname, &
81 & setparaccess = .false.)
82 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
85 vinfo( 2)=
'half-saturation for phytoplankton PO4 uptake'
86 vinfo( 3)=
'micromole_PO4 liter-1'
87 status=def_var(ng, model, piofile, piovar, pio_type, &
88 & 1, (/phydim/), aval, vinfo, ncname, &
89 & setparaccess = .false.)
90 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
93 vinfo( 2)=
'half-saturation for phytoplankton Fe uptake'
94 vinfo( 3)=
'micromole_Fe liter-1'
95 status=def_var(ng, model, piofile, piovar, pio_type, &
96 & 1, (/phydim/), aval, vinfo, ncname, &
97 & setparaccess = .false.)
98 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
100 vinfo( 1)=
'GtALG_max'
101 vinfo( 2)=
'maximum phytoplankton 24 hour growth rate'
103 status=def_var(ng, model, piofile, piovar, pio_type, &
104 & 1, (/phydim/), aval, vinfo, ncname, &
105 & setparaccess = .false.)
106 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
110 &
'phytoplankton temperature base for exponential response'
112 status=def_var(ng, model, piofile, piovar, pio_type, &
113 & 1, (/phydim/), aval, vinfo, ncname, &
114 & setparaccess = .false.)
115 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
118 vinfo( 2)=
'phytoplankton exponential temperature factor'
119 vinfo( 3)=
'Celsius-1'
120 status=def_var(ng, model, piofile, piovar, pio_type, &
121 & 1, (/phydim/), aval, vinfo, ncname, &
122 & setparaccess = .false.)
123 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
126 vinfo( 2)=
'nitrate uptake inhibition for NH4'
127 vinfo( 3)=
'micromole-1'
128 status=def_var(ng, model, piofile, piovar, pio_type, &
129 & 1, (/phydim/), aval, vinfo, ncname, &
130 & setparaccess = .false.)
131 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
133 vinfo( 1)=
'maxC2nALG'
134 vinfo( 2)=
'maximum phytoplankton C:N ratio'
135 vinfo( 3)=
'micromole_C micromole_N-1'
136 status=def_var(ng, model, piofile, piovar, pio_type, &
137 & 1, (/phydim/), aval, vinfo, ncname, &
138 & setparaccess = .false.)
139 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
141 vinfo( 1)=
'minC2nALG'
142 vinfo( 2)=
'balanced phytoplankton C:N ratio'
143 vinfo( 3)=
'micromole_C micromole_N-1'
144 status=def_var(ng, model, piofile, piovar, pio_type, &
145 & 1, (/phydim/), aval, vinfo, ncname, &
146 & setparaccess = .false.)
147 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
149 vinfo( 1)=
'C2nALGminABS'
150 vinfo( 2)=
'absolute minimum phytoplankton C:N ratio'
151 vinfo( 3)=
'micromole_C micromole_N-1'
152 status=def_var(ng, model, piofile, piovar, pio_type, &
153 & 1, (/phydim/), aval, vinfo, ncname, &
154 & setparaccess = .false.)
155 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
157 vinfo( 1)=
'maxC2SiALG'
158 vinfo( 2)=
'maximum phytoplankton C:Si ratio'
159 vinfo( 3)=
'micromole_C micromole_Si-1'
160 status=def_var(ng, model, piofile, piovar, pio_type, &
161 & 1, (/phydim/), aval, vinfo, ncname, &
162 & setparaccess = .false.)
163 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
165 vinfo( 1)=
'minC2SiALG'
166 vinfo( 2)=
'balanced phytoplankton C:Si ratio'
167 vinfo( 3)=
'micromole_C micromole_Si-1'
168 status=def_var(ng, model, piofile, piovar, pio_type, &
169 & 1, (/phydim/), aval, vinfo, ncname, &
170 & setparaccess = .false.)
171 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
173 vinfo( 1)=
'C2SiALGminABS'
174 vinfo( 2)=
'absolute minimum phytoplankton C:Si ratio'
175 vinfo( 3)=
'micromole_C micromole_Si-1'
176 status=def_var(ng, model, piofile, piovar, pio_type, &
177 & 1, (/phydim/), aval, vinfo, ncname, &
178 & setparaccess = .false.)
179 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
181 vinfo( 1)=
'maxC2pALG'
182 vinfo( 2)=
'maximum phytoplankton C:P ratio'
183 vinfo( 3)=
'micromole_C micromole_P-1'
184 status=def_var(ng, model, piofile, piovar, pio_type, &
185 & 1, (/phydim/), aval, vinfo, ncname, &
186 & setparaccess = .false.)
187 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
189 vinfo( 1)=
'minC2pALG'
190 vinfo( 2)=
'balanced phytoplankton C:P ratio'
191 vinfo( 3)=
'micromole_C micromole_P-1'
192 status=def_var(ng, model, piofile, piovar, pio_type, &
193 & 1, (/phydim/), aval, vinfo, ncname, &
194 & setparaccess = .false.)
195 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
197 vinfo( 1)=
'C2pALGminABS'
198 vinfo( 2)=
'absolute minimum phytoplankton C:P ratio'
199 vinfo( 3)=
'micromole_C micromole_P-1'
200 status=def_var(ng, model, piofile, piovar, pio_type, &
201 & 1, (/phydim/), aval, vinfo, ncname, &
202 & setparaccess = .false.)
203 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
205 vinfo( 1)=
'maxC2FeALG'
206 vinfo( 2)=
'maximum phytoplankton C:Fe ratio'
207 vinfo( 3)=
'micromole_C micromole_Fe-1'
208 status=def_var(ng, model, piofile, piovar, pio_type, &
209 & 1, (/phydim/), aval, vinfo, ncname, &
210 & setparaccess = .false.)
211 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
213 vinfo( 1)=
'minC2FeALG'
214 vinfo( 2)=
'balanced phytoplankton C:Fe ratio'
215 vinfo( 3)=
'micromole_C micromole_Fe-1'
216 status=def_var(ng, model, piofile, piovar, pio_type, &
217 & 1, (/phydim/), aval, vinfo, ncname, &
218 & setparaccess = .false.)
219 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
221 vinfo( 1)=
'C2FeALGminABS'
222 vinfo( 2)=
'absolute minimum phytoplankton C:Fe ratio'
223 vinfo( 3)=
'micromole_C micromole_Fe-1'
224 status=def_var(ng, model, piofile, piovar, pio_type, &
225 & 1, (/phydim/), aval, vinfo, ncname, &
226 & setparaccess = .false.)
227 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
230 vinfo( 2)=
'maximum quantum yield'
231 vinfo( 3)=
'micromole_C micromole_quanta-1'
232 status=def_var(ng, model, piofile, piovar, pio_type, &
233 & 1, (/phydim/), aval, vinfo, ncname, &
234 & setparaccess = .false.)
235 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
238 vinfo( 2)=
'compensation light level'
239 vinfo( 3)=
'micromole_quanta'
240 status=def_var(ng, model, piofile, piovar, pio_type, &
241 & 1, (/phydim/), aval, vinfo, ncname, &
242 & setparaccess = .false.)
243 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
246 vinfo( 2)=
'light level for onset of photoinhibition'
247 vinfo( 3)=
'micromole_quanta'
248 status=def_var(ng, model, piofile, piovar, pio_type, &
249 & 1, (/phydim/), aval, vinfo, ncname, &
250 & setparaccess = .false.)
251 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
253 vinfo( 1)=
'inhib_fac'
254 vinfo( 2)=
'exponential decay factor for light limited growth'
255 vinfo( 3)=
'micromole_quanta-1'
256 status=def_var(ng, model, piofile, piovar, pio_type, &
257 & 1, (/phydim/), aval, vinfo, ncname, &
258 & setparaccess = .false.)
259 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
261 vinfo( 1)=
'C2CHL_max'
262 vinfo( 2)=
'maximum lighted limited C:Chl ratio'
263 vinfo( 3)=
'microgram_C microgram_Chl-1'
264 status=def_var(ng, model, piofile, piovar, pio_type, &
265 & 1, (/phydim/), aval, vinfo, ncname, &
266 & setparaccess = .false.)
267 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
270 vinfo( 2)=
'rate of change in light limited C:CHL ratio'
271 vinfo( 3)=
'microgram_C microgram_Chl-1 micromole_quanta-1'
272 status=def_var(ng, model, piofile, piovar, pio_type, &
273 & 1, (/phydim/), aval, vinfo, ncname, &
274 & setparaccess = .false.)
275 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
278 vinfo( 2)=
'minimum lighted limited C:Chl ratio'
279 vinfo( 3)=
'microgram_C microgram_Chl-1'
280 status=def_var(ng, model, piofile, piovar, pio_type, &
281 & 1, (/phydim/), aval, vinfo, ncname, &
282 & setparaccess = .false.)
283 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
286 vinfo( 2)=
'rate of change in nutient limited C:Chl ratio'
287 vinfo( 3)=
'microgram_C microgram_Chl-1 micromole_N micromole_C-1'
288 status=def_var(ng, model, piofile, piovar, pio_type, &
289 & 1, (/phydim/), aval, vinfo, ncname, &
290 & setparaccess = .false.)
291 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
294 vinfo( 2)=
'minimum nutrient limited C:CHL ratio'
295 vinfo( 3)=
'microgram_C microgram_Chl-1'
296 status=def_var(ng, model, piofile, piovar, pio_type, &
297 & 1, (/phydim/), aval, vinfo, ncname, &
298 & setparaccess = .false.)
299 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
302 vinfo( 2)=
'rate of change in package effect'
303 vinfo( 3)=
'microgram_Chl microgram_C-1'
304 status=def_var(ng, model, piofile, piovar, pio_type, &
305 & 1, (/phydim/), aval, vinfo, ncname, &
306 & setparaccess = .false.)
307 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
310 vinfo( 2)=
'maximum package effect'
311 vinfo( 3)=
'microgram_Chl microgram_C-1'
312 status=def_var(ng, model, piofile, piovar, pio_type, &
313 & 1, (/phydim/), aval, vinfo, ncname, &
314 & setparaccess = .false.)
315 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
318 vinfo( 2)=
'rate of change in the Chl_b:CHL_a ratio'
320 &
'microgram_Chl_b microgram_Chl_a-1 microgram_Chl_a microgrma_C-1'
321 status=def_var(ng, model, piofile, piovar, pio_type, &
322 & 1, (/phydim/), aval, vinfo, ncname, &
323 & setparaccess = .false.)
324 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
327 vinfo( 2)=
'maximum Chl_b:Chl_a ratio'
328 vinfo( 3)=
'microgram_Chl_b microgram_Chl_a-1'
329 status=def_var(ng, model, piofile, piovar, pio_type, &
330 & 1, (/phydim/), aval, vinfo, ncname, &
331 & setparaccess = .false.)
332 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
335 vinfo( 2)=
'rate of change in the Chl_c:Chl_a ratio'
337 &
'microgram_Chl_c microgram_Chl_a-1 microgram_Chl_a microgrma_C-1'
338 status=def_var(ng, model, piofile, piovar, pio_type, &
339 & 1, (/phydim/), aval, vinfo, ncname, &
340 & setparaccess = .false.)
341 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
344 vinfo( 2)=
'maximum Chl_c:Chl_a ratio'
345 vinfo( 3)=
'microgram_Chl_c microgram_Chl_a-1'
346 status=def_var(ng, model, piofile, piovar, pio_type, &
347 & 1, (/phydim/), aval, vinfo, ncname, &
348 & setparaccess = .false.)
349 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
352 vinfo( 2)=
'rate of change in the PSC:Chl_a ratio'
354 &
'microgram_PSC microgram_Chl_a-1 microgram_Chl_a microgram_C-1'
355 status=def_var(ng, model, piofile, piovar, pio_type, &
356 & 1, (/phydim/), aval, vinfo, ncname, &
357 & setparaccess = .false.)
358 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
361 vinfo( 2)=
'maximum PSC:Chl_a ratio'
362 vinfo( 3)=
'microgram_PSC microgram_Chl_a-1'
363 status=def_var(ng, model, piofile, piovar, pio_type, &
364 & 1, (/phydim/), aval, vinfo, ncname, &
365 & setparaccess = .false.)
366 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
369 vinfo( 2)=
'rate of change in the PPC:Chl_a ratio'
371 &
'microgram_PPC microgram_Chl_a-1 microgram_Chl_a microgram_C-1'
372 status=def_var(ng, model, piofile, piovar, pio_type, &
373 & 1, (/phydim/), aval, vinfo, ncname, &
374 & setparaccess = .false.)
375 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
378 vinfo( 2)=
'maximum PPC:Chl_a ratio'
379 vinfo( 3)=
'microgram_Chl_c microgram_Chl_a-1'
380 status=def_var(ng, model, piofile, piovar, pio_type, &
381 & 1, (/phydim/), aval, vinfo, ncname, &
382 & setparaccess = .false.)
383 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
386 vinfo( 2)=
'rate of change in the LPUb:Chl_a ratio'
388 &
'microgram_LPUb microgram_Chl_a-1 microgram_Chl_a microgram_C-1'
389 status=def_var(ng, model, piofile, piovar, pio_type, &
390 & 1, (/phydim/), aval, vinfo, ncname, &
391 & setparaccess = .false.)
392 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
395 vinfo( 2)=
'Maximum LPUb:Chl_a ratio'
396 vinfo( 3)=
'microgram_HPUb microgram_Chl_a-1'
397 status=def_var(ng, model, piofile, piovar, pio_type, &
398 & 1, (/phydim/), aval, vinfo, ncname, &
399 & setparaccess = .false.)
400 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
403 vinfo( 2)=
'rate of change in the HPUb:Chl_a ratio'
405 &
'microgram_HPUb microgram_Chl_a-1 microgram_Chl_a microgram_C-1'
406 status=def_var(ng, model, piofile, piovar, pio_type, &
407 & 1, (/phydim/), aval, vinfo, ncname, &
408 & setparaccess = .false.)
409 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
412 vinfo( 2)=
'maximum HPUb:Chl_a ratio'
413 vinfo( 3)=
'microgram_HPUb microgram_Chl_a-1'
414 status=def_var(ng, model, piofile, piovar, pio_type, &
415 & 1, (/phydim/), aval, vinfo, ncname, &
416 & setparaccess = .false.)
417 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
420 vinfo( 2)=
'proportion of grazing stress apportioned to DOM'
421 status=def_var(ng, model, piofile, piovar, pio_type, &
422 & 1, (/phydim/), aval, vinfo, ncname, &
423 & setparaccess = .false.)
424 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
427 vinfo( 2)=
'proportion of grazing stress apportioned to fecal'
428 status=def_var(ng, model, piofile, piovar, pio_type, &
429 & 1, (/phydim/), aval, vinfo, ncname, &
430 & setparaccess = .false.)
431 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
434 vinfo( 2)=
'proportion of grazing stress that is recycled'
435 status=def_var(ng, model, piofile, piovar, pio_type, &
436 & 1, (/phydim/), aval, vinfo, ncname, &
437 & setparaccess = .false.)
438 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
441 vinfo( 2)=
'proportion of daily production lost to excretion'
442 status=def_var(ng, model, piofile, piovar, pio_type, &
443 & 1, (/phydim/), aval, vinfo, ncname, &
444 & setparaccess = .false.)
445 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
448 vinfo( 2)=
'phytoplankton sinking speed'
449 vinfo( 3)=
'meter day-1'
450 status=def_var(ng, model, piofile, piovar, pio_type, &
451 & 1, (/phydim/), aval, vinfo, ncname, &
452 & setparaccess = .false.)
453 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
456 vinfo( 2)=
'phytoplankton grazing parameter'
457 status=def_var(ng, model, piofile, piovar, pio_type, &
458 & 1, (/phydim/), aval, vinfo, ncname, &
459 & setparaccess = .false.)
460 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
462 vinfo( 1)=
'MinRefuge'
463 vinfo( 2)=
'refuge phytoplankton population'
464 vinfo( 3)=
'micromole_C liter-1'
465 status=def_var(ng, model, piofile, piovar, pio_type, &
466 & 1, (/phydim/), aval, vinfo, ncname, &
467 & setparaccess = .false.)
468 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
470 vinfo( 1)=
'RefugeDep'
471 vinfo( 2)=
'maximum refuge phytoplankton depth'
473 status=def_var(ng, model, piofile, piovar, pio_type, &
474 & 1, (/phydim/), aval, vinfo, ncname, &
475 & setparaccess = .false.)
476 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
479 vinfo( 2)=
'normalized volume factor'
480 status=def_var(ng, model, piofile, piovar, pio_type, &
481 & 1, (/phydim/), aval, vinfo, ncname, &
482 & setparaccess = .false.)
483 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
485 vinfo( 1)=
'Norm_Surf'
486 vinfo( 2)=
'normalized surface area factor'
487 status=def_var(ng, model, piofile, piovar, pio_type, &
488 & 1, (/phydim/), aval, vinfo, ncname, &
489 & setparaccess = .false.)
490 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
493 vinfo( 2)=
'half-saturation constant for DOP uptake'
494 vinfo( 3)=
'micromole_DOP liter-1'
495 status=def_var(ng, model, piofile, piovar, pio_type, &
496 & 1, (/phydim/), aval, vinfo, ncname, &
497 & setparaccess = .false.)
498 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
500 vinfo( 1)=
'C2pALKPHOS'
501 vinfo( 2)=
'C:P ratio where DOP uptake begins'
502 vinfo( 3)=
'micromole_C micromole_P-1'
503 status=def_var(ng, model, piofile, piovar, pio_type, &
504 & 1, (/phydim/), aval, vinfo, ncname, &
505 & setparaccess = .false.)
506 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
509 vinfo( 2)=
'Half Saturation Constant for DON uptake'
510 vinfo( 3)=
'micromole_DON liter-1'
511 status=def_var(ng, model, piofile, piovar, pio_type, &
512 & 1, (/phydim/), aval, vinfo, ncname, &
513 & setparaccess = .false.)
514 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
516 vinfo( 1)=
'C2nNupDON'
517 vinfo( 2)=
'C:N ratio where DON uptake begins'
518 vinfo( 3)=
'micromole_C micromole_N-1'
519 status=def_var(ng, model, piofile, piovar, pio_type, &
520 & 1, (/phydim/), aval, vinfo, ncname, &
521 & setparaccess = .false.)
522 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
525 vinfo( 2)=
'half-saturation constant for bacteria DOC uptake'
526 vinfo( 3)=
'micromole_DOC liter-1'
527 status=def_var(ng, model, piofile, piovar, pio_type, &
528 & 1, (/bacdim/), aval, vinfo, ncname, &
529 & setparaccess = .false.)
530 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
532 vinfo( 1)=
'GtBAC_max'
533 vinfo( 2)=
'maximum 24 hour bacterial growth rate'
535 status=def_var(ng, model, piofile, piovar, pio_type, &
536 & 1, (/bacdim/), aval, vinfo, ncname, &
537 & setparaccess = .false.)
538 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
541 vinfo( 2)=
'bacteria temperature base for exponential response'
543 status=def_var(ng, model, piofile, piovar, pio_type, &
544 & 1, (/bacdim/), aval, vinfo, ncname, &
545 & setparaccess = .false.)
546 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
549 vinfo( 2)=
'bacteria exponential temperature factor'
551 status=def_var(ng, model, piofile, piovar, pio_type, &
552 & 1, (/bacdim/), aval, vinfo, ncname, &
553 & setparaccess = .false.)
554 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
557 vinfo( 2)=
'carbon to nitrogen ratio of bacteria'
558 vinfo( 3)=
'micromole_C micromole_N-1'
559 status=def_var(ng, model, piofile, piovar, pio_type, &
560 & 1, (/0/), aval, vinfo, ncname, &
561 & setparaccess = .false.)
562 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
565 vinfo( 2)=
'carbon to phosphorus ratio of bacteria'
566 vinfo( 3)=
'micromole_C micromole_P-1'
567 status=def_var(ng, model, piofile, piovar, pio_type, &
568 & 1, (/0/), aval, vinfo, ncname, &
569 & setparaccess = .false.)
570 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
573 vinfo( 2)=
'carbon to iron ratio of bacteria'
574 vinfo( 3)=
'micromole_C micromole_Fe-1'
575 status=def_var(ng, model, piofile, piovar, pio_type, &
576 & 1, (/0/), aval, vinfo, ncname, &
577 & setparaccess = .false.)
578 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
582 &
'proportion of bacteria grazing stress apportioned to DOM'
583 status=def_var(ng, model, piofile, piovar, pio_type, &
584 & 1, (/0/), aval, vinfo, ncname, &
585 & setparaccess = .false.)
586 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
590 &
'proportion of bacteria grazing stress apportioned to fecal'
591 status=def_var(ng, model, piofile, piovar, pio_type, &
592 & 1, (/0/), aval, vinfo, ncname, &
593 & setparaccess = .false.)
594 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
597 vinfo( 2)=
'proportion of bacteria grazing stress recycled'
598 status=def_var(ng, model, piofile, piovar, pio_type, &
599 & 1, (/0/), aval, vinfo, ncname, &
600 & setparaccess = .false.)
601 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
605 &
'bacterial recalcitrant C excretion as proportion of uptake'
606 status=def_var(ng, model, piofile, piovar, pio_type, &
607 & 1, (/0/), aval, vinfo, ncname, &
608 & setparaccess = .false.)
609 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
612 vinfo( 2)=
'bacterial recalcitrant excretion carbon:nitrogen ratio'
613 vinfo( 3)=
'micromole_C micromole_N-1'
614 status=def_var(ng, model, piofile, piovar, pio_type, &
615 & 1, (/0/), aval, vinfo, ncname, &
616 & setparaccess = .false.)
617 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
620 vinfo( 2)=
'bacterial gross growth carbon efficiency'
621 status=def_var(ng, model, piofile, piovar, pio_type, &
622 & 1, (/0/), aval, vinfo, ncname, &
623 & setparaccess = .false.)
624 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
627 vinfo( 2)=
'maximum bacterial nitrification rate'
629 status=def_var(ng, model, piofile, piovar, pio_type, &
630 & 1, (/0/), aval, vinfo, ncname, &
631 & setparaccess = .false.)
632 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
635 vinfo( 2)=
'half-saturation constant for bacterial nitrification'
636 vinfo( 3)=
'micromole_NH4 liter-1'
637 status=def_var(ng, model, piofile, piovar, pio_type, &
638 & 1, (/0/), aval, vinfo, ncname, &
639 & setparaccess = .false.)
640 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
642 vinfo( 1)=
'cDOCfrac_c'
644 &
'colored fraction of DOC from phytoplankton and bacteria losses'
645 status=def_var(ng, model, piofile, piovar, pio_type, &
646 & 1, (/domdim/), aval, vinfo, ncname, &
647 & setparaccess = .false.)
648 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
650 vinfo( 1)=
'RtUVR_DIC'
651 vinfo( 2)=
'UV degradation of DOC into DIC at 410 nanometer'
652 vinfo( 3)=
'micromole meter-1 liter-1 hour-1'
653 status=def_var(ng, model, piofile, piovar, pio_type, &
654 & 1, (/0/), aval, vinfo, ncname, &
655 & setparaccess = .false.)
656 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
658 vinfo( 1)=
'RtUVR_DOC'
660 &
'UV degradation of DOC to colorless labile DOC at 410 nanometer'
661 vinfo( 3)=
'micromole meter-1 liter-1 hour-1'
662 status=def_var(ng, model, piofile, piovar, pio_type, &
663 & 1, (/0/), aval, vinfo, ncname, &
664 & setparaccess = .false.)
665 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
668 vinfo( 2)=
'fecal sinking flux'
669 vinfo( 3)=
'meter day-1'
670 status=def_var(ng, model, piofile, piovar, pio_type, &
671 & 1, (/fecdim/), aval, vinfo, ncname, &
672 & setparaccess = .false.)
673 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
677 &
'fecal regeneration temperature base for exponential response'
679 status=def_var(ng, model, piofile, piovar, pio_type, &
680 & 1, (/fecdim/), aval, vinfo, ncname, &
681 & setparaccess = .false.)
682 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
685 vinfo( 2)=
'fecal regeneration exponential temperature factor'
686 vinfo( 3)=
'Celsius-1'
687 status=def_var(ng, model, piofile, piovar, pio_type, &
688 & 1, (/fecdim/), aval, vinfo, ncname, &
689 & setparaccess = .false.)
690 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
693 vinfo( 2)=
'fecal carbon regeneration rate'
695 status=def_var(ng, model, piofile, piovar, pio_type, &
696 & 1, (/fecdim/), aval, vinfo, ncname, &
697 & setparaccess = .false.)
698 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
701 vinfo( 2)=
'fecal nitrogen regeneration rate'
703 status=def_var(ng, model, piofile, piovar, pio_type, &
704 & 1, (/fecdim/), aval, vinfo, ncname, &
705 & setparaccess = .false.)
706 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
709 vinfo( 2)=
'fecal silica regeneration rate'
711 status=def_var(ng, model, piofile, piovar, pio_type, &
712 & 1, (/fecdim/), aval, vinfo, ncname, &
713 & setparaccess = .false.)
714 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
717 vinfo( 2)=
'fecal phosphorus regeneration rate'
719 status=def_var(ng, model, piofile, piovar, pio_type, &
720 & 1, (/fecdim/), aval, vinfo, ncname, &
721 & setparaccess = .false.)
722 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
725 vinfo( 2)=
'fecal iron regeneration rate'
727 status=def_var(ng, model, piofile, piovar, pio_type, &
728 & 1, (/fecdim/), aval, vinfo, ncname)
729 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN