lines = f.readlines()
lines = map(string.strip, lines)
lines = map(string.split, lines)
names = [line[0] for line in lines]
name_sts = array([float(line[7]) / 30. for line in lines])
file = sys.argv[1]
filt = pat.search(file).group(1)
f = FITS.FITS(file)
xs = f['CRVAL1'] + (arange(1, f['NAXIS1'] + 1) - f['CRPIX1']) * f['CDELT1']
ys = f['CRVAL2'] + (arange(1, f['NAXIS2'] + 1) - f['CRPIX2']) * f['CDELT2']
N = f['NAXIS1'] * f['NAXIS2']
X = dstack(meshgrid(xs, ys)).reshape((N, 2))
fl = interp2D(X, ravel(f.data()))
efl = interp2D(X, ravel(FITS.qread(file.replace('mean', 'std'))))
x, y, v, efluxes, mesh, xoff, xscale, yoff, yscale, mean_flux = getdata(
filt, False)
x = x * xscale + xoff
y = y * yscale + yoff
ids = concatenate([nonzero(diff(y))[0], array([len(y) - 1])])
sts = y[ids]
N = len(sts)
n_col = int(ceil(sqrt(N)))
n_row = int(ceil(float(N) / n_col))
print N, n_col, n_row
