PROGRAM printmap2

c	Read and print map written by FSFOUR2

	CHARACTER*80	filnam
	INTEGER		itle(20), is(2,3,24)
	REAL		para(6), ts(3,24), map(500)

	WRITE (6, 1000) 'Enter map file name:'
1000	FORMAT (a)
	READ (5, 1000) filnam
	OPEN (unit=1, form='UNFORMATTED', file=filnam, status='OLD', readonly)

	WRITE (6, 1000) 'Enter print file name:'
	READ (5, 1000) filnam
	OPEN (unit=2, form='FORMATTED', file=filnam, status='UNKNOWN',
	1	carriagecontrol='LIST')

	READ (1) itle, noset, para, nsym, ncent, latt, npic

	WRITE (2, 1001) itle
1001	FORMAT (20a4)
	WRITE (2, 1002) para
1002	FORMAT (6f10.4)
	WRITE (2, 1003) noset, nsym, ncent, latt, npic
1003	FORMAT (10i8)

	DO i = 1,nsym
	  READ (1) (ts(j,i), (is(k,j,i), k=1,2), j=1,3)
	  WRITE (2, 1004) (ts(j,i), (is(k,j,i), k=1,2), j=1,3)
1004	  FORMAT (3(f10.4, 2i5, 3x))
	endDO

	READ (1) nx, nz, ny, scale, norm
	WRITE (2, 1005) nx, ny, nz, scale, norn
1005	FORMAT (3i8, e15.7, i8) 	

	DO iy = 1,ny
	  DO iz = 1,nz
	    READ (1) (map(i), i=1,nx)
	    WRITE (2, 1006) (map(i), i=1,nx)
1006	    FORMAT (10f12.6)
	  endDO
	endDO

	CLOSE (unit=1)
	CLOSE (unit=2)
	end