integer bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)

     &      + g_objf_state_final(i,j,bi,bj,k+0*Nr)**2/stdDevtheta**2

     &      + g_objf_state_final(i,j,bi,bj,k+1*Nr)**2/stdDevsalt**2

     &      + g_objf_state_final(i,j,bi,bj,k+2*Nr)**2/stdDevuvel**2

     &      + g_objf_state_final(i,j,bi,bj,k+3*Nr)**2/stdDevvvel**2

     &      + g_objf_state_final(i,j,bi,bj,1+4*Nr)**2/stdDevetan**2

      DO bj=myByLo(myThid),myByHi(myThid)

      do bj = jtlo, jthi

                tmpfld3d(i,j,k,bi,bj) = tmpfld3d(i,j,k,bi,bj)

     &                + advvel(i,j,k,bi,bj)

      do bj = jtlo, jthi

                tmpfld3d(i,j,k,bi,bj) = tmpfld3d(i,j,k,bi,bj)

     &                + aduvel(i,j,k,bi,bj)

      do bj = jtlo, jthi

                tmpfld3d(i,j,k,bi,bj) = tmpfld3d(i,j,k,bi,bj)

     &                + adsalt(i,j,k,bi,bj)*fac

      do bj = jtlo, jthi

                tmpfld3d(i,j,k,bi,bj) = tmpfld3d(i,j,k,bi,bj)

     &                + adtheta(i,j,k,bi,bj)*fac

      integer bi, bj

      DO bj = jtlo, jthi

            adtheta(i,j,k,bi,bj) = 

     &       g_objf_state_final(i,j,bi,bj,k)

            adsalt(i,j,k,bi,bj) = 

     &       g_objf_state_final(i,j,bi,bj,1*Nr+k)

            aduvel(i,j,k,bi,bj) = 

     &       g_objf_state_final(i,j,bi,bj,2*Nr+k)

            advvel(i,j,k,bi,bj) = 

     &       g_objf_state_final(i,j,bi,bj,3*Nr+k)

          adetan(i,j,bi,bj) = 

     &       g_objf_state_final(i,j,bi,bj,4*Nr+1)

      do bj = jtlo, jthi

              tmpfld2d(i,j,bi,bj) = tmpfld2d(i,j,bi,bj)

     &              + adetan(i,j,bi,bj)

      integer bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)

            adobjf_state_final(i,j,bi,bj,k) = 

     &       g_objf_state_final(i,j,bi,bj,k)

           adobjf_state_final(i,j,bi,bj,4*Nr+1) = 

     &      g_objf_state_final(i,j,bi,bj,4*Nr+1)

         ReCountX(j,bj) = CountX(j,bj)

         IF (yC(1,j,bi,bj) .LE. -62.0) THEN

            ReCountX(j,bj) = 90.

         ELSE IF (yC(1,j,bi,bj) .EQ. 74.0) THEN

            ReCountX(j,bj) = CountX(j,bj) + 9.0

         ELSE IF (yC(1,j,bi,bj) .EQ. 70.0) THEN

            ReCountX(j,bj) = CountX(j,bj) + 8.0

         ELSE IF (yC(1,j,bi,bj) .EQ. 66.0) THEN

            ReCountX(j,bj) = CountX(j,bj) + 5.0

         ELSE IF (yC(1,j,bi,bj) .EQ. 62.0) THEN

            ReCountX(j,bj) = CountX(j,bj) + 1.0

          S(i,j,bj) = sin(yC(i,j,bi,bj)*deg2rad)

          P2(i,j,bj) = 0.5*(3.*S(i,j,bj)**2 - 1.)

          P4(i,j,bj) = 0.12*(35.*S(i,j,bj)**4 - 30.*S(i,j,bj)**2 + 3.)

          IF (yC(i,j,bi,bj) .LT. 0.) THEN

          SW(j,bj) = 0.25*Q0*(1 + Q2*P2(i,j,bj))*

     <         (1 - A0 - A2*P2(i,j,bj) - A4*P4(i,j,bj) )

          T(j,bj) = T0(no_so) + T2(no_so)*P2(i,j,bj)

          LW(j,bj) = LW0 + LW1*(T(j,bj)-t_mlt)

             LW(j,bj) = LW(j,bj) - 

         Qnet(i,j,bi,bj) = -1.0*( SW(j,bj) - LW(j,bj) - 

     <        0.000728e4      - 0.00678e4*S(i,j,bj) + 

     <        0.0955e4*S(i,j,bj)**2 + 0.0769e4*S(i,j,bj)**3 - 

     <        0.8508e4*S(i,j,bj)**4 - 0.3581e4*S(i,j,bj)**5 + 

     <        2.9240e4*S(i,j,bj)**6 + 0.8311e4*S(i,j,bj)**7 -   

     <        4.9548e4*S(i,j,bj)**8 - 0.8808e4*S(i,j,bj)**9 + 

     <        4.0644e4*S(i,j,bj)**10 +0.3409e4*S(i,j,bj)**11 - 

     <        1.2893e4*S(i,j,bj)**12 )

          IF (ReCountX(j,bj) .GT. 0.) THEN

             Qnet(i,j,bi,bj) = (-90./ReCountX(j,bj))*

     <            ( SW(j,bj) - LW(j,bj) - 

     <            Hd35(no_so)*(3.054e1 - 3.763e1*S(i,j,bj) + 

     <        1.892e2*S(i,j,bj)**2 + 3.041e2*S(i,j,bj)**3 - 

     <        1.540e3*S(i,j,bj)**4 - 9.586e2*S(i,j,bj)**5 + 

     <        2.939e3*S(i,j,bj)**6 + 1.219e3*S(i,j,bj)**7 -   

     <        2.550e3*S(i,j,bj)**8 - 5.396e2*S(i,j,bj)**9 + 

     <        8.119e2*S(i,j,bj)**10)

             Qnet(i,j,bi,bj) = 0.

          EmPmR(i,j,bi,bj) = -1.e-3*Fw35(no_so)

     <    *(-0.8454e5*S(i,j,bj)**14 + 0.5367e5*S(i,j,bj)**13 

     <    +3.3173e5*S(i,j,bj)**12 - 1.8965e5*S(i,j,bj)**11 

     <    -5.1701e5*S(i,j,bj)**10

     <    +2.6240e5*S(i,j,bj)**9 + 4.077e5*S(i,j,bj)**8 

     <    -1.791e5*S(i,j,bj)**7

     <    -1.7231e5*S(i,j,bj)**6 + 0.6229e5*S(i,j,bj)**5 

     <    +0.3824e5*S(i,j,bj)**4

     <    -0.1017e5*S(i,j,bj)**3 - 0.0387e5*S(i,j,bj)**2 

     <    +0.00562e5*S(i,j,bj)  + 0.0007743e5)

          IF (yC(i,j,bi,bj) .LT. -40.) THEN

           EmPmR(i,j,bi,bj) = -1.e-3*(Fw35(no_so)*

     <            (-6.5 + 35.3 + 71.7*S(i,j,bj)  

     <           - 1336.3*S(i,j,bj)**2 - 425.8*S(i,j,bj)**3 

     <           + 5434.8*S(i,j,bj)**4 + 707.9*S(i,j,bj)**5 

     <           - 6987.7*S(i,j,bj)**6 - 360.4*S(i,j,bj)**7 

     <           + 2855.0*S(i,j,bj)**8)

           IF (xC(i,j,bi,bj) .GT. 284. 

     <      .OR. xC(i,j,bi,bj) .LT. 28.) THEN

              EmPmR(i,j,bi,bj) = -1.e-3*(Fw35(no_so)*

     <             (-6.5 -2.878 + 3.157e2*S(i,j,bj) -

     <             2.388e3*S(i,j,bj)**2 - 4.101e3*S(i,j,bj)**3 + 

     <             1.963e4*S(i,j,bj)**4 + 1.534e4*S(i,j,bj)**5 - 

     <             6.556e4*S(i,j,bj)**6 - 2.478e4*S(i,j,bj)**7 + 

     <             1.083e5*S(i,j,bj)**8 + 1.85e4*S(i,j,bj)**9 - 

     <             8.703e4*S(i,j,bj)**10 - 5.276e3*S(i,j,bj)**11 + 

     <             2.703e4*S(i,j,bj)**12)

              EmPmR(i,j,bi,bj) = -1.e-3*(Fw35(no_so)

     <             *(-6.5 +51.89 + 4.916e2*S(i,j,bj) - 

     <             1.041e3*S(i,j,bj)**2 - 7.546e3*S(i,j,bj)**3 + 

     <             2.335e3*S(i,j,bj)**4 + 3.449e4*S(i,j,bj)**5 + 

     <             6.702e3*S(i,j,bj)**6 - 6.601e4*S(i,j,bj)**7 - 

     <             2.594e4*S(i,j,bj)**8 + 5.652e4*S(i,j,bj)**9 + 

     <             2.738e4*S(i,j,bj)**10 - 1.795e4*S(i,j,bj)**11 - 

     <             9.486e3*S(i,j,bj)**12)

      INTEGER bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)

          act2 = bj - myByLo(myThid)

          S(i,j,bj) = 0.0

          P2(i,j,bj) = 0.0

          P4(i,j,bj) = 0.0

         SW(j,bj) = 0.0

         LW(j,bj) = 0.0

         Hd(j,bj) = 0.0

         Fw(j,bj) = 0.0

         T(j,bj) = 0.0

         ReCountX(j,bj) = 0.0

     I             bi, bj, iMin, iMax, jMin, jMax,

      INTEGER bi,bj

          surfaceForcingU(i,j,bi,bj) = 

     &      fu(i,j,bi,bj)*horiVertRatio*recip_rhoConst

     &           + winPert(i,j,bi,bj)

     &            *drF(kSurface)*hFacW(i,j,kSurface,bi,bj)

          surfaceForcingV(i,j,bi,bj) = 

     &      fv(i,j,bi,bj)*horiVertRatio*recip_rhoConst

          surfaceForcingT(i,j,bi,bj) = 

     &      -Qnet(i,j,bi,bj)*recip_Cp*horiVertRatio*recip_rhoConst

     &      -lambdaThetaZonRelax*maskC(i,j,kSurface,bi,bj)*

     &         (theta(i,j,kSurface,bi,bj)-ZonalMeanSST(j,bj))

     &        *drF(kSurface)*hFacC(i,j,kSurface,bi,bj)

          surfaceForcingS(i,j,bi,bj) = 

     &      EmPmR(i,j,bi,bj)*convertFW2Salt*convertEmP2rUnit

     &      +Run(i,j,bi,bj)*scale_runoff

     &      -lambdaSaltClimRelax(i,j,bi,bj)*maskC(i,j,kSurface,bi,bj)

     &        *(salt(i,j,kSurface,bi,bj)-SSS(i,j,bi,bj))

     &        *drF(kSurface)*hFacC(i,j,kSurface,bi,bj)

      INTEGER i, j, k, bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)

         ZonalMeanSST(j,bj) = 0.0

         CountX(j,bj) = 0.0

      DO bj=myByLo(myThid),myByHi(myThid)

         CountX_tile(j,bi,bj) = 0.0

         ZonalMean_tile(j,bi,bj) = 0.0

          ZonalMean_tile(j,bi,bj) = ZonalMean_tile(j,bi,bj) +

     &           theta(i,j,k,bi,bj)

          CountX_tile(j,bi,bj) = CountX_tile(j,bi,bj) +  

     &         maskC(i,j,k,bi,bj)

         ZonalMeanSST(j,bj) = ZonalMeanSST(j,bj) + 

     &    ZonalMean_tile(j,bi,bj)

         CountX(j,bj) = CountX(j,bj) + CountX_tile(j,bi,bj)

      DO bj=myByLo(myThid),myByHi(myThid)

        _GLOBAL_SUM_R8( CountX(j,bj), myThid )

        _GLOBAL_SUM_R8( ZonalMeanSST(j,bj), myThid )

      DO bj=myByLo(myThid),myByHi(myThid)

        IF ( CountX(j,bj) .GT. 0.0) THEN

          ZonalMeanSST(j,bj) = ZonalMeanSST(j,bj)/CountX(j,bj)

         winPert(i,j,bi,bj) = maskW(i,j,k,bi,bj)*

     &        (1.0/1.66)*(0.75*winPert(i,j,bi,bj) + 

      INTEGER i, j, bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)

      INTEGER bi,bj,i,j

       DO bj = myByLo(myThid), myByHi(myThid)

           distY = (yC(i,j,bi,bj)-(yC0))/lY

           fu(i,j,bi,bj) = tauX

       DO bj = myByLo(myThid), myByHi(myThid)

           fv(i,j,bi,bj) = 0.0

      INTEGER bi,bj,i,j,k

      DO bj = myByLo(myThid), myByHi(myThid)

          Run   (i,j,bi,bj) = 0.

          winPert(i,j,bi,bj) = 0.

      DO bj = myByLo(myThid), myByHi(myThid)

     I       bi, bj, iMin, iMax, jMin, jMax,

# define _recip_hFacS(i,j,k,bi,bj) recip_hFacS(i,j,k,bi,bj)*maskS(i,j,k,bi,bj)

          ELSEIF ( yC(i,j,bi,bj) .GE. lat(2) .AND. 

     &       yC(i,j,bi,bj) .LT. lat(3)) THEN

             TmlN_tile(bi,bj) = TmlN_tile(bi,bj) + 

     &           theta(i,j,k,bi,bj)

            CountmlN_tile(bi,bj) = CountmlN_tile(bi,bj) + 

     &           maskC(i,j,k,bi,bj)

        TmlS = TmlS + TmlS_tile(bi,bj)

        CountmlS = CountmlS + CountmlS_tile(bi,bj)

        TmlN = TmlN + TmlN_tile(bi,bj)

        CountmlN = CountmlN + CountmlN_tile(bi,bj)

        TtS = TtS + TtS_tile(bi,bj)

        CounttS = CounttS + CounttS_tile(bi,bj)

        TtN = TtN + TtN_tile(bi,bj)

        CounttN = CounttN + CounttN_tile(bi,bj)

      INTEGER i, j, k, bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)

          TmlS_tile(bi,bj) = 0.

          TmlN_tile(bi,bj) = 0.

          TtS_tile(bi,bj) = 0.

          TtN_tile(bi,bj) = 0.

          CountmlS_tile(bi,bj) = 0.

          CountmlN_tile(bi,bj) = 0.

          CounttS_tile(bi,bj) = 0.

          CounttN_tile(bi,bj) = 0.

          IF ( yC(i,j,bi,bj) .GE. -lat(3) .AND. 

     &         yC(i,j,bi,bj) .LE. -lat(2)) THEN

            TmlS_tile(bi,bj) = TmlS_tile(bi,bj) + 

     &           theta(i,j,k,bi,bj)

            CountmlS_tile(bi,bj) = CountmlS_tile(bi,bj) + 

     &           maskC(i,j,k,bi,bj)

          ELSEIF ( yC(i,j,bi,bj) .GT. -lat(2) .AND. 

     &       yC(i,j,bi,bj) .LT. -lat(1)) THEN 

            TtS_tile(bi,bj) = TtS_tile(bi,bj) + 

     &           theta(i,j,k,bi,bj)

            CounttS_tile(bi,bj) = CounttS_tile(bi,bj) + 

     &           maskC(i,j,k,bi,bj)

          ELSEIF ( yC(i,j,bi,bj) .GT. lat(1) .AND. 

     &       yC(i,j,bi,bj) .LT. lat(2)) THEN

            TtN_tile(bi,bj) = TtN_tile(bi,bj) + 

     &           theta(i,j,k,bi,bj)

            CounttN_tile(bi,bj) = CounttN_tile(bi,bj) + 

     &           maskC(i,j,k,bi,bj)

          if ( ABS(uwind(i,j,bi,bj)) .GT. 100. .AND.

     &         maskW(i,j,1,bi,bj) .NE. 0. ) then

     &          bi, bj, i, j, myiter, uwind(i,j,bi,bj)

          if ( ABS(vwind(i,j,bi,bj)) .GT. 100. .AND.

     &         maskS(i,j,1,bi,bj) .NE. 0. ) then

     &          bi, bj, i, j, myiter, vwind(i,j,bi,bj)

          if ( ( wspeed(i,j,bi,bj) .LT. 0. .OR.

     &         wspeed(i,j,bi,bj) .GT. 100. ) .AND.

     &         maskS(i,j,1,bi,bj) .NE. 0. ) then

     &          bi, bj, i, j, myiter, wspeed(i,j,bi,bj)

          if ( (atemp(i,j,bi,bj) .LT. 183 .OR.

     &         atemp(i,j,bi,bj) .GT. 343 ) .AND.

     &