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 nemuro ecosystem model input parameters in
output netcdf **
10** files. it is included in routine
"def_info.F". **
12************************************************************************
19 vinfo( 2)=
'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
26 vinfo( 2)=
'light attenuation by seawater'
28 status=def_var(ng, model, piofile, piovar, pio_type, &
29 & 1, (/0/), aval, vinfo, ncname, &
30 & setparaccess = .false.)
31 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
34 vinfo( 2)=
'light attenuation due to small phytoplankton, '// &
35 &
'self-shading coefficient'
36 vinfo( 3)=
'meter2 millimole_N-1'
37 status=def_var(ng, model, piofile, piovar, pio_type, &
38 & 1, (/0/), aval, vinfo, ncname, &
39 & setparaccess = .false.)
40 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
43 vinfo( 2)=
'light attenuation due to large phytoplankton, '// &
44 &
'self-shading coefficient'
45 vinfo( 3)=
'meter2 millimole_N-1'
46 status=def_var(ng, model, piofile, piovar, pio_type, &
47 & 1, (/0/), aval, vinfo, ncname, &
48 & setparaccess = .false.)
49 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
52 vinfo( 2)=
'photosynthetically available radiation fraction'
53 status=def_var(ng, model, piofile, piovar, pio_type, &
54 & 1, (/0/), aval, vinfo, ncname, &
55 & setparaccess = .false.)
56 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
59 vinfo( 2)=
'small phytoplankton initial slope of P-I curve'
60 vinfo( 3)=
'meter2 watt-1 day-1'
61 status=def_var(ng, model, piofile, piovar, pio_type, &
62 & 1, (/0/), aval, vinfo, ncname, &
63 & setparaccess = .false.)
64 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
67 vinfo( 2)=
'Large phytoplankton initial slope of P-I curve'
68 vinfo( 3)=
'meter2 watt-1 day-1'
69 status=def_var(ng, model, piofile, piovar, pio_type, &
70 & 1, (/0/), aval, vinfo, ncname, &
71 & setparaccess = .false.)
72 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
75 vinfo( 2)=
'small phytoplankton photoinhibition coefficient'
76 vinfo( 3)=
'meter2 watt-1 day-1'
77 status=def_var(ng, model, piofile, piovar, pio_type, &
78 & 1, (/0/), aval, vinfo, ncname, &
79 & setparaccess = .false.)
80 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
83 vinfo( 2)=
'large phytoplankton photoinhibition coefficient'
84 vinfo( 3)=
'meter2 watt-1 day-1'
85 status=def_var(ng, model, piofile, piovar, pio_type, &
86 & 1, (/0/), aval, vinfo, ncname, &
87 & setparaccess = .false.)
88 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
91 vinfo( 2)=
'small phytoplankton maximum photosynthetic rate '// &
94 status=def_var(ng, model, piofile, piovar, pio_type, &
95 & 1, (/0/), aval, vinfo, ncname, &
96 & setparaccess = .false.)
97 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
100 vinfo( 2)=
'large phytoplankton maximum photosynthetic rate '// &
103 status=def_var(ng, model, piofile, piovar, pio_type, &
104 & 1, (/0/), aval, vinfo, ncname, &
105 & setparaccess = .false.)
106 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
109 vinfo( 2)=
'small phytoplankton half-saturation constant '// &
111 vinfo( 3)=
'millimole_N meter-3'
112 status=def_var(ng, model, piofile, piovar, pio_type, &
113 & 1, (/0/), aval, vinfo, ncname, &
114 & setparaccess = .false.)
115 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
118 vinfo( 2)=
'large phytoplankton half saturation constant '// &
120 vinfo( 3)=
'millimole_N meter-3'
121 status=def_var(ng, model, piofile, piovar, pio_type, &
122 & 1, (/0/), aval, vinfo, ncname, &
123 & setparaccess = .false.)
124 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
128 'small phytoplankton half saturation constant for ammonium'
129 vinfo( 3)=
'millimole_N meter-3'
130 status=def_var(ng, model, piofile, piovar, pio_type, &
131 & 1, (/0/), aval, vinfo, ncname, &
132 & setparaccess = .false.)
133 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
136 vinfo( 2)=
'large phytoplankton half saturation constant '// &
138 vinfo( 3)=
'millimole_N meter-3'
139 status=def_var(ng, model, piofile, piovar, pio_type, &
140 & 1, (/0/), aval, vinfo, ncname, &
141 & setparaccess = .false.)
142 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
145 vinfo( 2)=
'large phytoplankton half saturation constant '// &
147 vinfo( 3)=
'millimole_Si meter-3'
148 status=def_var(ng, model, piofile, piovar, pio_type, &
149 & 1, (/0/), aval, vinfo, ncname, &
150 & setparaccess = .false.)
151 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
154 vinfo( 2)=
'small phytoplankton ammonium inhibition coefficient'
155 vinfo( 3)=
'meter3 millimole_N-1'
156 status=def_var(ng, model, piofile, piovar, pio_type, &
157 & 1, (/0/), aval, vinfo, ncname, &
158 & setparaccess = .false.)
159 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
162 vinfo( 2)=
'small phytoplankton ammonium inhibition coefficient'
163 vinfo( 3)=
'meter3 millimole_N-1'
164 status=def_var(ng, model, piofile, piovar, pio_type, &
165 & 1, (/0/), aval, vinfo, ncname, &
166 & setparaccess = .false.)
167 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
170 vinfo( 2)=
'small phytoplankton temperature coefficient for '// &
171 &
'photosynthetic rate'
172 vinfo( 3)=
'Celsius-1'
173 status=def_var(ng, model, piofile, piovar, pio_type, &
174 & 1, (/0/), aval, vinfo, ncname, &
175 & setparaccess = .false.)
176 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
179 vinfo( 2)=
'large phytoplankton temperature coefficient for '// &
180 &
'photosynthetic rate'
181 vinfo( 3)=
'Celsius-1'
182 status=def_var(ng, model, piofile, piovar, pio_type, &
183 & 1, (/0/), aval, vinfo, ncname, &
184 & setparaccess = .false.)
185 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
188 vinfo( 2)=
'small phytoplankton respiration rate at 0 Celsius'
190 status=def_var(ng, model, piofile, piovar, pio_type, &
191 & 1, (/0/), aval, vinfo, ncname, &
192 & setparaccess = .false.)
193 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
196 vinfo( 2)=
'large phytoplankton respiration rate at 0 Celsius'
198 status=def_var(ng, model, piofile, piovar, pio_type, &
199 & 1, (/0/), aval, vinfo, ncname, &
200 & setparaccess = .false.)
201 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
204 vinfo( 2)=
'small phytoplankton temperature coefficient '// &
206 vinfo( 3)=
'Celsius-1'
207 status=def_var(ng, model, piofile, piovar, pio_type, &
208 & 1, (/0/), aval, vinfo, ncname, &
209 & setparaccess = .false.)
210 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
213 vinfo( 2)=
'large phytoplankton temperature coefficient '// &
215 vinfo( 3)=
'Celsius-1'
216 status=def_var(ng, model, piofile, piovar, pio_type, &
217 & 1, (/0/), aval, vinfo, ncname, &
218 & setparaccess = .false.)
219 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
222 vinfo( 2)=
'small phytoplankton ratio of extracellular '// &
223 &
'excretion to photosynthesis'
224 status=def_var(ng, model, piofile, piovar, pio_type, &
225 & 1, (/0/), aval, vinfo, ncname, &
226 & setparaccess = .false.)
227 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
230 vinfo( 2)=
'large phytoplankton ratio of extracellular '// &
231 &
'excretion to photosynthesis'
232 status=def_var(ng, model, piofile, piovar, pio_type, &
233 & 1, (/0/), aval, vinfo, ncname, &
234 & setparaccess = .false.)
235 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
238 vinfo( 2)=
'small phytoplankton mortality rate at 0 Celsius'
239 vinfo( 3)=
'meter3 millimole_N-1 day-1'
240 status=def_var(ng, model, piofile, piovar, pio_type, &
241 & 1, (/0/), aval, vinfo, ncname, &
242 & setparaccess = .false.)
243 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
246 vinfo( 2)=
'large phytoplankton mortality rate at 0 Celsius'
247 vinfo( 3)=
'meter3 millimole_N-1 day-1'
248 status=def_var(ng, model, piofile, piovar, pio_type, &
249 & 1, (/0/), aval, vinfo, ncname, &
250 & setparaccess = .false.)
251 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
254 vinfo( 2)=
'small phytoplankton temperature coefficient for '// &
256 vinfo( 3)=
'Celsius-1'
257 status=def_var(ng, model, piofile, piovar, pio_type, &
258 & 1, (/0/), aval, vinfo, ncname, &
259 & setparaccess = .false.)
260 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
263 vinfo( 2)=
'large phytoplankton temperature coefficient for '// &
265 vinfo( 3)=
'Celsius-1'
266 status=def_var(ng, model, piofile, piovar, pio_type, &
267 & 1, (/0/), aval, vinfo, ncname, &
268 & setparaccess = .false.)
269 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
272 vinfo( 2)=
'small zooplankton maximum grazing rate on '// &
273 &
'small phytoplankton at 0 Celsius'
275 status=def_var(ng, model, piofile, piovar, pio_type, &
276 & 1, (/0/), aval, vinfo, ncname, &
277 & setparaccess = .false.)
278 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
281 vinfo( 2)=
'large zooplankton maximum grazing rate on '// &
282 &
'small phytoplankton at 0 Celsius'
284 status=def_var(ng, model, piofile, piovar, pio_type, &
285 & 1, (/0/), aval, vinfo, ncname, &
286 & setparaccess = .false.)
287 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
290 vinfo( 2)=
'large zooplankton maximum grazing rate on '// &
291 &
'large phytoplankton at 0 Celsius'
293 status=def_var(ng, model, piofile, piovar, pio_type, &
294 & 1, (/0/), aval, vinfo, ncname, &
295 & setparaccess = .false.)
296 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
299 vinfo( 2)=
'large zooplankton maximum grazing rate on '// &
300 &
'small zooplankton at 0 Celsius'
302 status=def_var(ng, model, piofile, piovar, pio_type, &
303 & 1, (/0/), aval, vinfo, ncname, &
304 & setparaccess = .false.)
305 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
308 vinfo( 2)=
'predator zooplankton maximum grazing rate on '// &
309 &
'large phytoplankton at 0 Celsius'
311 status=def_var(ng, model, piofile, piovar, pio_type, &
312 & 1, (/0/), aval, vinfo, ncname, &
313 & setparaccess = .false.)
314 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
317 vinfo( 2)=
'predator zooplankton maximum grazing rate on '// &
318 &
'small zooplankton at 0 Celsius'
320 status=def_var(ng, model, piofile, piovar, pio_type, &
321 & 1, (/0/), aval, vinfo, ncname, &
322 & setparaccess = .false.)
323 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
326 vinfo( 2)=
'predator zooplankton maximum grazing rate on '// &
327 &
'large zooplankton at 0 Celsius'
329 status=def_var(ng, model, piofile, piovar, pio_type, &
330 & 1, (/0/), aval, vinfo, ncname, &
331 & setparaccess = .false.)
332 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
335 vinfo( 2)=
'small zooplankton temperature coefficient '// &
337 vinfo( 3)=
'Celsius-1'
338 status=def_var(ng, model, piofile, piovar, pio_type, &
339 & 1, (/0/), aval, vinfo, ncname, &
340 & setparaccess = .false.)
341 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
344 vinfo( 2)=
'large zooplankton temperature coefficient '// &
346 vinfo( 3)=
'Celsius-1'
347 status=def_var(ng, model, piofile, piovar, pio_type, &
348 & 1, (/0/), aval, vinfo, ncname, &
349 & setparaccess = .false.)
350 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
353 vinfo( 2)=
'predator zooplankton temperature coefficient '// &
355 vinfo( 3)=
'Celsius-1'
356 status=def_var(ng, model, piofile, piovar, pio_type, &
357 & 1, (/0/), aval, vinfo, ncname, &
358 & setparaccess = .false.)
359 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
362 vinfo( 2)=
'small zooplankton Ivlev constant'
363 vinfo( 3)=
'meter3 millimole_N'
364 status=def_var(ng, model, piofile, piovar, pio_type, &
365 & 1, (/0/), aval, vinfo, ncname, &
366 & setparaccess = .false.)
367 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
370 vinfo( 2)=
'large zooplankton Ivlev constant'
371 vinfo( 3)=
'meter3 millimole_N'
372 status=def_var(ng, model, piofile, piovar, pio_type, &
373 & 1, (/0/), aval, vinfo, ncname, &
374 & setparaccess = .false.)
375 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
378 vinfo( 2)=
'predator zooplankton Ivlev constant'
379 vinfo( 3)=
'meter3 millimole_N'
380 status=def_var(ng, model, piofile, piovar, pio_type, &
381 & 1, (/0/), aval, vinfo, ncname, &
382 & setparaccess = .false.)
383 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
385#ifdef HOLLING_GRAZING
387 vinfo( 2)=
'small zooplankton squared half-saturation '// &
388 &
'coefficient for ingestion on small phytoplankton'
389 vinfo( 3)=
'millimole_N2 meter-6'
390 status=def_var(ng, model, piofile, piovar, pio_type, &
391 & 1, (/0/), aval, vinfo, ncname, &
392 & setparaccess = .false.)
393 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
396 vinfo( 2)=
'large zooplankton squared half-saturation '// &
397 &
'coefficient for ingestion on small phytoplankton'
398 vinfo( 3)=
'millimole_N2 meter-6'
399 status=def_var(ng, model, piofile, piovar, pio_type, &
400 & 1, (/0/), aval, vinfo, ncname, &
401 & setparaccess = .false.)
402 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
405 vinfo( 2)=
'large zooplankton squared half-saturation '// &
406 &
'coefficient for ingestion on large phytoplankton'
407 vinfo( 3)=
'millimole_N2 meter-6'
408 status=def_var(ng, model, piofile, piovar, pio_type, &
409 & 1, (/0/), aval, vinfo, ncname, &
410 & setparaccess = .false.)
411 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
414 vinfo( 2)=
'large zooplankton squared half-saturation '// &
415 &
'coefficient for ingestion on small zooplankton'
416 vinfo( 3)=
'millimole_N2 meter-6'
417 status=def_var(ng, model, piofile, piovar, pio_type, &
418 & 1, (/0/), aval, vinfo, ncname, &
419 & setparaccess = .false.)
420 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
423 vinfo( 2)=
'predator zooplankton squared half-saturation '// &
424 &
'coefficient for ingestion on large phytoplankton'
425 vinfo( 3)=
'millimole_N2 meter-6'
426 status=def_var(ng, model, piofile, piovar, pio_type, &
427 & 1, (/0/), aval, vinfo, ncname, &
428 & setparaccess = .false.)
429 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
432 vinfo( 2)=
'predator zooplankton squared half-saturation '// &
433 &
'coefficient for ingestion on small zooplankton'
434 vinfo( 3)=
'millimole_N2 meter-6'
435 status=def_var(ng, model, piofile, piovar, pio_type, &
436 & 1, (/0/), aval, vinfo, ncname, &
437 & setparaccess = .false.)
438 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
441 vinfo( 2)=
'predator zooplankton squared half-saturation '// &
442 &
'coefficient for ingestion on large zooplankton'
443 vinfo( 3)=
'millimole_N2 meter-6'
444 status=def_var(ng, model, piofile, piovar, pio_type, &
445 & 1, (/0/), aval, vinfo, ncname, &
446 & setparaccess = .false.)
447 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
451 vinfo( 1)=
'PS2ZSstar'
452 vinfo( 2)=
'small zooplankton threshold value for grazing '// &
453 &
'on small phytoplankton'
454 vinfo( 3)=
'millimole_N meter-3'
455 status=def_var(ng, model, piofile, piovar, pio_type, &
456 & 1, (/0/), aval, vinfo, ncname, &
457 & setparaccess = .false.)
458 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
460 vinfo( 1)=
'PS2ZLstar'
461 vinfo( 2)=
'large zooplankton threshold value for grazing '// &
462 &
'on small phytoplankton'
463 vinfo( 3)=
'millimole_N meter-3'
464 status=def_var(ng, model, piofile, piovar, pio_type, &
465 & 1, (/0/), aval, vinfo, ncname, &
466 & setparaccess = .false.)
467 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
469 vinfo( 1)=
'PL2ZLstar'
470 vinfo( 2)=
'large zooplankton threshold value for grazing '// &
471 &
'on large phytoplankton'
472 vinfo( 3)=
'millimole_N meter-3'
473 status=def_var(ng, model, piofile, piovar, pio_type, &
474 & 1, (/0/), aval, vinfo, ncname, &
475 & setparaccess = .false.)
476 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
478 vinfo( 1)=
'ZS2ZLstar'
479 vinfo( 2)=
'large zooplankton threshold value for grazing '// &
480 &
'on small zooplankton'
481 vinfo( 3)=
'millimole_N meter-3'
482 status=def_var(ng, model, piofile, piovar, pio_type, &
483 & 1, (/0/), aval, vinfo, ncname, &
484 & setparaccess = .false.)
485 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
487 vinfo( 1)=
'PL2ZPstar'
488 vinfo( 2)=
'predator zooplankton threshold value for grazing '// &
489 &
'on large phytoplankton'
490 vinfo( 3)=
'millimole_N meter-3'
491 status=def_var(ng, model, piofile, piovar, pio_type, &
492 & 1, (/0/), aval, vinfo, ncname, &
493 & setparaccess = .false.)
494 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
496 vinfo( 1)=
'ZS2ZPstar'
497 vinfo( 2)=
'predator zooplankton threshold value for grazing '// &
498 &
'on small zooplankton'
499 vinfo( 3)=
'millimole_N meter-3'
500 status=def_var(ng, model, piofile, piovar, pio_type, &
501 & 1, (/0/), aval, vinfo, ncname, &
502 & setparaccess = .false.)
503 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
505 vinfo( 1)=
'ZL2ZPstar'
506 vinfo( 2)=
'predator zooplankton threshold value for grazing '// &
507 &
'on large zooplankton'
508 vinfo( 3)=
'millimole_N meter-3'
509 status=def_var(ng, model, piofile, piovar, pio_type, &
510 & 1, (/0/), aval, vinfo, ncname, &
511 & setparaccess = .false.)
512 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
516 vinfo( 2)=
'predator zooplankton grazing inhibition coefficient'// &
517 &
'on large phytoplankton'
518 vinfo( 3)=
'meter3 millimole_N-1'
519 status=def_var(ng, model, piofile, piovar, pio_type, &
520 & 1, (/0/), aval, vinfo, ncname, &
521 & setparaccess = .false.)
522 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
525 vinfo( 2)=
'predator zooplankton grazing inhibition coefficient'// &
526 &
'on small zooplankton'
527 vinfo( 3)=
'meter3 millimole_N-1'
528 status=def_var(ng, model, piofile, piovar, pio_type, &
529 & 1, (/0/), aval, vinfo, ncname, &
530 & setparaccess = .false.)
531 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
534 vinfo( 2)=
'small zooplankton mortality rate at 0 Celsius'
535 vinfo( 3)=
'meter3 millimole_N-1 day-1'
536 status=def_var(ng, model, piofile, piovar, pio_type, &
537 & 1, (/0/), aval, vinfo, ncname, &
538 & setparaccess = .false.)
539 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
542 vinfo( 2)=
'large zooplankton mortality rate at 0 Celsius'
543 vinfo( 3)=
'meter3 millimole_N-1 day-1'
544 status=def_var(ng, model, piofile, piovar, pio_type, &
545 & 1, (/0/), aval, vinfo, ncname, &
546 & setparaccess = .false.)
547 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
550 vinfo( 2)=
'predator zooplankton mortality rate at 0 Celsius'
551 vinfo( 3)=
'meter3 millimole_N-1 day-1'
552 status=def_var(ng, model, piofile, piovar, pio_type, &
553 & 1, (/0/), aval, vinfo, ncname, &
554 & setparaccess = .false.)
555 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
558 vinfo( 2)=
'small zooplankton assimilation efficiency'
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)=
'large zooplankton assimilation efficiency'
566 status=def_var(ng, model, piofile, piovar, pio_type, &
567 & 1, (/0/), aval, vinfo, ncname, &
568 & setparaccess = .false.)
569 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
572 vinfo( 2)=
'predator zooplankton assimilation efficiency'
573 status=def_var(ng, model, piofile, piovar, pio_type, &
574 & 1, (/0/), aval, vinfo, ncname, &
575 & setparaccess = .false.)
576 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
579 vinfo( 2)=
'small zooplankton growth efficiency'
580 status=def_var(ng, model, piofile, piovar, pio_type, &
581 & 1, (/0/), aval, vinfo, ncname, &
582 & setparaccess = .false.)
583 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
586 vinfo( 2)=
'large zooplankton growth efficiency'
587 status=def_var(ng, model, piofile, piovar, pio_type, &
588 & 1, (/0/), aval, vinfo, ncname, &
589 & setparaccess = .false.)
590 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
593 vinfo( 2)=
'predator zooplankton growth efficiency'
594 status=def_var(ng, model, piofile, piovar, pio_type, &
595 & 1, (/0/), aval, vinfo, ncname, &
596 & setparaccess = .false.)
597 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
600 vinfo( 2)=
'NH4 nitrification rate at 0 Celsius'
602 status=def_var(ng, model, piofile, piovar, pio_type, &
603 & 1, (/0/), aval, vinfo, ncname, &
604 & setparaccess = .false.)
605 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
608 vinfo( 2)=
'PON decomposition to NH4 rate at 0 Celsius'
610 status=def_var(ng, model, piofile, piovar, pio_type, &
611 & 1, (/0/), aval, vinfo, ncname, &
612 & setparaccess = .false.)
613 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
616 vinfo( 2)=
'PON decomposition to DON rate at 0 Celsius'
618 status=def_var(ng, model, piofile, piovar, pio_type, &
619 & 1, (/0/), aval, vinfo, ncname, &
620 & setparaccess = .false.)
621 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
624 vinfo( 2)=
'DON decomposition to NH4 rate at 0 Celsius'
626 status=def_var(ng, model, piofile, piovar, pio_type, &
627 & 1, (/0/), aval, vinfo, ncname, &
628 & setparaccess = .false.)
629 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
632 vinfo( 2)=
'opal decomposition to silicate rate at 0 Celsius'
634 status=def_var(ng, model, piofile, piovar, pio_type, &
635 & 1, (/0/), aval, vinfo, ncname, &
636 & setparaccess = .false.)
637 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
640 vinfo( 2)=
'temperature coefficient for NH4 nitrification'
641 vinfo( 3)=
'Celsius-1'
642 status=def_var(ng, model, piofile, piovar, pio_type, &
643 & 1, (/0/), aval, vinfo, ncname, &
644 & setparaccess = .false.)
645 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
648 vinfo( 2)=
'temperature coefficient for PON decomposition to DON'
649 vinfo( 3)=
'Celsius-1'
650 status=def_var(ng, model, piofile, piovar, pio_type, &
651 & 1, (/0/), aval, vinfo, ncname, &
652 & setparaccess = .false.)
653 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
656 vinfo( 2)=
'temperature coefficient for PON decomposition to NH4'
657 vinfo( 3)=
'Celsius-1'
658 status=def_var(ng, model, piofile, piovar, pio_type, &
659 & 1, (/0/), aval, vinfo, ncname, &
660 & setparaccess = .false.)
661 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
664 vinfo( 2)=
'temperature coefficient for DON decomposition to NH4'
665 vinfo( 3)=
'Celsius-1'
666 status=def_var(ng, model, piofile, piovar, pio_type, &
667 & 1, (/0/), aval, vinfo, ncname, &
668 & setparaccess = .false.)
669 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
672 vinfo( 2)=
'temperature coefficient for opal decomposition to '// &
674 vinfo( 3)=
'Celsius-1'
675 status=def_var(ng, model, piofile, piovar, pio_type, &
676 & 1, (/0/), aval, vinfo, ncname, &
677 & setparaccess = .false.)
678 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
681 vinfo( 2)=
'Si:N ratio'
682 vinfo( 3)=
'millimole_Si millimole_N-1'
683 status=def_var(ng, model, piofile, piovar, pio_type, &
684 & 1, (/0/), aval, vinfo, ncname, &
685 & setparaccess = .false.)
686 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
689 vinfo( 2)=
'PON setting velocity'
690 vinfo( 3)=
'meter day-1'
691 status=def_var(ng, model, piofile, piovar, pio_type, &
692 & 1, (/0/), aval, vinfo, ncname, &
693 & setparaccess = .false.)
694 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
697 vinfo( 2)=
'opal setting velocity'
698 vinfo( 3)=
'meter day-1'
699 status=def_var(ng, model, piofile, piovar, pio_type, &
700 & 1, (/0/), aval, vinfo, ncname, &
701 & setparaccess = .false.)
702 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN