introduced optional argument for interspersed coupling

full/coupler_main.F90

@@ -71,7 +71,7 @@
 !! The three components of coupler: @ref coupler_main , flux_exchange_mod, and surface_flux_mod
 !! are configured through three namelists
 !! * \ref coupler_config "coupler_nml"
-!! * \ref flux_exchange_conf "flux_exchange_nml" 
+!! * \ref flux_exchange_conf "flux_exchange_nml"
 !! * \ref surface_flux_config "surface_flux_nml"
 !!
 !!
@@ -880,6 +880,7 @@ program coupler_main
           if (do_ice .AND. Ice%fast_ice_pe) then
             if (ice_npes .NE. atmos_npes)call fms_mpp_set_current_pelist(Ice%fast_pelist)
             call update_ice_model_fast( Atmos_ice_boundary, Ice )
+            !MJH if (do_debug_concurrent_ice) then; chksum(Ice%FIA%bmelt)
           endif
           if (ice_npes .NE. atmos_npes) call fms_mpp_set_current_pelist(Atm%pelist)
           call fms_mpp_clock_end(newClockf)
@@ -1054,7 +1055,8 @@ program coupler_main
       if (combined_ice_and_ocean) then
         call flux_ice_to_ocean_stocks(Ice)
         call update_slow_ice_and_ocean(ice_ocean_driver_CS, Ice, Ocean_state, Ocean, &
-                      Ice_ocean_boundary, Time_ocean, Time_step_cpld )
+             Ice_ocean_boundary, Time_ocean, Time_step_cpld, OIB=ocean_ice_boundary )
+        !MJH
       else
       if (do_chksum) call ocean_chksum('update_ocean_model-', nc, Ocean, Ice_ocean_boundary)
       ! update_ocean_model since fluxes don't change here