w1=taucl(i,j,k,1)*(aiw(1,ib)+(aiw(2,ib)+(aiw(3,ib)

           ww=(w1+w2+w3)/taux

           g1=w1*(aig(1,ib)+(aig(2,ib)+(aig(3,ib)

           gg=(g1+g2+g3)/(w1+w2+w3)

          w1=-8.0

     *                 w1,p1,dwe,dpe,h11,h12,h13,trant)

     *                 w1,p1,dwe,dpe,h21,h22,h23,trant)

     *                 w1,p1,dwe,dpe,h81,h82,h83,trant)

          w1=-4.0

     *                w1,p1,dwe,dpe,c1,c2,c3,trant)

          w1=-6.0

     *                w1,p1,dwe,dpe,o1,o2,o3,trant)

      subroutine TABLUP(k1,k2,m,n,np,nx,nh,nt,sabs,spre,stem,w1,p1,

      _RL w1,p1,dwe,dpe

        we=(log10(x1)-w1)/dwe

        we=max(we,w1-2.*dwe)

      _RL w1,w2,w3,g1,g2,g3,ww,gg,ff,taux,reff1,reff2

      _RL  taux,reff1,reff2,w1,w2,g1,g2

           w1=(1.-(aia(ib,1)+(aia(ib,2)+

           ssaclt(i,j)=(w1+w2)/taux

           g1=(aig(ib,1)+(aig(ib,2)+aig(ib,3)*reff1)*reff1)*w1

           asyclt(i,j)=(g1+g2)/(w1+w2)

                CALL FLT_TRILINEAR(itx,jty,ktz,w1,wVel,4,bi,bj,myThid)

                  w1 = w1 + PORT_RAND_NORM()*flt_noise

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

     &                          + flt_deltaT*w1*scalez

      _RL ww, w1, ktz, kz, scalez

                CALL FLT_TRILINEAR(ix,jy,kz,w1,wVel,4,bi,bj,myThid)

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

                ktz1=kz+0.5*flt_deltaT*w1*scalez

     &                          * (w1 + 2*w2 + 2*w3 + w4) * scalez

      _RL w1, w2, w3, w4, ktz1, ktz2, ktz3, kz, scalez