c
c reef model building program
c
c............. code deleted here .....................
c          the velocity at depth is determined by line segments 
c
           nline = 6
c
           bm(1) =  150.0
           bm(2) =  400.0
           bm(3) =  700.0
           bm(4) =  900.0
           bm(5) =  1300.0
           bm(6) =  1600.0
c
           am(1) = +0.02
           am(2) = -0.03
           am(3) = -0.02
           am(4) = -0.02
           am(5) = -0.03
           am(6) = -0.00
c
           vz(1) = 2000.0
           vz(2) = 2600.0
           vz(3) = 3200.0
           vz(4) = 2000.0
           vz(5) = 3600.0
           vz(6) = 3900.0
c
           dx  = 5.0
           dz  = 5.0
c
c          reef circle
c
           l_reef = 3
           xc = (nx-1)*dx*0.5
           zc = 800.0
           rc = 320.0
           salt = 4500.0
c
c          create initial water layer
c
           do i=1,nx
            do j=1,nz
            src(i,j) = 0.0
            vv(i,j) = 1500.0
            enddo
           enddo
c
c         now place a circle for the reef 
c
           ilmax = nx
           irmax = 1
c
           do i=1,nx
            do j=1,nz
            x = (i-1)*dx
            z = (j-1)*dz
            r = sqrt((x-xc)**2+(z-zc)**2)
c
c           if point is inside circle it's the reef
c
           if( r .lt.  rc ) then
               if( i .lt. ilmax ) ilmax = i
               if( i .gt. irmax ) irmax = i
           endif
           enddo
          enddo
c
c          now layer by layer make velocities under each
c             line 
c
c          each new layer over writes that below
c
           do k=1,nline
            do i=1,nx
             do j=1,nz
             x = (i-1)*dx
             z = (j-1)*dz
             zl = am(k)*x + bm(k)
c
             if( z .gt. zl ) then
              vmult = 1.0
              if(l_reef .eq. k) then
               if( i .gt. ilmax .and. i .lt. irmax ) vmult = 1.4
               if( i .gt. irmax ) vmult = 1.2
              endif
               vv(i,j) = vz(k)*vmult
             endif
            enddo
           enddo
          enddo
c
c.............. code deleted here ..................