I             maskFld3d, boundsVec, myThid

      _RL boundsVec(5)

      x0=boundsVec(1)

      x1=boundsVec(2)

      x2=boundsVec(3)

      x3=boundsVec(4)

      _RL boundsVec(5)

      x0=boundsVec(1)

      x1=boundsVec(2)

      x2=boundsVec(3)

      x3=boundsVec(4)

     I             maskFld3d,boundsVec,myThid

     I             boundsVec,myThid)

      _RL boundsVec(5)

      if (boundsVec(1).LT.boundsVec(4)) then

            IF (fieldCur(i,j,k,bi,bj).gt.boundsVec(4)) then

             fieldCur(i,j,k,bi,bj)=boundsVec(4)

            IF (fieldCur(i,j,k,bi,bj).lt.boundsVec(1)) then

             fieldCur(i,j,k,bi,bj)=boundsVec(1)

     I             boundsVec,myThid)

      _RL boundsVec(5)

      if (boundsVec(1).LT.boundsVec(4)) then

            IF (fieldCur(i,j,bi,bj).gt.boundsVec(4)) then

             fieldCur(i,j,bi,bj)=boundsVec(4)

            IF (fieldCur(i,j,bi,bj).lt.boundsVec(1)) then

             fieldCur(i,j,bi,bj)=boundsVec(1)

      boundsVec(1)=-2.

      boundsVec(2)=-1.9

      boundsVec(3)=39.

      boundsVec(4)=40.

      boundsVec(5)=5.

     & xx_theta_dummy, boundsVec, theta, maskc, paramSmooth,

      boundsVec(1)=29.

      boundsVec(2)=29.5

      boundsVec(3)=40.5

      boundsVec(4)=41.

      boundsVec(5)=5.

     & xx_salt_dummy, boundsVec, salt, maskc, paramSmooth, 

      boundsVec(1)=0.

      boundsVec(2)=0.

      boundsVec(3)=0.

      boundsVec(4)=0.

      boundsVec(5)=0.      

     & xx_uvel_dummy, boundsVec, uvel, maskw, paramSmooth, 

     & xx_vvel_dummy, boundsVec, vvel, masks, paramSmooth, 

      boundsVec(1)=0.

      boundsVec(2)=0.

      boundsVec(3)=0.

      boundsVec(4)=0.

      boundsVec(5)=0.

     & xx_etan_dummy, boundsVec, etan, maskc, paramSmooth, 

      boundsVec(1)=0.

      boundsVec(2)=0.

      boundsVec(3)=0.

      boundsVec(4)=0.

      boundsVec(5)=0.

     & xx_relaxsst_dummy, boundsVec, 

      boundsVec(1)=0.

      boundsVec(2)=0.

      boundsVec(3)=0.

      boundsVec(4)=0.

      boundsVec(5)=0.

     & xx_relaxsss_dummy, boundsVec, 

      boundsVec(1)=0.

      boundsVec(2)=1. _d -4

      boundsVec(3)=0.9 _d -2

      boundsVec(4)=1. _d -2

      boundsVec(5)=0.

     & xx_bottomdrag_dummy, boundsVec, bottomdragfld, maskc, 

      boundsVec(1)=-0.4

      boundsVec(2)=-0.3

      boundsVec(3)=0.3

      boundsVec(4)=0.4

      boundsVec(5)=0.

     & xx_edtaux_dummy, boundsVec, eddyPsiX, maskw, paramSmooth,

     & xx_edtauy_dummy, boundsVec, eddyPsiY, masks, paramSmooth, 

      boundsVec(1)=1. _d -6

      boundsVec(2)=2. _d -6

      boundsVec(3)=0.9 _d -3

      boundsVec(4)=1. _d -3

      boundsVec(5)=0.

     & xx_diffkr_dummy, boundsVec, diffkr, maskc, paramSmooth, 

      boundsVec(1)=1. _d 1

      boundsVec(2)=2. _d 1

      boundsVec(3)=0.9 _d 4

      boundsVec(4)=1. _d 4

      boundsVec(5)=0.

     & xx_kapgm_dummy, boundsVec, kapgm, maskc, paramSmooth, 

      boundsVec(1)=3. _d 2

      boundsVec(2)=4. _d 2

      boundsVec(3)=2.5 _d 3

      boundsVec(4)=3.0 _d 3

      boundsVec(5)=0.

     & xx_kapredi_dummy, boundsVec, kapredi, maskc, paramSmooth, 

      _RL boundsVec(5)

      if ( boundsVec(5).GT.0.) then

       tmpMax=boundsVec(5)/sqrt(wFld3d(i,j,k,bi,bj))

      CALL CTRL_BOUND_3D(paramFld3d,maskFld3d,boundsVec,myThid)

     & boundsVec, paramFld2d, maskFld3d, paramSmooth, mythid )

      _RL boundsVec(5),tmpMax,xxDummyCur

      if ( boundsVec(5).GT.0.) then

       tmpMax=boundsVec(5)/sqrt(wFld2d(i,j,bi,bj))

      CALL CTRL_BOUND_2D(paramFld2d,maskFld3d,boundsVec,myThid)

      _RL boundsVec(5),tmpMax,xxDummyCur

     & boundsVec, paramFld3d, maskFld3d, paramSmooth, mythid )