reads,'FOC.d0h',a,nx,ny
p1=a(53:152,357:456)
p2=a(53:152,205:304)
p3=a(53:152,53:152)
p4=a(205:304,357:456)
p5=a(205:304,205:304)
p6=a(357:456,357:456)
p7=a(357:456,205:304)
x=indgen(100)#replicate(1.,100)
y=replicate(1.,100)#indgen(100)
x1=total(p1*x)/total(p1) & y1=total(p1*y)/total(p1)
x2=total(p2*x)/total(p2) & y2=total(p2*y)/total(p2)
x3=total(p3*x)/total(p3) & y3=total(p3*y)/total(p3)
x4=total(p4*x)/total(p4) & y4=total(p4*y)/total(p4)
x5=total(p5*x)/total(p5) & y5=total(p5*y)/total(p5)
x6=total(p6*x)/total(p6) & y6=total(p6*y)/total(p6)
x7=total(p7*x)/total(p7) & y7=total(p7*y)/total(p7)
print,x1,x2,x3,x4,x5,x6,x7
print,y1,y2,y3,y4,y5,y6,y7
p1=shift(p1,52-x1,52-y1)
p2=shift(p2,52-x2,52-y2)
p3=shift(p3,52-x3,52-y3)
p4=shift(p4,52-x4,52-y4)
p5=shift(p5,52-x5,52-y5)
p6=shift(p6,52-x6,52-y6)
p7=shift(p7,52-x7,52-y7)
x1=total(p1*x)/total(p1) & y1=total(p1*y)/total(p1)
x2=total(p2*x)/total(p2) & y2=total(p2*y)/total(p2)
x3=total(p3*x)/total(p3) & y3=total(p3*y)/total(p3)
x4=total(p4*x)/total(p4) & y4=total(p4*y)/total(p4)
x5=total(p5*x)/total(p5) & y5=total(p5*y)/total(p5)
x6=total(p6*x)/total(p6) & y6=total(p6*y)/total(p6)
x7=total(p7*x)/total(p7) & y7=total(p7*y)/total(p7)
print,x1,x2,x3,x4,x5,x6,x7
print,y1,y2,y3,y4,y5,y6,y7
ptot=p1+p2+p3+p4+p5+p6+p7
tvscl,congrid(ptot,512,512)
stop
end