ROMS
Loading...
Searching...
No Matches
ecosim_def_pio.h
Go to the documentation of this file.
1/*
2** git $id$
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************************************************************************
8** **
9** defines ecosim bio-optical model input parameters in output **
10** netcdf files. it is included in routine "def_info.F". **
11** **
12************************************************************************
13*/
14
15!
16! Define EcoSim bio-optical model parameters.
17!
18 vinfo( 1)='BioIter'
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
24
25 vinfo( 1)='RtUVR_flag'
26 vinfo( 2)='switch to calculate CDOC UV photolysis.'
27 vinfo( 9)='.FALSE.'
28 vinfo(10)='.TRUE.'
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
33
34 vinfo( 1)='NFIX_flag'
35 vinfo( 2)='switch to calculate temperature based N fixation'
36 vinfo( 9)='.FALSE.'
37 vinfo(10)='.TRUE.'
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
42
43 vinfo( 1)='Regen_flag'
44 vinfo( 2)='switch to calculate fecal matter regeneration'
45 vinfo( 9)='.FALSE.'
46 vinfo(10)='.TRUE.'
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
51
52 vinfo( 1)='light'
53 vinfo( 2)='center wavelength of spectral irradiance bands'
54 vinfo( 3)='nanometer'
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
59
60 vinfo( 1)='HsNO3'
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
67
68 vinfo( 1)='HsNH4'
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
75
76 vinfo( 1)='HsSiO'
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
83
84 vinfo( 1)='HsPO4'
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
91
92 vinfo( 1)='HsFe'
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
99
100 vinfo( 1)='GtALG_max'
101 vinfo( 2)='maximum phytoplankton 24 hour growth rate'
102 vinfo( 3)='day-1'
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
107
108 vinfo( 1)='PhyTbase'
109 vinfo( 2)= &
110 & 'phytoplankton temperature base for exponential response'
111 vinfo( 3)='Celsius'
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
116
117 vinfo( 1)='PhyTfac'
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
124
125 vinfo( 1)='BET_'
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
132
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
140
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
148
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
156
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
164
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
172
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
180
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
188
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
196
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
204
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
212
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
220
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
228
229 vinfo( 1)='qu_yld'
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
236
237 vinfo( 1)='E0_comp'
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
244
245 vinfo( 1)='E0_inhib'
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
252
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
260
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
268
269 vinfo( 1)='mxC2Cl'
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
276
277 vinfo( 1)='b_C2Cl'
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
284
285 vinfo( 1)='mxC2Cn'
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
292
293 vinfo( 1)='b_C2Cn'
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
300
301 vinfo( 1)='mxPacEff'
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
308
309 vinfo( 1)='b_PacEff'
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
316
317 vinfo( 1)='mxChlB'
318 vinfo( 2)='rate of change in the Chl_b:CHL_a ratio'
319 vinfo( 3)= &
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
325
326 vinfo( 1)='b_ChlB'
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
333
334 vinfo( 1)='mxChlC'
335 vinfo( 2)='rate of change in the Chl_c:Chl_a ratio'
336 vinfo( 3)= &
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
342
343 vinfo( 1)='b_ChlC'
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
350
351 vinfo( 1)='mxPSC'
352 vinfo( 2)='rate of change in the PSC:Chl_a ratio'
353 vinfo( 3)= &
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
359
360 vinfo( 1)='b_PSC'
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
367
368 vinfo( 1)='mxPPC'
369 vinfo( 2)='rate of change in the PPC:Chl_a ratio'
370 vinfo( 3)= &
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
376
377 vinfo( 1)='b_PPC'
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
384
385 vinfo( 1)='mxLPUb'
386 vinfo( 2)='rate of change in the LPUb:Chl_a ratio'
387 vinfo( 3)= &
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
393
394 vinfo( 1)='b_LPUb'
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
401
402 vinfo( 1)='mxHPUb'
403 vinfo( 2)='rate of change in the HPUb:Chl_a ratio'
404 vinfo( 3)= &
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
410
411 vinfo( 1)='b_HPUb'
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
418
419 vinfo( 1)='FecDOC'
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
425
426 vinfo( 1)='FecPEL'
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
432
433 vinfo( 1)='FecCYC'
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
439
440 vinfo( 1)='ExALG'
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
446
447 vinfo( 1)='WS'
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
454
455 vinfo( 1)='HsGRZ'
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
461
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
469
470 vinfo( 1)='RefugeDep'
471 vinfo( 2)='maximum refuge phytoplankton depth'
472 vinfo( 3)='meter'
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
477
478 vinfo( 1)='Norm_Vol'
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
484
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
491
492 vinfo( 1)='HsDOP'
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
499
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
507
508 vinfo( 1)='HsDON'
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
515
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
523
524 vinfo( 1)='HsDOC_ba'
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
531
532 vinfo( 1)='GtBAC_max'
533 vinfo( 2)='maximum 24 hour bacterial growth rate'
534 vinfo( 3)='day-1'
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
539
540 vinfo( 1)='BacTbase'
541 vinfo( 2)='bacteria temperature base for exponential response'
542 vinfo( 3)='Celsius'
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
547
548 vinfo( 1)='BacTfac'
549 vinfo( 2)='bacteria exponential temperature factor'
550 vinfo( 3)='Celsius'
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
555
556 vinfo( 1)='C2nBAC'
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
563
564 vinfo( 1)='C2pBAC'
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
571
572 vinfo( 1)='C2FeBAC'
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
579
580 vinfo( 1)='BacDOC'
581 vinfo( 2)= &
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
587
588 vinfo( 1)='BacPEL'
589 vinfo( 2)= &
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
595
596 vinfo( 1)='BacCYC'
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
602
603 vinfo( 1)='ExBAC_c'
604 vinfo( 2)= &
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
610
611 vinfo( 1)='ExBacC2N'
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
618
619 vinfo( 1)='Bac_Ceff'
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
625
626 vinfo( 1)='RtNIT'
627 vinfo( 2)='maximum bacterial nitrification rate'
628 vinfo( 3)='day-1'
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
633
634 vinfo( 1)='HsNIT'
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
641
642 vinfo( 1)='cDOCfrac_c'
643 vinfo( 2)= &
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
649
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
657
658 vinfo( 1)='RtUVR_DOC'
659 vinfo( 2)= &
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
666
667 vinfo( 1)='WF'
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
674
675 vinfo( 1)='RegTbase'
676 vinfo( 2)= &
677 & 'fecal regeneration temperature base for exponential response'
678 vinfo( 3)='Celsius'
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
683
684 vinfo( 1)='RegTfac'
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
691
692 vinfo( 1)='RegCR'
693 vinfo( 2)='fecal carbon regeneration rate'
694 vinfo( 3)='day-1'
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
699
700 vinfo( 1)='RegNR'
701 vinfo( 2)='fecal nitrogen regeneration rate'
702 vinfo( 3)='day-1'
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
707
708 vinfo( 1)='RegSR'
709 vinfo( 2)='fecal silica regeneration rate'
710 vinfo( 3)='day-1'
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
715
716 vinfo( 1)='RegPR'
717 vinfo( 2)='fecal phosphorus regeneration rate'
718 vinfo( 3)='day-1'
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
723
724 vinfo( 1)='RegFR'
725 vinfo( 2)='fecal iron regeneration rate'
726 vinfo( 3)='day-1'
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
subroutine output(ng)
Definition output.F:4