diff --git a/components/mpas-seaice/bld/namelist_files/namelist_definition_mpassi.xml b/components/mpas-seaice/bld/namelist_files/namelist_definition_mpassi.xml
index efbaef07cb9d..ba6ae4be4bb1 100644
--- a/components/mpas-seaice/bld/namelist_files/namelist_definition_mpassi.xml
+++ b/components/mpas-seaice/bld/namelist_files/namelist_definition_mpassi.xml
@@ -1077,7 +1077,7 @@ Default: Defined in namelist_defaults.xml
 
 <entry id="config_use_humics" type="logical"
 	category="biogeochemistry" group="biogeochemistry">
-Use the humic matter tracer
+Use the humic (refractory dissolved organic matter) tracer
 
 Valid values: true or false
 Default: Defined in namelist_defaults.xml
@@ -1949,7 +1949,7 @@ Default: Defined in namelist_defaults.xml
 
 <entry id="config_mobility_type_humics" type="real"
 	category="biogeochemistry" group="biogeochemistry">
-Transport type of humics
+Transport type of humics (refractory dissolved organic matter)
 
 Valid values: -1 = entirely in the mobile phase; 0 = retention dominated; 1 = release dominated; 0.5 = equal but rapid exchange; 2 = equal but slow exchange
 Default: Defined in namelist_defaults.xml
diff --git a/components/mpas-seaice/cime_config/buildnml b/components/mpas-seaice/cime_config/buildnml
index 471fbf9233f2..29b017f6d8ba 100755
--- a/components/mpas-seaice/cime_config/buildnml
+++ b/components/mpas-seaice/cime_config/buildnml
@@ -934,22 +934,41 @@ def buildnml(case, caseroot, compname):
                 lines.append('    <var name="bergLatentHeatFlux"/>')
 
             if ice_bgc == 'ice_bgc':
-                lines.append('    <var name="totalVerticalBiologyIce"/>')
                 lines.append('    <var name="totalChlorophyll"/>')
                 lines.append('    <var name="primaryProduction"/>')
                 lines.append('    <var name="bioShortwaveFlux"/>')
                 lines.append('    <var name="netBrineHeight"/>')
                 lines.append('    <var name="newlyFormedIce"/>')
                 lines.append('    <var name="snowThicknessChange"/>')
-                lines.append('    <var name="verticalAlgaeIceCell"/>')
+                lines.append('    <var name="totalVerticalDiatomIce"/>')
+                lines.append('    <var name="totalVerticalSmallPlanktonIce"/>')
+                lines.append('    <var name="totalVerticalAlgaeCarbonIce"/>')
+                lines.append('    <var name="totalVerticalNitrateIce"/>')
+                lines.append('    <var name="totalVerticalAmmoniumIce"/>')
+                lines.append('    <var name="totalVerticalSilicateIce"/>')
+                lines.append('    <var name="totalVerticalProteinsIce"/>')
+                lines.append('    <var name="totalVerticalDOCLabileIce"/>')
+                lines.append('    <var name="totalVerticalHumicsIce"/>')
+                lines.append('    <var name="totalVerticalDissolvedIronIce"/>')
+                lines.append('    <var name="totalVerticalParticulateIronIce"/>')
+                lines.append('    <var name="totalVerticalBCIce"/>')
+                lines.append('    <var name="totalVerticalBCSnow"/>')
+                lines.append('    <var name="totalVerticalDustIce"/>')
+                lines.append('    <var name="totalVerticalDustSnow"/>')
+                lines.append('    <var name="totalVerticalDissolvedIronSnow"/>')
+                lines.append('    <var name="verticalAlgaeTotalCarbonIceCell"/>')
+                lines.append('    <var name="verticalDOCLabileIceCell"/>')
+                lines.append('    <var name="verticalDONIceCell"/>')
+                lines.append('    <var name="verticalDICIceCell"/>')
                 lines.append('    <var name="verticalNitrateIceCell"/>')
                 lines.append('    <var name="verticalSilicateIceCell"/>')
                 lines.append('    <var name="verticalAmmoniumIceCell"/>')
-                lines.append('    <var name="verticalDMSIceCell"/>')
-                lines.append('    <var name="verticalDMSPdIceCell"/>')
-                lines.append('    <var name="verticalNonreactiveIceCell"/>')
                 lines.append('    <var name="verticalHumicsIceCell"/>')
                 lines.append('    <var name="verticalDissolvedIronIceCell"/>')
+                lines.append('    <var name="verticalBCTotalIceCell"/>')
+                lines.append('    <var name="verticalDustTotalIceCell"/>')
+                lines.append('    <var name="verticalBCTotalSnowCell"/>')
+                lines.append('    <var name="verticalDustTotalSnowCell"/>')
                 lines.append('    <var name="iceEnthalpyCell"/>')
                 lines.append('    <var name="iceSalinityCell"/>')
                 lines.append('    <var name="oceanBioFluxes"/>')
@@ -1009,15 +1028,35 @@ def buildnml(case, caseroot, compname):
                 lines.append('    <var name="snowAlbedoCell"/>')
                 lines.append('    <var name="bareIceAlbedoCell"/>')
                 lines.append('    <var name="pondAlbedoCell"/>')
-                lines.append('    <var name="verticalAlgaeIceCell"/>')
+                lines.append('    <var name="totalVerticalDiatomIce"/>')
+                lines.append('    <var name="totalVerticalSmallPlanktonIce"/>')
+                lines.append('    <var name="totalVerticalAlgaeCarbonIce"/>')
+                lines.append('    <var name="totalVerticalNitrateIce"/>')
+                lines.append('    <var name="totalVerticalAmmoniumIce"/>')
+                lines.append('    <var name="totalVerticalSilicateIce"/>')
+                lines.append('    <var name="totalVerticalProteinsIce"/>')
+                lines.append('    <var name="totalVerticalDOCLabileIce"/>')
+                lines.append('    <var name="totalVerticalHumicsIce"/>')
+                lines.append('    <var name="totalVerticalDissolvedIronIce"/>')
+                lines.append('    <var name="totalVerticalParticulateIronIce"/>')
+                lines.append('    <var name="totalVerticalBCIce"/>')
+                lines.append('    <var name="totalVerticalBCSnow"/>')
+                lines.append('    <var name="totalVerticalDustIce"/>')
+                lines.append('    <var name="totalVerticalDustSnow"/>')
+                lines.append('    <var name="totalVerticalDissolvedIronSnow"/>')
+                lines.append('    <var name="verticalAlgaeTotalCarbonIceCell"/>')
+                lines.append('    <var name="verticalDOCLabileIceCell"/>')
+                lines.append('    <var name="verticalDONIceCell"/>')
+                lines.append('    <var name="verticalDICIceCell"/>')
                 lines.append('    <var name="verticalNitrateIceCell"/>')
                 lines.append('    <var name="verticalSilicateIceCell"/>')
                 lines.append('    <var name="verticalAmmoniumIceCell"/>')
-                lines.append('    <var name="verticalDMSIceCell"/>')
-                lines.append('    <var name="verticalDMSPdIceCell"/>')
-                lines.append('    <var name="verticalNonreactiveIceCell"/>')
                 lines.append('    <var name="verticalHumicsIceCell"/>')
                 lines.append('    <var name="verticalDissolvedIronIceCell"/>')
+                lines.append('    <var name="verticalBCTotalIceCell"/>')
+                lines.append('    <var name="verticalDustTotalIceCell"/>')
+                lines.append('    <var name="verticalBCTotalSnowCell"/>')
+                lines.append('    <var name="verticalDustTotalSnowCell"/>')
                 lines.append('    <var name="iceEnthalpyCell"/>')
                 lines.append('    <var name="iceSalinityCell"/>')
                 lines.append('    <var name="oceanBioFluxes"/>')
diff --git a/components/mpas-seaice/src/Registry.xml b/components/mpas-seaice/src/Registry.xml
index 7c842ca2f3c4..78830161c62d 100644
--- a/components/mpas-seaice/src/Registry.xml
+++ b/components/mpas-seaice/src/Registry.xml
@@ -866,7 +866,7 @@
 			icepack_name="tr_bgc_PON"
 		/>
 		<nml_option name="config_use_humics" type="logical" default_value="false" units="unitless"
-			description="Use the humic matter tracer"
+			description="Use the humic (refractory dissolved organic matter) tracer"
 			possible_values="true or false"
 			icepack_name="tr_bgc_hum"
 		/>
@@ -948,77 +948,77 @@
 			possible_values="positive real number"
 			icepack_name="frazil_scav"
 		/>
-		<nml_option name="config_ratio_Si_to_N_diatoms" type="real" default_value="1.80" units="mol silicate/mol nitrogen"
+		<nml_option name="config_ratio_Si_to_N_diatoms" type="real" default_value="1.80" units="mol Si (mol N)-1"
 			description="Sets the ratio of silicate to nitrogen in diatoms"
 			possible_values="positive real of order 1.8"
 			icepack_name="ratio_Si2N_diatoms"
 		/>
-		<nml_option name="config_ratio_Si_to_N_small_plankton" type="real" default_value="0.00" units="mol silicate/mol nitrogen"
+		<nml_option name="config_ratio_Si_to_N_small_plankton" type="real" default_value="0.00" units="mol Si (mol N)-1"
 			description="Sets the ratio of silicate to nitrogen in small plankton"
 			possible_values="0"
 			icepack_name="ratio_Si2N_sp"
 		/>
-		<nml_option name="config_ratio_Si_to_N_phaeocystis" type="real" default_value="0.00" units="mol silicate/mol nitrogen"
+		<nml_option name="config_ratio_Si_to_N_phaeocystis" type="real" default_value="0.00" units="mol Si (mol N)-1"
 			description="Sets the ratio of silicate to nitrogen in phaeocystis"
 			possible_values="0"
 			icepack_name="ratio_Si2N_phaeo"
 		/>
-		<nml_option name="config_ratio_S_to_N_diatoms" type="real" default_value="0.03" units="mol Sulfur/mol nitrogen"
+		<nml_option name="config_ratio_S_to_N_diatoms" type="real" default_value="0.03" units="mol S (mol N)-1"
 			description="Sets the ratio of sulfur to nitrogen in diatoms"
 			possible_values="order of 0.03"
 			icepack_name="ratio_S2N_diatoms"
 		/>
-		<nml_option name="config_ratio_S_to_N_small_plankton" type="real" default_value="0.03" units="mol sulfur/mol nitrogen"
+		<nml_option name="config_ratio_S_to_N_small_plankton" type="real" default_value="0.03" units="mol S (mol N)-1"
 			description="Sets the ratio of sulfur to nitrogen in small plankeon"
 			possible_values="order of 0.03"
 			icepack_name="ratio_S2N_sp"
 		/>
-		<nml_option name="config_ratio_S_to_N_phaeocystis" type="real" default_value="0.03" units="mol sulfur/mol nitrogen"
+		<nml_option name="config_ratio_S_to_N_phaeocystis" type="real" default_value="0.03" units="mol S (mol N)-1"
 			description="Sets the ratio of sulfur to nitrogen in phaeocystis"
 			possible_values="order of 0.03"
 			icepack_name="ratio_S2N_phaeo"
 		/>
-		<nml_option name="config_ratio_Fe_to_C_diatoms" type="real" default_value="0.0033" units="mmol iron/mol Carbon"
+		<nml_option name="config_ratio_Fe_to_C_diatoms" type="real" default_value="0.0033" units="mmol Fe (mol C)-1"
 			description="Sets the ratio of iron to carbon in diatoms"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2C_diatoms"
 		/>
-		<nml_option name="config_ratio_Fe_to_C_small_plankton" type="real" default_value="0.0033" units="mmol iron/mol carbon"
+		<nml_option name="config_ratio_Fe_to_C_small_plankton" type="real" default_value="0.0033" units="mmol Fe (mol C)-1"
 			description="Sets the ratio of iron to carbon in small plankton"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2C_sp"
 		/>
-		<nml_option name="config_ratio_Fe_to_C_phaeocystis" type="real" default_value="0.1" units="mmol iron/mol carbon"
+		<nml_option name="config_ratio_Fe_to_C_phaeocystis" type="real" default_value="0.1" units="mmol Fe (mol C)-1"
 			description="Sets the ratio of iron to carbon in phaeocystis"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2C_phaeo"
 		/>
-		<nml_option name="config_ratio_Fe_to_N_diatoms" type="real" default_value="0.023" units="mmol iron/mol nitrogen"
+		<nml_option name="config_ratio_Fe_to_N_diatoms" type="real" default_value="0.023" units="mmol Fe (mol N)-1"
 			description="Sets the ratio of iron to nitrogen in diatoms"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2N_diatoms"
 		/>
-		<nml_option name="config_ratio_Fe_to_N_small_plankton" type="real" default_value="0.023" units="mmol iron/mol nitrogen"
+		<nml_option name="config_ratio_Fe_to_N_small_plankton" type="real" default_value="0.023" units="mmol Fe (mol N)-1"
 			description="Sets the ratio of iron to nitrogen in small plankton"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2N_sp"
 		/>
-		<nml_option name="config_ratio_Fe_to_N_phaeocystis" type="real" default_value="0.7" units="mmol iron/mol nitrogen"
+		<nml_option name="config_ratio_Fe_to_N_phaeocystis" type="real" default_value="0.7" units="mmol Fe (mol N)-1"
 			description="Sets the ratio of iron to nitrogen in phaeocystis"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2N_phaeo"
 		/>
-		<nml_option name="config_ratio_Fe_to_DON" type="real" default_value="0.023" units="mmol iron/mol nitrogen"
+		<nml_option name="config_ratio_Fe_to_DON" type="real" default_value="0.023" units="mmol Fe (mol N)-1"
 			description="Sets the ratio of iron to dissolved organic nitrogen"
 			possible_values="true or false"
 			icepack_name="ratio_Fe2DON"
 		/>
-		<nml_option name="config_ratio_Fe_to_DOC_saccharids" type="real" default_value="0.1" units="mmol iron/mol carbon"
+		<nml_option name="config_ratio_Fe_to_DOC_saccharids" type="real" default_value="0.1" units="mmol Fe (mol C)-1"
 			description="Sets the ratio of iron to poly/saccharid carbon"
 			possible_values="positive real number"
 			icepack_name="ratio_Fe2DOC_s"
 		/>
