ROMS
Loading...
Searching...
No Matches
npzd_iron_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 npzd iron(fiechter, et al. 2009) ecosystem model input **
10** parameters in output netcdf files. it is included in routine **
11** "def_info.F". **
12** **
13************************************************************************
14*/
15
16!
17! Define NPZD iron (Fiechter, et al. 2009) ecosystem model parameters.
18!
19 vinfo( 1)='BioIter'
20 vinfo( 2)='number of iterations to achieve convergence'
21 status=def_var(ng, model, piofile, piovar, pio_int, &
22 & 1, (/0/), aval, vinfo, ncname, &
23 & setparaccess = .false.)
24 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
25
26 vinfo( 1)='PARfrac'
27 vinfo( 2)='photosynthetically available radiation fraction'
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)='AttSW'
34 vinfo( 2)='light attenuation due to sea water'
35 vinfo( 3)='meter-1'
36 status=def_var(ng, model, piofile, piovar, pio_type, &
37 & 1, (/0/), aval, vinfo, ncname, &
38 & setparaccess = .false.)
39 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
40
41 vinfo( 1)='AttPhy'
42 vinfo( 2)='light attenuation due to phytoplankton'
43 vinfo( 3)='meter2 millimole_N-1'
44 status=def_var(ng, model, piofile, piovar, pio_type, &
45 & 1, (/0/), aval, vinfo, ncname, &
46 & setparaccess = .false.)
47 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
48
49 vinfo( 1)='PhyIS'
50 vinfo( 2)='phytoplankton initial slope of the P-I curve'
51 vinfo( 3)='meter2 watt-1'
52 status=def_var(ng, model, piofile, piovar, pio_type, &
53 & 1, (/0/), aval, vinfo, ncname, &
54 & setparaccess = .false.)
55 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
56
57 vinfo( 1)='Vm_NO3'
58 vinfo( 2)='nitrate uptake rate'
59 vinfo( 3)='day-1'
60 status=def_var(ng, model, piofile, piovar, pio_type, &
61 & 1, (/0/), aval, vinfo, ncname, &
62 & setparaccess = .false.)
63 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
64
65 vinfo( 1)='PhyMRD'
66 vinfo( 2)='phytoplankton mortality rate to the detritus pool'
67 vinfo( 3)='day-1'
68 status=def_var(ng, model, piofile, piovar, pio_type, &
69 & 1, (/0/), aval, vinfo, ncname, &
70 & setparaccess = .false.)
71 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
72
73 vinfo( 1)='PhyMRN'
74 vinfo( 2)='phytoplankton mortality rate to the Nitrogen pool'
75 vinfo( 3)='day-1'
76 status=def_var(ng, model, piofile, piovar, pio_type, &
77 & 1, (/0/), aval, vinfo, ncname, &
78 & setparaccess = .false.)
79 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
80
81 vinfo( 1)='K_NO3'
82 vinfo( 2)='inverse half-saturation for phytoplankton NO3 uptake'
83 vinfo( 3)='meter3 millimole-1'
84 status=def_var(ng, model, piofile, piovar, pio_type, &
85 & 1, (/0/), aval, vinfo, ncname, &
86 & setparaccess = .false.)
87 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
88
89 vinfo( 1)='Ivlev'
90 vinfo( 2)='Ivlev constant for zooplankton grazing'
91 status=def_var(ng, model, piofile, piovar, pio_type, &
92 & 1, (/0/), aval, vinfo, ncname, &
93 & setparaccess = .false.)
94 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
95
96 vinfo( 1)='ZooGR'
97 vinfo( 2)='zooplankton grazing rate'
98 vinfo( 3)='day-1'
99 status=def_var(ng, model, piofile, piovar, pio_type, &
100 & 1, (/0/), aval, vinfo, ncname, &
101 & setparaccess = .false.)
102 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
103
104 vinfo( 1)='ZooEED'
105 vinfo( 2)='zooplankton excretion efficiency to detritus pool'
106 status=def_var(ng, model, piofile, piovar, pio_type, &
107 & 1, (/0/), aval, vinfo, ncname, &
108 & setparaccess = .false.)
109 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
110
111 vinfo( 1)='ZooEEN'
112 vinfo( 2)='zooplankton excretion efficiency to Nitrogen pool'
113 status=def_var(ng, model, piofile, piovar, pio_type, &
114 & 1, (/0/), aval, vinfo, ncname, &
115 & setparaccess = .false.)
116 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
117
118 vinfo( 1)='ZooMRD'
119 vinfo( 2)='zooplankton mortality rate to the detritus pool'
120 vinfo( 3)='day-1'
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)='ZooMRN'
127 vinfo( 2)='zootoplankton mortality rate to the Nitrogen pool'
128 vinfo( 3)='day-1'
129 status=def_var(ng, model, piofile, piovar, pio_type, &
130 & 1, (/0/), aval, vinfo, ncname, &
131 & setparaccess = .false.)
132 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
133
134 vinfo( 1)='DetRR'
135 vinfo( 2)='detritus remineralization rate'
136 vinfo( 3)='day-1'
137 status=def_var(ng, model, piofile, piovar, pio_type, &
138 & 1, (/0/), aval, vinfo, ncname, &
139 & setparaccess = .false.)
140 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
141
142 vinfo( 1)='wPhy'
143 vinfo( 2)='phytoplankton sinking rate'
144 vinfo( 3)='m day-1'
145 status=def_var(ng, model, piofile, piovar, pio_type, &
146 & 1, (/0/), aval, vinfo, ncname, &
147 & setparaccess = .false.)
148 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
149
150 vinfo( 1)='wDet'
151 vinfo( 2)='detrital sinking rate'
152 vinfo( 3)='m day-1'
153 status=def_var(ng, model, piofile, piovar, pio_type, &
154 & 1, (/0/), aval, vinfo, ncname, &
155 & setparaccess = .false.)
156 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
157
158#ifdef IRON_LIMIT
159 vinfo( 1)='T_Fe'
160 vinfo( 2)='iron uptake timescale'
161 vinfo( 3)='day'
162 status=def_var(ng, model, piofile, piovar, pio_type, &
163 & 1, (/0/), aval, vinfo, ncname, &
164 & setparaccess = .false.)
165 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
166
167 vinfo( 1)='A_Fe'
168 vinfo( 2)='empirical Fe:C power'
169 status=def_var(ng, model, piofile, piovar, pio_type, &
170 & 1, (/0/), aval, vinfo, ncname, &
171 & setparaccess = .false.)
172 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
173
174 vinfo( 1)='B_Fe'
175 vinfo( 2)='empirical Fe:C coefficient'
176 vinfo( 3)='mole_C-1'
177 status=def_var(ng, model, piofile, piovar, pio_type, &
178 & 1, (/0/), aval, vinfo, ncname, &
179 & setparaccess = .false.)
180 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
181
182 vinfo( 1)='K_FeC'
183 vinfo( 2)='Fe:C ratio'
184 vinfo( 3)='millimole_Fe mole_C-1'
185 status=def_var(ng, model, piofile, piovar, pio_type, &
186 & 1, (/0/), aval, vinfo, ncname, &
187 & setparaccess = .false.)
188 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
189
190 vinfo( 1)='FeRR'
191 vinfo( 2)='iron remineralization rate'
192 vinfo( 3)='day-1'
193 status=def_var(ng, model, piofile, piovar, pio_type, &
194 & 1, (/0/), aval, vinfo, ncname, &
195 & setparaccess = .false.)
196 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
197
198# ifdef IRON_RELAX
199 vinfo( 1)='Fe_Hmin'
200 vinfo( 2)='minimum depth for dissolved iron relaxation'
201 vinfo( 3)='meter'
202 status=def_var(ng, model, piofile, piovar, pio_type, &
203 & 1, (/0/), aval, vinfo, ncname, &
204 & setparaccess = .false.)
205 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
206
207 vinfo( 1)='Fe_NudgTime'
208 vinfo( 2)='nudging time for dissolved iron relaxation'
209 vinfo( 3)='day'
210 status=def_var(ng, model, piofile, piovar, pio_type, &
211 & 1, (/0/), aval, vinfo, ncname, &
212 & setparaccess = .false.)
213 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
214
215 vinfo( 1)='Fe_max'
216 vinfo( 2)='dissolved iron value for nudging'
217 vinfo( 3)='millimole_Fe meter-3'
218 status=def_var(ng, model, piofile, piovar, pio_type, &
219 & 1, (/0/), aval, vinfo, ncname, &
220 & setparaccess = .false.)
221 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
222# endif
223#endif
subroutine output(ng)
Definition output.F:4