Run wlmGALEX images through to get SB HXZ: 02/25/10 Get exposure times: (Check exp FITS images) NGA_WLM-fd-exp.fits 1420.71 NGA_WLM-nd-exp.fits 1291.93 ======> fuv: 1420.71 s; nuv: 1291.93 s 1)Cut out smaller images: cutout subimage size: 1600*1600 centered around the galaxy center of cutout: imcopy NGA_WLM-fd-int.fits[1132:2731,1113:2712] wlmfcut.fit imcopy NGA_WLM-nd-int.fits[1132:2731,1113:2712] wlmncut.fit -------------------------------------------- 2)Remove foreground stars and background galaxies: # imedit wlmncut.fit wlmncl.fit logfile=wlmncut.log imedit wlmfcut.fit wlmfcl.fit logfile=wlmfcut.log ========> wlmncl.fit; wlmfcl.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 wlmfcl.fit wlmfblkav.fit 10 10 blkavg wlmncl.fit wlmnblkav.fit 10 10 ========> wlmfblkav.fit; wlmnblkav.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 = 78.073717) initial isophote center X (y0 = 79.301524) initial isophote center Y (ellip0 = 0.19) initial ellipticity (pa0 = 70.6) initial position angle (degrees) (sma0 = 30.) initial semi-major axis lenght (minsma = 30.) minimum semi-major axis lenght (maxsma = 80.) 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 wlmfblkav.fit fsky_ell ellip wlmnblkav.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 14.0 arcmin annulus (13th annul in *_ell.lst) # The sky variations around the galaxy are -- 5.58E-5 for FUV, 2.17E-4 for NUV imarith wlmfcl.fit - 2.63E-4 wlmfms.fit imarith wlmncl.fit - 0.00291 wlmnms.fit =========> wlmnms.fit; wlmfms.fit #Add SKYMEAN keyword to these fits headers sethead wlmfms.fit skymean= 2.63E-4 sethead wlmnms.fit skymean= 0.00291 ------------------------------------------- #Geometrically transform the NUV and FUV images to match the pixel scale and orientation of # the original V image: #I used 4 field stars for geometry transformation center wlmncut.fit center wlmv.fits txdump wlmv.fits.ctr.1 xc,yc, yes > v.xy txdump wlmncut.fit.ctr.1 xc,yc yes > n.xy joinlines v.xy n.xy > vn.xy imheader wlmv.fits wlmv.fits[1009,991][real]: wlm geomap vn.xy vn.geomap 1 1009 1 991 fitgeo=general functio=polyn xxorder=2 yyorder=2 xyorder=2 xxterms=half =====> V image is 1009 * 991. rms in geomap is 0.014 and 0.010 for x and y geotran wlmfms.fit wlmfms_geo.fit vn.geomap vn.xy geotran wlmnms.fit wlmnms_geo.fit vn.geomap vn.xy geotran wlmfcut.fit wlmfcut_geo.fit vn.geomap vn.xy geotran wlmncut.fit wlmncut_geo.fit vn.geomap vn.xy =====> wlmfms_geo.fit; wlmnms_geo.fit; wlmfcut.fit; wlmncut_geo.fit ---------------------------------------------------------------------------------------- Ellipse photometry: ### #Use wlmncl.tab as inellip to do photometry # #. ellipse wlmfms_geo.fit wlmf_ell inellip=wlmncl.tab inter- region+ ellipse wlmnms_geo.fit wlmn_ell inellip=wlmncl.tab inter- region+ ====> wlmf_ell.tab; wlmn_ell.tab isoexam wlmf_ell.tab isoexam wlmn_ell.tab ====> wlmf_ell.ps; wlmn_ell.ps tprint wlmf_ell.tab > wlmf.phot tprint wlmn_ell.tab > wlmn.phot ====> wlmf.phot; wlmn.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. wlmf.phot wlmn.phot wlmuv.final82 9 1420.71 1291.93 0.000150314 0.00166317 0.44 1.0 0.068 1.134 204.22 118.44 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 =====> wlmuv.final82 ---------------------------------------------------------------------------------------- Plot: uv_82.sm sm input uv_82.sm =====> lgs3uv_82.eps -------------------------------------- Fit exponential disk: fields wlmuv.final82 2,5 lines="18-26" | polyfit stdin 1 > wlmn_.expfit This gives you mu_NUV = mu_NUV^0 + slope * R(arcmin) ../expfit.exe ### This calculates R_D for you. give it wlmn_.expfit mv expfit.out82 wlmn_.expout ===========> R_D= 1.393 +/- 0.072 arcmin R_D= 0.405 +/- 0.021 kpc mu_NUV_0= 23.780 +/- 0.198 y(100kpc)= 291.66 ----------------------------------- Color: ../avecolors.exe Final file=wlmuv.final82 Beg line for FUV-NUV(0=none)=18 End line for FUV-NUV(0=none)=26 =======>Ave FUV-NUV=0.281 +/- 0.192 fields wlmuv.final82 3,9 lines="18-21" | polyfit stdin 1 > wlmfmn_inner.expfit This calculates FUV-NUV = zeropt + slope * R(kpc) =====> zeropt = 0.07922376 +/- 0.01802845; slope = 0.008684075 +/- 0.02863643 fields wlmuv.final82 3,9 lines="21-26" | polyfit stdin 1 > wlmfmn_outer.expfit =====> zeropt = -0.1424424 +/- 0.08379252; slope = 0.3183489 +/- 0.04883752 y(100kpc)= 31.6924 # ../intersect.exe 1st file returned by POLYFIT: wlmfmn_inner.expfit 2nd file returned by POLYFIT: wlmfmn_outer.expfit ------> intersection x0 : 0.716 intersection x0 in Kpc: 0.208 intersection y0: 0.085 ==========> Rbr = 0.716 arcmin = 0.208 kpc uvbr = 0.085 ---------------- 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 wlmuv_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 =====> wlmallphot82.dat