do isect = 1, nsect

     &        'cost_trans_mer',isect,'.',optimcycle

     &     'ECCO_TRANS_MER_YLAT: section at lat: ', ylat(isect)

      mvsum(isect) = 0.0

        tv(isect) = 0.0

        sv(isect) = 0.0

        mv(isect) = 0.0

        mvmin(isect) = 0.0

        mvmax(isect) = 0.0

        kmin(isect) = 0

        kmax(isect) = 0

                   if (  yc(i,j,bi,bj) .ge. ylat(isect) .and.

     $                   yc(i,j,bi,bj) .lt. ylat(isect)+del_y .and.

     $                   xc(i,j,bi,bj) .ge. beglon(isect) .and.

     $                   xc(i,j,bi,bj) .le. endlon(isect) ) then

                       tv(isect) = tv(isect) + p5*(tbar(i,j,k,bi,bj)

                       sv(isect) = sv(isect) + p5*(sbar(i,j,k,bi,bj)

                       mv(isect) = mv(isect) + p5*(hFacC(i,j,k,bi,bj)

              mvlev(isect,k) = mv(isect)

        _GLOBAL_SUM_RL( tv(isect), mythid )

        _GLOBAL_SUM_RL( sv(isect), mythid )

        _GLOBAL_SUM_RL( mv(isect), mythid )

           _GLOBAL_SUM_RL( mvlev(isect,k), mythid )

        mvmin(isect) = mvlev(isect,1)

        mvmax(isect) = mvlev(isect,1)

           if ( mvlev(isect,k) .GT. mvlev(isect,k-1) ) then

              mvmax(isect) = mvlev(isect,k)

              kmax(isect) = k

           if ( mvlev(isect,k) .LT. mvlev(isect,k-1) ) then

              mvmin(isect) = mvlev(isect,k)

              kmin(isect) = k

        mvsum(isect) = mvsum(isect) + mv(isect)

     &       'ECCO_TRANS_MER ', isect, irec, kmin(isect), kmax(isect),

     &       tv(isect), sv(isect), mv(isect), mvmin(isect), mvmax(isect)

     &     'ECCO_TRANS_MER_SUM  mvsum ', isect, mvsum(isect)

     &       isect, irec, kmin(isect), kmax(isect),

     &       tv(isect), sv(isect), mv(isect), mvmin(isect), mvmax(isect)

      integer isect

      do isect = 1, nsect

     &        'cost_trans_zon',isect,'.',optimcycle

     &     'ECCO_TRANS_ZON_XLON: section at lon: ', xlon(isect)

      musum(isect) = 0.0

        tu(isect) = 0.0

        su(isect) = 0.0

        mu(isect) = 0.0

        mumin(isect) = 0.0

        mumax(isect) = 0.0

        kmin(isect) = 0

        kmax(isect) = 0

                   if (  xc(i,j,bi,bj) .ge. xlon(isect) .and.

     $                   xc(i,j,bi,bj) .lt. xlon(isect)+del_x .and.

     $                   yc(i,j,bi,bj) .ge. beglat(isect) .and.

     $                   yc(i,j,bi,bj) .le. endlat(isect) ) then

                       tu(isect) = tu(isect) + p5*(tbar(i,j,k,bi,bj)

                       su(isect) = su(isect) + p5*(sbar(i,j,k,bi,bj)

                       mu(isect) = mu(isect) + p5*(hFacC(i,j,k,bi,bj)

              mulev(isect,k) = mu(isect)

        _GLOBAL_SUM_RL( tu(isect), mythid )

        _GLOBAL_SUM_RL( su(isect), mythid )

        _GLOBAL_SUM_RL( mu(isect), mythid )

           _GLOBAL_SUM_RL( mulev(isect,k), mythid )

        mumin(isect) = mulev(isect,1)

        mumax(isect) = mulev(isect,1)

           if ( mulev(isect,k) .GT. mulev(isect,k-1) ) then

              mumax(isect) = mulev(isect,k)

              kmax(isect) = k

           if ( mulev(isect,k) .LT. mulev(isect,k-1) ) then

              mumin(isect) = mulev(isect,k)

              kmin(isect) = k

        musum(isect) = musum(isect) + mu(isect)

     &       'ECCO_TRANS_ZON ', isect, irec, kmin(isect), kmax(isect),

     &       tu(isect), su(isect), mu(isect), mumin(isect), mumax(isect)

     &     'ECCO_TRANS_ZON_SUM  musum ', isect, musum(isect)

     &       isect, irec, kmin(isect), kmax(isect),

     &       tu(isect), su(isect), mu(isect), mumin(isect), mumax(isect)

      integer isect