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=25)
      real xd(8),yd(8),zd(8),fd(8)
      integer nr(8)
      data zero / 0. /
*
      open(7,file='bolsurf.spl')
      pm=(NODES-1)/2
* ... Coordinates:
      do 10 k=1,NODES-1
      do 10 j=1,NODES-1
      do 10 i=1,NODES-1
* ... Brick values:
      do 15 m=1,8
      xd(m)=float(i+mod(m-1,2))-pm
      yd(m)=float(j+mod((m-1)/2,2))-pm
      zd(m)=float(k+mod((m-1)/4,2))-pm
      term=xd(m)**2+yd(m)**2+zd(m)**2
15    fd(m)=10.**2-term
* ... AVS field:
c     write(13,*) fd(1)
* ... Brick bounds:
      dmin=fd(1)
      dmax=fd(1)
      do 20 m=2,8
      dmin=amin1(dmin,fd(m))
      dmax=amax1(dmax,fd(m))
20    continue
      if(zero.lt.dmin) goto 11
      if(zero.gt.dmax) goto 11
* ... Isosurfaces:
      call brick(xd,yd,zd,fd,zero)
11    continue
10    continue
      stop 'OK'
      end