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 npzd iron(fiechter, et al. 2009) ecosystem model input **
10** parameters in
output netcdf files. it is included in routine **
13************************************************************************
20 vinfo( 2)=
'number of iterations to achieve convergence'
21 status=def_var(ng, model, ncid, varid, nf90_int, &
22 & 1, (/0/), aval, vinfo, ncname, &
23 & setparaccess = .false.)
24 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
27 vinfo( 2)=
'photosynthetically available radiation fraction'
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 due to sea water'
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)=
'light attenuation due to phytoplankton'
43 vinfo( 3)=
'meter2 millimole_N-1'
44 status=def_var(ng, model, ncid, varid, nf_type, &
45 & 1, (/0/), aval, vinfo, ncname, &
46 & setparaccess = .false.)
47 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
50 vinfo( 2)=
'phytoplankton initial slope of the P-I curve'
51 vinfo( 3)=
'meter2 watt-1'
52 status=def_var(ng, model, ncid, varid, nf_type, &
53 & 1, (/0/), aval, vinfo, ncname, &
54 & setparaccess = .false.)
55 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
58 vinfo( 2)=
'nitrate uptake rate'
60 status=def_var(ng, model, ncid, varid, nf_type, &
61 & 1, (/0/), aval, vinfo, ncname, &
62 & setparaccess = .false.)
63 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
66 vinfo( 2)=
'phytoplankton mortality rate to the detritus pool'
68 status=def_var(ng, model, ncid, varid, nf_type, &
69 & 1, (/0/), aval, vinfo, ncname, &
70 & setparaccess = .false.)
71 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
74 vinfo( 2)=
'phytoplankton mortality rate to the Nitrogen pool'
76 status=def_var(ng, model, ncid, varid, nf_type, &
77 & 1, (/0/), aval, vinfo, ncname, &
78 & setparaccess = .false.)
79 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
82 vinfo( 2)=
'inverse half-saturation for phytoplankton NO3 uptake'
83 vinfo( 3)=
'meter3 millimole-1'
84 status=def_var(ng, model, ncid, varid, nf_type, &
85 & 1, (/0/), aval, vinfo, ncname, &
86 & setparaccess = .false.)
87 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
90 vinfo( 2)=
'Ivlev constant for zooplankton grazing'
91 status=def_var(ng, model, ncid, varid, nf_type, &
92 & 1, (/0/), aval, vinfo, ncname, &
93 & setparaccess = .false.)
94 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
97 vinfo( 2)=
'zooplankton grazing rate'
99 status=def_var(ng, model, ncid, varid, nf_type, &
100 & 1, (/0/), aval, vinfo, ncname, &
101 & setparaccess = .false.)
102 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
105 vinfo( 2)=
'zooplankton excretion efficiency to detritus pool'
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)=
'zooplankton excretion efficiency to Nitrogen pool'
113 status=def_var(ng, model, ncid, varid, nf_type, &
114 & 1, (/0/), aval, vinfo, ncname, &
115 & setparaccess = .false.)
116 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
119 vinfo( 2)=
'zooplankton mortality rate to the detritus pool'
121 status=def_var(ng, model, ncid, varid, nf_type, &
122 & 1, (/0/), aval, vinfo, ncname, &
123 & setparaccess = .false.)
124 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
127 vinfo( 2)=
'zootoplankton mortality rate to the Nitrogen pool'
129 status=def_var(ng, model, ncid, varid, nf_type, &
130 & 1, (/0/), aval, vinfo, ncname, &
131 & setparaccess = .false.)
132 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
135 vinfo( 2)=
'detritus remineralization rate'
137 status=def_var(ng, model, ncid, varid, nf_type, &
138 & 1, (/0/), aval, vinfo, ncname, &
139 & setparaccess = .false.)
140 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
143 vinfo( 2)=
'phytoplankton sinking rate'
145 status=def_var(ng, model, ncid, varid, nf_type, &
146 & 1, (/0/), aval, vinfo, ncname, &
147 & setparaccess = .false.)
148 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
151 vinfo( 2)=
'detrital sinking rate'
153 status=def_var(ng, model, ncid, varid, nf_type, &
154 & 1, (/0/), aval, vinfo, ncname, &
155 & setparaccess = .false.)
156 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
160 vinfo( 2)=
'iron uptake timescale'
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)=
'empirical Fe:C power'
169 status=def_var(ng, model, ncid, varid, nf_type, &
170 & 1, (/0/), aval, vinfo, ncname, &
171 & setparaccess = .false.)
172 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
175 vinfo( 2)=
'empirical Fe:C coefficient'
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)=
'Fe:C ratio'
184 vinfo( 3)=
'millimole_Fe mole_C-1'
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)=
'iron remineralization rate'
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
200 vinfo( 2)=
'minimum depth for dissolved iron relaxation'
202 status=def_var(ng, model, ncid, varid, nf_type, &
203 & 1, (/0/), aval, vinfo, ncname, &
204 & setparaccess = .false.)
205 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
207 vinfo( 1)=
'Fe_NudgTime'
208 vinfo( 2)=
'nudging time for dissolved iron relaxation'
210 status=def_var(ng, model, ncid, varid, nf_type, &
211 & 1, (/0/), aval, vinfo, ncname, &
212 & setparaccess = .false.)
213 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
216 vinfo( 2)=
'dissolved iron value for nudging'
217 vinfo( 3)=
'millimole_Fe meter-3'
218 status=def_var(ng, model, ncid, varid, nf_type, &
219 & 1, (/0/), aval, vinfo, ncname, &
220 & setparaccess = .false.)
221 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN