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 fennel et al. (2006) ecosystem model input parameters in **
10**
output netcdf 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, ncid, varid, nf90_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, ncid, varid, nf_type, &
29 & 1, (/0/), aval, vinfo, ncname, &
30 & setparaccess = .false.)
31 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
34 vinfo( 2)=
'light attenuation by chlorophyll'
35 vinfo( 3)=
'meter-2 milligram_Chl-1'
36 status=def_var(ng, model, ncid, varid, nf_type, &
37 & 1, (/0/), aval, vinfo, ncname, &
38 & setparaccess = .false.)
39 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
42 vinfo( 2)=
'photosynthetically available radiation fraction'
43 status=def_var(ng, model, ncid, varid, nf_type, &
44 & 1, (/0/), aval, vinfo, ncname, &
45 & setparaccess = .false.)
46 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
49 vinfo( 2)=
'Eppley temperature-limited growth parameter'
50 status=def_var(ng, model, ncid, varid, nf_type, &
51 & 1, (/0/), aval, vinfo, ncname, &
52 & setparaccess = .false.)
53 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
56 vinfo( 2)=
'radiation threshold for nitrification'
57 vinfo( 3)=
'watt meter-2'
58 status=def_var(ng, model, ncid, varid, nf_type, &
59 & 1, (/0/), aval, vinfo, ncname, &
60 & setparaccess = .false.)
61 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
64 vinfo( 2)=
'half-saturation radiation for nitrification'
65 vinfo( 3)=
'watt meter-2'
66 status=def_var(ng, model, ncid, varid, nf_type, &
67 & 1, (/0/), aval, vinfo, ncname, &
68 & setparaccess = .false.)
69 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
72 vinfo( 2)=
'nitrification rate'
74 status=def_var(ng, model, ncid, varid, nf_type, &
75 & 1, (/0/), aval, vinfo, ncname, &
76 & setparaccess = .false.)
77 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
80 vinfo( 2)=
'inverse half-saturation for phytoplankton NO3 uptake'
81 vinfo( 3)=
'meter3 millimole_N-1'
82 status=def_var(ng, model, ncid, varid, nf_type, &
83 & 1, (/0/), aval, vinfo, ncname, &
84 & setparaccess = .false.)
85 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
88 vinfo( 2)=
'inverse half-saturation for phytoplankton NH4 uptake'
89 vinfo( 3)=
'meter3 millimole_N-1'
90 status=def_var(ng, model, ncid, varid, nf_type, &
91 & 1, (/0/), aval, vinfo, ncname, &
92 & setparaccess = .false.)
93 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
96 vinfo( 2)=
'inverse half-saturation for phytoplankton PO4 uptake'
97 vinfo( 3)=
'meter3 millimole_P-1'
98 status=def_var(ng, model, ncid, varid, nf_type, &
99 & 1, (/0/), aval, vinfo, ncname, &
100 & setparaccess = .false.)
101 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
104 vinfo( 2)=
'zooplankton half-saturation constant for ingestion'
106 status=def_var(ng, model, ncid, varid, nf_type, &
107 & 1, (/0/), aval, vinfo, ncname, &
108 & setparaccess = .false.)
109 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
112 vinfo( 2)=
'maximum chlorophyll to carbon ratio'
113 vinfo( 3)=
'milligram_chl milligram_carbon-1'
114 status=def_var(ng, model, ncid, varid, nf_type, &
115 & 1, (/0/), aval, vinfo, ncname, &
116 & setparaccess = .false.)
117 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
120 vinfo( 2)=
'minimum chlorophyll threshold'
121 vinfo( 3)=
'milligram_chl meter-3'
122 status=def_var(ng, model, ncid, varid, nf_type, &
123 & 1, (/0/), aval, vinfo, ncname, &
124 & setparaccess = .false.)
125 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
128 vinfo( 2)=
'phytoplankton Carbon:Nitrogen ratio'
129 vinfo( 3)=
'mole_C mole_N-1'
130 status=def_var(ng, model, ncid, varid, nf_type, &
131 & 1, (/0/), aval, vinfo, ncname, &
132 & setparaccess = .false.)
133 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
136 vinfo( 2)=
'phytoplankton Phosphate:Nitrogen ratio'
137 vinfo( 3)=
'mole_P mole_N-1'
138 status=def_var(ng, model, ncid, varid, nf_type, &
139 & 1, (/0/), aval, vinfo, ncname, &
140 & setparaccess = .false.)
141 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
144 vinfo( 2)=
'phytoplankton, NH4 inhibition parameter'
145 vinfo( 3)=
'millimole_N-1'
146 status=def_var(ng, model, ncid, varid, nf_type, &
147 & 1, (/0/), aval, vinfo, ncname, &
148 & setparaccess = .false.)
149 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
152 vinfo( 2)=
'phytoplankton, initial slope of P-I curve'
153 vinfo( 3)=
'milligram_C milligram_Chl-1 watt-1 meter2 day-1'
154 status=def_var(ng, model, ncid, varid, nf_type, &
155 & 1, (/0/), aval, vinfo, ncname, &
156 & setparaccess = .false.)
157 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
160 vinfo( 2)=
'minimum phytoplankton threshold'
161 vinfo( 3)=
'millimole_N meter-3'
162 status=def_var(ng, model, ncid, varid, nf_type, &
163 & 1, (/0/), aval, vinfo, ncname, &
164 & setparaccess = .false.)
165 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
168 vinfo( 2)=
'phytoplankton mortality rate'
170 status=def_var(ng, model, ncid, varid, nf_type, &
171 & 1, (/0/), aval, vinfo, ncname, &
172 & setparaccess = .false.)
173 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
176 vinfo( 2)=
'zooplankton mitrogen assimilation efficiency'
177 status=def_var(ng, model, ncid, varid, nf_type, &
178 & 1, (/0/), aval, vinfo, ncname, &
179 & setparaccess = .false.)
180 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
183 vinfo( 2)=
'rate for zooplankton basal metabolism'
185 status=def_var(ng, model, ncid, varid, nf_type, &
186 & 1, (/0/), aval, vinfo, ncname, &
187 & setparaccess = .false.)
188 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
191 vinfo( 2)=
'zooplankton Carbon:Nitrogen ratio'
192 vinfo( 3)=
'mole_C mole_N-1'
193 status=def_var(ng, model, ncid, varid, nf_type, &
194 & 1, (/0/), aval, vinfo, ncname, &
195 & setparaccess = .false.)
196 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
199 vinfo( 2)=
'zooplankton specific excretion rate'
201 status=def_var(ng, model, ncid, varid, nf_type, &
202 & 1, (/0/), aval, vinfo, ncname, &
203 & setparaccess = .false.)
204 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
207 vinfo( 2)=
'zooplankton maximum growth rate'
209 status=def_var(ng, model, ncid, varid, nf_type, &
210 & 1, (/0/), aval, vinfo, ncname, &
211 & setparaccess = .false.)
212 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
215 vinfo( 2)=
'minimum zooplankton threshold'
216 vinfo( 3)=
'millimole_N meter-3'
217 status=def_var(ng, model, ncid, varid, nf_type, &
218 & 1, (/0/), aval, vinfo, ncname, &
219 & setparaccess = .false.)
220 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
223 vinfo( 2)=
'zooplankton mortality rate'
225 status=def_var(ng, model, ncid, varid, nf_type, &
226 & 1, (/0/), aval, vinfo, ncname, &
227 & setparaccess = .false.)
228 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
231 vinfo( 2)=
'rate of large detritus nitrogen re-mineralization'
233 status=def_var(ng, model, ncid, varid, nf_type, &
234 & 1, (/0/), aval, vinfo, ncname, &
235 & setparaccess = .false.)
236 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
239 vinfo( 2)=
'rate of large detritus carbon re-mineralization'
241 status=def_var(ng, model, ncid, varid, nf_type, &
242 & 1, (/0/), aval, vinfo, ncname, &
243 & setparaccess = .false.)
244 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
247 vinfo( 2)=
'coagulation rate'
249 status=def_var(ng, model, ncid, varid, nf_type, &
250 & 1, (/0/), aval, vinfo, ncname, &
251 & setparaccess = .false.)
252 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
255 vinfo( 2)=
'remineralization rate for small detritus nitrogen'
257 status=def_var(ng, model, ncid, varid, nf_type, &
258 & 1, (/0/), aval, vinfo, ncname, &
259 & setparaccess = .false.)
260 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
263 vinfo( 2)=
'remineralization rate for small detritus carbon'
265 status=def_var(ng, model, ncid, varid, nf_type, &
266 & 1, (/0/), aval, vinfo, ncname, &
267 & setparaccess = .false.)
268 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
271 vinfo( 2)=
'remineralization rate for river detritus nitrogen'
273 status=def_var(ng, model, ncid, varid, nf_type, &
274 & 1, (/0/), aval, vinfo, ncname, &
275 & setparaccess = .false.)
276 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
279 vinfo( 2)=
'remineralization rate for river detritus carbon'
281 status=def_var(ng, model, ncid, varid, nf_type, &
282 & 1, (/0/), aval, vinfo, ncname, &
283 & setparaccess = .false.)
284 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
287 vinfo( 2)=
'vertical sinking velocity for phytoplankton'
288 vinfo( 3)=
'meter day-1'
289 status=def_var(ng, model, ncid, varid, nf_type, &
290 & 1, (/0/), aval, vinfo, ncname, &
291 & setparaccess = .false.)
292 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
295 vinfo( 2)=
'vertical sinking velocity for large detritus'
296 vinfo( 3)=
'meter day-1'
297 status=def_var(ng, model, ncid, varid, nf_type, &
298 & 1, (/0/), aval, vinfo, ncname, &
299 & setparaccess = .false.)
300 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
303 vinfo( 2)=
'vertical sinking velocity for small detritus'
304 vinfo( 3)=
'meter day-1'
305 status=def_var(ng, model, ncid, varid, nf_type, &
306 & 1, (/0/), aval, vinfo, ncname, &
307 & setparaccess = .false.)
308 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
311 vinfo( 2)=
'partial pressure of CO2 in the air'
312 vinfo( 3)=
'parts per million by volume'
313 status=def_var(ng, model, ncid, varid, nf_type, &
314 & 1, (/0/), aval, vinfo, ncname, &
315 & setparaccess = .false.)
316 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN