<p><b>dwj07@fsu.edu</b> 2012-09-12 08:38:29 -0600 (Wed, 12 Sep 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Fixing some issues with the CESM build.<br>
</p><hr noshade><pre><font color="gray">Copied: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/MPASO_CplIndices.F (from rev 2152, branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F)
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/MPASO_CplIndices.F                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/MPASO_CplIndices.F        2012-09-12 14:38:29 UTC (rev 2153)
@@ -0,0 +1,125 @@
+module MPASO_CplIndices
+  
+  use seq_flds_mod
+  use mct_mod
+
+  implicit none
+
+  SAVE
+  public                               ! By default make data private
+
+  ! ocn -&gt; drv
+
+  integer :: index_o2x_So_t      
+  integer :: index_o2x_So_u
+  integer :: index_o2x_So_v
+  integer :: index_o2x_So_s
+  integer :: index_o2x_So_dhdx
+  integer :: index_o2x_So_dhdy
+  integer :: index_o2x_Fioo_q
+  integer :: index_o2x_Faoo_fco2_ocn
+  integer :: index_o2x_Faoo_fdms_ocn
+
+  ! drv -&gt; ocn
+
+  integer :: index_x2o_Si_ifrac        ! fractional ice wrt ocean
+  integer :: index_x2o_So_duu10n       ! 10m wind speed squared           (m^2/s^2)
+  integer :: index_x2o_Sa_pslv         ! sea-level pressure               (Pa)
+  integer :: index_x2o_Sa_co2prog      ! bottom atm level prognostic CO2
+  integer :: index_x2o_Sa_co2diag      ! bottom atm level diagnostic CO2
+  integer :: index_x2o_Foxx_taux       ! zonal wind stress (taux)         (W/m2   )
+  integer :: index_x2o_Foxx_tauy       ! meridonal wind stress (tauy)     (W/m2   )
+  integer :: index_x2o_Foxx_swnet      ! net short-wave heat flux         (W/m2   )
+  integer :: index_x2o_Foxx_sen        ! sensible heat flux               (W/m2   )
+  integer :: index_x2o_Foxx_lat        
+  integer :: index_x2o_Foxx_lwup       ! longwave radiation (up)          (W/m2   )
+  integer :: index_x2o_Faxa_lwdn       ! longwave radiation (down)        (W/m2   )
+  integer :: index_x2o_Fioi_melth      ! heat flux from snow &amp; ice melt   (W/m2   )
+  integer :: index_x2o_Fioi_meltw      ! snow melt flux                   (kg/m2/s)
+  integer :: index_x2o_Fioi_salt       ! salt                             (kg(salt)/m2/s)
+  integer :: index_x2o_Foxx_evap       ! evaporation flux                 (kg/m2/s)
+  integer :: index_x2o_Faxa_prec         
+  integer :: index_x2o_Faxa_snow       ! water flux due to snow           (kg/m2/s)
+  integer :: index_x2o_Faxa_rain       ! water flux due to rain           (kg/m2/s)
+  integer :: index_x2o_Faxa_bcphidry   ! flux: Black   Carbon hydrophilic dry deposition
+  integer :: index_x2o_Faxa_bcphodry   ! flux: Black   Carbon hydrophobic dry deposition
+  integer :: index_x2o_Faxa_bcphiwet   ! flux: Black   Carbon hydrophilic wet deposition
+  integer :: index_x2o_Faxa_ocphidry   ! flux: Organic Carbon hydrophilic dry deposition
+  integer :: index_x2o_Faxa_ocphodry   ! flux: Organic Carbon hydrophobic dry deposition
+  integer :: index_x2o_Faxa_ocphiwet   ! flux: Organic Carbon hydrophilic dry deposition
+  integer :: index_x2o_Faxa_dstwet1    ! flux: Size 1 dust -- wet deposition
+  integer :: index_x2o_Faxa_dstwet2    ! flux: Size 2 dust -- wet deposition
+  integer :: index_x2o_Faxa_dstwet3    ! flux: Size 3 dust -- wet deposition
+  integer :: index_x2o_Faxa_dstwet4    ! flux: Size 4 dust -- wet deposition
+  integer :: index_x2o_Faxa_dstdry1    ! flux: Size 1 dust -- dry deposition
+  integer :: index_x2o_Faxa_dstdry2    ! flux: Size 2 dust -- dry deposition
+  integer :: index_x2o_Faxa_dstdry3    ! flux: Size 3 dust -- dry deposition
+  integer :: index_x2o_Faxa_dstdry4    ! flux: Size 4 dust -- dry deposition
+  integer :: index_x2o_Forr_roff       ! river runoff flux                (kg/m2/s)
+  integer :: index_x2o_Forr_ioff       ! ice runoff flux                  (kg/m2/s)
+
+contains
+
+  subroutine MPASO_CplIndicesSet( )
+
+    type(mct_aVect) :: o2x      ! temporary
+    type(mct_aVect) :: x2o      ! temporary
+
+    ! Determine attribute vector indices
+
+    ! create temporary attribute vectors
+    call mct_aVect_init(x2o, rList=seq_flds_x2o_fields, lsize=1)
+    call mct_aVect_init(o2x, rList=seq_flds_o2x_fields, lsize=1)
+
+    index_o2x_So_t          = mct_avect_indexra(o2x,'So_t')
+    index_o2x_So_u          = mct_avect_indexra(o2x,'So_u')
+    index_o2x_So_v          = mct_avect_indexra(o2x,'So_v')
+    index_o2x_So_s          = mct_avect_indexra(o2x,'So_s')
+    index_o2x_So_dhdx       = mct_avect_indexra(o2x,'So_dhdx')
+    index_o2x_So_dhdy       = mct_avect_indexra(o2x,'So_dhdy')
+    index_o2x_Fioo_q        = mct_avect_indexra(o2x,'Fioo_q')
+    index_o2x_Faoo_fco2_ocn = mct_avect_indexra(o2x,'Faoo_fco2_ocn',perrWith='quiet')
+    index_o2x_Faoo_fdms_ocn = mct_avect_indexra(o2x,'Faoo_fdms_ocn',perrWith='quiet')
+    index_x2o_Si_ifrac      = mct_avect_indexra(x2o,'Si_ifrac')
+    index_x2o_Sa_pslv       = mct_avect_indexra(x2o,'Sa_pslv')
+    index_x2o_So_duu10n     = mct_avect_indexra(x2o,'So_duu10n')
+
+    index_x2o_Foxx_tauy     = mct_avect_indexra(x2o,'Foxx_tauy')
+    index_x2o_Foxx_taux     = mct_avect_indexra(x2o,'Foxx_taux')
+    index_x2o_Foxx_swnet    = mct_avect_indexra(x2o,'Foxx_swnet')
+    index_x2o_Foxx_lat      = mct_avect_indexra(x2o,'Foxx_lat')
+    index_x2o_Foxx_sen      = mct_avect_indexra(x2o,'Foxx_sen')
+    index_x2o_Foxx_lwup     = mct_avect_indexra(x2o,'Foxx_lwup')
+    index_x2o_Faxa_lwdn     = mct_avect_indexra(x2o,'Faxa_lwdn')
+    index_x2o_Fioi_melth    = mct_avect_indexra(x2o,'Fioi_melth')   
+    index_x2o_Fioi_meltw    = mct_avect_indexra(x2o,'Fioi_meltw')
+    index_x2o_Fioi_salt     = mct_avect_indexra(x2o,'Fioi_salt')   
+    index_x2o_Faxa_prec     = mct_avect_indexra(x2o,'Faxa_prec')   
+    index_x2o_Faxa_snow     = mct_avect_indexra(x2o,'Faxa_snow')   
+    index_x2o_Faxa_rain     = mct_avect_indexra(x2o,'Faxa_rain')   
+    index_x2o_Foxx_evap     = mct_avect_indexra(x2o,'Foxx_evap')
+    index_x2o_Forr_roff     = mct_avect_indexra(x2o,'Forr_roff')
+    index_x2o_Forr_ioff     = mct_avect_indexra(x2o,'Forr_ioff')
+    index_x2o_Faxa_bcphidry = mct_avect_indexra(x2o,'Faxa_bcphidry')
+    index_x2o_Faxa_bcphodry = mct_avect_indexra(x2o,'Faxa_bcphodry')
+    index_x2o_Faxa_bcphiwet = mct_avect_indexra(x2o,'Faxa_bcphiwet')
+    index_x2o_Faxa_ocphidry = mct_avect_indexra(x2o,'Faxa_ocphidry')
+    index_x2o_Faxa_ocphodry = mct_avect_indexra(x2o,'Faxa_ocphodry')
+    index_x2o_Faxa_ocphiwet = mct_avect_indexra(x2o,'Faxa_ocphiwet')
+    index_x2o_Faxa_dstdry1  = mct_avect_indexra(x2o,'Faxa_dstdry1')
+    index_x2o_Faxa_dstdry2  = mct_avect_indexra(x2o,'Faxa_dstdry2')
+    index_x2o_Faxa_dstdry3  = mct_avect_indexra(x2o,'Faxa_dstdry3')
+    index_x2o_Faxa_dstdry4  = mct_avect_indexra(x2o,'Faxa_dstdry4')
+    index_x2o_Faxa_dstwet1  = mct_avect_indexra(x2o,'Faxa_dstwet1')
+    index_x2o_Faxa_dstwet2  = mct_avect_indexra(x2o,'Faxa_dstwet2')
+    index_x2o_Faxa_dstwet3  = mct_avect_indexra(x2o,'Faxa_dstwet3')
+    index_x2o_Faxa_dstwet4  = mct_avect_indexra(x2o,'Faxa_dstwet4')
+    index_x2o_Sa_co2prog    = mct_avect_indexra(x2o,'Sa_co2prog',perrWith='quiet')
+    index_x2o_Sa_co2diag    = mct_avect_indexra(x2o,'Sa_co2diag',perrWith='quiet')
+
+    call mct_aVect_clean(x2o)
+    call mct_aVect_clean(o2x)
+
+  end subroutine MPASO_CplIndicesSet
+
+end module MPASO_CplIndices

Modified: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/Makefile
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/Makefile        2012-09-11 22:19:34 UTC (rev 2152)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/Makefile        2012-09-12 14:38:29 UTC (rev 2153)
@@ -1,16 +1,21 @@
-all: cpl_mct cpl_share
-        ar ru libdrver.a cpl_mct/*.o
-        ar ru libdrver.a cpl_share/*.o
-        cp cpl_mct/*.o .
-        cp cpl_share/*.o .
+.SUFFIXES: .F .o
 
-cpl_mct:
-        ( cd cpl_mct; $(MAKE) all )
+OBJS = ocn_comp_mct.o \
+       ocn_communicator.o \
+       MPASO_CplIndices.o
 
-cpl_share:
-        (cd cpl_share; $(MAKE) all )
+all: $(OBJS)
 
+ocn_comp_mct.o: ocn_communicator.o
+
+ocn_communicator.o:
+
+MPASO_CplIndices.o:
+
 clean:
-        $(RM) *.o libdrver.a
-        ( cd cpl_mct; $(MAKE) clean )
-        ( cd cpl_share; $(MAKE) clean )
+        $(RM) *.o *.mod *.f90
+
+.F.o:
+        $(RM) $@ $*.mod
+        $(CPP) $(CPPFLAGS) $(CPPINCLUDES) $&lt; &gt; $*.f90
+        $(FC) $(FFLAGS) -c $*.f90 $(FCINCLUDES) -I../framework -I../core_$(CORE) -I../external/esmf_time_f90

Copied: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/ocn_communicator.F (from rev 2152, branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F)
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/ocn_communicator.F                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/ocn_communicator.F        2012-09-12 14:38:29 UTC (rev 2153)
@@ -0,0 +1,5 @@
+module ocn_communicator
+
+  integer, public :: mpi_communicator_ocn
+
+end module ocn_communicator

Copied: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/ocn_comp_mct.F (from rev 2152, branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F)
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/ocn_comp_mct.F                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/ocn_comp_mct.F        2012-09-12 14:38:29 UTC (rev 2153)
@@ -0,0 +1,1616 @@
+module ocn_comp_mct
+
+!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
+!BOP
+! !MODULE: ocn_comp_mct
+! !INTERFACE:
+
+! !DESCRIPTION:
+!  This is the main driver for the Model for Predication Across Scales Ocean Model (MPAS-O).
+!
+! !REVISION HISTORY:
+!  SVN:$Id:
+!
+! !USES:
+   use mct_mod
+   use esmf
+   use seq_flds_mod
+   use seq_cdata_mod
+   use seq_infodata_mod
+   use seq_timemgr_mod
+   use seq_comm_mct,      only : seq_comm_suffix, seq_comm_inst, seq_comm_name
+   use shr_file_mod 
+   use shr_cal_mod,       only : shr_cal_date2ymd
+   use shr_sys_mod
+   use perf_mod
+   use ocn_communicator,  only: mpi_communicator_ocn
+
+   use MPASO_CplIndices
+
+   use mpas_framework
+   use mpas_core
+   use mpas_kind_types
+!
+! !PUBLIC MEMBER FUNCTIONS:
+  implicit none
+  public :: ocn_init_mct
+  public :: ocn_run_mct
+  public :: ocn_final_mct
+  SAVE
+  private                              ! By default make data private
+
+!
+! ! PUBLIC DATA:
+!
+! !REVISION HISTORY:
+! Author: Mariana Vertenstein
+!
+!EOP
+! !PRIVATE MODULE FUNCTIONS:
+  private :: ocn_export_mct
+  private :: ocn_import_mct
+  private :: ocn_SetGSMap_mct
+  private :: ocn_domain_mct
+!
+! !PRIVATE MODULE VARIABLES
+
+  logical ::   &amp;
+       ldiag_cpl = .false.
+
+  integer, private ::   &amp;
+      cpl_write_restart,   &amp;! flag id for write restart
+      cpl_write_history,   &amp;! flag id for write history
+      cpl_write_tavg,      &amp;! flag id for write tavg      
+      cpl_diag_global,     &amp;! flag id for computing diagnostics
+      cpl_diag_transp       ! flag id for computing diagnostics
+
+  real (kind=RKIND),   &amp;
+      dimension(:,:,:,:), allocatable ::  &amp;
+      SBUFF_SUM           ! accumulated sum of send buffer quantities
+                          ! for averaging before being sent
+   real (kind=RKIND) ::  &amp;
+     tlast_coupled
+
+   integer  ::   &amp;
+      nsend, nrecv
+
+   character(len=StrKIND) :: &amp;
+      runtype         
+
+   type(seq_infodata_type), pointer :: &amp;
+      infodata   
+
+   !! MPAS-O Datatypes
+   type (dm_info), pointer :: dminfo
+   type (domain_type), pointer :: domain
+   type (io_output_object), save :: output_obj
+   integer :: output_frame
+
+!=======================================================================
+
+contains
+
+!***********************************************************************
+!BOP
+!
+! !IROUTINE: ocn_init_mct
+!
+! !INTERFACE:
+  subroutine ocn_init_mct( EClock, cdata_o, x2o_o, o2x_o, NLFilename )
+!
+! !DESCRIPTION:
+! Initialize POP 
+!
+! !INPUT/OUTPUT PARAMETERS:
+
+    type(ESMF_Clock)            , intent(in)    :: EClock
+    type(seq_cdata)             , intent(inout) :: cdata_o
+    type(mct_aVect)             , intent(inout) :: x2o_o, o2x_o
+    character(len=*), optional  , intent(in)    :: NLFilename ! Namelist filename
+!
+! !REVISION HISTORY:
+! Author: Mariana Vertenstein
+!EOP
+!-----------------------------------------------------------------------
+!
+!  local variables
+!
+!-----------------------------------------------------------------------
+
+    integer ::  &amp;
+       OCNID,       &amp;
+       mpicom_o,    &amp;
+       lsize,       &amp;
+       start_ymd,   &amp;
+       start_tod,   &amp;
+       start_year,  &amp;
+       start_day,   &amp;
+       start_month, &amp;
+       start_hour,  &amp;
+       iyear,       &amp;
+       ocn_cpl_dt,  &amp;
+       mpas_o_cpl_dt,  &amp;
+       shrlogunit,  &amp;  ! old values
+       shrloglev       ! old values
+
+    type(mct_gsMap), pointer :: &amp;
+       gsMap_o
+
+    type(mct_gGrid), pointer :: &amp;
+       dom_o
+
+    integer :: &amp;
+       errorCode         ! error code
+
+    integer :: &amp;
+       nThreads
+
+    real (kind=RKIND) ::  &amp;
+       precadj
+
+    integer :: iam,ierr 
+    character(len=StrKIND)  :: starttype          ! infodata start type
+
+    integer :: lbnum
+
+!-----------------------------------------------------------------------
+!
+!  set cdata pointers
+!
+!-----------------------------------------------------------------------
+
+   errorCode = 0
+
+   call seq_cdata_setptrs(cdata_o, ID=OCNID, mpicom=mpicom_o, &amp;
+        gsMap=gsMap_o, dom=dom_o, infodata=infodata)
+
+!  MPASO_MCT_OCNID =  OCNID
+!  MPASO_MCT_gsMap_o =&gt; gsMap_o
+!  MPASO_MCT_dom_o   =&gt; dom_o
+
+#if (defined _MEMTRACE)
+    call MPI_comm_rank(mpicom_o,iam,ierr)
+    if(iam == 0) then
+        lbnum=1
+        call memmon_dump_fort('memmon.out','ocn_init_mct:start::',lbnum) 
+    endif
+#endif
+
+
+    ! The following communicator module variable will be utilized in init_communicate that
+    ! is called by initial - this is done to make the code backwards compatible
+
+    mpi_communicator_ocn = mpicom_o
+
+!-----------------------------------------------------------------------
+!
+!  initialize the model run 
+!
+!-----------------------------------------------------------------------
+
+    call MPASO_CplIndicesSet()
+
+!   call seq_infodata_GetData( infodata, case_name=runid )
+   
+    call seq_infodata_GetData( infodata, start_type=starttype)
+
+
+    !!!DWJ CHECK RUN TYPES
+    if (     trim(starttype) == trim(seq_infodata_start_type_start)) then
+       runtype = &quot;initial&quot;
+    else if (trim(starttype) == trim(seq_infodata_start_type_cont) ) then
+       runtype = &quot;continue&quot;
+    else if (trim(starttype) == trim(seq_infodata_start_type_brnch)) then
+       runtype = &quot;branch&quot;
+    else
+       call mpas_dmpar_global_abort(' ocn_comp_mct ERROR: unknown starttype')
+    end if
+    !!!DWJ CHECK RUN TYPES
+
+
+    !TODO: check for consistency of pop runid and runtype with seq_infodata
+   
+!-----------------------------------------------------------------------
+!
+!  first initializaiton phase of mpas-o
+!  initialize mpas-o because grid information is needed for
+!  creation of GSMap_ocn.
+!  call mpas-o initialization routines
+!
+!-----------------------------------------------------------------------
+
+!  inst_name   = seq_comm_name(OCNID)
+!  inst_index  = seq_comm_inst(OCNID)
+!  inst_suffix = seq_comm_suffix(OCNID)
+
+   call t_startf('mpas-o_init')
+   call mpas_framework_init(dminfo, domain)
+
+!-----------------------------------------------------------------------
+!
+!  register non-standard incoming fields
+!
+!-----------------------------------------------------------------------
+
+!  if (index_x2o_Sa_co2prog &gt; 0) then
+!     call named_field_register('ATM_CO2_PROG', ATM_CO2_PROG_nf_ind)
+!  endif
+!  if (index_x2o_Sa_co2diag &gt; 0) then
+!     call named_field_register('ATM_CO2_DIAG', ATM_CO2_DIAG_nf_ind)
+!  endif
+!  call register_string('mpas-o_init_coupled')
+!  call flushm (stdout)
+
+!-----------------------------------------------------------------------
+!
+!  second initialization phase of mpas-o
+!
+!-----------------------------------------------------------------------
+
+   call mpas_input_state_for_domain(domain)
+
+!  call mpas_core_init(domain, timeStamp)
+
+!-----------------------------------------------------------------------
+!
+!  initialize time-stamp information
+!
+!-----------------------------------------------------------------------
+
+!  call ccsm_char_date_and_time
+
+   call t_stopf ('mpas-o_init')
+
+!----------------------------------------------------------------------------
+!
+! reset shr logging to my log file
+!
+!----------------------------------------------------------------------------
+
+    call shr_file_getLogUnit (shrlogunit)
+    call shr_file_getLogLevel(shrloglev)
+!   call shr_file_setLogUnit (stdout)
+   
+!-----------------------------------------------------------------------
+!
+!  check for consistency of mpas-o and sync clock initial time
+!
+!-----------------------------------------------------------------------
+
+   if (runtype == 'initial') then
+      call seq_timemgr_EClockGetData(EClock, &amp;
+           start_ymd=start_ymd, start_tod=start_tod)
+      call shr_cal_date2ymd(start_ymd,start_year,start_month,start_day)
+
+!     if (iyear0 /= start_year) then
+!        call mpas_dmpar_global_abort(' iyear0 does not match start_year')
+!     end if
+!     if (imonth0 /= start_month) then
+!        call mpas_dmpar_global_abort(' imonth0 does not match start_year')
+!     end if
+!#ifndef _HIRES 
+!      if (seconds_this_day /= start_tod) then
+!         call document ('ocn_init_mct', 'sec0     ', seconds_this_day)
+!         call document ('ocn_init_mct', 'start_tod ', start_tod)
+!         call exit_POP(sigAbort,' sec0 does not start_tod')
+!      end if
+!#endif
+   end if
+
+!-----------------------------------------------------------------------
+!
+!  initialize MCT attribute vectors and indices
+!
+!-----------------------------------------------------------------------
+
+   call t_startf ('mpas-o_mct_init')
+
+   call ocn_SetGSMap_mct( mpicom_o, OCNID, GSMap_o )
+   lsize = mct_gsMap_lsize(gsMap_o, mpicom_o)
+
+   ! Initialize mct ocn domain (needs ocn initialization info)
+   
+   call ocn_domain_mct( lsize, gsMap_o, dom_o )
+   
+   ! Inialize mct attribute vectors
+   
+   call mct_aVect_init(x2o_o, rList=seq_flds_x2o_fields, lsize=lsize)
+   call mct_aVect_zero(x2o_o)
+   
+   call mct_aVect_init(o2x_o, rList=seq_flds_o2x_fields, lsize=lsize) 
+   call mct_aVect_zero(o2x_o)
+   
+   nsend = mct_avect_nRattr(o2x_o)
+   nrecv = mct_avect_nRattr(x2o_o)
+!  allocate (SBUFF_SUM(nx_block,ny_block,max_blocks_clinic,nsend))
+
+!-----------------------------------------------------------------------
+!
+!  Initialize flags and shortwave absorption profile
+!  Note that these cpl_write_xxx flags have no freqency options
+!  set; therefore, they will retain a default value of .false.
+!  unless they are explicitly set .true.  at the appropriate times
+!
+!-----------------------------------------------------------------------
+
+!  call init_time_flag('cpl_write_restart',cpl_write_restart, owner = 'ocn_init_mct')
+!  call init_time_flag('cpl_write_history',cpl_write_history, owner = 'ocn_init_mct')
+!  call init_time_flag('cpl_write_tavg'   ,cpl_write_tavg,    owner = 'ocn_init_mct')
+!  call init_time_flag('cpl_diag_global'  ,cpl_diag_global,   owner = 'ocn_init_mct')
+!  call init_time_flag('cpl_diag_transp'  ,cpl_diag_transp,   owner = 'ocn_init_mct')
+
+!  lsmft_avail = .true.
+!  tlast_coupled = c0
+
+!-----------------------------------------------------------------------
+!
+!   initialize necessary  coupling info
+!
+!-----------------------------------------------------------------------
+
+    call seq_timemgr_EClockGetData(EClock, dtime=ocn_cpl_dt)
+!   mpas_o_cpl_dt = seconds_in_day / ncouple_per_day
+!   if (mpas_o_cpl_dt /= ocn_cpl_dt) then
+!      write(stdout,*)'mpas-o_cpl_dt= ',mpas_o_cpl_dt, &amp;
+!                    ' ocn_cpl_dt= ',ocn_cpl_dt   
+!      call mpas_dmpar_global_abort('ERROR mpas-o_cpl_dt and ocn_cpl_dt must be identical')
+!   end if
+
+!-----------------------------------------------------------------------
+!
+!  send intial state to driver
+!
+!-----------------------------------------------------------------------
+
+!  if ( lsend_precip_fact )  then
+!     precadj = precip_fact * 1.0e6_r8  
+!     call seq_infodata_PutData( infodata, precip_fact=precadj)
+!  end if
+
+   call ocn_export_mct(o2x_o, errorCode)  
+   if (errorCode /= 0) then
+      call mpas_dmpar_global_abort('ERROR in ocn_export_mct')
+   endif
+
+   call t_stopf ('mpas-o_mct_init')
+
+!  call seq_infodata_PutData( infodata, ocn_nx = nx_global , ocn_ny = ny_global)
+   call seq_infodata_PutData( infodata, ocn_prognostic=.true., ocnrof_prognostic=.true.)
+
+!----------------------------------------------------------------------------
+!
+! Reset shr logging to original values
+!
+!----------------------------------------------------------------------------
+
+    call shr_file_setLogUnit (shrlogunit)
+    call shr_file_setLogLevel(shrloglev)
+
+#if (defined _MEMTRACE)
+    if(iam  == 0) then
+!        write(6,*) 'ocn_init_mct:end::'
+        lbnum=1
+        call memmon_dump_fort('memmon.out','ocn_init_mct:end::',lbnum) 
+        call memmon_reset_addr()
+    endif
+#endif
+
+!-----------------------------------------------------------------------
+!
+!  document orbital parameters
+!
+!-----------------------------------------------------------------------
+
+!  if (registry_match('qsw_distrb_iopt_cosz')) then
+!  call seq_infodata_GetData(infodata, &amp;
+!     orb_eccen=orb_eccen, orb_mvelpp=orb_mvelpp, orb_lambm0=orb_lambm0, orb_obliqr=orb_obliqr)
+
+!    write(stdout,*) ' '
+!   endif
+
+!-----------------------------------------------------------------------
+!
+!  Now document all time flags, because this is the last step of pop2 
+!    initialization
+!
+!-----------------------------------------------------------------------
+
+!  call document_time_flags
+
+!-----------------------------------------------------------------------
+!
+!  output delimiter to log file
+!
+!-----------------------------------------------------------------------
+
+   !DWJ IO Units
+!  if (my_task == master_task) then
+!     write(stdout,blank_fmt)
+!     write(stdout,'(&quot; End of initialization&quot;)')
+!     write(stdout,blank_fmt)
+!     write(stdout,ndelim_fmt)
+!     call POP_IOUnitsFlush(POP_stdout)
+#ifdef CCSMCOUPLED
+!     call POP_IOUnitsFlush(stdout)
+#endif
+!  endif
+
+!-----------------------------------------------------------------------
+!EOC
+
+ end subroutine ocn_init_mct
+
+!***********************************************************************
+!BOP
+!
+! !IROUTINE: ocn_run_mct
+!
+! !INTERFACE:
+  subroutine ocn_run_mct( EClock, cdata_o, x2o_o, o2x_o)
+
+    use mpas_kind_types
+    use mpas_grid_types
+    use mpas_io_output
+    use mpas_timer
+
+    implicit none
+!
+! !DESCRIPTION:
+! Run POP for a coupling interval
+!
+! !INPUT/OUTPUT PARAMETERS:
+    type(ESMF_Clock)            , intent(in)    :: EClock
+    type(seq_cdata)             , intent(inout) :: cdata_o
+    type(mct_aVect)             , intent(inout) :: x2o_o
+    type(mct_aVect)             , intent(inout) :: o2x_o
+
+!!
+!! !REVISION HISTORY:
+!! Author: Mariana Vertenstein
+!!EOP
+!!-----------------------------------------------------------------------
+!!
+!!  local variables
+!!
+!!-----------------------------------------------------------------------
+!
+!    integer :: errorCode           ! error flag
+!
+!    integer :: itimestep
+!
+!    integer :: &amp;
+!         ymd, &amp;          ! POP2 current date (YYYYMMDD)
+!         tod, &amp;          ! POP2 current time of day (sec)
+!         ymd_sync, &amp;     ! Sync clock current date (YYYYMMDD)
+!         tod_sync, &amp;     ! Sync clcok current time of day (sec)
+!         shrlogunit,  &amp;  ! old values
+!         shrloglev       ! old values
+!
+!    character(len=char_len_long) :: &amp;
+!         fname
+!
+!    type(seq_infodata_type), pointer :: &amp;
+!       infodata   ! Input init object
+!
+!    real (kind=RKIND) ::  &amp;
+!         precadj
+!
+!    logical :: &amp;
+!         lcoupled,   &amp;  ! temporary
+!         rstwr,      &amp;  ! true =&gt; write restart at end of day
+!         first_time = .true.
+!
+!    character (char_len)  :: message
+!
+!    integer(int_kind) :: info_debug
+!
+!    integer :: lbnum
+!
+!#if (defined _MEMTRACE)
+!    if(my_task == 0 ) then
+!       lbnum=1
+!       call memmon_dump_fort('memmon.out',SubName//':start::',lbnum) 
+!    endif
+!#endif
+!
+!!-----------------------------------------------------------------------
+!!
+!!  start up the main timer
+!!
+!!-----------------------------------------------------------------------
+!
+!   call mpas_timer_start('total time')
+!
+!
+!!-----------------------------------------------------------------------
+!!
+!! reset shr logging to my log file
+!!
+!!----------------------------------------------------------------------------
+!
+!    errorCode = POP_Success
+!
+!    call shr_file_getLogUnit (shrlogunit)
+!    call shr_file_getLogLevel(shrloglev)
+!    call shr_file_setLogUnit (stdout)
+!
+!    call seq_cdata_setptrs(cdata_o, infodata=infodata)
+!
+!!----------------------------------------------------------------------------
+!!
+!! restart flag (rstwr) will assume only an eod restart for now
+!!
+!!----------------------------------------------------------------------------
+!
+!    call seq_infodata_GetData( infodata, info_debug=info_debug)
+!    if (info_debug &gt;= 2) then
+!       ldiag_cpl = .true. 
+!    else
+!       ldiag_cpl = .false.
+!    endif
+!
+!    rstwr = seq_timemgr_RestartAlarmIsOn(EClock)
+!
+!!-----------------------------------------------------------------------
+!!
+!!  advance the model in time over coupling interval
+!!  write restart dumps and archiving
+!!
+!!-----------------------------------------------------------------------
+!
+!    ! Note that all ocean time flags are evaluated each timestep in time_manager
+!    ! ocn_import_mct is analogous to pop_unpack_fluxes_from_coupler in cpl6 
+!    ! ocn_export_mct is analogous to prepare_send_to_coupler in cpl6
+!    ! tlast_coupled is set to zero at the end of ocn_export_mct
+!
+!    itimestep = 0
+!    timestep: do
+!       ! obtain import state from driver
+!       if (check_time_flag(cpl_ts) .or. nsteps_run == 0) then
+!          call ocn_import_mct(x2o_o, errorCode)   
+!
+!          if (errorCode /= POP_Success) then
+!             call POP_ErrorPrint(errorCode)
+!             call exit_POP(sigAbort, 'ERROR in step')
+!          endif
+!
+!          !DWJ Need to write on mpas side
+!!         call pop_set_coupled_forcing 
+!       end if
+!
+!
+!    end do timestep
+!
+!    advance: do 
+!
+!       call step(errorCode)
+!
+!       if (errorCode /= POP_Success) then
+!          call POP_ErrorPrint(errorCode)
+!          call exit_POP(sigAbort, 'ERROR in step')
+!       endif
+!
+!       if (check_KE(100.0_r8)) then
+!          !*** exit if energy is blowing
+!          call output_driver
+!          call exit_POP(sigAbort,'ERROR: k.e. &gt; 100 ')
+!       endif
+!       call output_driver
+!       
+!       ! return export state to driver
+!       call pop_sum_buffer()  
+!       if (check_time_flag(cpl_ts)) then
+!          call ocn_export_mct(o2x_o, errorCode)
+!          if (errorCode /= POP_Success) then
+!             call POP_ErrorPrint(errorCode)
+!             call exit_POP(sigAbort, 'ERROR in ocn_export_mct')
+!          endif
+!
+!          exit advance
+!       end if
+!       
+!    enddo advance
+!
+!    if ( lsend_precip_fact ) then
+!       precadj = precip_fact * 1.0e6_r8  
+!       call seq_infodata_PutData( infodata, precip_fact=precadj )
+!    end if
+!    
+!!--------------------------------------------------------------------
+!!
+!! check that internal clock is in sync with master clock
+!!
+!!--------------------------------------------------------------------
+!
+!    ymd = iyear*10000 + imonth*100 + iday
+!    tod = ihour*seconds_in_hour + iminute*seconds_in_minute + isecond
+!    if ( .not. seq_timemgr_EClockDateInSync( EClock, ymd, tod ) )then
+!       call seq_timemgr_EClockGetData( EClock, curr_ymd=ymd_sync, &amp;
+!          curr_tod=tod_sync )
+!       write(stdout,*)' pop2 ymd=',ymd     ,'  pop2 tod= ',tod
+!       write(stdout,*)' sync ymd=',ymd_sync,'  sync tod= ',tod_sync
+!       write(stdout,*)' Internal pop2 clock not in sync with Sync Clock'
+!       call shr_sys_abort( SubName// &amp;
+!          &quot;:: Internal pop2 clock not in sync with Sync Clock&quot;)
+!    end if
+!   
+!!----------------------------------------------------------------------------
+!!
+!! Reset shr logging to original values
+!!
+!!----------------------------------------------------------------------------
+!
+!   call shr_file_setLogUnit (shrlogunit)
+!   call shr_file_setLogLevel(shrloglev)
+!
+!   call mpas_timer_stop('total time')
+!
+!#if (defined _MEMTRACE)
+!    if(my_task == 0) then
+!       lbnum=1
+!       call memmon_dump_fort('memmon.out',SubName//':end::',lbnum) 
+!       call memmon_reset_addr()
+!    endif
+!#endif
+!-----------------------------------------------------------------------
+!EOC
+
+  end subroutine ocn_run_mct
+
+!***********************************************************************
+!BOP
+!
+! !IROUTINE: ocn_final_mct
+!
+! !INTERFACE:
+  subroutine ocn_final_mct( EClock, cdata_o, x2o_o, o2x_o)
+!
+! !DESCRIPTION:
+! Finalize POP
+!
+! !USES:
+! !ARGUMENTS:
+    type(ESMF_Clock)            , intent(in)    :: EClock
+    type(seq_cdata)             , intent(inout) :: cdata_o
+    type(mct_aVect)             , intent(inout) :: x2o_o
+    type(mct_aVect)             , intent(inout) :: o2x_o
+!
+! !REVISION HISTORY:
+! Author: Fei Liu
+!EOP
+!BOC
+!-----------------------------------------------------------------------
+!
+!  local variables
+!
+!-----------------------------------------------------------------------
+
+   integer :: &amp;
+      errorCode         ! error code
+
+!-----------------------------------------------------------------------
+
+    call mpas_output_state_finalize(output_obj, domain % dminfo)
+    call mpas_timer_write()
+    call mpas_core_finalize(domain)
+    call mpas_framework_finalize(dminfo, domain)
+
+  end subroutine ocn_final_mct
+
+!***********************************************************************
+!BOP
+!IROUTINE: ocn_SetGSMap_mct
+! !INTERFACE:
+
+ subroutine ocn_SetGSMap_mct( mpicom_ocn, OCNID, gsMap_ocn )
+
+   use mpas_dmpar
+
+! !DESCRIPTION:
+!  This routine mct global seg maps for the pop decomposition
+!
+! !REVISION HISTORY:
+!  same as module
+
+! !INPUT/OUTPUT PARAMETERS:
+
+    implicit none
+    integer        , intent(in)    :: mpicom_ocn
+    integer        , intent(in)    :: OCNID
+    type(mct_gsMap), intent(inout) :: gsMap_ocn
+
+!EOP
+!BOC
+!-----------------------------------------------------------------------
+!
+!  local variables
+!
+!-----------------------------------------------------------------------
+
+    integer,allocatable :: &amp;
+      gindex(:)
+
+    integer ::   &amp;
+      i,j, k, n, iblock, &amp;
+      lsize, gsize,   &amp;
+      ier
+
+    type (block_type), pointer :: block_ptr 
+
+!-----------------------------------------------------------------------
+!  Build the POP grid numbering for MCT
+!  NOTE:  Numbering scheme is: West to East and South to North starting
+!  at the south pole.  Should be the same as what's used in SCRIP
+!-----------------------------------------------------------------------
+
+    n = 0
+    block_ptr =&gt; domain % blocklist
+    do while(associated(block_ptr))
+      n = n + block_ptr % mesh % nCellsSolve
+      block_ptr =&gt; block_ptr % next
+    end do
+
+    lsize = n
+    call mpas_dmpar_sum_int(dminfo, lsize, gsize)
+    allocate(gindex(lsize),stat=ier)
+
+    n = 0
+    block_ptr =&gt; domain % blocklist
+    do while(associated(block_ptr))
+      do i = 1, block_ptr % mesh % nCellsSolve
+        n = n + 1
+        gindex(n) = block_ptr % mesh % indexToCellID % array(i)
+      end do
+      block_ptr =&gt; block_ptr % next
+    end do
+
+    call mct_gsMap_init( gsMap_ocn, gindex, mpicom_ocn, OCNID, lsize, gsize )
+
+    deallocate(gindex)
+
+!-----------------------------------------------------------------------
+!EOC
+
+  end subroutine ocn_SetGSMap_mct
+
+!***********************************************************************
+!BOP
+! !IROUTINE: ocn_domain_mct
+! !INTERFACE:
+
+ subroutine ocn_domain_mct( lsize, gsMap_o, dom_o )
+
+! !DESCRIPTION:
+!  This routine mct global seg maps for the pop decomposition
+!
+! !REVISION HISTORY:
+!  same as module
+!
+! !INPUT/OUTPUT PARAMETERS:
+
+    implicit none
+    integer        , intent(in)    :: lsize
+    type(mct_gsMap), intent(in)    :: gsMap_o
+    type(mct_ggrid), intent(inout) :: dom_o     
+
+!EOP
+!BOC
+!-----------------------------------------------------------------------
+!
+!  local variables
+!
+!-----------------------------------------------------------------------
+
+    integer, pointer :: idata(:)
+
+    real(kind=RKIND), pointer :: data(:)
+    real(kind=RKIND) :: r2d
+
+    integer :: i,j, k, n, ier
+
+    type (block_type), pointer :: block_ptr
+
+    r2d = 180.0/pii
+
+!-------------------------------------------------------------------
+!
+!  initialize mct domain type, lat/lon in degrees,
+!  area in radians^2, mask is 1 (ocean), 0 (non-ocean)
+!
+!-------------------------------------------------------------------
+
+    call mct_gGrid_init( GGrid=dom_o, CoordChars=trim(seq_flds_dom_coord), &amp;
+       OtherChars=trim(seq_flds_dom_other), lsize=lsize )
+    call mct_aVect_zero(dom_o%data)
+    allocate(data(lsize))
+
+!-------------------------------------------------------------------
+!
+! Determine global gridpoint number attribute, GlobGridNum, which is set automatically by MCT
+!
+!-------------------------------------------------------------------
+
+!   call mct_gsMap_orderedPoints(gsMap_o, my_task, idata)
+    call mct_gGrid_importIAttr(dom_o,'GlobGridNum',idata,lsize)
+
+!-------------------------------------------------------------------
+!
+! Determine domain (numbering scheme is: West to East and South to North to South pole)
+! Initialize attribute vector with special value
+!
+!-------------------------------------------------------------------
+
+    data(:) = -9999.0_R8 
+    call mct_gGrid_importRAttr(dom_o,&quot;lat&quot;  ,data,lsize) 
+    call mct_gGrid_importRAttr(dom_o,&quot;lon&quot;  ,data,lsize) 
+    call mct_gGrid_importRAttr(dom_o,&quot;area&quot; ,data,lsize) 
+    call mct_gGrid_importRAttr(dom_o,&quot;aream&quot;,data,lsize) 
+    data(:) = 0.0_R8     
+    call mct_gGrid_importRAttr(dom_o,&quot;mask&quot;,data,lsize) 
+    call mct_gGrid_importRAttr(dom_o,&quot;frac&quot;,data,lsize) 
+
+!-------------------------------------------------------------------
+!
+! Fill in correct values for domain components
+!
+!-------------------------------------------------------------------
+
+    n = 0
+    block_ptr =&gt; domain % blocklist
+    do while(associated(block_ptr))
+      do i = 1, block_ptr % mesh % nCellsSolve
+        n = n + 1
+        data(n) = block_ptr % mesh % lonCell % array(i) * r2d
+      end do
+      
+      block_ptr =&gt; block_ptr % next
+    end do
+    call mct_gGrid_importRattr(dom_o,&quot;lon&quot;,data,lsize) 
+
+    n = 0
+    block_ptr =&gt; domain % blocklist
+    do while(associated(block_ptr))
+      do i = 1, block_ptr % mesh % nCellsSolve
+        n = n + 1
+        data(n) = block_ptr % mesh % latCell % array(i) * r2d
+      end do
+      block_ptr =&gt; block_ptr % next
+    end do
+    call mct_gGrid_importRattr(dom_o,&quot;lat&quot;,data,lsize) 
+
+    n = 0
+    block_ptr =&gt; domain % blocklist
+    do while(associated(block_ptr))
+      do i = 1, block_ptr % mesh % nCellsSolve
+        n = n + 1
+        data(n) = block_ptr % mesh % areaCell % array(i)
+!       data(n) = block_ptr % mesh % areaCell % array(i) / (sphere_radius * sphere_radius)
+      end do
+      block_ptr =&gt; block_ptr % next
+    end do
+    call mct_gGrid_importRattr(dom_o,&quot;area&quot;,data,lsize) 
+
+    data(:) = 1.0_RKIND ! No land cells in MPAS-O, only Ocean cells
+    call mct_gGrid_importRattr(dom_o,&quot;mask&quot;,data,lsize) 
+    call mct_gGrid_importRattr(dom_o,&quot;frac&quot;,data,lsize) 
+
+    deallocate(data)
+    deallocate(idata)
+
+!-----------------------------------------------------------------------
+!EOC
+
+  end subroutine ocn_domain_mct
+
+
+!***********************************************************************
+!BOP
+! !IROUTINE: ocn_import_mct
+! !INTERFACE:
+
+ subroutine ocn_import_mct(x2o_o, errorCode)
+
+! !DESCRIPTION:
+!-----------------------------------------------------------------------
+!  This routine receives message from cpl7 driver
+!
+!    The following fields are always received from the coupler:
+! 
+!    o  taux   -- zonal wind stress (taux)                 (W/m2   )
+!    o  tauy   -- meridonal wind stress (tauy)             (W/m2   )
+!    o  snow   -- water flux due to snow                   (kg/m2/s)
+!    o  rain   -- water flux due to rain                   (kg/m2/s)
+!    o  evap   -- evaporation flux                         (kg/m2/s)
+!    o  meltw  -- snow melt flux                           (kg/m2/s)
+!    o  salt   -- salt                                     (kg(salt)/m2/s)
+!    o  swnet  -- net short-wave heat flux                 (W/m2   )
+!    o  sen    -- sensible heat flux                       (W/m2   )
+!    o  lwup   -- longwave radiation (up)                  (W/m2   )
+!    o  lwdn   -- longwave radiation (down)                (W/m2   )
+!    o  melth  -- heat flux from snow&amp;ice melt             (W/m2   )
+!    o  ifrac  -- ice fraction
+!    o  roff   -- river runoff flux                        (kg/m2/s)
+!    o  ioff   -- ice runoff flux                          (kg/m2/s)
+! 
+!    The following fields are sometimes received from the coupler,
+!      depending on model options:
+! 
+!    o  pslv   -- sea-level pressure                       (Pa)
+!    o  duu10n -- 10m wind speed squared                   (m^2/s^2)
+!    o  co2prog-- bottom atm level prognostic co2
+!    o  co2diag-- bottom atm level diagnostic co2
+! 
+!-----------------------------------------------------------------------
+!
+! !REVISION HISTORY:
+!  same as module
+
+! !INPUT/OUTPUT PARAMETERS:
+
+    type(mct_aVect)   , intent(inout) :: x2o_o
+
+! !OUTPUT PARAMETERS:
+
+   integer, intent(out) :: &amp;
+      errorCode              ! returned error code
+
+!EOP
+!BOC
+!-----------------------------------------------------------------------
+!
+!  local variables
+!
+!-----------------------------------------------------------------------
+
+!   character (len=StrKIND) ::   &amp;
+!      label,                 &amp;
+!      message
+! 
+!   integer ::  &amp;
+!      i,j,k,n,iblock
+!
+!   real (kind=RKIND), dimension(nx_block,ny_block) ::  &amp;
+!      WORKB
+!
+!   real (kind=RKIND), dimension(nx_block,ny_block,max_blocks_clinic) ::   &amp;
+!      WORK1, WORK2        ! local work space
+!
+!   real (r8) ::  &amp;
+!      m2percm2,  &amp;
+!      gsum
+!
+!   type (block) :: this_block ! local block info
+!
+!!-----------------------------------------------------------------------
+!!
+!!  zero out padded cells 
+!!
+!!-----------------------------------------------------------------------
+!
+!   errorCode = POP_Success
+!
+!   WORK1 = c0
+!   WORK2 = c0
+!
+!!-----------------------------------------------------------------------
+!!
+!!  unpack and distribute wind stress, then convert to correct units
+!!  and rotate components to local coordinates
+!!
+!!-----------------------------------------------------------------------
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!      this_block = get_block(blocks_clinic(iblock),iblock)
+!
+!      do j=this_block%jb,this_block%je
+!      do i=this_block%ib,this_block%ie
+!         n = n + 1
+!         WORK1(i,j,iblock) = x2o_o%rAttr(index_x2o_Foxx_taux,n)
+!         WORK2(i,j,iblock) = x2o_o%rAttr(index_x2o_Foxx_tauy,n)
+!      enddo
+!      enddo
+!   enddo ! iblock
+!
+!   !***
+!   !*** do NOT perform halo updates now, because vector updates must
+!   !***   be done after the rotation is completed.
+!   !***
+!
+!!-----------------------------------------------------------------------
+!!
+!!  rotate true zonal/meridional wind stress into local coordinates,
+!!  convert to dyne/cm**2, and shift SMFT to U grid
+!!
+!!  halo updates are performed in subroutine rotate_wind_stress, 
+!!  following the rotation
+!!
+!!-----------------------------------------------------------------------
+!
+!      call rotate_wind_stress(WORK1, WORK2)
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!      this_block = get_block(blocks_clinic(iblock),iblock)
+!
+!!-----------------------------------------------------------------------
+!!
+!!  unpack and distribute fresh water flux and salt flux
+!!
+!!  NOTE: if there are code changes associated with changing the names or
+!!        the number of fluxes received from the coupler, then subroutine
+!!        update_ghost_cells_coupler_fluxes will need to be modified also
+!!
+!!-----------------------------------------------------------------------
+!
+!
+!      do j=this_block%jb,this_block%je
+!      do i=this_block%ib,this_block%ie
+!         n = n + 1
+!         SNOW_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Faxa_snow,n)
+!         WORKB (i,j       ) = x2o_o%rAttr(index_x2o_Faxa_rain,n)
+!         EVAP_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Foxx_evap,n)
+!         MELT_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Fioi_meltw,n)
+!         ROFF_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Forr_roff,n)
+!         IOFF_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Forr_ioff,n)
+!         SALT_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Fioi_salt,n)
+!
+!         PREC_F(i,j,iblock) = WORKB(i,j) + SNOW_F(i,j,iblock)    ! rain + snow
+!
+!         WORKB(i,j        ) = x2o_o%rAttr(index_x2o_Foxx_swnet,n)
+!         SHF_QSW(i,j,iblock) = WORKB(i,j)*  &amp;
+!            RCALCT(i,j,iblock)*hflux_factor  !  convert from W/m**2
+!         SENH_F(i,j,iblock)  = x2o_o%rAttr(index_x2o_Foxx_sen,n)
+!         LWUP_F(i,j,iblock)  = x2o_o%rAttr(index_x2o_Foxx_lwup,n)
+!         LWDN_F(i,j,iblock)  = x2o_o%rAttr(index_x2o_Faxa_lwdn,n)
+!         MELTH_F(i,j,iblock) = x2o_o%rAttr(index_x2o_Fioi_melth,n)
+!
+!         WORKB(i,j       ) = x2o_o%rAttr(index_x2o_Si_ifrac,n)
+!         IFRAC(i,j,iblock) = WORKB(i,j) * RCALCT(i,j,iblock)
+!
+!         !***  converting from Pa to dynes/cm**2
+!         WORKB(i,j       ) = x2o_o%rAttr(index_x2o_Sa_pslv,n)
+!         ATM_PRESS(i,j,iblock) = c10 * WORKB(i,j) * RCALCT(i,j,iblock)
+!
+!         !***  converting from m**2/s**2 to cm**2/s**2
+!         WORKB(i,j       ) = x2o_o%rAttr(index_x2o_So_duu10n,n)
+!         U10_SQR(i,j,iblock) = cmperm * cmperm * WORKB(i,j) * RCALCT(i,j,iblock)
+!
+!      enddo
+!      enddo
+!
+!   enddo
+!
+!!-----------------------------------------------------------------------
+!!
+!!  incoming data quality control
+!!
+!!-----------------------------------------------------------------------
+!#ifdef CCSMCOUPLED
+!      if ( any(IOFF_F &lt; c0) ) then
+!        call shr_sys_abort ('Error: incoming IOFF_F is negative')
+!      endif
+!#endif
+!
+!
+!!-----------------------------------------------------------------------
+!!
+!!  update ghost cells for fluxes received from the coupler
+!!
+!!-----------------------------------------------------------------------
+!
+!   call update_ghost_cells_coupler_fluxes(errorCode)
+!
+!   if (errorCode /= POP_Success) then
+!      call POP_ErrorSet(errorCode, &amp;
+!         'ocn_import_mct: error in update_ghost_cells_coupler_fluxes')
+!      return
+!   endif
+!
+!!-----------------------------------------------------------------------
+!!
+!!  unpack atmospheric CO2
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (index_x2o_Sa_co2prog &gt; 0) then
+!      n = 0
+!      do iblock = 1, nblocks_clinic
+!         this_block = get_block(blocks_clinic(iblock),iblock)
+!
+!         do j=this_block%jb,this_block%je
+!         do i=this_block%ib,this_block%ie
+!            n = n + 1
+!            WORK1(i,j,iblock) = x2o_o%rAttr(index_x2o_Sa_co2prog,n)
+!         enddo
+!         enddo
+!      enddo
+!
+!      call POP_HaloUpdate(WORK1,POP_haloClinic,          &amp;
+!                       POP_gridHorzLocCenter,          &amp;
+!                       POP_fieldKindScalar, errorCode, &amp;
+!                       fillValue = 0.0_POP_r8)
+!
+!      if (errorCode /= POP_Success) then
+!         call POP_ErrorSet(errorCode, &amp;
+!            'ocn_import_mct: error updating PROG CO2 halo')
+!         return
+!      endif
+!
+!      call named_field_set(ATM_CO2_PROG_nf_ind, WORK1)
+!   endif
+!
+!   if (index_x2o_Sa_co2diag &gt; 0) then
+!      n = 0
+!      do iblock = 1, nblocks_clinic
+!         this_block = get_block(blocks_clinic(iblock),iblock)
+!
+!         do j=this_block%jb,this_block%je
+!         do i=this_block%ib,this_block%ie
+!            n = n + 1
+!            WORK1(i,j,iblock) = x2o_o%rAttr(index_x2o_Sa_co2diag,n)
+!         enddo
+!         enddo
+!      enddo
+!
+!      call POP_HaloUpdate(WORK1,POP_haloClinic,          &amp;
+!                       POP_gridHorzLocCenter,          &amp;
+!                       POP_fieldKindScalar, errorCode, &amp;
+!                       fillValue = 0.0_POP_r8)
+!
+!      if (errorCode /= POP_Success) then
+!         call POP_ErrorSet(errorCode, &amp;
+!            'ocn_import_mct: error updating DIAG CO2 halo')
+!         return
+!      endif
+!
+!      call named_field_set(ATM_CO2_DIAG_nf_ind, WORK1)
+!   endif
+! 
+!!-----------------------------------------------------------------------
+!!
+!!  receive orbital parameters
+!!
+!!-----------------------------------------------------------------------
+!
+!   call seq_infodata_GetData(infodata, &amp;
+!      orb_eccen=orb_eccen, orb_mvelpp=orb_mvelpp, orb_lambm0=orb_lambm0, orb_obliqr=orb_obliqr)
+!
+!
+!!-----------------------------------------------------------------------
+!!
+!!  diagnostics
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (ldiag_cpl) then
+!
+!     write(message,'(6a,1x,5a)')  &amp;
+!         ' Global averages of fluxes received from cpl at ',  &amp;
+!           cyear,'/',cmonth ,'/',cday,  chour,':',cminute,':',csecond
+!     call document ('pop_recv_from_coupler', trim(message))
+! 
+!     m2percm2  = mpercm*mpercm
+!     do k = 1,nrecv
+!
+!         n = 0
+!         do iblock = 1, nblocks_clinic
+!            this_block = get_block(blocks_clinic(iblock),iblock)
+!
+!            do j=this_block%jb,this_block%je
+!            do i=this_block%ib,this_block%ie
+!               n = n + 1
+!               WORK1(i,j,iblock) = x2o_o%rAttr(k,n)  ! mult. by TAREA in global_sum_prod
+!            enddo
+!            enddo
+!         enddo
+!
+!         gsum = global_sum_prod(WORK1 , TAREA, distrb_clinic, &amp;
+!                                 field_loc_center, RCALCT)*m2percm2
+!         if (my_task == master_task) then
+!            call seq_flds_getField(label,k,seq_flds_x2o_fields)
+!            write(stdout,1100)'ocn','recv', label ,gsum
+!            call shr_sys_flush(stdout)
+!         endif
+!      enddo
+!   endif
+!
+!
+!1100  format ('comm_diag ', a3, 1x, a4, 1x, a8, 1x, es26.19:, 1x, a6)
+
+!-----------------------------------------------------------------------
+!EOC
+
+ end subroutine ocn_import_mct
+!***********************************************************************
+!BOP
+! !IROUTINE: ocn_export_mct
+! !INTERFACE:
+
+ subroutine ocn_export_mct(o2x_o, errorCode)   
+
+! !DESCRIPTION:
+!  This routine calls the routines necessary to send pop fields to
+!  the CCSM cpl7 driver
+!
+! !REVISION HISTORY:
+!  same as module
+!
+! !INPUT/OUTPUT PARAMETERS:
+
+   type(mct_aVect)   , intent(inout) :: o2x_o
+
+! !OUTPUT PARAMETERS:
+
+   integer, intent(out) :: &amp;
+      errorCode              ! returned error code
+
+!EOP
+!BOC
+!!-----------------------------------------------------------------------
+!!
+!!  local variables
+!!
+!!-----------------------------------------------------------------------
+!
+!   integer (int_kind) :: n, iblock
+!           
+!   character (char_len)    :: label
+! 
+!   integer (int_kind) ::  &amp;
+!      i,j,k
+!
+!   real (r8), dimension(nx_block,ny_block) ::   &amp;
+!      WORK1, WORK2,      &amp;! local work space
+!      WORK3, WORK4
+!
+!   real (r8), dimension(nx_block,ny_block,max_blocks_clinic) ::   &amp;
+!        WORKA               ! local work space with full block dimension
+!
+!   real (r8) ::   &amp;
+!      m2percm2,   &amp;
+!      gsum
+!
+!   type (block) :: this_block ! local block info
+!
+!!-----------------------------------------------------------------------
+!
+!!-----------------------------------------------------------------------
+!!
+!!  initialize control buffer
+!!
+!!-----------------------------------------------------------------------
+!
+!   errorCode = POP_Success
+!
+!!-----------------------------------------------------------------------
+!!
+!!     interpolate onto T-grid points and rotate on T grid
+!!
+!!-----------------------------------------------------------------------
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!     this_block = get_block(blocks_clinic(iblock),iblock)
+!
+!     call ugrid_to_tgrid(WORK3,SBUFF_SUM(:,:,iblock,index_o2x_So_u),iblock)
+!     call ugrid_to_tgrid(WORK4,SBUFF_SUM(:,:,iblock,index_o2x_So_v),iblock)
+!
+!     WORK1 = (WORK3*cos(ANGLET(:,:,iblock))+WORK4*sin(-ANGLET(:,:,iblock)))  &amp;
+!            * mpercm/tlast_coupled
+!     WORK2 = (WORK4*cos(ANGLET(:,:,iblock))-WORK3*sin(-ANGLET(:,:,iblock)))  &amp;
+!            * mpercm/tlast_coupled
+!
+!     do j=this_block%jb,this_block%je
+!     do i=this_block%ib,this_block%ie
+!        n = n + 1
+!        o2x_o%rAttr(index_o2x_So_u,n) = WORK1(i,j)
+!        o2x_o%rAttr(index_o2x_So_v,n) = WORK2(i,j)
+!     enddo
+!     enddo
+!  enddo
+!
+!!-----------------------------------------------------------------------
+!!
+!!     convert and pack surface temperature
+!!
+!!-----------------------------------------------------------------------
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!      this_block = get_block(blocks_clinic(iblock),iblock)
+!      do j=this_block%jb,this_block%je
+!      do i=this_block%ib,this_block%ie
+!         n = n + 1
+!         o2x_o%rAttr(index_o2x_So_t,n) =   &amp;
+!             SBUFF_SUM(i,j,iblock,index_o2x_So_t)/tlast_coupled + T0_Kelvin
+!      enddo
+!      enddo
+!   enddo
+!
+!!-----------------------------------------------------------------------
+!!
+!!     convert and pack salinity
+!!
+!!-----------------------------------------------------------------------
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!      this_block = get_block(blocks_clinic(iblock),iblock)
+!      do j=this_block%jb,this_block%je
+!      do i=this_block%ib,this_block%ie
+!         n = n + 1
+!         o2x_o%rAttr(index_o2x_So_s,n) =   &amp;
+!             SBUFF_SUM(i,j,iblock,index_o2x_So_s)*salt_to_ppt/tlast_coupled
+!      enddo
+!      enddo
+!   enddo
+!
+!!-----------------------------------------------------------------------
+!!
+!!     interpolate onto T-grid points, then rotate on T grid
+!!
+!!-----------------------------------------------------------------------
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!      this_block = get_block(blocks_clinic(iblock),iblock)
+!      call ugrid_to_tgrid(WORK3,SBUFF_SUM(:,:,iblock,index_o2x_So_dhdx),iblock)
+!      call ugrid_to_tgrid(WORK4,SBUFF_SUM(:,:,iblock,index_o2x_So_dhdy),iblock)
+! 
+!      WORK1 = (WORK3*cos(ANGLET(:,:,iblock)) + WORK4*sin(-ANGLET(:,:,iblock)))  &amp;
+!              /grav/tlast_coupled
+!      WORK2 = (WORK4*cos(ANGLET(:,:,iblock)) - WORK3*sin(-ANGLET(:,:,iblock)))  &amp;
+!              /grav/tlast_coupled
+!
+!      do j=this_block%jb,this_block%je
+!      do i=this_block%ib,this_block%ie
+!         n = n + 1
+!         o2x_o%rAttr(index_o2x_So_dhdx,n) = WORK1(i,j)
+!         o2x_o%rAttr(index_o2x_So_dhdy,n) = WORK2(i,j)
+!      enddo
+!      enddo
+!   enddo
+!
+!!-----------------------------------------------------------------------
+!!
+!!     pack heat flux due to freezing/melting (W/m^2)
+!!     QFLUX computation and units conversion occurs in ice.F
+!!
+!!-----------------------------------------------------------------------
+!
+!   n = 0
+!   do iblock = 1, nblocks_clinic
+!      this_block = get_block(blocks_clinic(iblock),iblock)
+!      do j=this_block%jb,this_block%je
+!      do i=this_block%ib,this_block%ie
+!         n = n + 1
+!         o2x_o%rAttr(index_o2x_Fioo_q,n) = QFLUX(i,j,iblock)
+!      enddo
+!      enddo
+!   enddo
+!
+!   tlast_ice = c0
+!   AQICE     = c0
+!   QICE      = c0
+!
+!!-----------------------------------------------------------------------
+!!
+!!     pack co2 flux, if requested (kg CO2/m^2/s)
+!!     units conversion occurs where co2 flux is computed
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (index_o2x_Faoo_fco2_ocn &gt; 0) then
+!      n = 0
+!      do iblock = 1, nblocks_clinic
+!         this_block = get_block(blocks_clinic(iblock),iblock)
+!         do j=this_block%jb,this_block%je
+!         do i=this_block%ib,this_block%ie
+!            n = n + 1
+!            o2x_o%rAttr(index_o2x_Faoo_fco2_ocn,n) = &amp;
+!               SBUFF_SUM(i,j,iblock,index_o2x_Faoo_fco2_ocn)/tlast_coupled
+!         enddo
+!         enddo
+!      enddo
+!   endif
+!
+!!-----------------------------------------------------------------------
+!!
+!!     diagnostics
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (ldiag_cpl) then
+!      call ccsm_char_date_and_time
+!!DEBUG      write(message,'(6a,1x,5a)')' Global averages of fluxes sent to cpl at ', &amp;
+!!DEBUG           cyear,'/',cmonth, '/',cday,  chour,':',cminute,':',csecond
+!!DEBUG      call document ('pop_send_to_coupler', message)
+!      write(stdout,*)'pop_send_to_coupler'
+!
+!      m2percm2  = mpercm*mpercm
+!      do k = 1,nsend
+!        n = 0
+!        do iblock = 1, nblocks_clinic
+!           this_block = get_block(blocks_clinic(iblock),iblock)
+!           do j=this_block%jb,this_block%je
+!           do i=this_block%ib,this_block%ie
+!              n = n + 1
+!              WORKA(i,j,iblock) = o2x_o%rAttr(k,n)
+!           enddo
+!           enddo
+!        enddo
+!
+!        call POP_HaloUpdate(WORKA,POP_haloClinic,          &amp;
+!                       POP_gridHorzLocCenter,          &amp;
+!                       POP_fieldKindScalar, errorCode, &amp;
+!                       fillValue = 0.0_POP_r8)
+!       
+!         if (errorCode /= POP_Success) then
+!            call POP_ErrorSet(errorCode, &amp;
+!               'ocn_export_mct: error updating halo for state')
+!            return
+!         endif
+!
+!        gsum = global_sum_prod(WORKA , TAREA, distrb_clinic, &amp;
+!                                   field_loc_center, RCALCT)*m2percm2
+!        if (my_task == master_task) then
+!           call seq_flds_getField(label,k,seq_flds_o2x_fields)
+!           write(stdout,1100)'ocn','send', label ,gsum
+!        endif
+!      enddo ! k
+!      if (my_task == master_task) call shr_sys_flush(stdout)
+!   endif
+!
+!1100 format ('comm_diag ', a3, 1x, a4, 1x, a8, 1x, es26.19:, 1x, a6)
+!    
+!    tlast_coupled = c0
+
+!-----------------------------------------------------------------------
+!EOC
+
+ end subroutine ocn_export_mct
+
+!***********************************************************************
+
+!BOP
+! !IROUTINE: pop_sum_buffer
+! !INTERFACE:
+
+ subroutine pop_sum_buffer
+
+! !DESCRIPTION:
+!  This routine accumulates sums for averaging fields to
+!  be sent to the coupler
+!
+! !REVISION HISTORY:
+!  same as module
+! 
+!EOP
+!BOC
+
+!#ifdef CCSMCOUPLED
+!!-----------------------------------------------------------------------
+!!
+!!  local variables
+!!
+!!-----------------------------------------------------------------------
+!
+!   real (r8), dimension(nx_block,ny_block,max_blocks_clinic) ::  &amp;
+!      WORK                ! local work arrays
+!
+!   real (r8) ::   &amp;
+!      delt,             &amp; ! time interval since last step
+!      delt_last           ! time interval for previous step
+!
+!   integer (int_kind) :: &amp;
+!      iblock,           &amp; ! block index
+!      sflux_co2_nf_ind = 0! named field index of fco2
+!
+!   logical (log_kind) :: &amp;
+!      first = .true.      ! only true for first call
+!
+!   save first
+!
+!!-----------------------------------------------------------------------
+!!
+!!  zero buffer if this is the first time after a coupling interval
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (tlast_coupled == c0) SBUFF_SUM = c0
+!   WORK = c0
+!
+!!-----------------------------------------------------------------------
+!!
+!!  update time since last coupling
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (avg_ts .or. back_to_back) then
+!      delt = p5*dtt
+!   else
+!      delt =    dtt
+!   endif
+!   tlast_coupled = tlast_coupled + delt
+!
+!!-----------------------------------------------------------------------
+!!
+!!  allow for fco2 field to not be registered on first call
+!!     because init_forcing is called before init_passive_tracers
+!!  use weight from previous timestep because flux used here is that
+!!     computed during the previous timestep
+!!
+!!-----------------------------------------------------------------------
+!
+!   if (index_o2x_Faoo_fco2_ocn &gt; 0) then
+!      if (sflux_co2_nf_ind == 0) then
+!         call named_field_get_index('SFLUX_CO2', sflux_co2_nf_ind, &amp;
+!                                    exit_on_err=.not. first)
+!      endif
+!
+!      if (avg_ts .or. back_to_back) then
+!         delt_last = p5*dtt
+!      else
+!         delt_last =    dtt
+!      endif
+!   endif
+!
+!!-----------------------------------------------------------------------
+!!
+!!  accumulate sums of U,V,T,S and GRADP
+!!  accumulate sum of co2 flux, if requested
+!!     implicitly use zero flux if fco2 field not registered yet
+!!  ice formation flux is handled separately in ice routine
+!!
+!!-----------------------------------------------------------------------
+!
+!   !$OMP PARALLEL DO PRIVATE(iblock)
+!   do iblock = 1, nblocks_clinic
+!   SBUFF_SUM(:,:,iblock,index_o2x_So_u) =   &amp;
+!      SBUFF_SUM(:,:,iblock,index_o2x_So_u) + delt*  &amp;
+!                                   UVEL(:,:,1,curtime,iblock)
+!
+!   SBUFF_SUM(:,:,iblock,index_o2x_So_v) =   &amp;
+!      SBUFF_SUM(:,:,iblock,index_o2x_So_v) + delt*  &amp;
+!                                   VVEL(:,:,1,curtime,iblock)
+!
+!   SBUFF_SUM(:,:,iblock,index_o2x_So_t ) =   &amp;
+!      SBUFF_SUM(:,:,iblock,index_o2x_So_t ) + delt*  &amp;
+!                                   TRACER(:,:,1,1,curtime,iblock)
+!
+!   SBUFF_SUM(:,:,iblock,index_o2x_So_s ) =   &amp;
+!      SBUFF_SUM(:,:,iblock,index_o2x_So_s ) + delt*  &amp;
+!                                   TRACER(:,:,1,2,curtime,iblock)
+!
+!   SBUFF_SUM(:,:,iblock,index_o2x_So_dhdx) =   &amp;
+!      SBUFF_SUM(:,:,iblock,index_o2x_So_dhdx) + delt*  &amp;
+!                                   GRADPX(:,:,curtime,iblock)
+!
+!   SBUFF_SUM(:,:,iblock,index_o2x_So_dhdy) =   &amp;
+!      SBUFF_SUM(:,:,iblock,index_o2x_So_dhdy) + delt*  &amp;
+!                                   GRADPY(:,:,curtime,iblock)
+!
+!   if (index_o2x_Faoo_fco2_ocn &gt; 0 .and. sflux_co2_nf_ind &gt; 0) then
+!      call named_field_get(sflux_co2_nf_ind, iblock, WORK(:,:,iblock))
+!      SBUFF_SUM(:,:,iblock,index_o2x_Faoo_fco2_ocn) = &amp;
+!         SBUFF_SUM(:,:,iblock,index_o2x_Faoo_fco2_ocn) + delt_last*WORK(:,:,iblock)
+!   endif
+!
+!   enddo
+!   !$OMP END PARALLEL DO
+!
+!   first = .false.
+!
+!#endif
+!
+!-----------------------------------------------------------------------
+!EOC
+
+ end subroutine pop_sum_buffer
+ 
+!***********************************************************************
+
+end module ocn_comp_mct
+
+!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

</font>
</pre>