program MAIN
*             ====
* This code was developed & is CopyLefted by:
*
* Han de Bruijn; Production&Services "A little bit of Physics would be  (===)
* TUD Computing Centre; P.O. Box 354  NO Idleness in Mathematics"(HdB) @-O^O-@
* 2600 AJ  Delft; The Netherlands -- http://pubwww.tudelft.nl/~rcpshdb  #/_\#
* E-mail: Han.deBruijn@RC.TUDelft.NL  Fax: +31 15 278 3787. Tel: 2751.   ###
*
      parameter (NODES=1500,LINES=4500)
      real xx(NODES),yy(NODES),zz(NODES)
      real xd(3),yd(3),zd(3)
      integer nr(3)
      logical done(LINES)
*
* ... Coordinates:
      open(1,file='coor2d.dat',status='old')
      do 10 k=1,NODES
      read(1,*,end=11) xx(k),yy(k),zz(k)
      NN=k
   10 continue
      stop 'NODES too small'
   11 close(1)
*
* ... Search path:
      open(1,file='found2d.dat')
      do 15 k=1,LINES
15    done(k)=.false.
      do 45 k=1,LINES
      read(1,*,end=46) last
45    done(last)=.true.
      stop 'LINES too small'
46    close(1)
*
* ... Connectivity:
      open(1,file='topol2d.dat')
* ... PostScript plotfile:
      open(7,file='meshpath.ps')
      call plots
      write(7,'(a)') '0. setgray'
      do 20 k=1,LINES
      read(1,*,end=21) (nr(j),j=1,3)
      do 25 i=1,3
      ii=nr(i)+1
      xd(i)=xx(ii)
      yd(i)=yy(ii)
25    continue
* ... Mesh plotting:
      if(done(k)) then
      if(k.ne.last) write(7,'(a)') '0.7 setgray'
      call plot(xd(3),yd(3),'u')
      do 55 i=1,3
      call plot(xd(i),yd(i),'d')
55    continue
      write(7,'(a)') 'closepath fill'
      endif
      write(7,'(a)') '0. setgray'
      call plot(xd(3),yd(3),'u')
      do 35 i=1,3
      call plot(xd(i),yd(i),'d')
35    continue
20    continue
      stop 'LINES too small'
21    close(1)
*
* ... Ending:
      write(7,'(a)') 'stroke showpage'
      close(7)
      stop 'OK'
      end