<p><b>dwj07@fsu.edu</b> 2012-09-11 15:04:29 -0600 (Tue, 11 Sep 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Setting up CESM driver make system.<br>
        It's probably not working yet, but this is a good place to start.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/cesm_coupling/src/Makefile
===================================================================
--- branches/ocean_projects/cesm_coupling/src/Makefile        2012-09-11 20:55:28 UTC (rev 2150)
+++ branches/ocean_projects/cesm_coupling/src/Makefile        2012-09-11 21:04:29 UTC (rev 2151)
@@ -37,7 +37,7 @@
         ar ru lib$(CORE).a framework/*.o
         ar ru lib$(CORE).a operators/*.o
         ar ru lib$(CORE).a core_$(CORE)/*.o
-        ar ru lib$(CORE).a cesm_driver/*.o
+        ar ru lib$(CORE).a $(CORE)_cesm_driver/*.o
 else
         $(LINKER) $(LDFLAGS) -o $(CORE)_model.exe driver/*.o -L. -ldycore -lops -lframework $(LIBS) -I./external/esmf_time_f90 -L./external/esmf_time_f90 -lesmf_time
 endif
@@ -63,7 +63,7 @@
 
 ifeq &quot;$(CESM)&quot; &quot;true&quot;
 drver:  reg_includes externals frame ops dycore
-        ( cd cesm_driver; $(MAKE) CPPFLAGS=&quot;$(CPPFLAGS)&quot; CPPINCLUDES=&quot;$(CPPINCLUDES)&quot; all ) 
+        ( cd $(CORE)_cesm_driver; $(MAKE) CPPFLAGS=&quot;$(CPPFLAGS)&quot; CPPINCLUDES=&quot;$(CPPINCLUDES)&quot; all ) 
 else
 drver:  reg_includes externals frame ops dycore
         ( cd driver; $(MAKE) CPPFLAGS=&quot;$(CPPFLAGS)&quot; CPPINCLUDES=&quot;$(CPPINCLUDES)&quot; all ) 

Added: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/Makefile
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/Makefile                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/Makefile        2012-09-11 21:04:29 UTC (rev 2151)
@@ -0,0 +1,16 @@
+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 .
+
+cpl_mct:
+        ( cd cpl_mct; $(MAKE) all )
+
+cpl_share:
+        (cd cpl_share; $(MAKE) all )
+
+clean:
+        $(RM) *.o libdrver.a
+        ( cd cpl_mct; $(MAKE) clean )
+        ( cd cpl_share; $(MAKE) clean )

Added: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/Makefile
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/Makefile                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/Makefile        2012-09-11 21:04:29 UTC (rev 2151)
@@ -0,0 +1,18 @@
+.SUFFIXES: .F .o
+
+OBJS = ocn_comp_mct.o \
+       ocn_communicator.o
+
+all: $(OBJS)
+
+ocn_comp_mct.o:
+
+ocn_communicator.o:
+
+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/cpl_mct/ocn_communicator.F (from rev 2150, branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F90)
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F        2012-09-11 21:04:29 UTC (rev 2151)
@@ -0,0 +1,5 @@
+module ocn_communicator
+
+  integer, public :: mpi_communicator_ocn
+
+end module ocn_communicator

Deleted: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F90
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F90        2012-09-11 20:55:28 UTC (rev 2150)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_communicator.F90        2012-09-11 21:04:29 UTC (rev 2151)
@@ -1,5 +0,0 @@
-module ocn_communicator
-
-  integer, public :: mpi_communicator_ocn
-
-end module ocn_communicator

Copied: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F (from rev 2150, branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F90)
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F        2012-09-11 21:04:29 UTC (rev 2151)
@@ -0,0 +1,1633 @@
+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 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 (log_kind) ::   &amp;
+       ldiag_cpl = .false.
+
+  integer (int_kind), 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 (r8),   &amp;
+      dimension(:,:,:,:), allocatable ::  &amp;
+      SBUFF_SUM           ! accumulated sum of send buffer quantities
+                          ! for averaging before being sent
+   real (r8) ::  &amp;
+     tlast_coupled
+
+   integer (int_kind)  ::   &amp;
+      nsend, nrecv
+
+   character(char_len) :: &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(int_kind) ::  &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 (POP_i4) :: &amp;
+       errorCode         ! error code
+
+    integer (int_kind) :: &amp;
+       nThreads
+
+    real (r8) ::  &amp;
+       precadj
+
+    integer (int_kind) :: iam,ierr 
+    character(len=32)  :: 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
+       write(stdout,*) 'ocn_comp_mct ERROR: unknown starttype'
+       call exit_POP(sigAbort,' 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
+         if(master_task == my_task)   then
+            call document ('ocn_init_mct', 'iyear0     ', iyear0)
+            call document ('ocn_init_mct', 'start_year ', start_year)
+         endif
+         call exit_POP(sigAbort,' iyear0 does not match start_year')
+      end if
+      if (imonth0 /= start_month) then
+         if(master_task == my_task)   then
+            call document ('ocn_init_mct', 'imonth0     ', imonth0)
+            call document ('ocn_init_mct', 'start_month ', start_month)
+         endif
+         call exit_POP(sigAbort,' imonth0 does not match start_year')
+      end if
+      if (iday0 /= start_day) then
+         if(master_task == my_task)   then
+            call document ('ocn_init_mct', 'iday0     ', iday0)
+            call document ('ocn_init_mct', 'start_day ', start_day)
+         endif
+      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,*) ' '
+     call document ('ocn_import_mct', 'orb_eccen   ',  orb_eccen)
+     call document ('ocn_import_mct', 'orb_mvelpp  ',  orb_mvelpp)
+     call document ('ocn_import_mct', 'orb_lambm0  ',  orb_lambm0)
+     call document ('ocn_import_mct', 'orb_obliqr  ',  orb_obliqr)
+    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:
+    use POP_FinalMod
+! !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 (POP_i4) :: &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 )
+
+! !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 (int_kind) ::   &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
+    gsize = n
+    call mpas_dmpar_sum_int(dminfo, 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 (int_kind) :: i,j, k, n, ier
+
+    type (block_type) :: 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 (POP_i4), intent(out) :: &amp;
+      errorCode              ! returned error code
+
+!EOP
+!BOC
+!-----------------------------------------------------------------------
+!
+!  local variables
+!
+!-----------------------------------------------------------------------
+
+   character (char_len) ::   &amp;
+      label,                 &amp;
+      message
+ 
+   integer (int_kind) ::  &amp;
+      i,j,k,n,iblock
+
+   real (r8), dimension(nx_block,ny_block) ::  &amp;
+      WORKB
+
+   real (r8), 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 (POP_i4), 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
+
+!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Deleted: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F90
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F90        2012-09-11 20:55:28 UTC (rev 2150)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_mct/ocn_comp_mct.F90        2012-09-11 21:04:29 UTC (rev 2151)
@@ -1,1633 +0,0 @@
-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 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 (log_kind) ::   &amp;
-       ldiag_cpl = .false.
-
-  integer (int_kind), 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 (r8),   &amp;
-      dimension(:,:,:,:), allocatable ::  &amp;
-      SBUFF_SUM           ! accumulated sum of send buffer quantities
-                          ! for averaging before being sent
-   real (r8) ::  &amp;
-     tlast_coupled
-
-   integer (int_kind)  ::   &amp;
-      nsend, nrecv
-
-   character(char_len) :: &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(int_kind) ::  &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 (POP_i4) :: &amp;
-       errorCode         ! error code
-
-    integer (int_kind) :: &amp;
-       nThreads
-
-    real (r8) ::  &amp;
-       precadj
-
-    integer (int_kind) :: iam,ierr 
-    character(len=32)  :: 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
-       write(stdout,*) 'ocn_comp_mct ERROR: unknown starttype'
-       call exit_POP(sigAbort,' 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
-         if(master_task == my_task)   then
-            call document ('ocn_init_mct', 'iyear0     ', iyear0)
-            call document ('ocn_init_mct', 'start_year ', start_year)
-         endif
-         call exit_POP(sigAbort,' iyear0 does not match start_year')
-      end if
-      if (imonth0 /= start_month) then
-         if(master_task == my_task)   then
-            call document ('ocn_init_mct', 'imonth0     ', imonth0)
-            call document ('ocn_init_mct', 'start_month ', start_month)
-         endif
-         call exit_POP(sigAbort,' imonth0 does not match start_year')
-      end if
-      if (iday0 /= start_day) then
-         if(master_task == my_task)   then
-            call document ('ocn_init_mct', 'iday0     ', iday0)
-            call document ('ocn_init_mct', 'start_day ', start_day)
-         endif
-      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,*) ' '
-     call document ('ocn_import_mct', 'orb_eccen   ',  orb_eccen)
-     call document ('ocn_import_mct', 'orb_mvelpp  ',  orb_mvelpp)
-     call document ('ocn_import_mct', 'orb_lambm0  ',  orb_lambm0)
-     call document ('ocn_import_mct', 'orb_obliqr  ',  orb_obliqr)
-    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:
-    use POP_FinalMod
-! !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 (POP_i4) :: &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 )
-
-! !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 (int_kind) ::   &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
-    gsize = n
-    call mpas_dmpar_sum_int(dminfo, 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 (int_kind) :: i,j, k, n, ier
-
-    type (block_type) :: 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 (POP_i4), intent(out) :: &amp;
-      errorCode              ! returned error code
-
-!EOP
-!BOC
-!-----------------------------------------------------------------------
-!
-!  local variables
-!
-!-----------------------------------------------------------------------
-
-   character (char_len) ::   &amp;
-      label,                 &amp;
-      message
- 
-   integer (int_kind) ::  &amp;
-      i,j,k,n,iblock
-
-   real (r8), dimension(nx_block,ny_block) ::  &amp;
-      WORKB
-
-   real (r8), 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 (POP_i4), 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
-
-!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Copied: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F (from rev 2150, branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F90)
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F        2012-09-11 21:04:29 UTC (rev 2151)
@@ -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

Deleted: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F90
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F90        2012-09-11 20:55:28 UTC (rev 2150)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/MPASO_CplIndices.F90        2012-09-11 21:04:29 UTC (rev 2151)
@@ -1,125 +0,0 @@
-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

Added: branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/Makefile
===================================================================
--- branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/Makefile                                (rev 0)
+++ branches/ocean_projects/cesm_coupling/src/ocean_cesm_driver/cpl_share/Makefile        2012-09-11 21:04:29 UTC (rev 2151)
@@ -0,0 +1,15 @@
+.SUFFIXES: .F .o
+
+OBJS = MPASO_CplIndices.o
+
+all: $(OBJS)
+
+MPASO_CplIndices.o:
+
+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

</font>
</pre>