c
      parameter(nxmax=5557, nzmax=890)
c
      real vel(nxmax),
     +     den(nxmax),lam(nxmax),mu(nxmax),qp(nxmax),qs(nxmax)
c
      xmin=0.
      xmax=250.02
      zmin=0.
      zmax=40.005
      gsize=.045
c
      nx=nint((xmax-xmin)/gsize)+1
      nz=nint((zmax-zmin)/gsize)+1
c
      open(35, file='vp', form='unformatted')
      open(36, file='vp.ascii')
c
      open(41, file='density', form='unformatted', status='old')
      open(42, file='lambda', form='unformatted', status='old')
      open(43, file='mu', form='unformatted', status='old')
      open(44, file='qp', form='unformatted', status='old')
      open(45, file='qs', form='unformatted', status='old')
c
      do j=1,nz
         z=zmin+float(j-1)*gsize
c
         read(41) (den(i),i=1,nx)
         read(42) (lam(i),i=1,nx)
         read(43) (mu(i),i=1,nx)
         read(44) (qp(i),i=1,nx)
         read(45) (qs(i),i=1,nx)
c
         do i=1,nx
            x=xmin+float(i-1)*gsize
c
            vel(i)=((lam(i)+2.*mu(i))/den(i))**.5
            write(36,5) x,z,vel(i)
5           format(3f10.3)
         enddo
c
         write(35) (vel(i),i=1,nx)
c
      enddo
c
      stop
      end