-		<nml_option name="config_ratio_Fe_to_DOC_lipids" type="real" default_value="0.033" units="mmol iron/mol carbon"
+		<nml_option name="config_ratio_Fe_to_DOC_lipids" type="real" default_value="0.033" units="mmol Fe (mol C)-1"
 			description="Sets the ratio of iron to lipid carbon"
 			possible_values="positive real number"
 			icepack_name="ratio_Fe2DOC_l"
@@ -1038,12 +1038,12 @@
 			possible_values="zero or positive real number"
 			icepack_name="tau_max"
 		/>
-		<nml_option name="config_algal_maximum_velocity" type="real" default_value="0.0000001" units="m/s"
+		<nml_option name="config_algal_maximum_velocity" type="real" default_value="0.0000001" units="m s-1"
 			description="Maximum speed at which diatoms can move"
 			possible_values="zero or positive real number"
 			icepack_name="algal_vel"
 		/>
-		<nml_option name="config_ratio_Fe_to_dust" type="real" default_value="0.035" units="mg iron/g dust"
+		<nml_option name="config_ratio_Fe_to_dust" type="real" default_value="0.035" units="g Fe (g dust)-1"
 			description="Mass ratio of iron in dust"
 			possible_values="zero or positive real number"
 			icepack_name="R_dFe2dust"
@@ -1053,77 +1053,77 @@
 			possible_values="real number between and including zero and one"
 			icepack_name="dustFe_sol"
 		/>
-		<nml_option name="config_chla_absorptivity_of_diatoms" type="real" default_value="0.03" units="1/m/(mg chla/m^3)"
+		<nml_option name="config_chla_absorptivity_of_diatoms" type="real" default_value="0.03" units="m-1 (mg chla)-1 m-3"
 			description="Absorptivity of chlorophyll in diatoms"
 			possible_values="positive real"
 			icepack_name="chlabs_diatoms"
 		/>
-		<nml_option name="config_chla_absorptivity_of_small_plankton" type="real" default_value="0.01" units="1/m/(mg chla/m^3)"
+		<nml_option name="config_chla_absorptivity_of_small_plankton" type="real" default_value="0.01" units="m-1 (mg chla)-1 m-3"
 			description="Absorptivity of chlorophyll in small plankton"
 			possible_values="positive real"
 			icepack_name="chlabs_sp"
 		/>
-		<nml_option name="config_chla_absorptivity_of_phaeocystis" type="real" default_value="0.05" units="1/m/(mg chla/m^3)"
+		<nml_option name="config_chla_absorptivity_of_phaeocystis" type="real" default_value="0.05" units="m-1 (mg chla)-1 m-3"
 			description="Absorptivity of chlorophyll in phaeocystis"
 			possible_values="positive real"
 			icepack_name="chlabs_phaeo"
 		/>
-		<nml_option name="config_light_attenuation_diatoms" type="real" default_value="0.3" units="1/(W/m^2)"
+		<nml_option name="config_light_attenuation_diatoms" type="real" default_value="0.3" units="m2 W-1"
 			description="Controls the ability of diatoms to utilize light"
 			possible_values="positive real less than 1"
 			icepack_name="alpha2max_low_diatoms"
 		/>
-		<nml_option name="config_light_attenuation_small_plankton" type="real" default_value="0.2" units="1/(W/m^2)"
+		<nml_option name="config_light_attenuation_small_plankton" type="real" default_value="0.2" units="m2 W-1"
 			description="Controls the ability of small plankton to utilize light"
 			possible_values="positive real less than 1"
 			icepack_name="alpha2max_low_sp"
 		/>
-		<nml_option name="config_light_attenuation_phaeocystis" type="real" default_value="0.17" units="1/(W/m^2)"
+		<nml_option name="config_light_attenuation_phaeocystis" type="real" default_value="0.17" units="m2 W-1"
 			description="Controls the ability of phaeocystis to utilize light"
 			possible_values="positive real less than 1"
 			icepack_name="alpha2max_low_phaeo"
 		/>
-		<nml_option name="config_light_inhibition_diatoms" type="real" default_value="0.001" units="1/(W/m^2)"
+		<nml_option name="config_light_inhibition_diatoms" type="real" default_value="0.001" units="m2 W-1"
 			description="Controls the strength of light inhibition to growth for diatoms"
 			possible_values="positive real less than 1"
 			icepack_name="beta2max_diatoms"
 		/>
-		<nml_option name="config_light_inhibition_small_plankton" type="real" default_value="0.001" units="1/(W/m^2)"
+		<nml_option name="config_light_inhibition_small_plankton" type="real" default_value="0.001" units="m2 W-1"
 			description="Controls the strength of light inhibition to growth for small plankton"
 			possible_values="positive real less than 1"
 			icepack_name="beta2max_sp"
 		/>
-		<nml_option name="config_light_inhibition_phaeocystis" type="real" default_value="0.04" units="1/(W/m^2)"
+		<nml_option name="config_light_inhibition_phaeocystis" type="real" default_value="0.04" units="m2 W-1"
 			description="Controls the strength of light inhibition to growth for phaeocystis"
 			possible_values="positive real less than 1"
 			icepack_name="beta2max_phaeo"
 		/>
-		<nml_option name="config_maximum_growth_rate_diatoms" type="real" default_value="1.44" units="1/day"
+		<nml_option name="config_maximum_growth_rate_diatoms" type="real" default_value="1.44" units="day-1"
 			description="Maximum growth rate for diatoms"
 			possible_values="positive real"
 			icepack_name="mu_max_diatoms"
 		/>
-		<nml_option name="config_maximum_growth_rate_small_plankton" type="real" default_value="0.41" units="1/day"
+		<nml_option name="config_maximum_growth_rate_small_plankton" type="real" default_value="0.41" units="day-1"
 			description="Maximum growth rate for small plankton"
 			possible_values="positive real"
 			icepack_name="mu_max_sp"
 		/>
-		<nml_option name="config_maximum_growth_rate_phaeocystis" type="real" default_value="0.63" units="1/day"
+		<nml_option name="config_maximum_growth_rate_phaeocystis" type="real" default_value="0.63" units="day-1"
 			description="Maximum growth rate for phaeocystis"
 			possible_values="positive real"
 			icepack_name="mu_max_phaeo"
 		/>
-		<nml_option name="config_temperature_growth_diatoms" type="real" default_value="0.063" units="1/C"
+		<nml_option name="config_temperature_growth_diatoms" type="real" default_value="0.063" units="C-1"
 			description="Controls the temperature dependence of growth for diatoms"
 			possible_values="positive real less than 1"
 			icepack_name="grow_Tdep_diatoms"
 		/>
-		<nml_option name="config_temperature_growth_small_plankton" type="real" default_value="0.063" units="1/C"
+		<nml_option name="config_temperature_growth_small_plankton" type="real" default_value="0.063" units="C-1"
 			description="Controls the temperature dependence of growth for small plankton"
 			possible_values="positive real less than 1"
 			icepack_name="grow_Tdep_sp"
 		/>
-		<nml_option name="config_temperature_growth_phaeocystis" type="real" default_value="0.063" units="1/C"
+		<nml_option name="config_temperature_growth_phaeocystis" type="real" default_value="0.063" units="C-1"
 			description="Controls the temperature dependence of growth for phaeocystis"
 			possible_values="positive real less than 1"
 			icepack_name="grow_Tdep_phaeo"
@@ -1143,107 +1143,108 @@
 			possible_values="zero or positive real less than 1"
 			icepack_name="fr_graze_phaeo"
 		/>
-		<nml_option name="config_mortality_diatoms" type="real" default_value="0.007" units="1/day"
+		<nml_option name="config_mortality_diatoms" type="real" default_value="0.007" units="day-1"
 			description="Mortality rate of diatoms"
 			possible_values="positive real"
 			icepack_name="mort_pre_diatoms"
 		/>
-		<nml_option name="config_mortality_small_plankton" type="real" default_value="0.007" units="1/day"
+		<nml_option name="config_mortality_small_plankton" type="real" default_value="0.007" units="day-1"
 			description="Mortality rate of small plankton"
 			possible_values="positive real"
 			icepack_name="mort_pre_sp"
 		/>
-		<nml_option name="config_mortality_phaeocystis" type="real" default_value="0.007" units="1/day"
+		<nml_option name="config_mortality_phaeocystis" type="real" default_value="0.007" units="day-1"
 			description="Mortality rate of phaeocystis"
 			possible_values="positive real"
 			icepack_name="mort_pre_phaeo"
 		/>
-		<nml_option name="config_temperature_mortality_diatoms" type="real" default_value="0.03" units="1/C"
+		<nml_option name="config_temperature_mortality_diatoms" type="real" default_value="0.03" units="C-1"
 			description="Temperature dependence of diatom mortality rate"
 			possible_values="positive real"
 			icepack_name="mort_Tdep_diatoms"
 		/>
-		<nml_option name="config_temperature_mortality_small_plankton" type="real" default_value="0.03" units="1/C"
+		<nml_option name="config_temperature_mortality_small_plankton" type="real" default_value="0.03" units="C-1"
 			description="Temperature dependence of small plankton mortality rate"
 			possible_values="positive real"
 			icepack_name="mort_Tdep_sp"
 		/>
-		<nml_option name="config_temperature_mortality_phaeocystis" type="real" default_value="0.03" units="1/C"
+		<nml_option name="config_temperature_mortality_phaeocystis" type="real" default_value="0.03" units="C-1"
 			description="Temperature dependence of phaeocystis mortality rate"
 			possible_values="positive real"
 			icepack_name="mort_Tdep_phaeo"
 		/>
-		<nml_option name="config_exudation_diatoms" type="real" default_value="0.0" units="1/day"
+		<nml_option name="config_exudation_diatoms" type="real" default_value="0.0" units="day-1"
 			description="Rate diatoms exude DOC"
 			possible_values="zero or positive real"
 			icepack_name="k_exude_diatoms"
 		/>
-		<nml_option name="config_exudation_small_plankton" type="real" default_value="0.0" units="1/day"
+		<nml_option name="config_exudation_small_plankton" type="real" default_value="0.0" units="day-1"
 			description="Rate small plankton exude DOC"
 			possible_values="zero or positive real"
 			icepack_name="k_exude_sp"
 		/>
-		<nml_option name="config_exudation_phaeocystis" type="real" default_value="0.0" units="1/day"
+		<nml_option name="config_exudation_phaeocystis" type="real" default_value="0.0" units="day-1"
 			description="Rate phaeocystis exude DOC"
 			possible_values="zero or positive real"
 			icepack_name="k_exude_phaeo"
 		/>
-		<nml_option name="config_nitrate_saturation_diatoms" type="real" default_value="1.0" units="mmol nitrogen/m^3"
+		<nml_option name="config_nitrate_saturation_diatoms" type="real" default_value="1.0" units="mmol N m-3"
 			description="Nitrate half saturation for diatom growth"
 			possible_values="positive real"
 			icepack_name="K_Nit_diatoms"
 		/>
-		<nml_option name="config_nitrate_saturation_small_plankton" type="real" default_value="1.0" units="mmol nitrogen/m^3"
+		<nml_option
+		    name="config_nitrate_saturation_small_plankton" type="real" default_value="1.0" units="mmol N m-3"
 			description="Nitrate half saturation for small plankton growth"
 			possible_values="positive real"
 			icepack_name="K_Nit_sp"
 		/>
-		<nml_option name="config_nitrate_saturation_phaeocystis" type="real" default_value="1.0" units="mmol nitrogen/m^3"
+		<nml_option name="config_nitrate_saturation_phaeocystis" type="real" default_value="1.0" units="mmol N m-3"
 			description="Nitrate half saturation for phaeocystis growth"
 			possible_values="positive real"
 			icepack_name="K_Nit_phaeo"
 		/>
-		<nml_option name="config_ammonium_saturation_diatoms" type="real" default_value="0.3" units="mmol nitrogen/m^3"
+		<nml_option name="config_ammonium_saturation_diatoms" type="real" default_value="0.3" units="mmol N m-3"
 			description="Ammonium half saturation for diatom growth"
 			possible_values="positive real"
 			icepack_name="K_Am_diatoms"
 		/>
-		<nml_option name="config_ammonium_saturation_small_plankton" type="real" default_value="0.3" units="mmol nitrogen/m^3"
+		<nml_option name="config_ammonium_saturation_small_plankton" type="real" default_value="0.3" units="mmol N m-3"
 			description="Ammonium half saturation for small plankton growth"
 			possible_values="positive real"
 			icepack_name="K_Am_sp"
 		/>
-		<nml_option name="config_ammonium_saturation_phaeocystis" type="real" default_value="0.3" units="mmol nitrogen/m^3"
+		<nml_option name="config_ammonium_saturation_phaeocystis" type="real" default_value="0.3" units="mmol N m-3"
 			description="Ammonium half saturation for phaeocystis growth"
 			possible_values="positive real"
 			icepack_name="K_Am_phaeo"
 		/>
-		<nml_option name="config_silicate_saturation_diatoms" type="real" default_value="4.0" units="mmol Si/m^3"
+		<nml_option name="config_silicate_saturation_diatoms" type="real" default_value="4.0" units="mmol Si m-3"
 			description="Silicate half saturation for diatom growth"
 			possible_values="zero or positive real"
 			icepack_name="K_Sil_diatoms"
 		/>
-		<nml_option name="config_silicate_saturation_small_plankton" type="real" default_value="0.0" units="mmol Si/m^3"
+		<nml_option name="config_silicate_saturation_small_plankton" type="real" default_value="0.0" units="mmol Si m-3"
 			description="Silicate half saturation for small plankton growth"
 			possible_values="zero or positive real"
 			icepack_name="K_Sil_sp"
 		/>
-		<nml_option name="config_silicate_saturation_phaeocystis" type="real" default_value="0.0" units="mmol Si/m^3"
+		<nml_option name="config_silicate_saturation_phaeocystis" type="real" default_value="0.0" units="mmol Si m-3"
 			description="Silicate half saturation for phaeocystis growth"
 			possible_values="zero or positive real"
 			icepack_name="K_Sil_phaeo"
 		/>
