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