PROGRAM axis2

c	Make PDB file for axis

	REAL	v(3)/ 34.040005, 18.409086, 14.65247 /
	REAL	w(3) / 0.18032, 0.52823, -0.82973 /
	REAL	x(3)

	OPEN (unit=1, form='FORMATTED', status='new', file='axis.pdb',
	1	carriagecontrol='LIST')

	DO i = 901, 999
	  di = (i-950)*2.
	  DO j = 1,3
	    x(j) = v(j) + di*w(j)
	  endDO
	  WRITE (1, 1000) i-900, i, x
1000	  FORMAT ('ATOM', i7, '  CA  GLY', i6, 4x, 3f8.3, '  1.00 20.00')
	endDO

	CLOSE (unit=1)
	end