-		<nml_option name="config_iron_saturation_diatoms" type="real" default_value="1.0" units="umol iron/m^3"
+		<nml_option name="config_iron_saturation_diatoms" type="real" default_value="1.0" units="umol Fe m-3"
 			description="Iron half saturation for diatom growth"
 			possible_values="positive real"
 			icepack_name="K_Fe_diatoms"
 		/>
-		<nml_option name="config_iron_saturation_small_plankton" type="real" default_value="0.2" units="umol iron/m^3"
+		<nml_option name="config_iron_saturation_small_plankton" type="real" default_value="0.2" units="umol Fe m-3"
 			description="Iron half saturation for small plankton growth"
 			possible_values="positive real"
 			icepack_name="K_Fe_sp"
 		/>
-		<nml_option name="config_iron_saturation_phaeocystis" type="real" default_value="0.1" units="umol iron/m^3"
+		<nml_option name="config_iron_saturation_phaeocystis" type="real" default_value="0.1" units="umol Fe m-3"
 			description="Iron half saturation for phaeocystis growth"
 			possible_values="positive real"
 			icepack_name="K_Fe_phaeo"
@@ -1253,7 +1254,7 @@
 			possible_values="zero or positive real less than or equal to 1"
 			icepack_name="f_don_protein"
 		/>
-		<nml_option name="config_degredation_of_DON" type="real" default_value="0.2" units="1/day"
+		<nml_option name="config_degredation_of_DON" type="real" default_value="0.2" units="day-1"
 			description="Bacterial degredation of DON"
 			possible_values="positive real"
 			icepack_name="kn_bac_protein"
@@ -1283,12 +1284,12 @@
 			possible_values="zero or positive real less than or equal to 1"
 			icepack_name="f_exude_l"
 		/>
-		<nml_option name="config_remineralization_saccharids" type="real" default_value="0.03" units="1/day"
+		<nml_option name="config_remineralization_saccharids" type="real" default_value="0.03" units="day-1"
 			description="Bacterial degredaton rate of poly/saccarids"
 			possible_values="positive real"
 			icepack_name="k_bac_s"
 		/>
-		<nml_option name="config_remineralization_lipids" type="real" default_value="0.03" units="1/day"
+		<nml_option name="config_remineralization_lipids" type="real" default_value="0.03" units="day-1"
 			description="Bacterial degredation rate of lipids"
 			possible_values="positive real"
 			icepack_name="k_bac_l"
@@ -1328,12 +1329,12 @@
 			possible_values="zero or positive real less than or equal to 1"
 			icepack_name="fr_dFe"
 		/>
-		<nml_option name="config_nitrification_rate" type="real" default_value="0.046" units="1/day"
+		<nml_option name="config_nitrification_rate" type="real" default_value="0.046" units="day-1"
 			description="Nitrification rate"
 			possible_values="zero or positive real"
 			icepack_name="k_nitrif"
 		/>
-		<nml_option name="config_desorption_loss_particulate_iron" type="real" default_value="3065.0" units="day"
+		<nml_option name="config_desorption_loss_particulate_iron" type="real" default_value="3065.0" units="day-1"
 			description="Desorption rate of particulate iron to dissolved iron"
 			possible_values="positive real"
 			icepack_name="t_iron_conv"
@@ -1409,7 +1410,7 @@
 			icepack_name="dmspdtye"
 		/>
 		<nml_option name="config_mobility_type_humics" type="real" default_value="0.0" units="unitless"
-			description="Transport type of humics"
+			description="Transport type of humics (refractory dissolved organic matter)"
 			possible_values="-1 = entirely in the mobile phase; 0 = retention dominated; 1 = release dominated; 0.5 = equal but rapid exchange; 2 = equal but slow exchange"
 			icepack_name="humtypes"
 		/>
@@ -1473,32 +1474,32 @@
 			possible_values="-1 = entirely in the mobile phase; 0 = retention dominated; 1 = release dominated; 0.5 = equal but rapid exchange; 2 = equal but slow exchange"
 			icepack_name="zaerotype_dust4"
 		/>
-		<nml_option name="config_ratio_C_to_N_diatoms" type="real" default_value="7.0" units="mol carbon/mol nitrogen"
+		<nml_option name="config_ratio_C_to_N_diatoms" type="real" default_value="7.0" units="mol C (mol N)-1"
 			description="Algal carbon to nitrogen mole ratio for diatoms"
 			possible_values="positive real"
 			icepack_name="ratio_C2N_diatoms"
 		/>
-		<nml_option name="config_ratio_C_to_N_small_plankton" type="real" default_value="7.0" units="mol carbon/mol nitrogen"
+		<nml_option name="config_ratio_C_to_N_small_plankton" type="real" default_value="7.0" units="mol C (mol N)-1"
 			description="Algal carbon to nitrogen mole ratio for small plankton"
 			possible_values="positive real"
 			icepack_name="ratio_C2N_sp"
 		/>
-		<nml_option name="config_ratio_C_to_N_phaeocystis" type="real" default_value="7.0" units="mol carbon/mol nitrogen"
+		<nml_option name="config_ratio_C_to_N_phaeocystis" type="real" default_value="7.0" units="mol C (mol N)-1"
 			description="Algal carbon to nitrogen mole ratio for phaeocystis"
 			possible_values="positive real"
 			icepack_name="ratio_C2N_phaeo"
 		/>
-		<nml_option name="config_ratio_chla_to_N_diatoms" type="real" default_value="2.1" units="g chla/mol nitrogen"
+		<nml_option name="config_ratio_chla_to_N_diatoms" type="real" default_value="2.1" units="g chla (mol N)-1"
 			description="Algal chlorophyll to nitrogen ratio for diatoms"
 			possible_values="positive real"
 			icepack_name="ratio_chl2N_diatoms"
 		/>
-		<nml_option name="config_ratio_chla_to_N_small_plankton" type="real" default_value="1.1" units="g chla/mol nitrogen"
+		<nml_option name="config_ratio_chla_to_N_small_plankton" type="real" default_value="1.1" units="g chla (mol N)-1"
 			description="Algal chlorophyll to nitrogen ratio for small plankton"
 			possible_values="positive real"
 			icepack_name="ratio_chl2N_sp"
 		/>
-		<nml_option name="config_ratio_chla_to_N_phaeocystis" type="real" default_value="0.84" units="g chla/mol nitrogen"
+		<nml_option name="config_ratio_chla_to_N_phaeocystis" type="real" default_value="0.84" units="g chla (mol N)-1"
 			description="Algal chlorophyll to nitrogen ratio for phaeocystis"
 			possible_values="positive real"
 			icepack_name="ratio_chl2N_phaeo"
@@ -1518,7 +1519,7 @@
 			possible_values="positive real"
 			icepack_name="F_abs_chl_phaeo"
 		/>
-		<nml_option name="config_ratio_C_to_N_proteins" type="real" default_value="5.0" units="mol carbon/mol nitrogen"
+		<nml_option name="config_ratio_C_to_N_proteins" type="real" default_value="5.0" units="mol C (mol N)-1"
 			description="Ratio of carbon to nitrogen in proteins"
 			possible_values="positive real"
 			icepack_name="ratio_C2N_proteins"
@@ -1565,17 +1566,17 @@
 			possible_values="Any positive real number."
 			icepack_name="ahmax"
 		/>
-		<nml_option name="config_ice_shortwave_tuning_parameter" type="real" default_value="0.0" units=""
+		<nml_option name="config_ice_shortwave_tuning_parameter" type="real" default_value="0.0" units="unitless"
 			description="Sea ice tuning parameter; +1 for 1sig increase in albedo."
 			possible_values=""
 			icepack_name="R_ice"
 		/>
-		<nml_option name="config_pond_shortwave_tuning_parameter" type="real" default_value="0.0" units=""
+		<nml_option name="config_pond_shortwave_tuning_parameter" type="real" default_value="0.0" units="unitless"
 			description="Ponded ice tuning parameter; +1 for 1sig increase in albedo."
 			possible_values=""
 			icepack_name="R_pnd"
 		/>
-		<nml_option name="config_snow_shortwave_tuning_parameter" type="real" default_value="1.5" units=""
+		<nml_option name="config_snow_shortwave_tuning_parameter" type="real" default_value="1.5" units="unitless"
 			description="Snow tuning parameter; +1 for ~.01 change in broadband albedo."
 			possible_values=""
 			icepack_name="R_snw"
@@ -1590,7 +1591,7 @@
 			possible_values=""
 			icepack_name="rsnw_mlt"
 		/>
-		<nml_option name="config_algae_absorption_coefficient" type="real" default_value="0.6" units=""
+		<nml_option name="config_algae_absorption_coefficient" type="real" default_value="0.6" units="unitless"
 			description="Algae absorption coefficient for 0.5 m thick layer."
 			possible_values=""
 			icepack_name="kalg"
@@ -1600,7 +1601,7 @@
 			possible_values="true or false"
 			icepack_name="sw_redist"
 		/>
-		<nml_option name="config_shortwave_redistribution_fraction" type="real" default_value="0.9" units="none"
+		<nml_option name="config_shortwave_redistribution_fraction" type="real" default_value="0.9" units="unitless"
 			description="Fraction of shortwave moved from a layer near the melting temperature to the surface."
 			possible_values="fraction between 0 and 1"
 			icepack_name="sw_frac"
@@ -1628,27 +1629,28 @@
 			possible_values="true or false"
 			icepack_name="use_smliq_pnd"
 		/>
-		<nml_option name="config_new_snow_density" type="real" default_value="100.0" units="kg/m^3"
+		<nml_option name="config_new_snow_density" type="real" default_value="100.0" units="kg m-3"
 			description="new snow density"
 			possible_values="positive real"
 			icepack_name="rhosnew"
 		/>
-		<nml_option name="config_max_snow_density" type="real" default_value="450.0" units="kg/m^3"
+		<nml_option name="config_max_snow_density" type="real" default_value="450.0" units="kg m-3"
 			description="maximum snow density"
 			possible_values="positive real"
 			icepack_name="rhosmax"
 		/>
-		<nml_option name="config_minimum_wind_compaction" type="real" default_value="10.0" units="m/s"
+		<nml_option name="config_minimum_wind_compaction" type="real" default_value="10.0" units="m s-1"
 			description="minimum wind speed to compact snow"
 			possible_values="positive real"
 			icepack_name="windmin"
 		/>
-		<nml_option name="config_wind_compaction_factor" type="real" default_value="27.3" units="kg s/m^4"
+		<nml_option name="config_wind_compaction_factor"
+			    type="real" default_value="27.3" units="kg s m-4"
 			description="wind compaction factor"
 			possible_values="positive real"
 			icepack_name="drhosdwind"
 		/>
-		<nml_option name="config_snow_redistribution_factor" type="real" default_value="0.3" units="none"
+		<nml_option name="config_snow_redistribution_factor" type="real" default_value="0.3" units="unitless"
 			description="snow loss factor for wind redistribution"
 			possible_values="positive real less than or equal to 1"
 			icepack_name="snwlvlfac"
@@ -1658,7 +1660,8 @@
 			possible_values="positive real"
 			icepack_name="rsnw_tmax"
 		/>
-		<nml_option name="config_snow_thermal_conductivity" type="real" default_value="0.3" units="W/m/deg"
+		<nml_option name="config_snow_thermal_conductivity"
+			    type="real" default_value="0.3" units="W m-1 deg-1"
 			description="thermal conductivity of snow"
 			possible_values="positive real"
 			icepack_name="ksno"
@@ -2689,8 +2692,9 @@
 		<package name="pkgColumnTracerAerosols" description="Evolves aerosol tracers in ice and snow"/>
 		<package name="pkgColumnTracerEffectiveSnowDensity" description="Computes effective snow density based on snow compaction"/>
 		<package name="pkgColumnTracerSnowGrainRadius" description="Evolves snow grain radius for radiative transfer based on dry and wet metamorphism"/>
-		<package name="pkgTracerBrine" description=""/>
-		<package name="pkgTracerMobileFraction" description=""/>
+		<package name="pkgTracerBrine" description="brine fraction tracer (ratio of brine penetration height to ice thickness)"/>
+		<package name="pkgTracerMobileFraction"
+			 description="fraction of vertical bgc/aerosol tracer freely moving in the brine"/>
 		<package name="pkgTracerSkeletalAlgae" description=""/>
 		<package name="pkgTracerSkeletalNitrate" description=""/>
 		<package name="pkgTracerSkeletalCarbon" description=""/>
@@ -2701,17 +2705,17 @@
 		<package name="pkgTracerSkeletalHumics" description=""/>
 		<package name="pkgTracerSkeletalDON" description=""/>
 		<package name="pkgTracerSkeletalIron" description=""/>
-		<package name="pkgTracerVerticalAlgae" description=""/>
-		<package name="pkgTracerVerticalNitrate" description=""/>
-		<package name="pkgTracerVerticalCarbon" description=""/>
-		<package name="pkgTracerVerticalAmmonium" description=""/>
-		<package name="pkgTracerVerticalSilicate" description=""/>
-		<package name="pkgTracerVerticalDMS" description=""/>
-		<package name="pkgTracerVerticalNonreactive" description=""/>
-		<package name="pkgTracerVerticalHumics" description=""/>
-		<package name="pkgTracerVerticalDON" description=""/>
-		<package name="pkgTracerVerticalIron" description=""/>
-		<package name="pkgTracerZAerosols" description=""/>
+		<package name="pkgTracerVerticalAlgae" description="vertical algal nitrogen tracers"/>
+		<package name="pkgTracerVerticalNitrate" description="vertical nitrate tracer"/>
+		<package name="pkgTracerVerticalCarbon" description="vertical carbon tracers"/>
+		<package name="pkgTracerVerticalAmmonium" description="vertical ammonium tracer"/>
+		<package name="pkgTracerVerticalSilicate" description="vertical silicate tracer"/>
+		<package name="pkgTracerVerticalDMS" description="vertical sulfur tracers"/>
+		<package name="pkgTracerVerticalNonreactive" description="vertical nonreactive nitrate tracer"/>
+		<package name="pkgTracerVerticalHumics" description="vertical humics (refractory dissolved organic carbon) tracer"/>
+		<package name="pkgTracerVerticalDON" description="vertical dissolved organic nitrogen tracers"/>
+		<package name="pkgTracerVerticalIron" description="vertical iron tracers"/>
+		<package name="pkgTracerZAerosols" description="vertical aerosol tracers"/>
 		<package name="pkgTracerZSalinity" description=""/>
 
 		<!-- other column packages -->
