diff --git a/full/ice_ocean_flux_exchange.F90 b/full/ice_ocean_flux_exchange.F90
index 4131c3ba..881a9f47 100644
--- a/full/ice_ocean_flux_exchange.F90
+++ b/full/ice_ocean_flux_exchange.F90
@@ -136,6 +136,7 @@ subroutine ice_ocean_flux_exchange_init(Time, Ice, Ocean, Ocean_state, ice_ocean
     allocate( ice_ocean_boundary%lprec    (is:ie,js:je) ) ;         ice_ocean_boundary%lprec = 0.0
     allocate( ice_ocean_boundary%fprec    (is:ie,js:je) ) ;         ice_ocean_boundary%fprec = 0.0
     allocate( ice_ocean_boundary%runoff   (is:ie,js:je) ) ;         ice_ocean_boundary%runoff = 0.0
+    allocate( ice_ocean_boundary%runoff_carbon (is:ie,js:je) ) ;    ice_ocean_boundary%runoff_carbon = 0.0
     allocate( ice_ocean_boundary%calving  (is:ie,js:je) ) ;         ice_ocean_boundary%calving = 0.0
     allocate( ice_ocean_boundary%runoff_hflx   (is:ie,js:je) ) ;    ice_ocean_boundary%runoff_hflx = 0.0
     allocate( ice_ocean_boundary%calving_hflx  (is:ie,js:je) ) ;    ice_ocean_boundary%calving_hflx = 0.0
@@ -292,6 +293,9 @@ subroutine flux_ice_to_ocean ( Ice, Ocean, Ice_Ocean_Boundary )
     if(ASSOCIATED(Ice_Ocean_Boundary%runoff) ) call flux_ice_to_ocean_redistribute( Ice, Ocean, &
          Ice%runoff, Ice_Ocean_Boundary%runoff, Ice_Ocean_Boundary%xtype, do_area_weighted_flux )
 
+    if(ASSOCIATED(Ice_Ocean_Boundary%runoff_carbon) ) call flux_ice_to_ocean_redistribute( Ice, Ocean, &
+         Ice%runoff_carbon, Ice_Ocean_Boundary%runoff_carbon, Ice_Ocean_Boundary%xtype, do_area_weighted_flux )
+
     if(ASSOCIATED(Ice_Ocean_Boundary%calving) ) call flux_ice_to_ocean_redistribute( Ice, Ocean, &
          Ice%calving, Ice_Ocean_Boundary%calving, Ice_Ocean_Boundary%xtype, do_area_weighted_flux )
 
@@ -340,6 +344,7 @@ subroutine flux_ice_to_ocean_finish ( Time, Ice_Ocean_Boundary )
     call fms_data_override('OCN', 'lprec',     Ice_Ocean_Boundary%lprec    , Time )
     call fms_data_override('OCN', 'fprec',     Ice_Ocean_Boundary%fprec    , Time )
     call fms_data_override('OCN', 'runoff',    Ice_Ocean_Boundary%runoff   , Time )
+    call fms_data_override('OCN', 'runoff_carbon',    Ice_Ocean_Boundary%runoff_carbon   , Time )
     call fms_data_override('OCN', 'calving',   Ice_Ocean_Boundary%calving  , Time )
     call fms_data_override('OCN', 'runoff_hflx',    Ice_Ocean_Boundary%runoff_hflx   , Time )
     call fms_data_override('OCN', 'calving_hflx',   Ice_Ocean_Boundary%calving_hflx  , Time )
diff --git a/full/land_ice_flux_exchange.F90 b/full/land_ice_flux_exchange.F90
index 6240a4ad..9460cffe 100644
--- a/full/land_ice_flux_exchange.F90
+++ b/full/land_ice_flux_exchange.F90
@@ -77,11 +77,13 @@ subroutine land_ice_flux_exchange_init(Land, Ice, land_ice_boundary, Dt_cpl_in,
 
     !allocate land_ice_boundary
     allocate( land_ice_boundary%runoff(is:ie,js:je) )
+    allocate( land_ice_boundary%runoff_carbon(is:ie,js:je) )
     allocate( land_ice_boundary%calving(is:ie,js:je) )
     allocate( land_ice_boundary%runoff_hflx(is:ie,js:je) )
     allocate( land_ice_boundary%calving_hflx(is:ie,js:je) )
     ! initialize values for override experiments (mjh)
     land_ice_boundary%runoff=0.0
+    land_ice_boundary%runoff_carbon=0.0
     land_ice_boundary%calving=0.0
     land_ice_boundary%runoff_hflx=0.0
     land_ice_boundary%calving_hflx=0.0
@@ -109,10 +111,9 @@ subroutine flux_land_to_ice( Time, Land, Ice, Land_Ice_Boundary )
                                                                     !! fluxes passed from land to ice
 
     integer                         :: ier
-    real, dimension(n_xgrid_runoff) :: ex_runoff, ex_calving, ex_runoff_hflx, ex_calving_hflx
+    real, dimension(n_xgrid_runoff) :: ex_runoff, ex_calving, ex_runoff_hflx, ex_calving_hflx, ex_runoff_DOC
     real, dimension(size(Land_Ice_Boundary%runoff,1),size(Land_Ice_Boundary%runoff,2),1) :: ice_buf
 
-    !Balaji
     call fms_mpp_clock_begin(cplClock)
     call fms_mpp_clock_begin(fluxLandIceClock)
 
@@ -121,11 +122,14 @@ subroutine flux_land_to_ice( Time, Land, Ice, Land_Ice_Boundary )
 
     if (do_runoff) then
        call fms_xgrid_put_to_xgrid ( Land%discharge,      'LND', ex_runoff,  xmap_runoff)
+       call fms_xgrid_put_to_xgrid ( Land%discharge_DOC,  'LND', ex_runoff_DOC, xmap_runoff)
        call fms_xgrid_put_to_xgrid ( Land%discharge_snow, 'LND', ex_calving, xmap_runoff)
        call fms_xgrid_put_to_xgrid ( Land%discharge_heat,      'LND', ex_runoff_hflx,  xmap_runoff)
        call fms_xgrid_put_to_xgrid ( Land%discharge_snow_heat, 'LND', ex_calving_hflx, xmap_runoff)
        call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_runoff,  xmap_runoff)
        Land_Ice_Boundary%runoff = ice_buf(:,:,1);
+       call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_runoff_DOC, xmap_runoff)
+       Land_Ice_Boundary%runoff_carbon = ice_buf(:,:,1);
        call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_calving, xmap_runoff)
        Land_Ice_Boundary%calving = ice_buf(:,:,1);
        call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_runoff_hflx,  xmap_runoff)
@@ -134,6 +138,7 @@ subroutine flux_land_to_ice( Time, Land, Ice, Land_Ice_Boundary )
        Land_Ice_Boundary%calving_hflx = ice_buf(:,:,1);
        !Balaji
        call fms_data_override('ICE', 'runoff' , Land_Ice_Boundary%runoff , Time)
+       call fms_data_override('ICE', 'runoff_DOC' , Land_Ice_Boundary%runoff_carbon , Time)
        call fms_data_override('ICE', 'calving', Land_Ice_Boundary%calving, Time)
        call fms_data_override('ICE', 'runoff_hflx' , Land_Ice_Boundary%runoff_hflx , Time)
        call fms_data_override('ICE', 'calving_hflx', Land_Ice_Boundary%calving_hflx, Time)