Run wlmGALEX images through to get SB HXZ: 02/25/10 Get exposure times: (Check exp FITS images) NGA_NGC1569-fd-exp.fits 6938.41 NGA_NGC1569-nd-exp.fits 6057.33 ======> fuv: 6938.41 s; nuv: 6057.33 s 1)Cut out smaller images: cutout subimage size: 1600*1600 centered around the galaxy center of cutout: imcopy NGA_NGC1569-fd-int.fits[1531:2330,1511:2310] n1569fcut.fit imcopy NGA_NGC1569-nd-int.fits[1531:2330,1511:2310] n1569ncut.fit -------------------------------------------- 2)Remove foreground stars and background galaxies: # imedit n1569fcut.fit n1569fcl.fit cursor=n1569fcut.log imedit n1569ncut.fit n1569ncl.fit cursor=n1569ncut.log ========> n1569ncl.fit; n1569fcl.fit ---------------------------------------------------------- 3)Get the background now that the extraneous stuff has been removed: #Note: #The low counts rates in GALEX images leads to poisson, rather than gaussian noise. #which leads to the failure of using the usual sky-fitting routines in IRAF #Fit the sky of the original images usually failed. #So I selected many sky regions around the galaxy to get a mean sky level and the associated variance # #Here I first average the original images with 10*10 bin size ======= blkavg n1569fcl.fit n1569fblkav.fit 10 10 blkavg n1569ncl.fit n1569nblkav.fit 10 10 ========> n1569fblkav.fit; n1569nblkav.fit #Then use ELLIPSE to get mean pixel intensity and the associated rms in the elliptical annular sky regions around the galaxy. #I fixed the ELLIPSE parameters to that used to do surface photometry. # PACKAGE = isophote TASK = geompar (x0 = 40.33248) initial isophote center X (y0 = 40.188214) initial isophote center Y (ellip0 = 0.45) initial ellipticity (pa0 = -58.9) initial position angle (degrees) (sma0 = 5.) initial semi-major axis lenght (minsma = 5.) minimum semi-major axis lenght (maxsma = 37.) maximum semi-major axis lenght (step = 3.) sma step between successive ellipses (linear = yes) linear sma step ? (maxrit = INDEF) maximum sma lenght for iterative mode (recente= no) allows finding routine to re-center x0-y0 ? (xylearn= no) updates pset with new x0-y0 ? (physica= no) physical coordinate system ? (mode = al) ------- ellip n1569fblkav.fit fsky_ell ellip n1569nblkav.fit nsky_ell tprint fsky_ell.tab columns='sma,intens,pix_var,ndata' > fsky_ell.lst tprint nsky_ell.tab columns='sma,intens,pix_var,ndata' > nsky_ell.lst # I choose the surface photometry intensity on the 5.0 arcmin annulus (6th annul in *_ell.lst) # The sky variations around the galaxy are -- 5.28E-5 for FUV, 2.81E-4 for NUV imarith n1569fcl.fit - 8.67E-4 n1569fms.fit imarith n1569ncl.fit - 0.00472 n1569nms.fit =========> n1569nms.fit; n1569fms.fit #Add SKYMEAN keyword to these fits headers sethead n1569fms.fit skymean= 8.67E-4 sethead n1569nms.fit skymean= 0.00472 ------------------------------------------- #Geometrically transform the NUV and FUV images to match the pixel scale and orientation of # the original V image: #I used vn.geomap vn.xy used before by deidre geotran n1569fms.fit n1569fms_geo.fit vn.geomap vn.xy geotran n1569nms.fit n1569nms_geo.fit vn.geomap vn.xy geotran n1569fcut.fit n1569fcut_geo.fit vn.geomap vn.xy geotran n1569ncut.fit n1569ncut_geo.fit vn.geomap vn.xy =====> n1569fms_geo.fit; n1569nms_geo.fit; n1569fcut.fit; n1569ncut_geo.fit ---------------------------------------------------------------------------------------- Ellipse photometry: ### #Use n1569n.tab as inellip to do photometry, and ngc1569fms_geo.pl; ngc1569nms_geo.pl used before #rename ngc1569fms_geo.pl; ngc1569nms_geo.pl to n1569fms_geo.pl; n1569nms_geo.pl # #. ellipse n1569fms_geo.fit n1569f_ell inellip=n1569n.tab inter- region+ ellipse n1569nms_geo.fit n1569n_ell inellip=n1569n.tab inter- region+ ====> n1569f_ell.tab; n1569n_ell.tab isoexam n1569f_ell.tab isoexam n1569n_ell.tab ====> n1569f_ell.ps; n1569n_ell.ps tprint n1569f_ell.tab > n1569f.phot tprint n1569n_ell.tab > n1569n.phot ====> n1569f.phot; n1569n.phot ---------------------------------------------------------------------------------------- Create input_galex file: FUV photometry file: 0=no FUV NUV photometry file: 0=no NUV Final output file Number of ellipses Integration time --- FUV Integration time --- NUV Average sky (counts/s/pixel) --- FUV Average sky (counts/s/pixel) --- NUV b/a D(Mpc) E(B-V) total scale of V image in "/pixel R_25 (arcsec), 9999.9 if none R_D (arcsec), 9999.9 if none The original input_galex file should work fine, but change to whatever filenames you are using. n1569f.phot n1569n.phot n1569uv.final82 10 6938.41 6057.33 6.79728e-05 0.000370048 0.55 2.5 0.558 0.42 9999.9 23.22 The sky is determined by scaling the sky to the V-band pixel size. scale as (V-band pixel size/1.5)**2. =====> input_galex ---------------------------------------------------------------------------------------- Calibrate the photometry and determine surface brightnesses: ../annuli_galex82.exe =====> n1569uv.final82 ---------------------------------------------------------------------------------------- Plot: uv_82.sm sm input uv_82.sm =====> n15693uv_82.eps -------------------------------------- Fit exponential disk: fields n1569uv.final82 2,5 lines="19-28" | polyfit stdin 1 > n1569n_.expfit This gives you mu_NUV = mu_NUV^0 + slope * R(arcmin) ../expfit.exe ### This calculates R_D for you. give it n1569n_.expfit mv expfit.out82 n1569n_.expout ===========> R_D= 0.307 +/- 0.027 arcmin R_D= 0.224 +/- 0.019 kpc mu_NUV_0= 18.725 +/- 0.370 y(100kpc)= 504.52 ----------------------------------- Color: ../avecolors.exe Final file=n1569uv.final82 Beg line for FUV-NUV(0=none)=19 End line for FUV-NUV(0=none)=28 =======>Ave FUV-NUV=0.630 +/- 0.120 Fit line to the first two inner points 0.08 kpc 0.449 0.23 kpc 0.494 =====> 0 kpc, 0.425 0.38 kpc, 0.539 fields n1569uv.final82 3,9 lines="21-28" | polyfit stdin 1 > n1569fmn_outer.expfit =====> zeropt = 0.9098857 +/- 0.04366181; slope = -0.2629039 +/- 0.04449817 y(1kpc)= 0.646982 # ---------------- Revise the FUV-NUV lines in uv_82.sm update : line 107 , line 108 for NUV SB fitting plot. update : line 182 , line 183 for UV color plot with new Ave color Redraw n1569uv_82.eps ------------------------------------------ Combine with other photometry: Enter in galex2/photfiles_rev.list and cp to tempfiles.list and edit out all but galaxy that you are working on, or do them all at the end. ../collectv_lt.exe =====> n1569allphot82.dat