@@ -2885,11 +2889,13 @@
 		<var name="brineFraction"
 		     type="real"
 		     dimensions="ONE nCategories nCells Time"
+		     units="1"
 		     packages="pkgTracerBrine"
 		/>
 		<var name="mobileFraction"
 		     type="real"
 		     dimensions="nZBGCTracers nCategories nCells Time"
+		     units="1"
 		     packages="pkgTracerMobileFraction"
 		/>
 		<var name="skeletalAlgaeConc"
@@ -2965,226 +2971,271 @@
 		<var name="verticalAlgaeConc"
 		     type="real"
 		     dimensions="nAlgaeLayers nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalAlgae"
 		/>
 		<var name="verticalDOCConc"
 		     type="real"
 		     dimensions="nDOCLayers nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDICConc"
 		     type="real"
 		     dimensions="nDICLayers nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDONConc"
 		     type="real"
 		     dimensions="nDONLayers nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalDON"
 		/>
 		<var name="verticalNitrateConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalNitrate"
 		/>
 		<var name="verticalSilicateConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol Si m-3"
 		     packages="pkgTracerVerticalSilicate"
 		/>
 		<var name="verticalAmmoniumConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalAmmonium"
 		/>
 		<var name="verticalDMSConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPpConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPdConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalNonreactiveConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalNonreactive"
 		/>
 		<var name="verticalHumicsConc"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalHumics"
 		/>
 		<var name="verticalParticulateIronConc"
 		     type="real"
 		     dimensions="nParticulateIronLayers nCategories nCells Time"
+		     units="umol Fe m-3"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalDissolvedIronConc"
 		     type="real"
 		     dimensions="nDissolvedIronLayers nCategories nCells Time"
+		     units="umol Fe m-3"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalAlgaeSnow"
 		     type="real"
 		     dimensions="nAlgaeSnowLayers nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalAlgae"
 		/>
 		<var name="verticalDOCSnow"
 		     type="real"
 		     dimensions="nDOCSnowLayers nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDICSnow"
 		     type="real"
 		     dimensions="nDICSnowLayers nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDONSnow"
 		     type="real"
 		     dimensions="nDONSnowLayers nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalDON"
 		/>
 		<var name="verticalNitrateSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalNitrate"
 		/>
 		<var name="verticalSilicateSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol Si m-3"
 		     packages="pkgTracerVerticalSilicate"
 		/>
 		<var name="verticalAmmoniumSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalAmmonium"
 		/>
 		<var name="verticalDMSSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPpSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPdSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalNonreactiveSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalNonreactive"
 		/>
 		<var name="verticalHumicsSnow"
 		     type="real"
 		     dimensions="TWO nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalHumics"
 		/>
 		<var name="verticalParticulateIronSnow"
 		     type="real"
 		     dimensions="nPartIronSnowLayers nCategories nCells Time"
+		     units="umol Fe m-3"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalDissolvedIronSnow"
 		     type="real"
 		     dimensions="nDisIronSnowLayers nCategories nCells Time"
+		     units="umol Fe m-3"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalAlgaeIce"
 		     type="real"
 		     dimensions="nAlgaeIceLayers nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalAlgae"
 		/>
 		<var name="verticalDOCIce"
 		     type="real"
 		     dimensions="nDOCIceLayers nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDICIce"
 		     type="real"
 		     dimensions="nDICIceLayers nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDONIce"
 		     type="real"
 		     dimensions="nDONIceLayers nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalDON"
 		/>
 		<var name="verticalNitrateIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalNitrate"
 		/>
 		<var name="verticalSilicateIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol Si m-3"
 		     packages="pkgTracerVerticalSilicate"
 		/>
 		<var name="verticalAmmoniumIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalAmmonium"
 		/>
 		<var name="verticalDMSIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPpIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPdIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol S m-3"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalNonreactiveIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol N m-3"
 		     packages="pkgTracerVerticalNonreactive"
 		/>
 		<var name="verticalHumicsIce"
 		     type="real"
 		     dimensions="nBioLayersP1 nCategories nCells Time"
+		     units="mmol C m-3"
 		     packages="pkgTracerVerticalHumics"
 		/>
 		<var name="verticalParticulateIronIce"
 		     type="real"
 		     dimensions="nPartIronIceLayers nCategories nCells Time"
+		     units="umol Fe m-3"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalDissolvedIronIce"
 		     type="real"
 		     dimensions="nDisIronIceLayers nCategories nCells Time"
+		     units="umol Fe m-3"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalAerosolsConc"
 		     type="real"
 		     dimensions="nzAerosolsLayers nCategories nCells Time"
+		     units="kg m-3"
 		     packages="pkgTracerZAerosols"
 		/>
 		<var name="verticalAerosolsSnow"
 		     type="real"
 		     dimensions="nzAerosolsSnowLayers nCategories nCells Time"
+		     units="kg m-3"
 		     packages="pkgTracerZAerosols"
 		/>
 		<var name="verticalAerosolsIce"
 		     type="real"
 		     dimensions="nzAerosolsIceLayers nCategories nCells Time"
+		     units="kg m-3"
 		     packages="pkgTracerZAerosols"
 		/>
 		<var name="verticalSalinity"
@@ -3346,6 +3397,7 @@
 		     type="real"
 		     dimensions="nCells Time"
 		     description="Grid cell average ratio of brine height to ice thickness"
+		     units="1"
 		     packages="pkgTracerBrine"
 		     icepack_name="fbrine"
 		/>
@@ -3422,241 +3474,325 @@
 		<var name="verticalAlgaeConcCell"
 		     type="real"
 		     dimensions="nAlgaeLayers nCells Time"
+		     units="mmol N m-3"
+		     description="Average algal nitrogen concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalAlgae"
 		/>
 		<var name="verticalDOCConcCell"
 		     type="real"
 		     dimensions="nDOCLayers nCells Time"
+		     units="mmol C m-3"
+		     description="Average dissolved organic carbon (polysaccarids, lipids) concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDICConcCell"
 		     type="real"
 		     dimensions="nDICLayers nCells Time"
+		     units="mmol C m-3"
+		     description="Average dissolved inorganic carbon concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDONConcCell"
 		     type="real"
 		     dimensions="nDONLayers nCells Time"
+		     units="mmol N m-3"
+		     description="Average dissolved organic nitrogen (proteins) concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalDON"
 		/>
 		<var name="verticalNitrateConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol N m-3"
+		     description="Average nitrate concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalNitrate"
 		/>
 		<var name="verticalSilicateConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol Si m-3"
+		     description="Average silicate concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalSilicate"
 		/>
 		<var name="verticalAmmoniumConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol N m-3"
+		     description="Average ammonium concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalAmmonium"
 		/>
 		<var name="verticalDMSConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol S m-3"
+		     description="Average dimethyl sulfide concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPpConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol S m-3"
+		     description="Average particulate dimethylsulfoniopropionate concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPdConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol S m-3"
+		     description="Average dissolved dimethylsulfoniopropionate concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalNonreactiveConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol N m-3"
+		     description="Average nonreactive nitrate concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalNonreactive"
 		/>
 		<var name="verticalHumicsConcCell"
 		     type="real"
 		     dimensions="nBioLayersP3 nCells Time"
+		     units="mmol C m-3"
+		     description="Average refractory dissolved organic carbon concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalHumics"
 		/>
 		<var name="verticalParticulateIronConcCell"
 		     type="real"
 		     dimensions="nParticulateIronLayers nCells Time"
+		     units="umol Fe m-3"
+		     description="Average particulate iron concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalDissolvedIronConcCell"
 		     type="real"
 		     dimensions="nDissolvedIronLayers nCells Time"
+		     units="umol Fe m-3"
+		     description="Average dissolved iron concentration in snow and ice bio layers"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalAlgaeSnowCell"
 		     type="real"
 		     dimensions="nAlgaeSnowLayers nCells Time"
+		     units="mmol N m-3"
+		     description="Average algal nitrogen concentration in snow bio layers"
 		     packages="pkgTracerVerticalAlgae"
 		/>
 		<var name="verticalDOCSnowCell"
 		     type="real"
 		     dimensions="nDOCSnowLayers nCells Time"
+		     units="mmol C m-3"
+		     description="Average dissolved organic carbon (polysaccarids, lipids) concentration in snow bio layers"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDICSnowCell"
 		     type="real"
 		     dimensions="nDICSnowLayers nCells Time"
+		     units="mmol C m-3"
+		     description="Average dissolved inorganic carbon concentration in snow bio layers"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDONSnowCell"
 		     type="real"
 		     dimensions="nDONSnowLayers nCells Time"
+		     units="mmol N m-3"
+		     description="Average dissolved organic nitrogen (proteins) concentration in snow bio layers"
 		     packages="pkgTracerVerticalDON"
 		/>
 		<var name="verticalNitrateSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol N m-3"
+		     description="Average nitrate concentration in snow bio layers"
 		     packages="pkgTracerVerticalNitrate"
 		/>
 		<var name="verticalSilicateSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol Si m-3"
+		     description="Average silicate concentration in snow bio layers"
 		     packages="pkgTracerVerticalSilicate"
 		/>
 		<var name="verticalAmmoniumSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol N m-3"
+		     description="Average ammonium concentration in snow bio layers"
 		     packages="pkgTracerVerticalAmmonium"
 		/>
 		<var name="verticalDMSSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol S m-3"
+		     description="Average dimethyl sulfide concentration in snow bio layers"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPpSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol S m-3"
+		     description="Average particulate dimethylsulfoniopropionate concentration in snow bio layers"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPdSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol S m-3"
+		     description="Average dissolved dimethylsulfoniopropionate concentration in snow bio layers"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalNonreactiveSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol N m-3"
+		     description="Average nonreactive nitrate concentration in snow bio layers"
 		     packages="pkgTracerVerticalNonreactive"
 		/>
 		<var name="verticalHumicsSnowCell"
 		     type="real"
 		     dimensions="TWO nCells Time"
+		     units="mmol C m-3"
+		     description="Average refractory dissolved organic carbon concentration in snow bio layers"
 		     packages="pkgTracerVerticalHumics"
 		/>
 		<var name="verticalParticulateIronSnowCell"
 		     type="real"
 		     dimensions="nPartIronSnowLayers nCells Time"
+		     units="umol Fe m-3"
+		     description="Average particulate iron concentration in snow bio layers"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalDissolvedIronSnowCell"
 		     type="real"
 		     dimensions="nDisIronSnowLayers nCells Time"
+		     units="umol Fe m-3"
+		     description="Average dissolved iron concentration in snow bio layers"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalAlgaeIceCell"
 		     type="real"
 		     dimensions="nAlgaeIceLayers nCells Time"
-		     units="mmol m-3"
+		     units="mmol N m-3"
 		     description="Cell average algal nitrogen for each type and ice bio-grid layer"
 		     packages="pkgTracerVerticalAlgae"
 		/>
+		<var name="verticalAlgaeTotalCarbonIceCell"
+		     type="real"
+		     dimensions="nBioLayersP1 nCells Time"
+		     units="mmol C m-3"
+		     description="Cell average total algal carbon in ice bio-grid layer"
+		     packages="pkgTracerVerticalAlgae"
+		/>
 		<var name="verticalDOCIceCell"
 		     type="real"
 		     dimensions="nDOCIceLayers nCells Time"
-		     units="mmol m-3"
+		     units="mmol C m-3"
 		     description="Cell average dissolved organic carbon for each type and ice bio-grid layer"
 		     packages="pkgTracerVerticalCarbon"
 		/>
+		<var name="verticalDOCLabileIceCell"
+		     type="real"
+		     dimensions="nBioLayersP1 nCells Time"
+		     units="mmol C m-3"
+		     description="Cell average total semi-labile dissolved organic carbon (polysaccharids+lipids+proteins)in ice bio-grid layer"
+		     packages="pkgTracerVerticalCarbon"
+		/>
 		<var name="verticalDICIceCell"
 		     type="real"
 		     dimensions="nDICIceLayers nCells Time"
-		     units="mmol m-3"
+		     units="mmol C m-3"
 		     description="Cell average dissolved inorganic carbon for each type and ice bio-grid layer"
 		     packages="pkgTracerVerticalCarbon"
 		/>
 		<var name="verticalDONIceCell"
 		     type="real"
 		     dimensions="nDONIceLayers nCells Time"
-		     units="mmol m-3"
+		     units="mmol N m-3"
 		     description="Cell average dissolved organic nitrogen for each type and ice bio-grid layer"
 		     packages="pkgTracerVerticalDON"
 		/>
 		<var name="verticalNitrateIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol N m-3"
 		     description="Cell average nitrate for each ice bio-grid layer"
 		     packages="pkgTracerVerticalNitrate"
 		/>
 		<var name="verticalSilicateIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol Si m-3"
 		     description="Cell average silicate for each ice bio-grid layer"
 		     packages="pkgTracerVerticalSilicate"
 		/>
 		<var name="verticalAmmoniumIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol N m-3"
 		     description="Cell average ammonium for each ice bio-grid layer"
 		     packages="pkgTracerVerticalAmmonium"
 		/>
 		<var name="verticalDMSIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol S m-3"
 		     description="Cell average DMS for each ice bio-grid layer"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPpIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol S m-3"
 		     description="Cell average particulate DMSP for each ice bio-grid layer"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalDMSPdIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol S m-3"
 		     description="Cell average dissolved DMSP for each ice bio-grid layer"
 		     packages="pkgTracerVerticalDMS"
 		/>
 		<var name="verticalNonreactiveIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
+		     units="mmol N m-3"
 		     description="Cell average nonreactive nitrate for each ice bio-grid layer"
 		     packages="pkgTracerVerticalNonreactive"
 		/>
 		<var name="verticalHumicsIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
-		     units="mmol m-3"
-		     description="Cell average humic carbon for each ice bio-grid layer"
+		     units="mmol C m-3"
+		     description="Cell average refractory dissolved organic carbon for each ice bio-grid layer"
 		     packages="pkgTracerVerticalHumics"
 		/>
 		<var name="verticalParticulateIronIceCell"
 		     type="real"
 		     dimensions="nPartIronIceLayers nCells Time"
-		     units="umol m-3"
+		     units="umol Fe m-3"
 		     description="Cell average particulate iron for each type and ice bio-grid layer"
 		     packages="pkgTracerVerticalIron"
 		/>
 		<var name="verticalDissolvedIronIceCell"
 		     type="real"
 		     dimensions="nDisIronIceLayers nCells Time"
-		     units="umol m-3"
+		     units="umol Fe m-3"
 		     description="Cell average dissolved iron for each type and ice bio-grid layer"
 		     packages="pkgTracerVerticalIron"
 		/>
+		<var name="verticalBCTotalSnowCell"
+		     type="real"
+		     dimensions="TWO nCells Time"
+		     units="kg m-3"
+		     description="Cell average total BC in each snow bio-grid layer"
+		     packages="pkgTracerZAerosols"
+		/>
+		<var name="verticalDustTotalSnowCell"
+		     type="real"
+		     dimensions="TWO nCells Time"
+		     units="kg m-3"
+		     description="Cell average total dust in each snow bio-grid layer"
+		     packages="pkgTracerZAerosols"
+		/>
 		<var name="verticalAerosolsSnowCell"
 		     type="real"
 		     dimensions="nzAerosolsSnowLayers nCells Time"
@@ -3671,9 +3807,25 @@
 		     description="Cell average aerosol for each type and ice bio-grid layer"
 		     packages="pkgTracerZAerosols"
 		/>
+		<var name="verticalBCTotalIceCell"
+		     type="real"
+		     dimensions="nBioLayersP1 nCells Time"
+		     units="kg m-3"
+		     description="Cell average total black carbon in ice bio-grid layer"
+		     packages="pkgTracerZAerosols"
+		/>
+		<var name="verticalDustTotalIceCell"
+		     type="real"
+		     dimensions="nBioLayersP1 nCells Time"
+		     units="kg m-3"
+		     description="Cell average total dust in ice bio-grid layer"
+		     packages="pkgTracerZAerosols"
+		/>
 		<var name="verticalAerosolsConcCell"
 		     type="real"
 		     dimensions="nzAerosolsLayers nCells Time"
+		     units="kg m-3"
+		     description="Cell average aerosol for each type in snow and ice bio-grid layer"
 		     packages="pkgTracerZAerosols"
 		/>
 		<var name="verticalSalinityCell"
@@ -5304,6 +5456,8 @@
 		<var name="newlyFormedIce"
 		     type="integer"
 		     dimensions="nCategories nCells"
+		     units="1"
+		     description="1 if ice is newly formed in category cell, 0 otherwise"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
 		     />
 		<var name="rayleighCriteriaReal"
@@ -5313,35 +5467,43 @@
 		<var name="atmosBioFluxes"
 		     type="real"
 		     dimensions="nZBGCTracers nCells Time"
+		     units="mmol m-2 s-1"
+		     description="Atmospheric deposition flux (positive down) of vertical biology (umol for iron) and aerosols (kg for dust/BC)"
 		     />
 		<var name="atmosIceBioFluxes"
 		     type="real"
 		     dimensions="nZBGCTracers nCells Time"
+		     units="mmol m-2 s-1"
+		     description="Atmospheric deposition flux directly into the ice (positive down) of vertical biology (umol for iron) and aerosols (kg for dust/BC)"
 		     />
 		<var name="snowIceBioFluxes"
 		     type="real"
 		     dimensions="nZBGCTracers nCells Time"
+		     units="mmol m-2 s-1"
+		     description="Flux from the snow into the ice (positive down) of vertical biology (umol for iron) and aerosols (kg for dust/BC)"
 		     />
 		<var name="atmosBlackCarbonFlux"
 		     type="real"
 		     dimensions="maxBCType nCells Time"
+		     units="kg m-2 s-1"
+		     description="Atmospheric flux (positive down) of black carbon types into the sea ice"
 		     />
 		<var name="atmosDustFlux"
 		     type="real"
-		     units="kg m-2 s-1"
 		     dimensions="maxDustType nCells Time"
+		     units="kg m-2 s-1"
 		     description="Total flux of dry and wet atmospheric dust positive into the ice column"
 		     />
 		<var name="atmosDryDustFlux"
 		     type="real"
-		     units="kg m-2 s-1"
 		     dimensions="maxDustType nCells Time"
+		     units="kg m-2 s-1"
 		     description="Flux of dry atmospheric dust positive into the ice column"
 		     />
 		<var name="atmosWetDustFlux"
 		     type="real"
-		     units="kg m-2 s-1"
 		     dimensions="maxDustType nCells Time"
+		     units="kg m-2 s-1"
 		     description="Flux of wet atmospheric dust positive into the ice column"
 		     />
 		<var name="totalCarbonContentCell"
@@ -5376,9 +5538,7 @@
 		     type="real"
 		     dimensions="nCells Time"
 		     units="mmol m-2 s-1"
-		     description="Total ocean carbon flux (positive
-				  into the ocean) weighted bysea ice
-				  area fraction"
+		     description="Total ocean carbon flux (positive into the ocean) weighted bysea ice area fraction"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
 		/>
 		<var name="oceanBioFluxes"
@@ -5425,267 +5585,266 @@
 		<var name="oceanAlgaeConc"
 		     type="real"
 		     dimensions="maxAlgaeType nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean algal nitrogen concentration"
+		     units="mmol N m-3"
+		     description="Surface ocean algal nitrogen concentration"
 		/>
 		<var name="oceanDOCConc"
 		     type="real"
 		     dimensions="maxDOCType nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean dissolved organic carbon concentration"
+		     units="mmol C m-3"
+		     description="Surface ocean dissolved organic carbon concentration"
 		/>
 		<var name="oceanDICConc"
 		     type="real"
 		     dimensions="maxDICType nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean dissolved inorganic carbon concentration"
+		     units="mmol C m-3"
+		     description="Surface ocean dissolved inorganic carbon concentration"
 		/>
 		<var name="oceanDONConc"
 		     type="real"
 		     dimensions="maxDONType nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean dissolved organic nitrogen concentration"
+		     units="mmol N m-3"
+		     description="Surface ocean dissolved organic nitrogen concentration"
 		/>
 		<var name="oceanNitrateConc"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean nitrate concentration"
+		     units="mmol N m-3"
+		     description="Surface ocean nitrate concentration"
 		/>
 		<var name="oceanSilicateConc"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean silicate concentration"
+		     units="mmol Si m-3"
+		     description="Surface ocean silicate concentration"
 		/>
 		<var name="oceanAmmoniumConc"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean ammonium concentration"
+		     units="mmol N m-3"
+		     description="Surface ocean ammonium concentration"
 		/>
 		<var name="oceanDMSConc"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean dimethylsulfide concentration"
+		     units="mmol S m-3"
+		     description="Surface ocean dimethylsulfide concentration"
 		/>
 		<var name="oceanDMSPConc"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean DMSP concentration"
+		     units="mmol S m-3"
+		     description="Surface ocean DMSP concentration"
 		/>
 		<var name="oceanHumicsConc"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-3"
-                     description="Surface ocean refractory dissolved
-				  organic carbon concentration"
+		     units="mmol C m-3"
+		     description="Surface ocean refractory dissolved organic carbon concentration"
 		/>
 		<var name="oceanParticulateIronConc"
 		     type="real"
 		     dimensions="maxIronType nCells Time"
-		     units="umol m-3"
-                     description="Surface ocean particulate iron concentration"
+		     units="umol Fe m-3"
+		     description="Surface ocean particulate iron concentration"
 		/>
 		<var name="oceanDissolvedIronConc"
 		     type="real"
 		     dimensions="maxIronType nCells Time"
-		     units="umol m-3"
-                     description="Surface ocean dissolved iron concentration"
+		     units="umol Fe m-3"
+		     description="Surface ocean dissolved iron concentration"
 		/>
 		<var name="oceanZAerosolConc"
 		     type="real"
 		     dimensions="maxAerosolType nCells Time"
 		     units="kg m-3"
-                     description="Surface ocean aerosol concentration"
+		     description="Surface ocean aerosol concentration"
 		/>
 		<var name="oceanAlgaeFlux"
 		     type="real"
 		     dimensions="maxAlgaeType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Algal nitrogen flux to ocean from ice column"
+		     units="mmol N m-2 s-1"
+		     description="Algal nitrogen (diatoms, small plantkon, Phaeocystis) flux to ocean from ice column"
 		/>
 		<var name="oceanDOCFlux"
 		     type="real"
 		     dimensions="maxDOCType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved organic carbon flux to ocean from ice column"
+		     units="mmol C m-2 s-1"
+		     description="Dissolved organic carbon flux (polysaccharids, lipids) to ocean from ice column"
 		/>
 		<var name="oceanDICFlux"
 		     type="real"
 		     dimensions="maxDICType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved inorganic carbon flux to ocean from ice column"
+		     units="mmol C m-2 s-1"
+		     description="Dissolved inorganic carbon flux to ocean from ice column"
 		/>
 		<var name="oceanDONFlux"
 		     type="real"
 		     dimensions="maxDONType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved organic nitrogen flux to ocean from ice column"
+		     units="mmol N m-2 s-1"
+		     description="Dissolved organic nitrogen flux to ocean from ice column"
 		/>
 		<var name="oceanNitrateFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Nitrate flux to ocean from ice column"
+		     units="mmol N m-2 s-1"
+		     description="Nitrate flux to ocean from ice column"
 		/>
 		<var name="oceanSilicateFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Silicate flux to ocean from ice column"
+		     units="mmol Si m-2 s-1"
+		     description="Silicate flux to ocean from ice column"
 		/>
 		<var name="oceanAmmoniumFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Ammonium flux to ocean from ice column"
+		     units="mmol N m-2 s-1"
+		     description="Ammonium flux to ocean from ice column"
 		/>
 		<var name="oceanDMSFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dimethylsulfide flux to ocean from ice column"
+		     units="mmol S m-2 s-1"
+		     description="Dimethylsulfide flux to ocean from ice column"
 		/>
 		<var name="oceanDMSPpFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Particulate DMSP flux to ocean from ice column"
+		     units="mmol S m-2 s-1"
+		     description="Particulate DMSP flux to ocean from ice column"
 		/>
 		<var name="oceanDMSPdFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved DMSP flux to ocean from ice column"
+		     units="mmol S m-2 s-1"
+		     description="Dissolved DMSP flux to ocean from ice column"
 		/>
 		<var name="oceanHumicsFlux"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Refractory DOC flux to ocean from ice column"
+		     units="mmol C m-2 s-1"
+		     description="Refractory DOC flux to ocean from ice column"
 		/>
 		<var name="oceanParticulateIronFlux"
 		     type="real"
 		     dimensions="maxIronType nCells Time"
-		     units="umol m-2 s-1"
-                     description="Particulate iron flux to ocean from ice column"
+		     units="umol Fe m-2 s-1"
+		     description="Particulate iron flux to ocean from ice column"
 		/>
 		<var name="oceanDissolvedIronFlux"
 		     type="real"
 		     dimensions="maxIronType nCells Time"
-		     units="umol m-2 s-1"
-                     description="Dissolved iron flux to ocean from ice column"
+		     units="umol Fe m-2 s-1"
+		     description="Dissolved iron flux to ocean from ice column"
 		/>
 		<var name="oceanDustIronFlux"
 		     type="real"
 		     dimensions="nCells Time"
 		     units="kg m-2 s-1"
-                     description="Total iron dust flux to ocean from ice column"
+		     description="Total dust flux to ocean from ice column"
 		/>
 		<var name="oceanBlackCarbonFlux"
 		     type="real"
 		     dimensions="nCells Time"
 		     units="kg m-2 s-1"
-                     description="Total black carbon flux to ocean from ice column"
+		     description="Total black carbon flux to ocean from ice column"
 		/>
 		<var name="oceanAlgaeFluxArea"
 		     type="real"
 		     dimensions="maxAlgaeType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Algal nitrogen flux to ocean weighted by ice area fraction"
+		     units="mmol N m-2 s-1"
+		     description="Algal nitrogen flux (diatom, small plankton, Phaeocystis) to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDOCFluxArea"
 		     type="real"
 		     dimensions="maxDOCType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved organic carbon flux to ocean weighted by ice area fraction"
+		     units="mmol C m-2 s-1"
+		     description="Dissolved organic carbon flux (polysaccharid, lipid) to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDICFluxArea"
 		     type="real"
 		     dimensions="maxDICType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved inorganic carbon flux to ocean weighted by ice area fraction"
+		     units="mmol C m-2 s-1"
+		     description="Dissolved inorganic carbon flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDONFluxArea"
 		     type="real"
 		     dimensions="maxDONType nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved organic nitrogen flux to ocean weighted by ice area fraction"
+		     units="mmol N m-2 s-1"
+		     description="Dissolved organic nitrogen flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanNitrateFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Nitrate flux to ocean weighted by ice area fraction"
+		     units="mmol N m-2 s-1"
+		     description="Nitrate flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanSilicateFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Silicate flux to ocean weighted by ice area fraction"
+		     units="mmol Si m-2 s-1"
+		     description="Silicate flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanAmmoniumFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Ammonium flux to ocean weighted by ice area fraction"
+		     units="mmol N m-2 s-1"
+		     description="Ammonium flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDMSFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dimethylsulfide flux to ocean weighted by ice area fraction"
+		     units="mmol S m-2 s-1"
+		     description="Dimethylsulfide flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDMSPpFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Particulate DMSP flux to ocean weighted by ice area fraction"
+		     units="mmol S m-2 s-1"
+		     description="Particulate DMSP flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDMSPdFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Dissolved DMSP flux to ocean weighted by ice area fraction"
+		     units="mmol S m-2 s-1"
+		     description="Dissolved DMSP flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanHumicsFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mmol m-2 s-1"
-                     description="Refractory DOC flux to ocean weighted by ice area fraction"
+		     units="mmol C m-2 s-1"
+		     description="Refractory DOC flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanParticulateIronFluxArea"
 		     type="real"
 		     dimensions="maxIronType nCells Time"
-		     units="umol m-2 s-1"
-                     description="Particulate iron flux to ocean weighted by ice area fraction"
+		     units="umol Fe m-2 s-1"
+		     description="Particulate iron flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDissolvedIronFluxArea"
 		     type="real"
 		     dimensions="maxIronType nCells Time"
