CALL FLT_TRILINEAR(ix,jy,kz,vv,vVel,2,bi,bj,myThid)

                CALL FLT_BILINEAR(ix,jy,vv,vVel,kc,2,bi,bj,myThid)

                vv = vv + vv*(PORT_RAND(seed)-0.5)*flt_noise

              jty=jy+0.5*flt_deltaT*vv*scaley

      _RL uu, vv, u1, v1

              CALL FLT_BILINEAR  (ix,jy,vv,vVel,  kp,2,bi,bj,myThid)

              tmp(11) = vv

      _RL uu, vv, tt, ss, pp

               CALL FLT_BILINEAR  (ix,jy,vv,vVel,  k,2,bi,bj,myThid)

               tmp(9+k+1*Nr) = vv

      _RL uu,vv,tt,ss, pp

     U     vv,w,wk1,wk2,

     U     sol,i,its,vv,w,wk1,wk2,

      _RL vv(n,im+1,nSx,nSy), w(n,im,nSx,nSy)

         vv(j,1,bi,bj) = rhs(j,bi,bj) - wk2(j,bi,bj)

      call SEAICE_SCALPROD(n, im+1, 1, 1, vv, vv, ro, myThid)

         vv(j,1,bi,bj) = vv(j,1,bi,bj)*t

         wk1(k,bi,bj) = vv(k,i,bi,bj)

         vv(k,i1,bi,bj)=wk2(k,bi,bj)

       call SEAICE_SCALPROD(n, im+1, j, i1, vv, vv, t, myThid)

          vv(k,i1,bi,bj) = vv(k,i1,bi,bj) - t*vv(k,j,bi,bj)

      call SEAICE_SCALPROD(n, im+1, i1, i1, vv, vv, t, myThid)

          vv(k,i1,bi,bj) = vv(k,i1,bi,bj)*t

          vv(k,1,bi,bj) = vv(k,1,bi,bj) + t*vv(k,j,bi,bj)

      _RL vv(nVec,im+1,nSx,nSy), w(nVec,im,nSx,nSy)

      COMMON /FGMRES_RL/ sol, rhs, vv, w

      _RL   vv(nmax)

      call OPTIM_READDATA ( nn, ctrlname, lheaderonly, ff, vv)

      _RL            vv(nn)

      _RL            vv

      pointer (pvv,vv(1))

               vv(icvoffset+icvcomp) = cbuff(icvcomp)

      _RL     vv(nn)

      _RL            vv

      pointer (pvv,vv(1))

      _RL   vv(nmax)

     O                      vv

                              cbuff(icvcomp) = vv(icvoffset + icvcomp)

      _RL     vv(nn)

     I                       vv