I                maskFld3d, boundsVec, myThid )

      _RL boundsVec(5)

      x0=boundsVec(1)

      x1=boundsVec(2)

      x2=boundsVec(3)

      x3=boundsVec(4)

     I                maskFld3d, boundsVec, myThid )

      _RL boundsVec(5)

      x0=boundsVec(1)

      x1=boundsVec(2)

      x2=boundsVec(3)

      x3=boundsVec(4)

      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)

     I                maskFld3d, 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                maskFld3d, boundsVec, myThid )

      _RL boundsVec(5)

      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)=4. _d -4

      boundsVec(4)=5. _d -4

      boundsVec(5)=0.

     & xx_diffkr_dummy, boundsVec, diffKr, maskc, paramSmooth,

      boundsVec(1)=1. _d 2

      boundsVec(2)=2. _d 2

      boundsVec(3)=0.9 _d 4

      boundsVec(4)=1. _d 4

      boundsVec(5)=0.

     & xx_kapgm_dummy, boundsVec, kapGM, maskC, paramSmooth,

      boundsVec(1)=1. _d 2

      boundsVec(2)=2. _d 2

      boundsVec(3)=0.9 _d 4

      boundsVec(4)=1.0 _d 4

      boundsVec(5)=0.

     & xx_kapredi_dummy, boundsVec, kapRedi, maskC, paramSmooth,

      _RL boundsVec(5)

     & boundsVec, paramFld3d, maskFld3d, paramSmooth, mythid )

      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