-		     units="umol m-2 s-1"
-                     description="Dissolved iron flux to ocean weighted by ice area fraction"
+		     units="umol Fe m-2 s-1"
+		     description="Dissolved iron flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanDustIronFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
 		     units="kg m-2 s-1"
-                     description="Total iron dust flux to ocean weighted by ice area fraction"
+		     description="Total iron flux to ocean weighted by ice area fraction"
 		/>
 		<var name="oceanBlackCarbonFluxArea"
 		     type="real"
 		     dimensions="nCells Time"
 		     units="kg m-2 s-1"
-                     description="Total black carbon flux to ocean weighted by ice area fraction"
+		     description="Total black carbon flux to ocean weighted by ice area fraction"
 		/>
 		<var name="bioTracerShortwave"
 		     type="real"
 		     dimensions="nShortwaveBio nCategories nCells"
 		     units="mg m-3"
-                     description="Chlorophyll and aerosol
+		     description="Chlorophyll and aerosol
 				  concentrations in the snow and on
 				  the shortwave vertical grid (verticalShortwaveGrid)"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
@@ -5754,22 +5913,23 @@
 		<var name="brineTopChange"
 		     type="real"
 		     dimensions="nCategories nCells"
+		     units="m"
+		     description="change in the brine penetration height near the ice top (positive up)"
+		     icepack_name="dhbr_top"
 		     packages="pkgTracerBrine"
 		/>
 		<var name="brineBottomChange"
 		     type="real"
 		     dimensions="nCategories nCells"
+		     units="m"
+		     description="change in the brine penetration height at the ice bottom (positive up)"
+		     icepack_name="dhbr_bot"
 		     packages="pkgTracerBrine"
 		/>
 		<var name="zSalinityIceDensity"
 		     type="real"
 		     dimensions="nCategories nCells Time"
 		/>
-		<var name="peakVerticalAlgalConc"
-		     type="real"
-		     dimensions="nAlgae nCells Time"
-		     packages="pkgTracerVerticalAlgae"
-		/>
 		<var name="netSpecificAlgalGrowthRate"
 		     type="real"
 		     dimensions="nCells Time"
@@ -5780,6 +5940,9 @@
 		<var name="verticalNitrogenLosses"
 		     type="real"
 		     dimensions="nBioLayersP3 nCategories nCells Time"
+		     units="mmol N m-3"
+		     description="Nitrogen losses to the sea ice system from reaction terms"
+		     icepack_name="zoo"
 		/>
 		<var name="netNitrateUptake"
 		     type="real"
@@ -5818,6 +5981,174 @@
 		     units="mmol m-2"
 		     description="Tracers are ordered: diatom nitrogen, smallPlankton nitrogen, phaeocystis nitrogen, nitrate, polysaccarid carbon, lipid carbon, dissolved inorganic carbon,ammonium, silicate, DMSPp, DMSPd, DMS, Nonreactive nitrate, Protein nitrogen, dissolved iron in umol/m2, particulate iron in umol/m2, humic carbon, black carbon1 in kg/m2, black carbon2 in kg/m2,dust1 in kg/m2, dust2 in kg/m2, dust3 in kg/m2, dust4 in kg/m2"
 		/>
+		<var name="totalVerticalDiatomIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total diatom nitrogen in the ice column"
+		/>
+		<var name="totalVerticalSmallPlanktonIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total small plankton nitrogen in the ice column"
+		/>
+		<var name="totalVerticalPhaeocystisIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total Phaeocystis nitrogen in the ice column"
+		/>
+		<var name="totalVerticalAlgaeCarbonIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol C m-2"
+		     description="Total algal carbon in the ice column"
+		/>
+		<var name="totalVerticalNitrateIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total nitrate in the ice column"
+		/>
+		<var name="totalVerticalPolysaccharidsIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol C m-2"
+		     description="Total DOC (polysaccharids) in the ice column"
+		/>
+		<var name="totalVerticalLipidsIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol C m-2"
+		     description="Total DOC (lipids) in the ice column"
+		/>
+		<var name="totalVerticalDOCLabileIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol C m-2"
+		     description="Total semi-labile dissolved organic carbon (polysaccharids+lipids+proteins) in the ice column"
+		/>
+		<var name="totalVerticalDICIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol C m-2"
+		     description="Total dissolved inorganic carbon in the ice column"
+		/>
+		<var name="totalVerticalAmmoniumIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total ammonium in the ice column"
+		/>
+		<var name="totalVerticalSilicateIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol Si m-2"
+		     description="Total silicate in the ice column"
+		/>
+		<var name="totalVerticalDMSPpIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol S m-2"
+		     description="Total particulate DMSP in the ice column"
+		/>
+		<var name="totalVerticalDMSPdIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol S m-2"
+		     description="Total dissolved DMSP in the ice column"
+		/>
+		<var name="totalVerticalDMSIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol S m-2"
+		     description="Total DMS in the ice column"
+		/>
+		<var name="totalVerticalNonreactiveIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total nonreactive nitrate in the ice column"
+		/>
+		<var name="totalVerticalProteinsIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol N m-2"
+		     description="Total DON (protein) in the ice column"
+		/>
+		<var name="totalVerticalDissolvedIronIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="umol Fe m-2"
+		     description="Total dissolved iron in the ice column"
+		/>
+		<var name="totalVerticalParticulateIronIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="umol Fe m-2"
+		     description="Total particulate iron in the ice column"
+		/>
+		<var name="totalVerticalHumicsIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="mmol C m-2"
+		     description="Total humics in the ice column"
+		/>
+		<var name="totalVerticalBC1Ice"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total BC1 in the ice column"
+		/>
+		<var name="totalVerticalBC2Ice"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total BC2 in the ice column"
+		/>
+		<var name="totalVerticalBCIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total BC1+BC2 in the ice column"
+		/>
+		<var name="totalVerticalBC1Snow"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total BC1 in the snow"
+		/>
+		<var name="totalVerticalBC2Snow"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total BC2 in the snow"
+		/>
+		<var name="totalVerticalBCSnow"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total BC1+BC2 in the snow"
+		/>
+		<var name="totalVerticalDustIce"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total dust in the ice column"
+		/>
+		<var name="totalVerticalDustSnow"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="kg m-2"
+		     description="Total dust in the snow"
+		/>
+		<var name="totalVerticalDissolvedIronSnow"
+		     type="real"
+		     dimensions="nCells Time"
+		     units="umol Fe m-2"
+		     description="Total dissolved iron in the snow"
+		/>
 		<var name="totalVerticalSalinity"
 		     type="real"
 		     dimensions="nCells Time"
@@ -5826,8 +6157,9 @@
 		<var name="totalChlorophyll"
 		     type="real"
 		     dimensions="nCells Time"
-		     units="mg m-2"
+		     units="mg Chla m-2"
 		     description="Ice integrated chlorophyll concentration per cell"
+		     icepack_name="totalChla"
 		/>
 		<var name="totalNitrate"
 		     type="real"
@@ -5847,38 +6179,61 @@
 		<var name="biologyGrid"
 		     type="real"
 		     dimensions="nBioLayersP2"
+		     units="1"
+		     description="nondimensional biology vertical grid points"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
+		     icepack_name="bgrid"
 		/>
 		<var name="interfaceBiologyGrid"
 		     type="real"
 		     dimensions="nBioLayersP1"
+		     units="1"
+		     description="nondimensional biology vertical interface points"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
+		     icepack_name="igrid"
 		/>
 		<var name="verticalGrid"
 		     type="real"
 		     dimensions="nIceLayersP1"
+		     units="1"
+		     description="MPAS vertical grid"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
+		     icepack_name="cgrid"
 		/>
 		<var name="interfaceGrid"
 		     type="real"
 		     dimensions="nIceLayersP1"
+		     units="1"
+		     description="MPAS vertical interface points"
+		     icepack_name="icgrid"
 		/>
 		<var name="verticalShortwaveGrid"
 		     type="real"
 		     dimensions="nIceLayersP1"
+		     units="1"
+		     description="Vertical ice grid for the shortwave radiation scheme"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
+		     icepack_name="swgrid"
 		/>
 		<var name="carbonToNitrogenRatioAlgae"
 		     type="real"
 		     dimensions="maxAlgaeType"
+		     units="mol C mol N-1"
+		     description="Carbon to nitrogen ratios for the algal types"
+		     icepack_name="CToN"
 		/>
 		<var name="carbonToNitrogenRatioDON"
 		     type="real"
 		     dimensions="maxDONType"
+		     units="mol C mol N-1"
+		     description="Carbon to nitrogen ratio for DON"
+		     icepack_name="CToN_DON"
 		/>
 		<var name="DOCPoolFractions"
 		     type="real"
 		     dimensions="maxDOCType"
+		     units="1"
+		     description="Fraction of DOC in polysaccharids, lipids and proteins"
 		/>
 	</var_struct>
 
@@ -5887,6 +6242,7 @@
 		<var name="bgridPorosityIceCell"
 		     type="real"
 		     dimensions="nBioLayersP1 nCells Time"
+		     units="1"
 		     description="Cell ice volume weighted average brine porosity on ice interface biological grid"
 		     packages="pkgColumnBiogeochemistry;pkgColumnPackage"
 		/>
diff --git a/components/mpas-seaice/src/shared/mpas_seaice_icepack.F b/components/mpas-seaice/src/shared/mpas_seaice_icepack.F
index 36899b4b10d7..9dc6f06ec342 100644
--- a/components/mpas-seaice/src/shared/mpas_seaice_icepack.F
+++ b/components/mpas-seaice/src/shared/mpas_seaice_icepack.F
@@ -3923,7 +3923,11 @@ subroutine column_biogeochemistry(domain)
           abortFlag = icepack_warnings_aborted()
 
           call get_cice_tracer_array_category(block, ciceTracerObject, &
-               tracerArrayCategory, iCell, setGetPhysicsTracers, setGetBGCTracers)
+             tracerArrayCategory, iCell, setGetPhysicsTracers, setGetBGCTracers)
+
+          if (config_use_vertical_tracers) &
+             call define_total_biogeochemistry_array_cell(block, ciceTracerObject, totalVerticalBiologyIce(:,iCell), &
+             totalVerticalBiologySnow(:,iCell), iCell)
 
           if (checkCarbon) then
             call seaice_total_carbon_content_category(block,totalCarbonCatFinal(:,iCell),iceAreaCategory(1,:,:),iceVolumeCategory(1,:,:),iCell)
@@ -8648,10 +8652,13 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
          nzAerosols
 
     type(MPAS_pool_type), pointer :: &
-         tracers_aggregate
+         tracers_aggregate, &
+         biogeochemistry
 
     real(kind=RKIND), dimension(:), pointer :: &
-         brineFractionCell
+         brineFractionCell, &
+         carbonToNitrogenRatioAlgae, &
+         carbonToNitrogenRatioDON
 
     real(kind=RKIND), dimension(:,:), pointer :: &
          skeletalAlgaeConcCell, &
@@ -8699,14 +8706,21 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
          verticalParticulateIronIceCell, &
          verticalDissolvedIronIceCell, &
          verticalAerosolsIceCell, &
-         verticalAerosolsSnowCell
+         verticalAerosolsSnowCell, &
+         verticalDOCLabileIceCell, &
+         verticalAlgaeTotalCarbonIceCell, &
+         verticalBCTotalIceCell, &
+         verticalDustTotalIceCell, &
+         verticalBCTotalSnowCell, &
+         verticalDustTotalSnowCell
 
     integer :: &
          iBioTracers, &
          iBioCount, &
          iLayers, &
          iIceCount, &
-         iSnowCount
+         iSnowCount, &
+         nAeroType
 
     call MPAS_pool_get_config(block % configs, "config_use_skeletal_biochemistry", config_use_skeletal_biochemistry)
     call MPAS_pool_get_config(block % configs, "config_use_vertical_biochemistry", config_use_vertical_biochemistry)
@@ -8736,6 +8750,7 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
     call MPAS_pool_get_dimension(block % dimensions, "nDissolvedIron", nDissolvedIron)
 
     call MPAS_pool_get_subpool(block % structs, "tracers_aggregate", tracers_aggregate)
+    call MPAS_pool_get_subpool(block % structs, "biogeochemistry", biogeochemistry)
 
     call MPAS_pool_get_array(tracers_aggregate, "skeletalAlgaeConcCell", skeletalAlgaeConcCell)
     call MPAS_pool_get_array(tracers_aggregate, "skeletalDOCConcCell", skeletalDOCConcCell)
@@ -8784,6 +8799,14 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
     call MPAS_pool_get_array(tracers_aggregate, "verticalAerosolsSnowCell", verticalAerosolsSnowCell)
     call MPAS_pool_get_array(tracers_aggregate, "verticalSalinityCell", verticalSalinityCell)
     call MPAS_pool_get_array(tracers_aggregate, "brineFractionCell", brineFractionCell)
+    call MPAS_pool_get_array(tracers_aggregate, "verticalAlgaeTotalCarbonIceCell", verticalAlgaeTotalCarbonIceCell)
+    call MPAS_pool_get_array(tracers_aggregate, "verticalDOCLabileIceCell", verticalDOCLabileIceCell)
+    call MPAS_pool_get_array(tracers_aggregate, "verticalBCTotalIceCell", verticalBCTotalIceCell)
+    call MPAS_pool_get_array(tracers_aggregate, "verticalDustTotalIceCell", verticalDustTotalIceCell)
+    call MPAS_pool_get_array(tracers_aggregate, "verticalDustTotalSnowCell", verticalDustTotalSnowCell)
+    call MPAS_pool_get_array(tracers_aggregate, "verticalBCTotalSnowCell", verticalBCTotalSnowCell)
+    call MPAS_pool_get_array(biogeochemistry, "carbonToNitrogenRatioAlgae", carbonToNitrogenRatioAlgae)
+    call MPAS_pool_get_array(biogeochemistry, "carbonToNitrogenRatioDON", carbonToNitrogenRatioDON)
 
     ! biogeochemistry
     ! brine height fraction
@@ -8866,12 +8889,15 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
           ! algal nitrogen
           do iBioTracers = 1, nAlgae
              iIceCount = (iBioTracers-1)*nBioLayersP1
+             verticalAlgaeTotalCarbonIceCell(:,iCell) = 0.0_RKIND
 
              do iLayers = 1,nBioLayersP1
                 iBiocount = iBiocount + 1
                 verticalAlgaeConcCell(iBioCount,iCell) = &
                      tracerArrayCell(tracerObject % index_algaeConc(iBioTracers)+iLayers-1)
                 verticalAlgaeIceCell(iLayers+iIceCount,iCell) = verticalAlgaeConcCell(iBioCount,iCell)
+                verticalAlgaeTotalCarbonIceCell(iLayers,iCell) = verticalAlgaeTotalCarbonIceCell(iLayers,iCell) + &
+                   carbonToNitrogenRatioAlgae(iBioTracers) * verticalAlgaeConcCell(iBioCount,iCell)
              enddo
              do iLayers = nBioLayersP1+1,nBioLayersP3
                 iBiocount = iBiocount + 1
@@ -8894,6 +8920,7 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
 
        if (config_use_carbon) then
           iBioCount = 0
+          verticalDOCLabileIceCell(:,iCell) = 0.0_RKIND
 
           ! DOC
           do iBioTracers = 1, nDOC
@@ -8904,6 +8931,8 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
                 verticalDOCConcCell(iBioCount,iCell) = &
                      tracerArrayCell(tracerObject % index_DOCConc(iBioTracers) + iLayers-1)
                 verticalDOCIceCell(iLayers+iIceCount,iCell) = verticalDOCConcCell(iBioCount,iCell)
+                verticalDOCLabileIceCell(iLayers, iCell) = verticalDOCLabileIceCell(iLayers, iCell) + &
+                   verticalDOCConcCell(iBioCount,iCell)
              enddo
              do iLayers = nBioLayersP1+1,nBioLayersP3
                 iBioCount = iBioCount + 1
@@ -8942,6 +8971,8 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
                 verticalDONConcCell(iBioCount,iCell) = &
                      tracerArrayCell(tracerObject % index_DONConc(iBioTracers) + iLayers-1)
                 verticalDONIceCell(iLayers+iIceCount,iCell) = verticalDONConcCell(iBioCount,iCell)
+                verticalDOCLabileIceCell(iLayers, iCell) = verticalDOCLabileIceCell(iLayers, iCell) + &
+                   verticalDONConcCell(iBioCount,iCell) * carbonToNitrogenRatioDON(iBioTracers)
              enddo
              do iLayers = nBioLayersP1+1,nBioLayersP3
                 iBioCount = iBioCount + 1
@@ -9058,15 +9089,25 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
        ! black carbon and dust aerosols
        if (config_use_zaerosols) then
           iBioCount = 0
+          verticalDustTotalIceCell(:,iCell) = 0.0_RKIND
+          verticalBCTotalIceCell(:,iCell) = 0.0_RKIND
+          verticalDustTotalSnowCell(:,iCell) = 0.0_RKIND
+          verticalBCTotalSnowCell(:,iCell) = 0.0_RKIND
+
           do iBioTracers = 1, nzAerosols
              iIceCount = (iBioTracers-1)*nBioLayersP1
              iSnowCount = (iBioTracers-1)*2
-
+             nAeroType = 1
+             if (iBioTracers .gt. 2) nAeroType = 0
              do iLayers = 1,nBioLayersP1
                 iBioCount = iBioCount + 1
                 verticalAerosolsConcCell(iBioCount,iCell) = &
                      tracerArrayCell(tracerObject % index_verticalAerosolsConc(iBioTracers)+iLayers-1)
                 verticalAerosolsIceCell(iLayers+iIceCount,iCell) = verticalAerosolsConcCell(iBioCount,iCell)
+                verticalBCTotalIceCell(iLayers, iCell) = verticalBCTotalIceCell(iLayers, iCell) + nAeroType * &
+                      verticalAerosolsConcCell(iBioCount,iCell)
+                verticalDustTotalIceCell(iLayers, iCell) = verticalDustTotalIceCell(iLayers, iCell) + (1-nAeroType) * &
+                      verticalAerosolsConcCell(iBioCount,iCell)
              enddo
              do iLayers = nBioLayersP1+1,nBioLayersP3
                 iBioCount = iBioCount + 1
@@ -9074,6 +9115,10 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
                      tracerArrayCell(tracerObject % index_verticalAerosolsConc(iBioTracers)+iLayers-1)
                 verticalAerosolsSnowCell(iLayers-nBioLayersP1+iSnowCount,iCell) = &
                      verticalAerosolsConcCell(iBioCount,iCell)
+                verticalBCTotalSnowCell(iLayers-nBioLayersP1, iCell) = verticalBCTotalSnowCell(iLayers-nBioLayersP1, iCell) + &
+                      nAeroType * verticalAerosolsConcCell(iBioCount,iCell)
+                verticalDustTotalSnowCell(iLayers-nBioLayersP1,iCell) = verticalDustTotalSnowCell(iLayers-nBioLayersP1,iCell) + &
+                      (1-nAeroType) * verticalAerosolsConcCell(iBioCount,iCell)
              enddo
           enddo
        endif
@@ -9088,6 +9133,274 @@ subroutine get_cice_biogeochemistry_tracer_array_cell(block, tracerObject, trace
 
   end subroutine get_cice_biogeochemistry_tracer_array_cell
 
+!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
+!
+!  define_total_biogeochemistry_array_cell
+!
+!> \brief
+!> \author Nicole Jeffery, LANL
+!> \date 16rd September 2024
+!> \details
+!>
+!   Defines total ice and snow arrays for biogeochemistry and aerosols diagnostics
+!
+!-----------------------------------------------------------------------
+
+  subroutine define_total_biogeochemistry_array_cell(block, tracerObject, totalVerticalBioIceCell, &
+         totalVerticalBioSnowCell, iCell)
+
+    type(block_type), intent(in) :: &
+         block
+
+    type(ciceTracerObjectType), intent(in) :: &
+         tracerObject
+
+    real(kind=RKIND), dimension(:), intent(in):: &
+         totalVerticalBioIceCell, &
+         totalVerticalBioSnowCell
+
+    integer, intent(in) :: &
+         iCell
+
+    logical, pointer :: &
+         config_use_vertical_biochemistry, &
+         config_use_nitrate, &
+         config_use_carbon, &
+         config_use_ammonium, &
+         config_use_silicate, &
+         config_use_DMS, &
+         config_use_nonreactive, &
+         config_use_humics, &
+         config_use_DON, &
+         config_use_iron, &
+         config_use_zaerosols
+
+    integer, pointer :: &
+         nAlgae, &
+         nDOC, &
+         nDIC, &
+         nDON, &
+         nParticulateIron, &
+         nDissolvedIron, &
+         nzAerosols
+
+    type(MPAS_pool_type), pointer :: &
+         biogeochemistry
+
+    real(kind=RKIND), dimension(:), pointer :: &
+         totalVerticalDiatomIce, &
+         totalVerticalSmallPlanktonIce, &
+         totalVerticalPhaeocystisIce, &
+         totalVerticalNitrateIce, &
+         totalVerticalPolysaccharidsIce, &
+         totalVerticalLipidsIce, &
+         totalVerticalDICIce, &
+         totalVerticalAmmoniumIce, &
+         totalVerticalSilicateIce, &
+         totalVerticalDMSPpIce, &
+         totalVerticalDMSPdIce, &
+         totalVerticalDMSIce, &
+         totalVerticalNonreactiveIce, &
+         totalVerticalHumicsIce, &
+         totalVerticalProteinsIce, &
+         totalVerticalDissolvedIronIce, &
+         totalVerticalParticulateIronIce, &
+         totalVerticalDustIce, &
+         totalVerticalDustSnow, &
+         totalVerticalBC1Ice, &
+         totalVerticalBC1Snow, &
+         totalVerticalBC2Ice, &
+         totalVerticalBC2Snow, &
+         totalVerticalBCIce, &
+         totalVerticalBCSnow, &
+         totalVerticalDissolvedIronSnow, &
+         totalVerticalAlgaeCarbonIce, &
+         totalVerticalDOCLabileIce, &
+         carbonToNitrogenRatioAlgae, &
+         carbonToNitrogenRatioDON
+
+    call MPAS_pool_get_config(block % configs, "config_use_vertical_biochemistry", config_use_vertical_biochemistry)
+    call MPAS_pool_get_config(block % configs, "config_use_nitrate", config_use_nitrate)
+    call MPAS_pool_get_config(block % configs, "config_use_carbon", config_use_carbon)
+    call MPAS_pool_get_config(block % configs, "config_use_ammonium",config_use_ammonium)
+    call MPAS_pool_get_config(block % configs, "config_use_silicate",config_use_silicate)
+    call MPAS_pool_get_config(block % configs, "config_use_DMS",config_use_DMS)
+    call MPAS_pool_get_config(block % configs, "config_use_nonreactive",config_use_nonreactive)
+    call MPAS_pool_get_config(block % configs, "config_use_humics",config_use_humics)
+    call MPAS_pool_get_config(block % configs, "config_use_DON",config_use_DON)
+    call MPAS_pool_get_config(block % configs, "config_use_iron",config_use_iron)
+    call MPAS_pool_get_config(block % configs, "config_use_zaerosols",config_use_zaerosols)
+
+    call MPAS_pool_get_dimension(block % dimensions, "nzAerosols", nzAerosols)
+    call MPAS_pool_get_dimension(block % dimensions, "nAlgae", nAlgae)
+    call MPAS_pool_get_dimension(block % dimensions, "nDOC", nDOC)
+    call MPAS_pool_get_dimension(block % dimensions, "nDIC", nDIC)
+    call MPAS_pool_get_dimension(block % dimensions, "nDON", nDON)
+    call MPAS_pool_get_dimension(block % dimensions, "nParticulateIron", nParticulateIron)
+    call MPAS_pool_get_dimension(block % dimensions, "nDissolvedIron", nDissolvedIron)
+
+    call MPAS_pool_get_subpool(block % structs, "biogeochemistry", biogeochemistry)
+
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDiatomIce", totalVerticalDiatomIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalSmallPlanktonIce", totalVerticalSmallPlanktonIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalPhaeocystisIce", totalVerticalPhaeocystisIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalPolysaccharidsIce", totalVerticalPolysaccharidsIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalLipidsIce", totalVerticalLipidsIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDICIce", totalVerticalDICIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalProteinsIce", totalVerticalProteinsIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalNitrateIce", totalVerticalNitrateIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalSilicateIce", totalVerticalSilicateIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalAmmoniumIce", totalVerticalAmmoniumIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDMSIce", totalVerticalDMSIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDMSPpIce", totalVerticalDMSPpIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDMSPdIce", totalVerticalDMSPdIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalNonreactiveIce", totalVerticalNonreactiveIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalHumicsIce", totalVerticalHumicsIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalParticulateIronIce", totalVerticalParticulateIronIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDissolvedIronIce", totalVerticalDissolvedIronIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDissolvedIronSnow", totalVerticalDissolvedIronSnow)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDustIce", totalVerticalDustIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDustSnow", totalVerticalDustSnow)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalBC1Ice", totalVerticalBC1Ice)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalBC1Snow", totalVerticalBC1Snow)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalBC2Ice", totalVerticalBC2Ice)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalBC2Snow", totalVerticalBC2Snow)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalBCIce", totalVerticalBCIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalBCSnow", totalVerticalBCSnow)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalAlgaeCarbonIce", totalVerticalAlgaeCarbonIce)
+    call MPAS_pool_get_array(biogeochemistry, "totalVerticalDOCLabileIce", totalVerticalDOCLabileIce)
+    call MPAS_pool_get_array(biogeochemistry, "carbonToNitrogenRatioAlgae", carbonToNitrogenRatioAlgae)
+    call MPAS_pool_get_array(biogeochemistry, "carbonToNitrogenRatioDON", carbonToNitrogenRatioDON)
+
+    ! biogeochemistry
+
+    if (config_use_vertical_biochemistry) then
+
+       ! algal nitrogen
+       totalVerticalDiatomIce(iCell) = totalVerticalBioIceCell(tracerObject % index_algaeConcLayer(1))
+       totalVerticalAlgaeCarbonIce(iCell) = totalVerticalDiatomIce(iCell) * carbonToNitrogenRatioAlgae(1)
+
+       SELECT CASE (nAlgae)
+       CASE (2)
+          totalVerticalSmallPlanktonIce(iCell) = totalVerticalBioIceCell(tracerObject % index_algaeConcLayer(2))
+          totalVerticalAlgaeCarbonIce(iCell) = totalVerticalAlgaeCarbonIce(iCell) + &
+             totalVerticalSmallPlanktonIce(iCell) * carbonToNitrogenRatioAlgae(2)
+       CASE (3)
+          totalVerticalSmallPlanktonIce(iCell) = totalVerticalBioIceCell(tracerObject % index_algaeConcLayer(2))
+          totalVerticalPhaeocystisIce(iCell) = totalVerticalBioIceCell(tracerObject % index_algaeConcLayer(3))
+          totalVerticalAlgaeCarbonIce(iCell) = totalVerticalAlgaeCarbonIce(iCell) + &
+             totalVerticalSmallPlanktonIce(iCell) * carbonToNitrogenRatioAlgae(2) + &
+             totalVerticalPhaeocystisIce(iCell) * carbonToNitrogenRatioAlgae(3)
+       END SELECT
+    end if
+
+    ! nitrate
+    if (config_use_nitrate) &
+      totalVerticalNitrateIce(iCell) = totalVerticalBioIceCell(tracerObject % index_nitrateConcLayer)
+
+    if (config_use_carbon) then
+
+       ! DOC
+       totalVerticalDOCLabileIce(iCell) = 0.0_RKIND
+       SELECT CASE (nDOC)
+       CASE (1)
+          totalVerticalPolysaccharidsIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DOCConcLayer(1))
+          totalVerticalDOCLabileIce(iCell) = totalVerticalPolysaccharidsIce(iCell)
+       CASE (2)
+          totalVerticalPolysaccharidsIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DOCConcLayer(1))
+          totalVerticalLipidsIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DOCConcLayer(1))
+          totalVerticalDOCLabileIce(iCell) = totalVerticalPolysaccharidsIce(iCell) + totalVerticalLipidsIce(iCell)
+       END SELECT
+
+       ! DIC
+       SELECT CASE (nDIC)
+       CASE (1)
+          totalVerticalDICIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DICConcLayer(1))
+       END SELECT
+    end if
+
+    if (config_use_DON) then
+
+       ! DON
+       SELECT CASE (nDON)
+       CASE (1)
+          totalVerticalProteinsIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DONConcLayer(1))
+          totalVerticalDOCLabileIce(iCell) = totalVerticalDOCLabileIce(iCell) + &
+             totalVerticalProteinsIce(iCell)* carbonToNitrogenRatioDON(1)
+       END SELECT
+    end if
+
+    ! ammonium
+    if (config_use_ammonium) &
+      totalVerticalAmmoniumIce(iCell) = totalVerticalBioIceCell(tracerObject % index_ammoniumConcLayer)
+
+    ! silicate
+    if (config_use_silicate) &
+      totalVerticalSilicateIce(iCell) = totalVerticalBioIceCell(tracerObject % index_silicateConcLayer)
+
+    ! DMS, DMSPp, DMSPd
+    if (config_use_DMS) then
+       totalVerticalDMSIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DMSConcLayer)
+       totalVerticalDMSPpIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DMSPpConcLayer)
+       totalVerticalDMSPdIce(iCell) = totalVerticalBioIceCell(tracerObject % index_DMSPdConcLayer)
+    end if
+
+    ! nonreactive
+    if (config_use_nonreactive) &
+      totalVerticalNonreactiveIce(iCell) = totalVerticalBioIceCell(tracerObject % index_nonreactiveConcLayer)
+
+    ! humic material (refractory DOC)
+    if (config_use_humics) then
+      totalVerticalHumicsIce(iCell) = totalVerticalBioIceCell(tracerObject % index_humicsConcLayer)
+    end if
+
+    if (config_use_iron) then
+
+      ! ParticulateIron
+      SELECT CASE (nParticulateIron)
+      CASE (1)
+         totalVerticalParticulateIronIce(iCell) = totalVerticalBioIceCell(tracerObject % index_particulateIronConcLayer(1))
+      END SELECT
+
+      ! DissolvedIron
+      SELECT CASE (nDissolvedIron)
+      CASE (1)
+         totalVerticalDissolvedIronIce(iCell) = totalVerticalBioIceCell(tracerObject % index_dissolvedIronConcLayer(1))
+         totalVerticalDissolvedIronSnow(iCell) = totalVerticalBioSnowCell(tracerObject % index_dissolvedIronConcLayer(1))
+      END SELECT
+
+    end if
+
+    ! black carbon and dust aerosols
+    if (config_use_zaerosols) then
+
+      SELECT CASE (nzAerosols)
+      CASE (1)
+         totalVerticalBC1Ice(iCell) = totalVerticalBioIceCell(tracerObject % index_verticalAerosolsConcLayer(1))
+         totalVerticalBC1Snow(iCell) = totalVerticalBioSnowCell(tracerObject % index_verticalAerosolsConcLayer(1))
+         totalVerticalBCIce(iCell) = totalVerticalBC1Ice(iCell)
+         totalVerticalBCSnow(iCell) = totalVerticalBC1Snow(iCell)
+      CASE (2)
+         totalVerticalBC1Ice(iCell) = totalVerticalBioIceCell(tracerObject % index_verticalAerosolsConcLayer(1))
+         totalVerticalBC1Snow(iCell) = totalVerticalBioSnowCell(tracerObject % index_verticalAerosolsConcLayer(1))
+         totalVerticalBC2Ice(iCell) = totalVerticalBioIceCell(tracerObject % index_verticalAerosolsConcLayer(2))
+         totalVerticalBC2Snow(iCell) = totalVerticalBioSnowCell(tracerObject % index_verticalAerosolsConcLayer(2))
+         totalVerticalBCIce(iCell) = totalVerticalBC1Ice(iCell) + totalVerticalBC2Ice(iCell)
+         totalVerticalBCSnow(iCell) = totalVerticalBC1Snow(iCell) + totalVerticalBC2Snow(iCell)
+      CASE (3)
+         totalVerticalBC1Ice(iCell) = totalVerticalBioIceCell(tracerObject % index_verticalAerosolsConcLayer(1))
+         totalVerticalBC1Snow(iCell) = totalVerticalBioSnowCell(tracerObject % index_verticalAerosolsConcLayer(1))
+         totalVerticalBC2Ice(iCell) = totalVerticalBioIceCell(tracerObject % index_verticalAerosolsConcLayer(2))
+         totalVerticalBC2Snow(iCell) = totalVerticalBioSnowCell(tracerObject % index_verticalAerosolsConcLayer(2))
+         totalVerticalDustIce(iCell) = totalVerticalBioIceCell(tracerObject % index_verticalAerosolsConcLayer(3))
+         totalVerticalDustSnow(iCell) = totalVerticalBioSnowCell(tracerObject % index_verticalAerosolsConcLayer(3))
+         totalVerticalBCIce(iCell) = totalVerticalBC1Ice(iCell) + totalVerticalBC2Ice(iCell)
+         totalVerticalBCSnow(iCell) = totalVerticalBC1Snow(iCell) + totalVerticalBC2Snow(iCell)
+      END SELECT
+   end if
+
+  end subroutine define_total_biogeochemistry_array_cell
+
 !-----------------------------------------------------------------------
 ! Init CICE parameters
 !-----------------------------------------------------------------------
