from astropy.io import fits
import numpy as np
from astropy.table import Table
rcsed = fits.open('/data1/chil/work/uv-to-nir-catalog/local/lsb_psg/psg_rcsed.fits')[1].data
galfile_hdul = fits.open('/data1/Data/Binospec/coma3/coma3_c138-20180407/obj_abs_slits_lin_preproc_full.fits.gz')
objid = 587741722823622833

phot_mags = fits.open('/data1/chil/work/uv-to-nir-catalog/local/lsb_psg/833_bb_mags.fits')[1].data
phot_mags_err = fits.open('/data1/chil/work/uv-to-nir-catalog/local/lsb_psg/833_bb_mags_err.fits')[1].data

totmag = np.array([phot_mags['FUV'],phot_mags['NUV'],phot_mags['u_CFHT'],phot_mags['g'],phot_mags['r'],phot_mags['corrmag_z'],phot_mags['i_CFHT'],phot_mags['J_CFHT'],phot_mags['Ks_CFHT']])
toterr = np.array([phot_mags_err['FUV'],phot_mags_err['NUV'],phot_mags_err['u_CFHT'],phot_mags_err['g'],phot_mags_err['r'],phot_mags_err['i_CFHT'],phot_mags_err['J_CFHT'],phot_mags_err['Ks_CFHT']])
#totkc = np.array([rcsed_ent['kcorr_fuv'],rcsed_ent['kcorr_nuv'],rcsed_ent['kcorr_u'],rcsed_ent['kcorr_g'],rcsed_ent['kcorr_r'],rcsed_ent['kcorr_i'],rcsed_ent['kcorr_z'],rcsed_ent['kcorr_Y'],rcsed_ent['kcorr_J'],rcsed_ent['kcorr_H'],rcsed_ent['kcorr_K'],kcrj,kcrj,kcrj])
print totmag

for i,hdu in enumerate(galfile_hdul):
	if ((i>0)) :
		shp = galfile_hdul[i].data.shape
		print shp
		data = np.zeros((shp[0], len(totmag)))
		data[:] = totmag if i==86
		galfile_hdul[i].data = data

galfile_hdul.writeto('photfile_psg_slit.fits')