diff --git a/unchanged/marbl_init_mod.F90 b/unchanged/marbl_init_mod.F90 index 862911e..fc504b7 100644 --- a/unchanged/marbl_init_mod.F90 +++ b/unchanged/marbl_init_mod.F90 @@ -895,6 +895,14 @@ subroutine marbl_init_interior_tendency_forcing_fields(& call interior_tendency_forcings(id)%set_rank(num_elements, 1, marbl_status_log, dim1 = num_levels) end if + ! DOC fluid Flux + if (id .eq. ind%docfluidy_id) then + found = .true. + interior_tendency_forcings(id)%metadata%varname = 'DOC fluid Flux' + interior_tendency_forcings(id)%metadata%field_units = '1/yr' + call interior_tendency_forcings(id)%set_rank(num_elements, 1, marbl_status_log, dim1 = num_levels) + end if + ! O2 Consumption Scale Factor if (id .eq. ind%o2_consumption_scalef_id) then found = .true.