@@ -14689,7 +15002,8 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
 
     type(MPAS_pool_type), pointer :: &
          biogeochemistryPool, &
-         diagnostics_biogeochemistryPool
+         diagnostics_biogeochemistryPool, &
+         tracers_aggregatePool
 
     ! biogeochemistry
     real(kind=RKIND), dimension(:), pointer :: &
@@ -14699,7 +15013,35 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
 !         zSalinityFlux, &        !echmod deprecate
 !         zSalinityGDFlux, &      !echmod deprecate
          totalChlorophyll, &
-         totalCarbonContentCell
+         totalCarbonContentCell, &
+         totalVerticalDiatomIce, &
+         totalVerticalSmallPlanktonIce, &
+         totalVerticalPhaeocystisIce, &
+         totalVerticalNitrateIce, &
+         totalVerticalPolysaccharidsIce, &
+         totalVerticalLipidsIce, &
+         totalVerticalDICIce, &
+         totalVerticalAmmoniumIce, &
+         totalVerticalSilicateIce, &
+         totalVerticalDMSPpIce, &
+         totalVerticalDMSPdIce, &
+         totalVerticalDMSIce, &
+         totalVerticalNonreactiveIce, &
+         totalVerticalHumicsIce, &
+         totalVerticalProteinsIce, &
+         totalVerticalDissolvedIronIce, &
+         totalVerticalParticulateIronIce, &
+         totalVerticalDustIce, &
+         totalVerticalDustSnow, &
+         totalVerticalBC1Ice, &
+         totalVerticalBC1Snow, &
+         totalVerticalBC2Ice, &
+         totalVerticalBC2Snow, &
+         totalVerticalBCIce, &
+         totalVerticalBCSnow, &
+         totalVerticalAlgaeCarbonIce, &
+         totalVerticalDOCLabileIce, &
+         totalVerticalDissolvedIronSnow
 
     real(kind=RKIND), dimension(:,:), pointer :: &
          oceanBioFluxes, &
@@ -14709,7 +15051,12 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
          totalVerticalBiologySnow, &
          bgridPorosityIceCell, &
          bgridSalinityIceCell, &
-         bgridTemperatureIceCell
+         bgridTemperatureIceCell, &
+         verticalBCTotalIceCell, &
+         verticalDustTotalIceCell, &
+         verticalBCTotalSnowCell, &
+         verticalDustTotalSnowCell
+
 
     real(kind=RKIND), dimension(:,:,:), pointer :: &
          bioTracerShortwave
@@ -14739,6 +15086,7 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
 
              call MPAS_pool_get_subpool(block % structs, "biogeochemistry", biogeochemistryPool)
              call MPAS_pool_get_subpool(block % structs, "diagnostics_biogeochemistry", diagnostics_biogeochemistryPool)
+             call MPAS_pool_get_subpool(block % structs, "tracers_aggregate", tracers_aggregatePool)
 
              if (config_use_vertical_tracers) then
                 call MPAS_pool_get_array(biogeochemistryPool, "primaryProduction", primaryProduction)
@@ -14747,13 +15095,12 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
                 call MPAS_pool_get_array(diagnostics_biogeochemistryPool, "bgridSalinityIceCell", bgridSalinityIceCell)
                 call MPAS_pool_get_array(diagnostics_biogeochemistryPool, "bgridPorosityIceCell", bgridPorosityIceCell)
                 call MPAS_pool_get_array(diagnostics_biogeochemistryPool, "bgridTemperatureIceCell", bgridTemperatureIceCell)
-
-                primaryProduction          = 0.0_RKIND
-                totalChlorophyll           = 0.0_RKIND
-                netSpecificAlgalGrowthRate = 0.0_RKIND
-                bgridSalinityIceCell       = 0.0_RKIND
-                bgridPorosityIceCell       = 0.0_RKIND
-                bgridTemperatureIceCell    = 0.0_RKIND
+                primaryProduction             = 0.0_RKIND
+                totalChlorophyll              = 0.0_RKIND
+                netSpecificAlgalGrowthRate    = 0.0_RKIND
+                bgridSalinityIceCell          = 0.0_RKIND
+                bgridPorosityIceCell          = 0.0_RKIND
+                bgridTemperatureIceCell       = 0.0_RKIND
 
              end if
 
@@ -14764,6 +15111,38 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
              call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBiologyIce", totalVerticalBiologyIce)
              call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBiologySnow", totalVerticalBiologySnow)
              call MPAS_pool_get_array(biogeochemistryPool, "totalCarbonContentCell", totalCarbonContentCell)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDiatomIce", totalVerticalDiatomIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalSmallPlanktonIce", totalVerticalSmallPlanktonIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalPhaeocystisIce", totalVerticalPhaeocystisIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalPolysaccharidsIce", totalVerticalPolysaccharidsIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalLipidsIce", totalVerticalLipidsIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDICIce", totalVerticalDICIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalProteinsIce", totalVerticalProteinsIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalNitrateIce", totalVerticalNitrateIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalSilicateIce", totalVerticalSilicateIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalAmmoniumIce", totalVerticalAmmoniumIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDMSIce", totalVerticalDMSIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDMSPpIce", totalVerticalDMSPpIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDMSPdIce", totalVerticalDMSPdIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalNonreactiveIce", totalVerticalNonreactiveIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalHumicsIce", totalVerticalHumicsIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalParticulateIronIce", totalVerticalParticulateIronIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDissolvedIronIce", totalVerticalDissolvedIronIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDissolvedIronSnow", totalVerticalDissolvedIronSnow)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDustIce", totalVerticalDustIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDustSnow", totalVerticalDustSnow)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBC1Ice", totalVerticalBC1Ice)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBC1Snow", totalVerticalBC1Snow)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBC2Ice", totalVerticalBC2Ice)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBC2Snow", totalVerticalBC2Snow)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBCIce", totalVerticalBCIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalBCSnow", totalVerticalBCSnow)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalAlgaeCarbonIce", totalVerticalAlgaeCarbonIce)
+             call MPAS_pool_get_array(biogeochemistryPool, "totalVerticalDOCLabileIce", totalVerticalDOCLabileIce)
+             call MPAS_pool_get_array(tracers_aggregatePool, "verticalBCTotalIceCell", verticalBCTotalIceCell)
+             call MPAS_pool_get_array(tracers_aggregatePool, "verticalDustTotalIceCell", verticalDustTotalIceCell)
+             call MPAS_pool_get_array(tracers_aggregatePool, "verticalDustTotalSnowCell", verticalDustTotalSnowCell)
+             call MPAS_pool_get_array(tracers_aggregatePool, "verticalBCTotalSnowCell", verticalBCTotalSnowCell)
 
 
              netBrineHeight             = 0.0_RKIND
@@ -14773,6 +15152,38 @@ subroutine seaice_icepack_reinitialize_diagnostics_bgc(domain)
              totalVerticalBiologyIce    = 0.0_RKIND
              totalVerticalBiologySnow   = 0.0_RKIND
              totalCarbonContentCell     = 0.0_RKIND
+             totalVerticalDiatomIce       = 0.0_RKIND
+             totalVerticalSmallPlanktonIce = 0.0_RKIND
+             totalVerticalPhaeocystisIce   = 0.0_RKIND
+             totalVerticalPolysaccharidsIce = 0.0_RKIND
+             totalVerticalLipidsIce        = 0.0_RKIND
+             totalVerticalDICIce           = 0.0_RKIND
+             totalVerticalProteinsIce      = 0.0_RKIND
+             totalVerticalNitrateIce       = 0.0_RKIND
+             totalVerticalSilicateIce      = 0.0_RKIND
+             totalVerticalAmmoniumIce      = 0.0_RKIND
+             totalVerticalDMSIce           = 0.0_RKIND
+             totalVerticalDMSPpIce         = 0.0_RKIND
+             totalVerticalDMSPdIce         = 0.0_RKIND
+             totalVerticalNonreactiveIce   = 0.0_RKIND
+             totalVerticalHumicsIce        = 0.0_RKIND
+             totalVerticalParticulateIronIce = 0.0_RKIND
+             totalVerticalDissolvedIronIce    = 0.0_RKIND
+             totalVerticalDissolvedIronSnow   = 0.0_RKIND
+             totalVerticalBC1Ice              = 0.0_RKIND
+             totalVerticalBC2Ice           = 0.0_RKIND
+             totalVerticalDustIce          = 0.0_RKIND
+             totalVerticalBC1Snow          = 0.0_RKIND
+             totalVerticalBC2Snow          = 0.0_RKIND
+             totalVerticalBCSnow           = 0.0_RKIND
+             totalVerticalBCIce            = 0.0_RKIND
+             totalVerticalDustSnow         = 0.0_RKIND
+             totalVerticalAlgaeCarbonIce   = 0.0_RKIND
+             totalVerticalDOCLabileIce     = 0.0_RKIND
+             verticalBCTotalIceCell        = 0.0_RKIND
+             verticalDustTotalIceCell      = 0.0_RKIND
+             verticalBCTotalSnowCell       = 0.0_RKIND
+             verticalDustTotalSnowCell     = 0.0_RKIND
 
           endif