real*8 x(20),y(100,100),z(100,100)
	real*8 pi,dx,dalpha,theta,alpha
	open(8,file='sphere.dat',access='sequential',status='unknown')
	pi=acos(-1.)
	n=19
	m=8
	dalpha=pi/2./(m-1)
	dx=pi/(n-1)
	do 10 i=1,n
	theta=pi-dx*(i-1)
	x(i)=cos(theta)
	do 20 j=1,m
	alpha=pi+dalpha*(j-1)
	y(i,j)=sin(theta)*cos(alpha)
	z(i,j)=sin(theta)*sin(alpha)
20	continue
10	continue
	mm=1
	write(8,*) n
	do 30 i=1,n
	if(i.eq.1.or.i.eq.n) write(8,*) x(i),mm
	if(i.ne.1.and.i.ne.n) write(8,*) x(i),m
	do 40 j=1,m
	write(8,*) y(i,j),z(i,j)
	if(i.eq.1.or.i.eq.n) go to 30
40	continue
30	continue
	stop
	end