global x1 y1 z1
global creatingclines

fid = fopen('setprob.data');
data = fscanf(fid,'%g',[2,1]);
fclose(fid);
xc1 = data(1);
yc1 = data(2);
Rsphere = 1;
r = Rsphere * cos(yc1);
z1 = Rsphere * sin(yc1);
x1 = r * cos(xc1);
y1 = r * sin(xc1);
