read "itall.txt":

lvl:=4:
phatlvl:=8:
numpsis:=(3^(lvl+1)-3)/2:
path:=cat("../data/","level",lvl,"/"):
T:=gasket(phatlvl):


linever:=our_line_vertices(phatlvl,[2],2,4):
for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i)):
	#ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
	#final_points:=SG_plot_line(T,8,12,linever,10):
	pointplot(SG_plot_line(T,8,12,linever,10),color=blue,style=LINE,title="blah");
end do;
print("done printing psis");