ROMS
Loading...
Searching...
No Matches
nemuro_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 nemuro ecosystem model input parameters in output netcdf **
10** files. it is included in routine "def_info.F". **
11** **
12************************************************************************
13*/
14
15!
16! Define Nemuro ecosystem model parameters.
17!
18 vinfo( 1)='BioIter'
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
24
25 vinfo( 1)='AttSW'
26 vinfo( 2)='light attenuation by seawater'
27 vinfo( 3)='meter-1'
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
32
33 vinfo( 1)='AttPS'
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
41
42 vinfo( 1)='AttPL'
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
50
51 vinfo( 1)='PARfrac'
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
57
58 vinfo( 1)='AlphaPS'
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
65
66 vinfo( 1)='AlphaPL'
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
73
74 vinfo( 1)='BetaPS'
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
81
82 vinfo( 1)='BetaPL'
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
89
90 vinfo( 1)='VmaxS'
91 vinfo( 2)='small phytoplankton maximum photosynthetic rate '// &
92 & 'at 0 Celsius'
93 vinfo( 3)='day-1'
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
98
99 vinfo( 1)='VmaxL'
100 vinfo( 2)='large phytoplankton maximum photosynthetic rate '// &
101 & 'at 0 Celsius'
102 vinfo( 3)='day-1'
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
107
108 vinfo( 1)='K_NO3S'
109 vinfo( 2)='small phytoplankton half-saturation constant '// &
110 & 'for nitrate'
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
116
117 vinfo( 1)='K_NO3L'
118 vinfo( 2)='large phytoplankton half saturation constant '// &
119 & 'for nitrate'
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
125
126 vinfo( 1)='K_NH4S'
127 vinfo( 2)= &
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
134
135 vinfo( 1)='K_NH4L'
136 vinfo( 2)='large phytoplankton half saturation constant '// &
137 & 'for ammonium'
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
143
144 vinfo( 1)='K_SiL'
145 vinfo( 2)='large phytoplankton half saturation constant '// &
146 & 'for silicate'
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
152
153 vinfo( 1)='PusaiS'
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
160
161 vinfo( 1)='PusaiL'
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
168
169 vinfo( 1)='KGppS'
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
177
178 vinfo( 1)='KGppL'
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
186
187 vinfo( 1)='ResPS0'
188 vinfo( 2)='small phytoplankton respiration rate at 0 Celsius'
189 vinfo( 3)='day-1'
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
194
195 vinfo( 1)='ResPL0'
196 vinfo( 2)='large phytoplankton respiration rate at 0 Celsius'
197 vinfo( 3)='day-1'
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
202
203 vinfo( 1)='KResPS'
204 vinfo( 2)='small phytoplankton temperature coefficient '// &
205 & 'for respiration'
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
211
212 vinfo( 1)='KResPL'
213 vinfo( 2)='large phytoplankton temperature coefficient '// &
214 & 'for respiration'
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
220
221 vinfo( 1)='GammaS'
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
228
229 vinfo( 1)='GammaL'
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
236
237 vinfo( 1)='MorPS0'
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
244
245 vinfo( 1)='MorPL0'
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
252
253 vinfo( 1)='KMorPS'
254 vinfo( 2)='small phytoplankton temperature coefficient for '// &
255 & 'mortality'
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
261
262 vinfo( 1)='KMorPL'
263 vinfo( 2)='large phytoplankton temperature coefficient for '// &
264 & 'mortality'
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
270
271 vinfo( 1)='GRmaxSps'
272 vinfo( 2)='small zooplankton maximum grazing rate on '// &
273 & 'small phytoplankton at 0 Celsius'
274 vinfo( 3)='day-1'
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
279
280 vinfo( 1)='GRmaxLps'
281 vinfo( 2)='large zooplankton maximum grazing rate on '// &
282 & 'small phytoplankton at 0 Celsius'
283 vinfo( 3)='day-1'
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
288
289 vinfo( 1)='GRmaxLpl'
290 vinfo( 2)='large zooplankton maximum grazing rate on '// &
291 & 'large phytoplankton at 0 Celsius'
292 vinfo( 3)='day-1'
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
297
298 vinfo( 1)='GRmaxLzs'
299 vinfo( 2)='large zooplankton maximum grazing rate on '// &
300 & 'small zooplankton at 0 Celsius'
301 vinfo( 3)='day-1'
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
306
307 vinfo( 1)='GRmaxPpl'
308 vinfo( 2)='predator zooplankton maximum grazing rate on '// &
309 & 'large phytoplankton at 0 Celsius'
310 vinfo( 3)='day-1'
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
315
316 vinfo( 1)='GRmaxPzs'
317 vinfo( 2)='predator zooplankton maximum grazing rate on '// &
318 & 'small zooplankton at 0 Celsius'
319 vinfo( 3)='day-1'
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
324
325 vinfo( 1)='GRmaxPzl'
326 vinfo( 2)='predator zooplankton maximum grazing rate on '// &
327 & 'large zooplankton at 0 Celsius'
328 vinfo( 3)='day-1'
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
333
334 vinfo( 1)='KGraS'
335 vinfo( 2)='small zooplankton temperature coefficient '// &
336 & 'for grazing'
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
342
343 vinfo( 1)='KGraL'
344 vinfo( 2)='large zooplankton temperature coefficient '// &
345 & 'for grazing'
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
351
352 vinfo( 1)='KGraP'
353 vinfo( 2)='predator zooplankton temperature coefficient '// &
354 & 'for grazing'
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
360
361 vinfo( 1)='LamS'
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
368
369 vinfo( 1)='LamL'
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
376
377 vinfo( 1)='LamP'
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
384
385#ifdef HOLLING_GRAZING
386 vinfo( 1)='K_PS2ZS'
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
394
395 vinfo( 1)='K_PS2ZL'
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
403
404 vinfo( 1)='K_PL2ZL'
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
412
413 vinfo( 1)='K_ZS2ZL'
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
421
422 vinfo( 1)='K_PL2ZP'
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
430
431 vinfo( 1)='K_ZS2ZP'
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
439
440 vinfo( 1)='K_ZL2ZP'
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
448
449#else
450
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
459
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
468
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
477
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
486
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
495
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
504
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
513#endif
514
515 vinfo( 1)='PusaiPL'
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
523
524 vinfo( 1)='PusaiZS'
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
532
533 vinfo( 1)='MorZS0'
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
540
541 vinfo( 1)='MorZL0'
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
548
549 vinfo( 1)='MorZP0'
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
556
557 vinfo( 1)='AlphaZS'
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
563
564 vinfo( 1)='AlphaZL'
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
570
571 vinfo( 1)='AlphaZP'
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
577
578 vinfo( 1)='BetaZS'
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
584
585 vinfo( 1)='BetaZL'
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
591
592 vinfo( 1)='BetaZP'
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
598
599 vinfo( 1)='Nit0'
600 vinfo( 2)='NH4 nitrification rate at 0 Celsius'
601 vinfo( 3)='day-1'
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
606
607 vinfo( 1)='VP2N0'
608 vinfo( 2)='PON decomposition to NH4 rate at 0 Celsius'
609 vinfo( 3)='day-1'
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
614
615 vinfo( 1)='VP2D0'
616 vinfo( 2)='PON decomposition to DON rate at 0 Celsius'
617 vinfo( 3)='day-1'
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
622
623 vinfo( 1)='VD2N0'
624 vinfo( 2)='DON decomposition to NH4 rate at 0 Celsius'
625 vinfo( 3)='day-1'
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
630
631 vinfo( 1)='VO2S0'
632 vinfo( 2)='opal decomposition to silicate rate at 0 Celsius'
633 vinfo( 3)='day-1'
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
638
639 vinfo( 1)='KNit'
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
646
647 vinfo( 1)='KP2D'
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
654
655 vinfo( 1)='KP2N'
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
662
663 vinfo( 1)='KD2N'
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
670
671 vinfo( 1)='KO2S'
672 vinfo( 2)='temperature coefficient for opal decomposition to '// &
673 & 'silicate'
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
679
680 vinfo( 1)='RSiN'
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
687
688 vinfo( 1)='setVPON'
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
695
696 vinfo( 1)='setVOpal'
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
subroutine output(ng)
Definition output.F:4