Commit a239436a authored by LUSTIG Peter's avatar LUSTIG Peter

wcs is now written to output file

parent 8e9816af
......@@ -116,10 +116,10 @@ flux = np.geomspace(1, 10, 3) * u.mJy
flux = [0]
flux = np.linspace(1, 3, 2)*u.mJy
flux = np.array([10]) * u.mJy
nsim = 6
nsim = 4
min_detection_threshold = 3
nsources = 2
nsources = 3
nsources = 0
outdir = Path("montecarlo_results/")
outdir = Path("testdir")
ncores = 2
......@@ -250,7 +250,7 @@ if nsources != 0:
fname = timeprefix + fname
outfile = outdir / fname
phdu = fits.PrimaryHDU()
phdu = fits.PrimaryHDU(header=jackknifes[0].wcs.to_header())
phdu.header['influx'] = '{}'.format(_flux)
phdu.header['nsim'] = nsim
phdu.header['sourcespersim'] = nsources
......@@ -316,7 +316,7 @@ if nsources == 0:
fname = timeprefix + fname
outfile = outdir / fname
phdu = fits.PrimaryHDU()
phdu = fits.PrimaryHDU(header=jackknifes[0].wcs.to_header())
phdu.header['influx'] = 0
phdu.header['nsim'] = nsim
phdu.header['sourcespersim'] = 0
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment