c      graphics test pattern program
c      file name  /home/g3/ohta/sample/test52.f  main52.f
c      f77 filename.f   a.out   lpr -c -s filename.ps
c      f77 -c filename.f   f77 fname1.o fname2.o   a.out
c
       parameter (ipci=2,ity=1,isz=1,mx=9,my=9)
       parameter (x1=1.0,y1=2.0,xl=17.0,yl=23.0)
       parameter (nx=62,ny=02,n1=nx*ny)
       dimension u(n1)
       character fname*40/"plot52f.ps"/
c
       open(10,file=fname,access="sequential")
       call xyopen
       no=1
       call plots(ipci,ity,isz,no)
       call factor(1.00)
c
       do 300 jj=1,my
       do 300 ii=1,mx
       ico=(ii-1)+10*(jj-1)
       num=ii+mx*(jj-1)
       i2=mod(num-1,mx)+1
       i3=mod((num-i2)/mx,my)+1
       dx=xl/float(mx)
       dy=yl/float(my)
c      x0=x1+dx*float(i2-1)
c      y0=y1+dy*float(my-i3)
       x0=x1+dx*float(ii-1)
       y0=y1+dy*float(my-jj)
c
       do 20 j=1,ny
       do 20 i=1,nx
       i1=i+nx*(j-1)
       u(i1)=float(i)/float(nx)
   20  continue
       call imagc1(nx,ny,x0,y0,dx,dy,u,ico)
c
  300  continue
       call plote
       stop
       end