Line data Source code
1 : !=======================================================================
2 : !
3 : ! Biogeochemistry driver
4 : !
5 : ! authors: Nicole Jeffery, LANL
6 : ! Scott Elliot, LANL
7 : ! Elizabeth C. Hunke, LANL
8 : !
9 : module icepack_zbgc
10 :
11 : use icepack_kinds
12 : use icepack_parameters, only: c0, c1, c2, p001, p1, p5, puny
13 : use icepack_parameters, only: depressT, rhosi, min_salin, salt_loss
14 : use icepack_parameters, only: fr_resp, algal_vel, R_dFe2dust, dustFe_sol, T_max
15 : use icepack_parameters, only: op_dep_min, fr_graze_s, fr_graze_e, fr_mort2min, fr_dFe
16 : use icepack_parameters, only: k_nitrif, t_iron_conv, max_loss, max_dfe_doc1
17 : use icepack_parameters, only: fr_resp_s, y_sk_DMS, t_sk_conv, t_sk_ox
18 : use icepack_parameters, only: scale_bgc, ktherm, skl_bgc, solve_zsal
19 : use icepack_parameters, only: z_tracers, fsal, conserv_check
20 :
21 : use icepack_tracers, only: nt_sice, nt_bgc_S, bio_index
22 : use icepack_tracers, only: tr_brine, nt_fbri, nt_qice, nt_Tsfc
23 : use icepack_tracers, only: nt_zbgc_frac
24 : use icepack_tracers, only: bio_index_o, bio_index
25 :
26 : use icepack_zbgc_shared, only: zbgc_init_frac
27 : use icepack_zbgc_shared, only: zbgc_frac_init
28 : use icepack_zbgc_shared, only: bgc_tracer_type, remap_zbgc
29 : use icepack_zbgc_shared, only: regrid_stationary
30 : use icepack_zbgc_shared, only: R_S2N, R_Si2N, R_Fe2C, R_Fe2N, R_Fe2DON, R_Fe2DOC
31 : use icepack_zbgc_shared, only: chlabs, alpha2max_low, beta2max
32 : use icepack_zbgc_shared, only: mu_max, grow_Tdep, fr_graze
33 : use icepack_zbgc_shared, only: mort_pre, mort_Tdep, k_exude
34 : use icepack_zbgc_shared, only: K_Nit, K_Am, K_Sil, K_Fe
35 : use icepack_zbgc_shared, only: f_don, kn_bac, f_don_Am, f_doc
36 : use icepack_zbgc_shared, only: f_exude, k_bac
37 : use icepack_zbgc_shared, only: tau_ret, tau_rel
38 : use icepack_zbgc_shared, only: R_C2N, R_CHL2N, f_abs_chl, R_C2N_DON
39 :
40 : use icepack_warnings, only: warnstr, icepack_warnings_add
41 : use icepack_warnings, only: icepack_warnings_setabort, icepack_warnings_aborted
42 :
43 : use icepack_brine, only: preflushing_changes, compute_microS_mushy
44 : use icepack_brine, only: update_hbrine, compute_microS
45 : use icepack_algae, only: zbio, sklbio
46 : use icepack_therm_shared, only: calculate_Tin_from_qin
47 : use icepack_itd, only: column_sum, column_conservation_check
48 : use icepack_zsalinity, only: zsalinity
49 :
50 : implicit none
51 :
52 : private
53 : public :: add_new_ice_bgc, &
54 : lateral_melt_bgc, &
55 : icepack_init_bgc, &
56 : icepack_init_zbgc, &
57 : icepack_biogeochemistry, &
58 : icepack_load_ocean_bio_array, &
59 : icepack_init_ocean_bio
60 :
61 : !=======================================================================
62 :
63 : contains
64 :
65 : !=======================================================================
66 :
67 : ! Adjust biogeochemical tracers when new frazil ice forms
68 :
69 97529376 : subroutine add_new_ice_bgc (dt, nblyr, &
70 : ncat, nilyr, nltrcr, &
71 97529376 : bgrid, cgrid, igrid, &
72 195058752 : aicen_init, vicen_init, vi0_init, &
73 97529376 : aicen, vicen, vsnon1, &
74 : vi0new, &
75 97529376 : ntrcr, trcrn, nbtrcr, &
76 97529376 : sss, ocean_bio, flux_bio, &
77 : hsurp)
78 :
79 : integer (kind=int_kind), intent(in) :: &
80 : nblyr , & ! number of bio layers
81 : ncat , & ! number of thickness categories
82 : nilyr , & ! number of ice layers
83 : nltrcr, & ! number of zbgc tracers
84 : nbtrcr , & ! number of biology tracers
85 : ntrcr ! number of tracers in use
86 :
87 : real (kind=dbl_kind), dimension (nblyr+2), intent(in) :: &
88 : bgrid ! biology nondimensional vertical grid points
89 :
90 : real (kind=dbl_kind), dimension (nblyr+1), intent(in) :: &
91 : igrid ! biology vertical interface points
92 :
93 : real (kind=dbl_kind), dimension (nilyr+1), intent(in) :: &
94 : cgrid ! CICE vertical coordinate
95 :
96 : real (kind=dbl_kind), intent(in) :: &
97 : dt ! time step (s)
98 :
99 : real (kind=dbl_kind), dimension (:), &
100 : intent(in) :: &
101 : aicen_init , & ! initial concentration of ice
102 : vicen_init , & ! intiial volume per unit area of ice (m)
103 : aicen , & ! concentration of ice
104 : vicen ! volume per unit area of ice (m)
105 :
106 : real (kind=dbl_kind), intent(in) :: &
107 : vsnon1 ! category 1 snow volume per unit area (m)
108 :
109 : real (kind=dbl_kind), dimension (:,:), &
110 : intent(inout) :: &
111 : trcrn ! ice tracers
112 :
113 : real (kind=dbl_kind), intent(in) :: &
114 : sss !sea surface salinity (ppt)
115 :
116 : real (kind=dbl_kind), intent(in) :: &
117 : vi0_init , & ! volume of new ice added to cat 1 (intial)
118 : vi0new ! volume of new ice added to cat 1
119 :
120 : real (kind=dbl_kind), intent(in) :: &
121 : hsurp ! thickness of new ice added to each cat
122 :
123 : real (kind=dbl_kind), dimension (:), &
124 : intent(inout) :: &
125 : flux_bio ! tracer flux to ocean from biology (mmol/m^2/s)
126 :
127 : real (kind=dbl_kind), dimension (:), &
128 : intent(in) :: &
129 : ocean_bio ! ocean concentration of biological tracer
130 :
131 : ! local
132 :
133 : integer (kind=int_kind) :: &
134 : location , & ! 1 (add frazil to bottom), 0 (add frazil throughout)
135 : n , & ! ice category index
136 : k ! ice layer index
137 :
138 : real (kind=dbl_kind) :: &
139 384288 : vbri1 , & ! starting volume of existing brine
140 384288 : vbri_init , & ! brine volume summed over categories
141 384288 : vbri_final ! brine volume summed over categories
142 :
143 : real (kind=dbl_kind) :: &
144 384288 : vsurp , & ! volume of new ice added to each cat
145 384288 : vtmp ! total volume of new and old ice
146 :
147 : real (kind=dbl_kind), dimension (ncat) :: &
148 99835104 : vbrin ! trcrn(nt_fbri,n)*vicen(n)
149 :
150 : real (kind=dbl_kind) :: &
151 384288 : vice_new ! vicen_init + vsurp
152 :
153 : real (kind=dbl_kind) :: &
154 384288 : Tmlts ! melting temperature (oC)
155 :
156 : character (len=char_len) :: &
157 : fieldid ! field identifier
158 :
159 : character(len=*),parameter :: subname='(add_new_ice_bgc)'
160 :
161 : !-----------------------------------------------------------------
162 : ! brine
163 : !-----------------------------------------------------------------
164 585176256 : vbrin(:) = c0
165 585176256 : do n = 1, ncat
166 487646880 : vbrin(n) = vicen_init(n)
167 585176256 : if (tr_brine) vbrin(n) = trcrn(nt_fbri,n)*vicen_init(n)
168 : enddo
169 :
170 97529376 : call column_sum (ncat, vbrin, vbri_init)
171 97529376 : if (icepack_warnings_aborted(subname)) return
172 :
173 97529376 : vbri_init = vbri_init + vi0_init
174 1804869888 : do k = 1, nbtrcr
175 6725040 : flux_bio(k) = flux_bio(k) &
176 1804869888 : - vi0_init/dt*ocean_bio(k)*zbgc_init_frac(k)
177 : enddo
178 : !-----------------------------------------------------------------
179 : ! Distribute bgc in new ice volume among all ice categories by
180 : ! increasing ice thickness, leaving ice area unchanged.
181 : !-----------------------------------------------------------------
182 :
183 : ! Diffuse_bio handles concentration changes from ice growth/melt
184 : ! ice area does not change
185 : ! add salt to the bottom , location = 1
186 :
187 97529376 : vsurp = c0
188 97529376 : vtmp = c0
189 :
190 585176256 : do n = 1,ncat
191 :
192 585176256 : if (hsurp > c0) then
193 :
194 36835090 : vtmp = vbrin(n)
195 36835090 : vsurp = hsurp * aicen_init(n)
196 36835090 : vbrin(n) = vbrin(n) + vsurp
197 36835090 : vice_new = vicen_init(n) + vsurp
198 36835090 : if (tr_brine .and. vicen(n) > c0) then
199 36770030 : trcrn(nt_fbri,n) = vbrin(n)/vicen(n)
200 65060 : elseif (tr_brine .and. vicen(n) <= c0) then
201 65060 : trcrn(nt_fbri,n) = c1
202 : endif
203 :
204 36835090 : if (nltrcr > 0) then
205 18468505 : location = 1
206 : call adjust_tracer_profile(nbtrcr, dt, ntrcr, &
207 50960 : aicen_init(n), &
208 50960 : vbrin(n), &
209 : vice_new, &
210 0 : trcrn(:,n), &
211 : vtmp, &
212 : vsurp, sss, &
213 : nilyr, nblyr, &
214 : solve_zsal, bgrid, &
215 : cgrid, &
216 0 : ocean_bio, igrid, &
217 18468505 : location)
218 18468505 : if (icepack_warnings_aborted(subname)) return
219 : endif ! nltrcr
220 : endif ! hsurp > 0
221 : enddo ! n
222 :
223 : !-----------------------------------------------------------------
224 : ! Combine bgc in new ice grown in open water with category 1 ice.
225 : !-----------------------------------------------------------------
226 :
227 97529376 : if (vi0new > c0) then
228 :
229 9623004 : vbri1 = vbrin(1)
230 9623004 : vbrin(1) = vbrin(1) + vi0new
231 9623004 : if (tr_brine .and. vicen(1) > c0) then
232 9623004 : trcrn(nt_fbri,1) = vbrin(1)/vicen(1)
233 0 : elseif (tr_brine .and. vicen(1) <= c0) then
234 0 : trcrn(nt_fbri,1) = c1
235 : endif
236 :
237 : ! Diffuse_bio handles concentration changes from ice growth/melt
238 : ! ice area changes
239 : ! add salt throughout, location = 0
240 :
241 9623004 : if (nltrcr > 0) then
242 4829547 : location = 0
243 : call adjust_tracer_profile(nbtrcr, dt, ntrcr, &
244 18045 : aicen(1), &
245 18045 : vbrin(1), &
246 18045 : vicen(1), &
247 0 : trcrn(:,1), &
248 : vbri1, &
249 : vi0new, sss, &
250 : nilyr, nblyr, &
251 : solve_zsal, bgrid, &
252 : cgrid, &
253 0 : ocean_bio, igrid, &
254 4829547 : location)
255 4829547 : if (icepack_warnings_aborted(subname)) return
256 :
257 4829547 : if (solve_zsal .and. vsnon1 .le. c0) then
258 0 : Tmlts = -trcrn(nt_sice,1)*depressT
259 0 : trcrn(nt_Tsfc,1) = calculate_Tin_from_qin(trcrn(nt_qice,1),Tmlts)
260 : endif ! solve_zsal
261 : endif ! nltrcr > 0
262 : endif ! vi0new > 0
263 :
264 97529376 : if (tr_brine .and. conserv_check) then
265 0 : call column_sum (ncat, vbrin, vbri_final)
266 0 : if (icepack_warnings_aborted(subname)) return
267 :
268 0 : fieldid = subname//':vbrin'
269 : call column_conservation_check (fieldid, &
270 : vbri_init, vbri_final, &
271 0 : puny)
272 0 : if (icepack_warnings_aborted(subname)) return
273 :
274 : endif ! conserv_check
275 :
276 : end subroutine add_new_ice_bgc
277 :
278 : !=======================================================================
279 :
280 : ! When sea ice melts laterally, flux bgc to ocean
281 :
282 1480757 : subroutine lateral_melt_bgc (dt, &
283 : ncat, nblyr, &
284 1480757 : rside, vicen, &
285 1480757 : trcrn, fzsal, &
286 1480757 : flux_bio, nbltrcr)
287 :
288 : integer (kind=int_kind), intent(in) :: &
289 : ncat , & ! number of thickness categories
290 : nblyr , & ! number of bio layers
291 : nbltrcr ! number of biology tracers
292 :
293 : real (kind=dbl_kind), intent(in) :: &
294 : dt ! time step (s)
295 :
296 : real (kind=dbl_kind), dimension(:), intent(in) :: &
297 : vicen ! volume per unit area of ice (m)
298 :
299 : real (kind=dbl_kind), dimension (:,:), intent(in) :: &
300 : trcrn ! tracer array
301 :
302 : real (kind=dbl_kind), intent(in) :: &
303 : rside ! fraction of ice that melts laterally
304 :
305 : real (kind=dbl_kind), intent(inout) :: &
306 : fzsal ! salt flux from layer Salinity (kg/m^2/s)
307 :
308 : real (kind=dbl_kind), dimension(:), intent(inout) :: &
309 : flux_bio ! biology tracer flux from layer bgc (mmol/m^2/s)
310 :
311 : ! local variables
312 :
313 : integer (kind=int_kind) :: &
314 : k , & ! layer index
315 : m , & !
316 : n ! category index
317 :
318 : real (kind=dbl_kind) :: &
319 2665 : zspace ! bio grid spacing
320 :
321 : character(len=*),parameter :: subname='(lateral_melt_bgc)'
322 :
323 1480757 : zspace = c1/(real(nblyr,kind=dbl_kind))
324 :
325 1480757 : if (solve_zsal) then
326 0 : do n = 1, ncat
327 0 : do k = 1,nblyr
328 0 : fzsal = fzsal + rhosi*trcrn(nt_fbri,n) &
329 0 : * vicen(n)*p001*zspace*trcrn(nt_bgc_S+k-1,n) &
330 0 : * rside/dt
331 : enddo
332 : enddo
333 : endif
334 :
335 29615140 : do m = 1, nbltrcr
336 170287055 : do n = 1, ncat
337 1294181618 : do k = 1, nblyr+1
338 2025400 : flux_bio(m) = flux_bio(m) + trcrn(nt_fbri,n) &
339 4050800 : * vicen(n)*zspace*trcrn(bio_index(m)+k-1,n) &
340 1268072635 : * rside/dt
341 : enddo
342 : enddo
343 : enddo
344 :
345 1480757 : end subroutine lateral_melt_bgc
346 :
347 : !=======================================================================
348 : !
349 : ! Add new ice tracers to the ice bottom and adjust the vertical profile
350 : !
351 : ! author: Nicole Jeffery, LANL
352 :
353 23298052 : subroutine adjust_tracer_profile (nbtrcr, dt, ntrcr, &
354 : aicen, vbrin, &
355 23298052 : vicen, trcrn, &
356 : vtmp, &
357 : vsurp, sss, &
358 : nilyr, nblyr, &
359 23298052 : solve_zsal, bgrid, &
360 46596104 : cgrid, ocean_bio, &
361 23298052 : igrid, location)
362 :
363 : integer (kind=int_kind), intent(in) :: &
364 : location , & ! 1 (add frazil to bottom), 0 (add frazil throughout)
365 : ntrcr , & ! number of tracers in use
366 : nilyr , & ! number of ice layers
367 : nbtrcr , & ! number of biology tracers
368 : nblyr ! number of biology layers
369 :
370 : real (kind=dbl_kind), intent(in) :: &
371 : dt ! timestep (s)
372 :
373 : real (kind=dbl_kind), intent(in) :: &
374 : aicen , & ! concentration of ice
375 : vicen , & ! volume of ice
376 : sss , & ! ocean salinity (ppt)
377 : ! hsurp , & ! flags new ice added to each cat
378 : vsurp , & ! volume of new ice added to each cat
379 : vtmp ! total volume of new and old ice
380 :
381 : real (kind=dbl_kind), dimension (nbtrcr), intent(in) :: &
382 : ocean_bio
383 :
384 : real (kind=dbl_kind), intent(in) :: &
385 : vbrin ! fbri*volume per unit area of ice (m)
386 :
387 : logical (kind=log_kind), intent(in) :: &
388 : solve_zsal
389 :
390 : real (kind=dbl_kind), dimension (nblyr+1), intent(in) :: &
391 : igrid ! zbio grid
392 :
393 : real (kind=dbl_kind), dimension (nblyr+2), intent(in) :: &
394 : bgrid ! zsal grid
395 :
396 : real (kind=dbl_kind), dimension (nilyr+1), intent(in) :: &
397 : cgrid ! CICE grid
398 :
399 : real (kind=dbl_kind), dimension (ntrcr), &
400 : intent(inout) :: &
401 : trcrn ! ice tracers
402 :
403 : ! local variables
404 :
405 : real (kind=dbl_kind), dimension (ntrcr+2) :: &
406 62812279 : trtmp0, & ! temporary, remapped tracers
407 62812279 : trtmp ! temporary, remapped tracers
408 :
409 : real (kind=dbl_kind) :: &
410 69005 : hin , & ! ice height
411 69005 : hinS_new, & ! brine height
412 69005 : temp_S
413 :
414 : integer (kind=int_kind) :: &
415 : k, m
416 :
417 : real (kind=dbl_kind), dimension (nblyr+1) :: &
418 47217149 : C_stationary ! stationary bulk concentration*h (mmol/m^2)
419 :
420 : real (kind=dbl_kind), dimension (nblyr) :: &
421 23919097 : S_stationary ! stationary bulk concentration*h (ppt*m)
422 :
423 : real(kind=dbl_kind) :: &
424 69005 : top_conc , & ! salinity or bgc ocean concentration of frazil
425 69005 : fluxb , & ! needed for regrid (set to zero here)
426 69005 : hbri_old , & ! previous timestep brine height
427 69005 : hbri ! brine height
428 :
429 : character(len=*),parameter :: subname='(adjust_tracer_profile)'
430 :
431 5521638324 : trtmp0(:) = c0
432 5521638324 : trtmp(:) = c0
433 23298052 : fluxb = c0
434 :
435 23298052 : if (location == 1 .and. vbrin > c0) then ! add frazil to bottom
436 :
437 18435899 : hbri = vbrin
438 18435899 : hbri_old = vtmp
439 18435899 : if (solve_zsal) then
440 0 : top_conc = sss * salt_loss
441 0 : do k = 1, nblyr
442 0 : S_stationary(k) = trcrn(nt_bgc_S+k-1)* hbri_old
443 : enddo
444 : call regrid_stationary (S_stationary, hbri_old, &
445 : hbri, dt, &
446 : ntrcr, &
447 : nblyr-1, top_conc, &
448 0 : bgrid(2:nblyr+1), fluxb )
449 0 : if (icepack_warnings_aborted(subname)) return
450 0 : do k = 1, nblyr
451 0 : trcrn(nt_bgc_S+k-1) = S_stationary(k)/hbri
452 0 : trtmp0(nt_sice+k-1) = trcrn(nt_bgc_S+k-1)
453 : enddo
454 : endif ! solve_zsal
455 :
456 368717980 : do m = 1, nbtrcr
457 350282081 : top_conc = ocean_bio(m)*zbgc_init_frac(m)
458 3152538729 : do k = 1, nblyr+1
459 3152538729 : C_stationary(k) = trcrn(bio_index(m) + k-1)* hbri_old
460 : enddo !k
461 : call regrid_stationary (C_stationary, hbri_old, &
462 : hbri, dt, &
463 : ntrcr, &
464 : nblyr, top_conc, &
465 350282081 : igrid, fluxb )
466 350282081 : if (icepack_warnings_aborted(subname)) return
467 3170974628 : do k = 1, nblyr+1
468 3152538729 : trcrn(bio_index(m) + k-1) = C_stationary(k)/hbri
469 : enddo !k
470 : enddo !m
471 :
472 18435899 : if (solve_zsal) then
473 0 : if (aicen > c0) then
474 0 : hinS_new = vbrin/aicen
475 0 : hin = vicen/aicen
476 : else
477 0 : hinS_new = c0
478 0 : hin = c0
479 : endif ! aicen
480 0 : temp_S = min_salin ! bio to cice
481 : call remap_zbgc(nilyr, &
482 : nt_sice, &
483 0 : trtmp0(1:ntrcr), trtmp, &
484 : 1, nblyr, &
485 : hin, hinS_new, &
486 0 : cgrid(2:nilyr+1), &
487 0 : bgrid(2:nblyr+1), temp_S )
488 0 : if (icepack_warnings_aborted(subname)) return
489 18435899 : do k = 1, nilyr
490 0 : trcrn(nt_sice+k-1) = trtmp(nt_sice+k-1)
491 : enddo ! k
492 : endif ! solve_zsal
493 :
494 4862153 : elseif (vbrin > c0) then ! add frazil throughout location == 0 .and.
495 :
496 43465923 : do k = 1, nblyr+1
497 38636376 : if (solve_zsal .and. k < nblyr + 1) then
498 0 : trcrn(nt_bgc_S+k-1) = (trcrn(nt_bgc_S+k-1) * vtmp &
499 0 : + sss*salt_loss * vsurp) / vbrin
500 0 : trtmp0(nt_sice+k-1) = trcrn(nt_bgc_S+k-1)
501 : endif ! solve_zsal
502 777557067 : do m = 1, nbtrcr
503 8228520 : trcrn(bio_index(m) + k-1) = (trcrn(bio_index(m) + k-1) * vtmp &
504 778213200 : + ocean_bio(m)*zbgc_init_frac(m) * vsurp) / vbrin
505 : enddo
506 : enddo
507 :
508 4829547 : if (solve_zsal) then
509 0 : if (aicen > c0) then
510 0 : hinS_new = vbrin/aicen
511 0 : hin = vicen/aicen
512 : else
513 0 : hinS_new = c0
514 0 : hin = c0
515 : endif !aicen
516 0 : temp_S = min_salin ! bio to cice
517 : call remap_zbgc(nilyr, &
518 : nt_sice, &
519 0 : trtmp0(1:ntrcr), trtmp, &
520 : 1, nblyr, &
521 : hin, hinS_new, &
522 0 : cgrid(2:nilyr+1), &
523 0 : bgrid(2:nblyr+1),temp_S )
524 0 : if (icepack_warnings_aborted(subname)) return
525 0 : do k = 1, nilyr
526 0 : trcrn(nt_sice+k-1) = trtmp(nt_sice+k-1)
527 : enddo !k
528 : endif ! solve_zsal
529 :
530 : endif ! location
531 :
532 : end subroutine adjust_tracer_profile
533 :
534 : !=======================================================================
535 : !autodocument_start icepack_init_bgc
536 : !
537 :
538 429185 : subroutine icepack_init_bgc(ncat, nblyr, nilyr, ntrcr_o, &
539 429185 : cgrid, igrid, ntrcr, nbtrcr, &
540 429185 : sicen, trcrn, sss, ocean_bio_all)
541 :
542 : integer (kind=int_kind), intent(in) :: &
543 : ncat , & ! number of thickness categories
544 : nilyr , & ! number of ice layers
545 : nblyr , & ! number of bio layers
546 : ntrcr_o,& ! number of tracers not including bgc
547 : ntrcr , & ! number of tracers in use
548 : nbtrcr ! number of bio tracers in use
549 :
550 : real (kind=dbl_kind), dimension (nblyr+1), intent(inout) :: &
551 : igrid ! biology vertical interface points
552 :
553 : real (kind=dbl_kind), dimension (nilyr+1), intent(inout) :: &
554 : cgrid ! CICE vertical coordinate
555 :
556 : real (kind=dbl_kind), dimension(nilyr, ncat), intent(in) :: &
557 : sicen ! salinity on the cice grid
558 :
559 : real (kind=dbl_kind), dimension (:,:), intent(inout) :: &
560 : trcrn ! subset of tracer array (only bgc)
561 :
562 : real (kind=dbl_kind), intent(in) :: &
563 : sss ! sea surface salinity (ppt)
564 :
565 : real (kind=dbl_kind), dimension (:), intent(inout) :: &
566 : ocean_bio_all ! fixed order, all values even for tracers false
567 :
568 : !autodocument_end
569 :
570 : ! local variables
571 :
572 : integer (kind=int_kind) :: &
573 : k , & ! vertical index
574 : n , & ! category index
575 : mm ! bio tracer index
576 :
577 : real (kind=dbl_kind), dimension (ntrcr+2) :: &
578 3654420 : trtmp ! temporary, remapped tracers
579 :
580 : character(len=*),parameter :: subname='(icepack_init_bgc)'
581 :
582 : !-----------------------------------------------------------------------------
583 : ! Skeletal Layer Model
584 : ! All bgc tracers are Bulk quantities in units of mmol or mg per m^3
585 : ! The skeletal layer model assumes a constant
586 : ! layer depth (sk_l) and porosity (phi_sk)
587 : !-----------------------------------------------------------------------------
588 429185 : if (skl_bgc) then
589 :
590 1219260 : do n = 1,ncat
591 17476060 : do mm = 1,nbtrcr
592 : ! bulk concentration (mmol or mg per m^3, or 10^-3 mmol/m^3)
593 17272850 : trcrn(bio_index(mm)-ntrcr_o, n) = ocean_bio_all(bio_index_o(mm))
594 : enddo ! nbtrcr
595 : enddo ! n
596 :
597 : !-----------------------------------------------------------------------------
598 : ! zbgc Model
599 : ! All bgc tracers are Bulk quantities in units of mmol or mg per m^3
600 : ! The vertical layer model uses prognosed porosity and layer depth
601 : !-----------------------------------------------------------------------------
602 :
603 : else ! not skl_bgc
604 :
605 225975 : if (scale_bgc .and. solve_zsal) then ! bulk concentration (mmol or mg per m^3)
606 0 : do n = 1,ncat
607 0 : do mm = 1,nbtrcr
608 0 : do k = 2, nblyr
609 0 : trcrn(bio_index(mm)+k-1-ntrcr_o,n) = &
610 0 : (p5*(trcrn(nt_bgc_S+k-1-ntrcr_o,n)+ trcrn(nt_bgc_S+k-2-ntrcr_o,n)) &
611 0 : / sss*ocean_bio_all(bio_index_o(mm)))
612 : enddo !k
613 0 : trcrn(nt_zbgc_frac-1+mm-ntrcr_o,n) = zbgc_frac_init(mm)
614 0 : trcrn(bio_index(mm)-ntrcr_o,n) = (trcrn(nt_bgc_S-ntrcr_o,n) &
615 0 : / sss*ocean_bio_all(bio_index_o(mm)))
616 0 : trcrn(bio_index(mm)+nblyr-ntrcr_o,n) = (trcrn(nt_bgc_S+nblyr-1-ntrcr_o,n) &
617 0 : / sss*ocean_bio_all(bio_index_o(mm)))
618 0 : trcrn(bio_index(mm)+nblyr+1-ntrcr_o:bio_index(mm)+nblyr+2-ntrcr_o,n) = c0 ! snow
619 : enddo ! mm
620 : enddo ! n
621 :
622 225975 : elseif (scale_bgc .and. ktherm == 2) then
623 53556075 : trtmp(:) = c0
624 1581825 : do n = 1,ncat
625 : call remap_zbgc(nilyr, &
626 : 1, &
627 0 : sicen(:,n), trtmp, &
628 : 0, nblyr+1, &
629 : c1, c1, &
630 0 : cgrid(2:nilyr+1), &
631 0 : igrid(1:nblyr+1), &
632 1129875 : sicen(1,n) )
633 1129875 : if (icepack_warnings_aborted(subname)) return
634 :
635 22823475 : do mm = 1,nbtrcr
636 194338500 : do k = 1, nblyr + 1
637 17411600 : trcrn(bio_index(mm)+k-1-ntrcr_o,n) = &
638 180446800 : (trtmp(k)/sss*ocean_bio_all(bio_index_o(mm)))
639 536690625 : trcrn(bio_index(mm)+nblyr+1-ntrcr_o:bio_index(mm)+nblyr+2-ntrcr_o,n) = c0 ! snow
640 : enddo ! k
641 : enddo ! mm
642 : enddo ! n
643 :
644 0 : elseif (nbtrcr > 0 .and. nt_fbri > 0) then ! not scale_bgc
645 :
646 0 : do n = 1,ncat
647 0 : do mm = 1,nbtrcr
648 0 : do k = 1, nblyr+1
649 0 : trcrn(bio_index(mm)+k-1-ntrcr_o,n) = ocean_bio_all(bio_index_o(mm)) &
650 0 : * zbgc_init_frac(mm)
651 0 : trcrn(bio_index(mm)+nblyr+1-ntrcr_o:bio_index(mm)+nblyr+2-ntrcr_o,n) = c0 ! snow
652 : enddo ! k
653 0 : trcrn(nt_zbgc_frac-1+mm-ntrcr_o,n) = zbgc_frac_init(mm)
654 : enddo ! mm
655 : enddo ! n
656 :
657 : endif ! scale_bgc
658 : endif ! skl_bgc
659 :
660 : end subroutine icepack_init_bgc
661 :
662 : !=======================================================================
663 : !autodocument_start icepack_init_zbgc
664 : !
665 :
666 5600 : subroutine icepack_init_zbgc ( &
667 11200 : R_Si2N_in, R_S2N_in, R_Fe2C_in, R_Fe2N_in, R_C2N_in, R_C2N_DON_in, &
668 11200 : R_chl2N_in, F_abs_chl_in, R_Fe2DON_in, R_Fe2DOC_in, chlabs_in, &
669 11200 : alpha2max_low_in, beta2max_in, mu_max_in, fr_graze_in, mort_pre_in, &
670 5600 : mort_Tdep_in, k_exude_in, K_Nit_in, K_Am_in, K_sil_in, K_Fe_in, &
671 5600 : f_don_in, kn_bac_in, f_don_Am_in, f_doc_in, f_exude_in, k_bac_in, &
672 11200 : grow_Tdep_in, zbgc_frac_init_in, &
673 5600 : zbgc_init_frac_in, tau_ret_in, tau_rel_in, bgc_tracer_type_in, &
674 : fr_resp_in, algal_vel_in, R_dFe2dust_in, dustFe_sol_in, T_max_in, &
675 : op_dep_min_in, fr_graze_s_in, fr_graze_e_in, fr_mort2min_in, fr_dFe_in, &
676 : k_nitrif_in, t_iron_conv_in, max_loss_in, max_dfe_doc1_in, &
677 : fr_resp_s_in, y_sk_DMS_in, t_sk_conv_in, t_sk_ox_in, fsal_in)
678 :
679 : real (kind=dbl_kind), optional :: R_C2N_in(:) ! algal C to N (mole/mole)
680 : real (kind=dbl_kind), optional :: R_chl2N_in(:) ! 3 algal chlorophyll to N (mg/mmol)
681 : real (kind=dbl_kind), optional :: F_abs_chl_in(:) ! to scale absorption in Dedd
682 : real (kind=dbl_kind), optional :: R_C2N_DON_in(:) ! increase compare to algal R_Fe2C
683 : real (kind=dbl_kind), optional :: R_Si2N_in(:) ! algal Sil to N (mole/mole)
684 : real (kind=dbl_kind), optional :: R_S2N_in(:) ! algal S to N (mole/mole)
685 : real (kind=dbl_kind), optional :: R_Fe2C_in(:) ! algal Fe to carbon (umol/mmol)
686 : real (kind=dbl_kind), optional :: R_Fe2N_in(:) ! algal Fe to N (umol/mmol)
687 : real (kind=dbl_kind), optional :: R_Fe2DON_in(:) ! Fe to N of DON (nmol/umol)
688 : real (kind=dbl_kind), optional :: R_Fe2DOC_in(:) ! Fe to C of DOC (nmol/umol)
689 :
690 : real (kind=dbl_kind), optional :: fr_resp_in ! frac of algal growth lost due to respiration
691 : real (kind=dbl_kind), optional :: algal_vel_in ! 0.5 cm/d(m/s) Lavoie 2005 1.5 cm/day
692 : real (kind=dbl_kind), optional :: R_dFe2dust_in ! g/g (3.5% content) Tagliabue 2009
693 : real (kind=dbl_kind), optional :: dustFe_sol_in ! solubility fraction
694 : real (kind=dbl_kind), optional :: T_max_in ! maximum temperature (C)
695 : real (kind=dbl_kind), optional :: op_dep_min_in ! Light attenuates for optical depths exceeding min
696 : real (kind=dbl_kind), optional :: fr_graze_s_in ! fraction of grazing spilled or slopped
697 : real (kind=dbl_kind), optional :: fr_graze_e_in ! fraction of assimilation excreted
698 : real (kind=dbl_kind), optional :: fr_mort2min_in ! fractionation of mortality to Am
699 : real (kind=dbl_kind), optional :: fr_dFe_in ! fraction of remineralized nitrogen
700 : ! (in units of algal iron)
701 : real (kind=dbl_kind), optional :: k_nitrif_in ! nitrification rate (1/day)
702 : real (kind=dbl_kind), optional :: t_iron_conv_in ! desorption loss pFe to dFe (day)
703 : real (kind=dbl_kind), optional :: max_loss_in ! restrict uptake to % of remaining value
704 : real (kind=dbl_kind), optional :: max_dfe_doc1_in ! max ratio of dFe to saccharides in the ice (nM Fe/muM C)
705 : real (kind=dbl_kind), optional :: fr_resp_s_in ! DMSPd fraction of respiration loss as DMSPd
706 : real (kind=dbl_kind), optional :: y_sk_DMS_in ! fraction conversion given high yield
707 : real (kind=dbl_kind), optional :: t_sk_conv_in ! Stefels conversion time (d)
708 : real (kind=dbl_kind), optional :: t_sk_ox_in ! DMS oxidation time (d)
709 : real (kind=dbl_kind), optional :: fsal_in ! salinity limitation factor (1)
710 :
711 : real (kind=dbl_kind), optional :: chlabs_in(:) ! chla absorption 1/m/(mg/m^3)
712 : real (kind=dbl_kind), optional :: alpha2max_low_in(:) ! light limitation (1/(W/m^2))
713 : real (kind=dbl_kind), optional :: beta2max_in(:) ! light inhibition (1/(W/m^2))
714 : real (kind=dbl_kind), optional :: mu_max_in(:) ! maximum growth rate (1/d)
715 : real (kind=dbl_kind), optional :: grow_Tdep_in(:) ! T dependence of growth (1/C)
716 : real (kind=dbl_kind), optional :: fr_graze_in(:) ! fraction of algae grazed
717 : real (kind=dbl_kind), optional :: mort_pre_in(:) ! mortality (1/day)
718 : real (kind=dbl_kind), optional :: mort_Tdep_in(:) ! T dependence of mortality (1/C)
719 : real (kind=dbl_kind), optional :: k_exude_in(:) ! algal carbon exudation rate (1/d)
720 : real (kind=dbl_kind), optional :: K_Nit_in(:) ! nitrate half saturation (mmol/m^3)
721 : real (kind=dbl_kind), optional :: K_Am_in(:) ! ammonium half saturation (mmol/m^3)
722 : real (kind=dbl_kind), optional :: K_Sil_in(:) ! silicon half saturation (mmol/m^3)
723 : real (kind=dbl_kind), optional :: K_Fe_in(:) ! iron half saturation or micromol/m^3
724 : real (kind=dbl_kind), optional :: f_don_in(:) ! fraction of spilled grazing to DON
725 : real (kind=dbl_kind), optional :: kn_bac_in(:) ! Bacterial degredation of DON (1/d)
726 : real (kind=dbl_kind), optional :: f_don_Am_in(:) ! fraction of remineralized DON to Am
727 : real (kind=dbl_kind), optional :: f_doc_in(:) ! fraction of mort_N that goes to each doc pool
728 : real (kind=dbl_kind), optional :: f_exude_in(:) ! fraction of exuded carbon to each DOC pool
729 : real (kind=dbl_kind), optional :: k_bac_in(:) ! Bacterial degredation of DOC (1/d)
730 :
731 : real (kind=dbl_kind), optional :: zbgc_frac_init_in(:) ! initializes mobile fraction
732 : real (kind=dbl_kind), optional :: bgc_tracer_type_in(:) ! described tracer in mobile or stationary phases
733 : real (kind=dbl_kind), optional :: zbgc_init_frac_in(:) ! fraction of ocean tracer concentration in new ice
734 : real (kind=dbl_kind), optional :: tau_ret_in(:) ! retention timescale (s), mobile to stationary phase
735 : real (kind=dbl_kind), optional :: tau_rel_in(:) ! release timescale (s), stationary to mobile phase
736 :
737 : !autodocument_end
738 :
739 : character(len=*),parameter :: subname='(icepack_init_zbgc)'
740 :
741 : !--------
742 :
743 14000 : if (present(R_C2N_in)) R_C2N(:) = R_C2N_in(:)
744 14000 : if (present(R_chl2N_in)) R_chl2N(:) = R_chl2N_in(:)
745 14000 : if (present(F_abs_chl_in)) F_abs_chl(:) = F_abs_chl_in(:)
746 8400 : if (present(R_C2N_DON_in)) R_C2N_DON(:) = R_C2N_DON_in(:)
747 14000 : if (present(R_Si2N_in)) R_Si2N(:) = R_Si2N_in(:)
748 14000 : if (present(R_S2N_in)) R_S2N(:) = R_S2N_in(:)
749 14000 : if (present(R_Fe2C_in)) R_Fe2C(:) = R_Fe2C_in(:)
750 14000 : if (present(R_Fe2N_in)) R_Fe2N(:) = R_Fe2N_in(:)
751 8400 : if (present(R_Fe2DON_in)) R_Fe2DON(:) = R_Fe2DON_in(:)
752 14000 : if (present(R_Fe2DOC_in)) R_Fe2DOC(:) = R_Fe2DOC_in(:)
753 :
754 5600 : if (present(fr_resp_in)) fr_resp = fr_resp_in
755 5600 : if (present(algal_vel_in)) algal_vel = algal_vel_in
756 5600 : if (present(R_dFe2dust_in)) R_dFe2dust = R_dFe2dust_in
757 5600 : if (present(dustFe_sol_in)) dustFe_sol = dustFe_sol_in
758 5600 : if (present(T_max_in)) T_max = T_max_in
759 5600 : if (present(op_dep_min_in)) op_dep_min = op_dep_min_in
760 5600 : if (present(fr_graze_s_in)) fr_graze_s = fr_graze_s_in
761 5600 : if (present(fr_graze_e_in)) fr_graze_e = fr_graze_e_in
762 5600 : if (present(fr_mort2min_in)) fr_mort2min = fr_mort2min_in
763 5600 : if (present(fr_dFe_in)) fr_dFe = fr_dFe_in
764 5600 : if (present(k_nitrif_in)) k_nitrif = k_nitrif_in
765 5600 : if (present(t_iron_conv_in)) t_iron_conv = t_iron_conv_in
766 5600 : if (present(max_loss_in)) max_loss = max_loss_in
767 5600 : if (present(max_dfe_doc1_in)) max_dfe_doc1 = max_dfe_doc1_in
768 5600 : if (present(fr_resp_s_in)) fr_resp_s = fr_resp_s_in
769 5600 : if (present(y_sk_DMS_in)) y_sk_DMS = y_sk_DMS_in
770 5600 : if (present(t_sk_conv_in)) t_sk_conv = t_sk_conv_in
771 5600 : if (present(t_sk_ox_in)) t_sk_ox = t_sk_ox_in
772 5600 : if (present(fsal_in)) fsal = fsal_in
773 :
774 14000 : if (present(chlabs_in)) chlabs(:) = chlabs_in(:)
775 14000 : if (present(alpha2max_low_in)) alpha2max_low(:) = alpha2max_low_in(:)
776 14000 : if (present(beta2max_in)) beta2max(:) = beta2max_in(:)
777 14000 : if (present(mu_max_in)) mu_max(:) = mu_max_in(:)
778 14000 : if (present(grow_Tdep_in)) grow_Tdep(:) = grow_Tdep_in(:)
779 14000 : if (present(fr_graze_in)) fr_graze(:) = fr_graze_in(:)
780 14000 : if (present(mort_pre_in)) mort_pre(:) = mort_pre_in(:)
781 14000 : if (present(mort_Tdep_in)) mort_Tdep(:) = mort_Tdep_in(:)
782 14000 : if (present(k_exude_in)) k_exude(:) = k_exude_in(:)
783 14000 : if (present(K_Nit_in)) K_Nit(:) = K_Nit_in(:)
784 14000 : if (present(K_Am_in)) K_Am(:) = K_Am_in(:)
785 14000 : if (present(K_Sil_in)) K_Sil(:) = K_Sil_in(:)
786 14000 : if (present(K_Fe_in)) K_Fe(:) = K_Fe_in(:)
787 8400 : if (present(f_don_in)) f_don(:) = f_don_in(:)
788 8400 : if (present(kn_bac_in)) kn_bac(:) = kn_bac_in(:)
789 8400 : if (present(f_don_Am_in)) f_don_Am(:) = f_don_Am_in(:)
790 14000 : if (present(f_doc_in)) f_doc(:) = f_doc_in(:)
791 14000 : if (present(f_exude_in)) f_exude(:) = f_exude_in(:)
792 14000 : if (present(k_bac_in)) k_bac(:) = k_bac_in(:)
793 :
794 86800 : if (present(zbgc_frac_init_in)) zbgc_frac_init(:) = zbgc_frac_init_in(:)
795 86800 : if (present(bgc_tracer_type_in)) bgc_tracer_type(:) = bgc_tracer_type_in(:)
796 86800 : if (present(zbgc_init_frac_in)) zbgc_init_frac(:) = zbgc_init_frac_in(:)
797 86800 : if (present(tau_ret_in)) tau_ret(:) = tau_ret_in(:)
798 86800 : if (present(tau_rel_in)) tau_rel(:) = tau_rel_in(:)
799 :
800 :
801 5600 : end subroutine icepack_init_zbgc
802 :
803 : !=======================================================================
804 : !autodocument_start icepack_biogeochemistry
805 : !
806 :
807 139062360 : subroutine icepack_biogeochemistry(dt, &
808 : ntrcr, nbtrcr, &
809 417187080 : upNO, upNH, iDi, iki, zfswin, &
810 139062360 : zsal_tot, darcy_V, grow_net, &
811 278124720 : PP_net, hbri,dhbr_bot, dhbr_top, Zoo,&
812 278124720 : fbio_snoice, fbio_atmice, ocean_bio, &
813 556249440 : first_ice, fswpenln, bphi, bTiz, ice_bio_net, &
814 139062360 : snow_bio_net, fswthrun, Rayleigh_criteria, &
815 139062360 : sice_rho, fzsal, fzsal_g, &
816 278124720 : bgrid, igrid, icgrid, cgrid, &
817 : nblyr, nilyr, nslyr, n_algae, n_zaero, ncat, &
818 : n_doc, n_dic, n_don, n_fed, n_fep, &
819 139062360 : meltbn, melttn, congeln, snoicen, &
820 139062360 : sst, sss, fsnow, meltsn, &
821 417187080 : hin_old, flux_bio, flux_bio_atm, &
822 417187080 : aicen_init, vicen_init, aicen, vicen, vsnon, &
823 278124720 : aice0, trcrn, vsnon_init, skl_bgc)
824 :
825 : real (kind=dbl_kind), intent(in) :: &
826 : dt ! time step
827 :
828 : integer (kind=int_kind), intent(in) :: &
829 : ncat, &
830 : nilyr, &
831 : nslyr, &
832 : nblyr, &
833 : ntrcr, &
834 : nbtrcr, &
835 : n_algae, n_zaero, &
836 : n_doc, n_dic, n_don, n_fed, n_fep
837 :
838 : real (kind=dbl_kind), dimension (:), intent(inout) :: &
839 : bgrid , & ! biology nondimensional vertical grid points
840 : igrid , & ! biology vertical interface points
841 : cgrid , & ! CICE vertical coordinate
842 : icgrid , & ! interface grid for CICE (shortwave variable)
843 : ocean_bio , & ! contains all the ocean bgc tracer concentrations
844 : fbio_snoice , & ! fluxes from snow to ice
845 : fbio_atmice , & ! fluxes from atm to ice
846 : dhbr_top , & ! brine top change
847 : dhbr_bot , & ! brine bottom change
848 : darcy_V , & ! darcy velocity positive up (m/s)
849 : hin_old , & ! old ice thickness
850 : sice_rho , & ! avg sea ice density (kg/m^3)
851 : ice_bio_net , & ! depth integrated tracer (mmol/m^2)
852 : snow_bio_net , & ! depth integrated snow tracer (mmol/m^2)
853 : flux_bio ! all bio fluxes to ocean
854 :
855 : logical (kind=log_kind), dimension (:), intent(inout) :: &
856 : first_ice ! distinguishes ice that disappears (e.g. melts)
857 : ! and reappears (e.g. transport) in a grid cell
858 : ! during a single time step from ice that was
859 : ! there the entire time step (true until ice forms)
860 :
861 : real (kind=dbl_kind), dimension (:,:), intent(inout) :: &
862 : Zoo , & ! N losses accumulated in timestep (ie. zooplankton/bacteria)
863 : ! mmol/m^3
864 : bphi , & ! porosity of layers
865 : bTiz , & ! layer temperatures interpolated on bio grid (C)
866 : zfswin , & ! Shortwave flux into layers interpolated on bio grid (W/m^2)
867 : iDi , & ! igrid Diffusivity (m^2/s)
868 : iki , & ! Ice permeability (m^2)
869 : trcrn ! tracers
870 :
871 : real (kind=dbl_kind), intent(inout) :: &
872 : grow_net , & ! Specific growth rate (/s) per grid cell
873 : PP_net , & ! Total production (mg C/m^2/s) per grid cell
874 : hbri , & ! brine height, area-averaged for comparison with hi (m)
875 : zsal_tot , & ! Total ice salinity in per grid cell (g/m^2)
876 : fzsal , & ! Total flux of salt to ocean at time step for conservation
877 : fzsal_g , & ! Total gravity drainage flux
878 : upNO , & ! nitrate uptake rate (mmol/m^2/d) times aice
879 : upNH ! ammonium uptake rate (mmol/m^2/d) times aice
880 :
881 : logical (kind=log_kind), intent(inout) :: &
882 : Rayleigh_criteria ! .true. means Ra_c was reached
883 :
884 : real (kind=dbl_kind), dimension (:,:), intent(in) :: &
885 : fswpenln ! visible SW entering ice layers (W m-2)
886 :
887 : real (kind=dbl_kind), dimension (:), intent(in) :: &
888 : fswthrun , & ! SW through ice to ocean (W/m^2)
889 : meltsn , & ! snow melt in category n (m)
890 : melttn , & ! top melt in category n (m)
891 : meltbn , & ! bottom melt in category n (m)
892 : congeln , & ! congelation ice formation in category n (m)
893 : snoicen , & ! snow-ice formation in category n (m)
894 : flux_bio_atm, & ! all bio fluxes to ice from atmosphere
895 : aicen_init , & ! initial ice concentration, for linear ITD
896 : vicen_init , & ! initial ice volume (m), for linear ITD
897 : vsnon_init , & ! initial snow volume (m), for aerosol
898 : aicen , & ! concentration of ice
899 : vicen , & ! volume per unit area of ice (m)
900 : vsnon ! volume per unit area of snow (m)
901 :
902 : real (kind=dbl_kind), intent(in) :: &
903 : aice0 , & ! open water area fraction
904 : sss , & ! sea surface salinity (ppt)
905 : sst , & ! sea surface temperature (C)
906 : fsnow ! snowfall rate (kg/m^2 s)
907 :
908 : logical (kind=log_kind), intent(in) :: &
909 : skl_bgc ! if true, solve skeletal biochemistry
910 :
911 : !autodocument_end
912 :
913 : ! local variables
914 :
915 : integer (kind=int_kind) :: &
916 : n, mm ! thickness category index
917 :
918 : real (kind=dbl_kind) :: &
919 553320 : hin , & ! new ice thickness
920 553320 : hsn , & ! snow thickness (m)
921 553320 : hbr_old , & ! old brine thickness before growh/melt
922 553320 : dhice , & ! change due to sublimation/condensation (m)
923 553320 : kavg , & ! average ice permeability (m^2)
924 553320 : bphi_o , & ! surface ice porosity
925 553320 : hbrin , & ! brine height
926 553320 : dh_direct ! surface flooding or runoff
927 :
928 : real (kind=dbl_kind), dimension (nblyr+2) :: &
929 : ! Defined on Bio Grid points
930 280880880 : bSin , & ! salinity on the bio grid (ppt)
931 280880880 : brine_sal , & ! brine salinity (ppt)
932 281987520 : brine_rho ! brine_density (kg/m^3)
933 :
934 : real (kind=dbl_kind), dimension (nblyr+1) :: &
935 : ! Defined on Bio Grid interfaces
936 280327560 : iphin , & ! porosity
937 280327560 : ibrine_sal , & ! brine salinity (ppt)
938 280327560 : ibrine_rho , & ! brine_density (kg/m^3)
939 142371840 : iTin ! Temperature on the interface grid (oC)
940 :
941 : real (kind=dbl_kind) :: &
942 553320 : sloss ! brine flux contribution from surface runoff (g/m^2)
943 :
944 : ! for bgc sk
945 : real (kind=dbl_kind) :: &
946 553320 : dh_bot_chl , & ! Chlorophyll may or may not flush
947 553320 : dh_top_chl , & ! Chlorophyll may or may not flush
948 553320 : darcy_V_chl
949 :
950 : character(len=*),parameter :: subname='(icepack_biogeochemistry)'
951 :
952 :
953 834374160 : do n = 1, ncat
954 :
955 : !-----------------------------------------------------------------
956 : ! initialize
957 : !-----------------------------------------------------------------
958 695311800 : hin_old(n) = c0
959 695311800 : if (aicen_init(n) > puny) then
960 160682 : hin_old(n) = vicen_init(n) &
961 58697512 : / aicen_init(n)
962 : else
963 636614288 : first_ice(n) = .true.
964 636614288 : if (tr_brine) trcrn(nt_fbri,n) = c1
965 636614288 : if (z_tracers) then
966 6391854060 : do mm = 1,nbtrcr
967 6391854060 : trcrn(nt_zbgc_frac-1+mm,n) = zbgc_frac_init(mm)
968 : enddo
969 : endif
970 636614288 : if (n == 1) Rayleigh_criteria = .false.
971 636614288 : if (solve_zsal) trcrn(nt_bgc_S:nt_bgc_S+nblyr-1,n) = c0
972 : endif
973 :
974 834374160 : if (aicen(n) > puny) then
975 :
976 58670776 : dh_top_chl = c0
977 58670776 : dh_bot_chl = c0
978 58670776 : darcy_V_chl= c0
979 411177178 : bSin(:) = c0
980 58670776 : hsn = c0
981 58670776 : hin = c0
982 58670776 : hbrin = c0
983 58670776 : kavg = c0
984 58670776 : bphi_o = c0
985 58670776 : sloss = c0
986 :
987 : !-----------------------------------------------------------------
988 : ! brine dynamics
989 : !-----------------------------------------------------------------
990 :
991 58670776 : dhbr_top(n) = c0
992 58670776 : dhbr_bot(n) = c0
993 :
994 58670776 : if (tr_brine) then
995 58670776 : if (trcrn(nt_fbri,n) .le. c0) trcrn(nt_fbri,n) = c1
996 :
997 58670776 : dhice = c0
998 160582 : call preflushing_changes (aicen (n), &
999 160582 : vicen (n), vsnon (n), &
1000 160582 : meltbn (n), melttn (n), &
1001 160582 : congeln (n), snoicen(n), &
1002 160582 : hin_old (n), dhice, &
1003 160582 : trcrn(nt_fbri,n), &
1004 160582 : dhbr_top(n), dhbr_bot(n), &
1005 : hbr_old, hin, &
1006 58670776 : hsn, first_ice(n) )
1007 58670776 : if (icepack_warnings_aborted(subname)) return
1008 :
1009 58670776 : if (solve_zsal) then
1010 :
1011 : call compute_microS (n, nilyr, nblyr, &
1012 0 : bgrid, cgrid, igrid, &
1013 0 : trcrn(1:ntrcr,n), hin_old(n), hbr_old, &
1014 0 : sss, sst, bTiz(:,n), &
1015 0 : iTin, bphi(:,n), kavg, &
1016 : bphi_o, Rayleigh_criteria, &
1017 0 : first_ice(n), bSin, brine_sal, &
1018 : brine_rho, iphin, ibrine_rho, &
1019 0 : ibrine_sal, sice_rho(n), sloss)
1020 0 : if (icepack_warnings_aborted(subname)) return
1021 : else
1022 :
1023 : ! Requires the average ice permeability = kavg(:)
1024 : ! and the surface ice porosity = zphi_o(:)
1025 : ! computed in "compute_microS" or from "thermosaline_vertical"
1026 :
1027 352506402 : iDi(:,n) = c0
1028 :
1029 : call compute_microS_mushy (nilyr, nblyr, &
1030 0 : bgrid, cgrid, igrid, &
1031 160582 : trcrn(:,n), hin_old(n), hbr_old, &
1032 0 : sss, sst, bTiz(:,n), &
1033 0 : iTin(:), bphi(:,n), kavg, &
1034 : bphi_o, bSin(:), &
1035 : brine_sal(:), brine_rho(:), iphin(:), &
1036 160582 : ibrine_rho(:), ibrine_sal(:), sice_rho(n), &
1037 58991940 : iDi(:,n) )
1038 58670776 : if (icepack_warnings_aborted(subname)) return
1039 :
1040 : endif ! solve_zsal
1041 :
1042 160582 : call update_hbrine (melttn(n), &
1043 160582 : meltsn (n), dt, &
1044 : hin, hsn, &
1045 160582 : hin_old (n), hbrin, &
1046 :
1047 : hbr_old, &
1048 160582 : trcrn(nt_fbri,n), &
1049 160582 : dhbr_top(n), dhbr_bot(n), &
1050 : dh_top_chl, dh_bot_chl, &
1051 : kavg, bphi_o, &
1052 160582 : darcy_V (n), darcy_V_chl, &
1053 160582 : bphi(2,n), aice0, &
1054 58670776 : dh_direct)
1055 58670776 : if (icepack_warnings_aborted(subname)) return
1056 :
1057 58670776 : hbri = hbri + hbrin * aicen(n)
1058 :
1059 58670776 : if (solve_zsal) then
1060 :
1061 : call zsalinity (n, dt, &
1062 0 : nilyr, bgrid, &
1063 0 : cgrid, igrid, &
1064 0 : trcrn(nt_bgc_S:nt_bgc_S+nblyr-1,n), &
1065 0 : trcrn(nt_qice:nt_qice+nilyr-1,n), &
1066 0 : trcrn(nt_sice:nt_sice+nilyr-1,n), &
1067 0 : ntrcr, trcrn(nt_fbri,n), &
1068 0 : bSin, bTiz(:,n), &
1069 0 : bphi(:,n), iphin, &
1070 0 : iki(:,n), hbr_old, &
1071 : hbrin, hin, &
1072 0 : hin_old(n), iDi(:,n), &
1073 0 : darcy_V(n), brine_sal, &
1074 : brine_rho, ibrine_sal, &
1075 : ibrine_rho, dh_direct, &
1076 : Rayleigh_criteria, &
1077 0 : first_ice(n), sss, &
1078 0 : sst, dhbr_top(n), &
1079 0 : dhbr_bot(n), &
1080 : fzsal, fzsal_g, &
1081 : bphi_o, nblyr, &
1082 0 : vicen(n), aicen_init(n), &
1083 0 : zsal_tot)
1084 0 : if (icepack_warnings_aborted(subname)) return
1085 :
1086 : endif ! solve_zsal
1087 :
1088 : endif ! tr_brine
1089 :
1090 : !-----------------------------------------------------------------
1091 : ! biogeochemistry
1092 : !-----------------------------------------------------------------
1093 :
1094 58670776 : if (z_tracers) then
1095 :
1096 : call zbio (dt, nblyr, &
1097 : nslyr, nilyr, &
1098 80291 : melttn(n), &
1099 80291 : meltsn(n), meltbn (n), &
1100 80291 : congeln(n), snoicen(n), &
1101 : nbtrcr, fsnow, &
1102 0 : ntrcr, trcrn(1:ntrcr,n), &
1103 80291 : bio_index(1:nbtrcr), aicen_init(n), &
1104 80291 : vicen_init(n), vsnon_init(n), &
1105 80291 : vicen(n), vsnon(n), &
1106 80291 : aicen(n), flux_bio_atm(1:nbtrcr), &
1107 : n, n_algae, &
1108 : n_doc, n_dic, &
1109 : n_don, &
1110 : n_fed, n_fep, &
1111 80291 : n_zaero, first_ice(n), &
1112 80291 : hin_old(n), ocean_bio(1:nbtrcr), &
1113 0 : bphi(:,n), iphin, &
1114 0 : iDi(:,n), &
1115 0 : fswpenln(:,n), &
1116 80291 : dhbr_top(n), dhbr_bot(n), &
1117 0 : zfswin(:,n), &
1118 : hbrin, hbr_old, &
1119 : ! darcy_V(n), darcy_V_chl, &
1120 80291 : darcy_V(n), &
1121 0 : bgrid, &
1122 0 : igrid, icgrid, &
1123 : bphi_o, &
1124 : iTin, &
1125 0 : Zoo(:,n), &
1126 0 : flux_bio(1:nbtrcr), dh_direct, &
1127 : upNO, upNH, &
1128 0 : fbio_snoice, fbio_atmice, &
1129 0 : PP_net, ice_bio_net (1:nbtrcr), &
1130 29736843 : snow_bio_net(1:nbtrcr),grow_net )
1131 29415679 : if (icepack_warnings_aborted(subname)) return
1132 :
1133 29255097 : elseif (skl_bgc) then
1134 :
1135 : call sklbio (dt, ntrcr, &
1136 : nbtrcr, n_algae, &
1137 : n_doc, &
1138 : n_dic, n_don, &
1139 : n_fed, n_fep, &
1140 0 : flux_bio (1:nbtrcr), ocean_bio(1:nbtrcr), &
1141 80291 : aicen (n), &
1142 80291 : meltbn (n), congeln (n), &
1143 80291 : fswthrun (n), first_ice(n), &
1144 0 : trcrn (1:ntrcr,n), &
1145 : PP_net, upNO, &
1146 29335388 : upNH, grow_net )
1147 29255097 : if (icepack_warnings_aborted(subname)) return
1148 :
1149 : endif ! skl_bgc
1150 :
1151 58670776 : first_ice(n) = .false.
1152 :
1153 : endif ! aicen > puny
1154 : enddo ! ncat
1155 :
1156 : end subroutine icepack_biogeochemistry
1157 :
1158 : !=======================================================================
1159 : !autodocument_start icepack_load_ocean_bio_array
1160 : ! basic initialization for ocean_bio_all
1161 :
1162 139874765 : subroutine icepack_load_ocean_bio_array(max_nbtrcr, &
1163 : max_algae, max_don, max_doc, max_dic, max_aero, max_fe, &
1164 139874765 : nit, amm, sil, dmsp, dms, algalN, &
1165 279749530 : doc, don, dic, fed, fep, zaeros, ocean_bio_all, hum)
1166 :
1167 : integer (kind=int_kind), intent(in) :: &
1168 : max_algae , & ! maximum number of algal types
1169 : max_dic , & ! maximum number of dissolved inorganic carbon types
1170 : max_doc , & ! maximum number of dissolved organic carbon types
1171 : max_don , & ! maximum number of dissolved organic nitrogen types
1172 : max_fe , & ! maximum number of iron types
1173 : max_aero , & ! maximum number of aerosols
1174 : max_nbtrcr ! maximum number of bio tracers
1175 :
1176 : real (kind=dbl_kind), intent(in) :: &
1177 : nit , & ! ocean nitrate (mmol/m^3)
1178 : amm , & ! ammonia/um (mmol/m^3)
1179 : sil , & ! silicate (mmol/m^3)
1180 : dmsp , & ! dmsp (mmol/m^3)
1181 : dms , & ! dms (mmol/m^3)
1182 : hum ! humic material (mmol/m^3)
1183 :
1184 : real (kind=dbl_kind), dimension (max_algae), intent(in) :: &
1185 : algalN ! ocean algal nitrogen (mmol/m^3) (diatoms, phaeo, pico)
1186 :
1187 : real (kind=dbl_kind), dimension (max_doc), intent(in) :: &
1188 : doc ! ocean doc (mmol/m^3) (proteins, EPS, lipid)
1189 :
1190 : real (kind=dbl_kind), dimension (max_don), intent(in) :: &
1191 : don ! ocean don (mmol/m^3)
1192 :
1193 : real (kind=dbl_kind), dimension (max_dic), intent(in) :: &
1194 : dic ! ocean dic (mmol/m^3)
1195 :
1196 : real (kind=dbl_kind), dimension (max_fe), intent(in) :: &
1197 : fed, fep ! ocean disolved and particulate fe (nM)
1198 :
1199 : real (kind=dbl_kind), dimension (max_aero), intent(in) :: &
1200 : zaeros ! ocean aerosols (mmol/m^3)
1201 :
1202 : real (kind=dbl_kind), dimension (max_nbtrcr), intent(inout) :: &
1203 : ocean_bio_all ! fixed order, all values even for tracers false
1204 :
1205 : !autodocument_end
1206 :
1207 : ! local variables
1208 :
1209 : integer (kind=int_kind) :: &
1210 : k, ks ! tracer indices
1211 :
1212 : character(len=*),parameter :: subname='(icepack_load_ocean_bio_array)'
1213 :
1214 4196242950 : ocean_bio_all(:) = c0
1215 :
1216 559499060 : do k = 1, max_algae
1217 419624295 : ocean_bio_all(k) = algalN(k) ! N
1218 419624295 : ks = max_algae + max_doc + max_dic + 1
1219 559499060 : ocean_bio_all(ks + k) = R_chl2N(k)*algalN(k)!chl
1220 : enddo
1221 :
1222 139874765 : ks = max_algae + 1
1223 559499060 : do k = 1, max_doc
1224 559499060 : ocean_bio_all(ks + k) = doc(k) ! doc
1225 : enddo
1226 139874765 : ks = ks + max_doc
1227 279749530 : do k = 1, max_dic
1228 279749530 : ocean_bio_all(ks + k) = dic(k) ! dic
1229 : enddo
1230 :
1231 139874765 : ks = 2*max_algae + max_doc + max_dic + 7
1232 279749530 : do k = 1, max_don
1233 279749530 : ocean_bio_all(ks + k) = don(k) ! don
1234 : enddo
1235 :
1236 139874765 : ks = max_algae + 1
1237 139874765 : ocean_bio_all(ks) = nit ! nit
1238 :
1239 139874765 : ks = 2*max_algae + max_doc + 2 + max_dic
1240 139874765 : ocean_bio_all(ks) = amm ! Am
1241 139874765 : ks = ks + 1
1242 139874765 : ocean_bio_all(ks) = sil ! Sil
1243 139874765 : ks = ks + 1
1244 576375 : ocean_bio_all(ks) = R_S2N(1)*algalN(1) & ! DMSPp
1245 576375 : + R_S2N(2)*algalN(2) &
1246 139874765 : + R_S2N(3)*algalN(3)
1247 139874765 : ks = ks + 1
1248 139874765 : ocean_bio_all(ks) = dmsp ! DMSPd
1249 139874765 : ks = ks + 1
1250 139874765 : ocean_bio_all(ks) = dms ! DMS
1251 139874765 : ks = ks + 1
1252 139874765 : ocean_bio_all(ks) = nit ! PON
1253 139874765 : ks = 2*max_algae + max_doc + 7 + max_dic + max_don
1254 419624295 : do k = 1, max_fe
1255 419624295 : ocean_bio_all(ks + k) = fed(k) ! fed
1256 : enddo
1257 139874765 : ks = ks + max_fe
1258 419624295 : do k = 1, max_fe
1259 419624295 : ocean_bio_all(ks + k) = fep(k) ! fep
1260 : enddo
1261 139874765 : ks = ks + max_fe
1262 979123355 : do k = 1, max_aero
1263 979123355 : ocean_bio_all(ks+k) = zaeros(k) ! zaero
1264 : enddo
1265 139874765 : ks = ks + max_aero + 1
1266 139874765 : ocean_bio_all(ks) = hum ! humics
1267 :
1268 139874765 : end subroutine icepack_load_ocean_bio_array
1269 :
1270 : !=======================================================================
1271 : !autodocument_start icepack_init_ocean_bio
1272 : ! Initialize ocean concentration
1273 :
1274 429185 : subroutine icepack_init_ocean_bio (amm, dmsp, dms, algalN, doc, dic, don, &
1275 429185 : fed, fep, hum, nit, sil, zaeros, max_dic, max_don, max_fe, max_aero,&
1276 429185 : CToN, CToN_DON)
1277 :
1278 : integer (kind=int_kind), intent(in) :: &
1279 : max_dic, &
1280 : max_don, &
1281 : max_fe, &
1282 : max_aero
1283 :
1284 : real (kind=dbl_kind), intent(out):: &
1285 : amm , & ! ammonium
1286 : dmsp , & ! DMSPp
1287 : dms , & ! DMS
1288 : hum , & ! humic material
1289 : nit , & ! nitrate
1290 : sil ! silicate
1291 :
1292 : real (kind=dbl_kind), dimension(:), intent(out):: &
1293 : algalN , & ! algae
1294 : doc , & ! DOC
1295 : dic , & ! DIC
1296 : don , & ! DON
1297 : fed , & ! Dissolved Iron
1298 : fep , & ! Particulate Iron
1299 : zaeros ! BC and dust
1300 :
1301 : real (kind=dbl_kind), dimension(:), intent(inout), optional :: &
1302 : CToN , & ! carbon to nitrogen ratio for algae
1303 : CToN_DON ! nitrogen to carbon ratio for proteins
1304 :
1305 : !autodocument_end
1306 :
1307 : ! local variables
1308 :
1309 : integer (kind=int_kind) :: &
1310 : k
1311 :
1312 : character(len=*),parameter :: subname='(icepack_init_ocean_bio)'
1313 :
1314 429185 : if (present(CToN)) then
1315 0 : CToN(1) = R_C2N(1)
1316 0 : CToN(2) = R_C2N(2)
1317 0 : CToN(3) = R_C2N(3)
1318 : endif
1319 :
1320 429185 : if (present(CToN_DON)) then
1321 0 : CToN_DON(1) = R_C2N_DON(1)
1322 : endif
1323 :
1324 429185 : amm = c1 ! ISPOL < 1 mmol/m^3
1325 429185 : dmsp = p1
1326 429185 : dms = p1
1327 429185 : algalN(1) = c1 !0.0026_dbl_kind ! ISPOL, Lannuzel 2013(pennate)
1328 429185 : algalN(2) = 0.0057_dbl_kind ! ISPOL, Lannuzel 2013(small plankton)
1329 429185 : algalN(3) = 0.0027_dbl_kind ! ISPOL, Lannuzel 2013(Phaeocystis)
1330 : ! 0.024_dbl_kind ! 5% of 1 mgchl/m^3
1331 429185 : doc(1) = 16.2_dbl_kind ! 18% saccharides
1332 429185 : doc(2) = 9.0_dbl_kind ! lipids
1333 429185 : doc(3) = c1 !
1334 858370 : do k = 1, max_dic
1335 858370 : dic(k) = c1
1336 : enddo
1337 858370 : do k = 1, max_don
1338 858370 : don(k) = 12.9_dbl_kind
1339 : ! 64.3_dbl_kind ! 72% Total DOC~90 mmolC/m^3 ISPOL with N:C of 0.2
1340 : enddo
1341 : !ki = 1
1342 : !if (trim(fe_data_type) == 'clim') ki = 2
1343 1287555 : do k = 1, max_fe ! ki, max_fe
1344 858370 : fed(k) = 0.4_dbl_kind ! c1 (nM) Lannuzel2007 DFe,
1345 : ! range 0.14-2.6 (nM) van der Merwe 2011
1346 : ! Tagliabue 2012 (0.4 nM)
1347 1287555 : fep(k) = c2 ! (nM) van der Merwe 2011
1348 : ! (0.6 to 2.9 nM ocean)
1349 : enddo
1350 429185 : hum = c1 ! mmol C/m^3
1351 429185 : nit = 12.0_dbl_kind
1352 429185 : sil = 25.0_dbl_kind
1353 3004295 : do k = 1, max_aero
1354 3004295 : zaeros(k) = c0
1355 : enddo
1356 :
1357 :
1358 429185 : end subroutine icepack_init_ocean_bio
1359 :
1360 : !=======================================================================
1361 :
1362 : end module icepack_zbgc
1363 :
1364 : !=======================================================================
|
