stilts tmatchn multimode=pairs nin=3 matcher=exact values1='objid' values2='sdss_objid' values3='objid' in1=chil_z_lt_003_iz_1.fit in2=galex_sdss_ukidss_lt_003_chil.xml.gz in3=UKIDSS_DR8_chil_z_lt_003_iz_1.fits ofmt=fits fixcols=dups join1=always ocmd='keepcols "objid_1 ra_1a dec_1a z zerr zconf petror50_r modelmag* *petromag* fuv_mag fuv_magerr nuv_mag nuv_magerr extinction_* aY aJ aH aK e_bv specClass"' > galaxies_11bands_lt003.fits

#stilts tpipe in=galaxies_11bands_lt003.fits ofmt=fits cmd='addcol corrmag_fuv "fuv_mag-extinction_g*8.29/3.793"; addcol corrmag_nuv "nuv_mag-extinction_g*8.87/3.793"; addcol corrmag_u "modelmag_u-extinction_u"; addcol corrmag_g "modelmag_g-extinction_g"; addcol corrmag_r "modelmag_r-extinction_r"; addcol corrmag_i "modelmag_i-extinction_i"; addcol corrmag_z "modelmag_z-extinction_z"; addcol corrmag_Y "yPetroMag-extinction_g/3.132+0.634"; addcol corrmag_J "j_1PetroMag-extinction_g/4.267+0.938"; addcol corrmag_H "hPetroMag-extinction_g/6.562+1.379"; addcol corrmag_K "kPetroMag-extinction_g/10.536+1.900"; ' > galaxies_11bands_lt003_extcorr.fits

stilts tpipe in=galaxies_11bands_lt003.fits ofmt=fits cmd='select "z>=0.007 && specClass==2"' cmd='addcol corrmag_fuv "fuv_mag-extinction_g*8.29/3.793"; addcol corrmag_nuv "nuv_mag-extinction_g*8.87/3.793"; addcol corrmag_u "modelmag_u-extinction_u"; addcol corrmag_g "modelmag_g-extinction_g"; addcol corrmag_r "modelmag_r-extinction_r"; addcol corrmag_i "modelmag_i-extinction_i"; addcol corrmag_z "modelmag_z-extinction_z"; addcol corrmag_Y "yPetroMag-extinction_g/3.132+0.634"; addcol corrmag_J "j_1PetroMag-extinction_g/4.267+0.938"; addcol corrmag_H "hPetroMag-extinction_g/6.562+1.379"; addcol corrmag_K "kPetroMag-extinction_g/10.536+1.900"; ' > galaxies_11bands_0007_003_extcorr.fits

stilts tpipe in=galaxies_11bands_0007_003_extcorr.fits ofmt=fits cmd='addcol FUV_u "corrmag_fuv-corrmag_u"; addcol NUV_g "corrmag_nuv-corrmag_g"; addcol u_r "corrmag_u-corrmag_r"; addcol g_r "corrmag_g-corrmag_r"; addcol g_i "corrmag_g-corrmag_i"; addcol r_z "corrmag_r-corrmag_z"; addcol Y_H "corrmag_Y-corrmag_H"; addcol J_K "corrmag_J-corrmag_K"; addcol H_K "corrmag_H-corrmag_H"; addcol kcorr_FUV "0-1.097750*pow(z, 1)*pow(FUV_u, 0)-0.383494*pow(z, 1)*pow(FUV_u, 1)+0.757580*pow(z, 1)*pow(FUV_u, 2)-0.144215*pow(z, 1)*pow(FUV_u, 3)-5.164474*pow(z, 2)*pow(FUV_u, 0)+11.595245*pow(z, 2)*pow(FUV_u, 1)-3.901377*pow(z, 2)*pow(FUV_u, 2)+0.328657*pow(z, 2)*pow(FUV_u, 3)+43.577339*pow(z, 3)*pow(FUV_u, 0)-41.071255*pow(z, 3)*pow(FUV_u, 1)+5.497657*pow(z, 3)*pow(FUV_u, 2)-89.706520*pow(z, 4)*pow(FUV_u, 0)+39.843071*pow(z, 4)*pow(FUV_u, 1)+55.914093*pow(z, 5)*pow(FUV_u, 0)"; addcol kcorr_NUV "0+2.517001*pow(z, 1)*pow(NUV_g, 0)-2.153217*pow(z, 1)*pow(NUV_g, 1)+0.613728*pow(z, 1)*pow(NUV_g, 2)-0.012136*pow(z, 1)*pow(NUV_g, 3)-22.935999*pow(z, 2)*pow(NUV_g, 0)+5.627578*pow(z, 2)*pow(NUV_g, 1)-0.798637*pow(z, 2)*pow(NUV_g, 2)-0.088559*pow(z, 2)*pow(NUV_g, 3)+90.253288*pow(z, 3)*pow(NUV_g, 0)+0.180280*pow(z, 3)*pow(NUV_g, 1)+1.208637*pow(z, 3)*pow(NUV_g, 2)-201.441422*pow(z, 4)*pow(NUV_g, 0)-9.612353*pow(z, 4)*pow(NUV_g, 1)+181.680908*pow(z, 5)*pow(NUV_g, 0)"; addcol kcorr_u "0+9.186063*pow(z, 1)*pow(u_r, 0)-6.961030*pow(z, 1)*pow(u_r, 1)+3.131392*pow(z, 1)*pow(u_r, 2)-0.397524*pow(z, 1)*pow(u_r, 3)-94.418167*pow(z, 2)*pow(u_r, 0)+34.886288*pow(z, 2)*pow(u_r, 1)-8.957370*pow(z, 2)*pow(u_r, 2)+1.090704*pow(z, 2)*pow(u_r, 3)+314.600128*pow(z, 3)*pow(u_r, 0)-28.369019*pow(z, 3)*pow(u_r, 1)-1.226428*pow(z, 3)*pow(u_r, 2)-505.852478*pow(z, 4)*pow(u_r, 0)+24.577915*pow(z, 4)*pow(u_r, 1)+278.303589*pow(z, 5)*pow(u_r, 0)"; addcol kcorr_g "0-3.671975*pow(z, 1)*pow(g_r, 0)+10.915098*pow(z, 1)*pow(g_r, 1)-6.483314*pow(z, 1)*pow(g_r, 2)-1.932091*pow(z, 1)*pow(g_r, 3)+74.797501*pow(z, 2)*pow(g_r, 0)-156.327652*pow(z, 2)*pow(g_r, 1)+251.265198*pow(z, 2)*pow(g_r, 2)-41.482807*pow(z, 2)*pow(g_r, 3)-705.632324*pow(z, 3)*pow(g_r, 0)+43.148262*pow(z, 3)*pow(g_r, 1)-1060.333496*pow(z, 3)*pow(g_r, 2)+219.644760*pow(z, 3)*pow(g_r, 3)+4637.382812*pow(z, 4)*pow(g_r, 0)+2550.372314*pow(z, 4)*pow(g_r, 1)+1242.749512*pow(z, 4)*pow(g_r, 2)-284.129944*pow(z, 4)*pow(g_r, 3)-16852.767578*pow(z, 5)*pow(g_r, 0)-5959.362305*pow(z, 5)*pow(g_r, 1)+48.727806*pow(z, 5)*pow(g_r, 2)+28531.974609*pow(z, 6)*pow(g_r, 0)+3114.860840*pow(z, 6)*pow(g_r, 1)-17274.367188*pow(z, 7)*pow(g_r, 0)"; addcol kcorr_r "0+1.691848*pow(z, 1)*pow(g_r, 0)-1.736085*pow(z, 1)*pow(g_r, 1)+4.401005*pow(z, 1)*pow(g_r, 2)-3.943418*pow(z, 1)*pow(g_r, 3)-21.236748*pow(z, 2)*pow(g_r, 0)+1.545018*pow(z, 2)*pow(g_r, 1)+29.598732*pow(z, 2)*pow(g_r, 2)+5.267979*pow(z, 2)*pow(g_r, 3)+63.996235*pow(z, 3)*pow(g_r, 0)-132.826462*pow(z, 3)*pow(g_r, 1)-59.209694*pow(z, 3)*pow(g_r, 2)+77.660744*pow(z, 4)*pow(g_r, 0)+251.961243*pow(z, 4)*pow(g_r, 1)-265.911682*pow(z, 5)*pow(g_r, 0)"; addcol kcorr_i "0-2.163365*pow(z, 1)*pow(g_i, 0)+2.893258*pow(z, 1)*pow(g_i, 1)+2.392071*pow(z, 1)*pow(g_i, 2)-2.110186*pow(z, 1)*pow(g_i, 3)-8.363442*pow(z, 2)*pow(g_i, 0)-30.825390*pow(z, 2)*pow(g_i, 1)+20.451029*pow(z, 2)*pow(g_i, 2)+4.460981*pow(z, 2)*pow(g_i, 3)+107.313011*pow(z, 3)*pow(g_i, 0)-34.375408*pow(z, 3)*pow(g_i, 1)-55.385445*pow(z, 3)*pow(g_i, 2)-108.093277*pow(z, 4)*pow(g_i, 0)+201.405807*pow(z, 4)*pow(g_i, 1)-168.219467*pow(z, 5)*pow(g_i, 0)"; addcol kcorr_z "0+0.413674*pow(z, 1)*pow(r_z, 0)-4.128458*pow(z, 1)*pow(r_z, 1)+13.264598*pow(z, 1)*pow(r_z, 2)-9.057702*pow(z, 1)*pow(r_z, 3)-7.995974*pow(z, 2)*pow(r_z, 0)-6.467990*pow(z, 2)*pow(r_z, 1)-16.459751*pow(z, 2)*pow(r_z, 2)+16.344189*pow(z, 2)*pow(r_z, 3)+45.932457*pow(z, 3)*pow(r_z, 0)+84.022560*pow(z, 3)*pow(r_z, 1)-11.926839*pow(z, 3)*pow(r_z, 2)-189.230835*pow(z, 4)*pow(r_z, 0)-113.353912*pow(z, 4)*pow(r_z, 1)+246.450256*pow(z, 5)*pow(r_z, 0)"; addcol kcorr_Y "0-3.137063*pow(z, 1)*pow(Y_H, 0)+8.748210*pow(z, 1)*pow(Y_H, 1)+2.936402*pow(z, 1)*pow(Y_H, 2)-11.170932*pow(z, 1)*pow(Y_H, 3)+23.819769*pow(z, 2)*pow(Y_H, 0)-59.077156*pow(z, 2)*pow(Y_H, 1)+18.851744*pow(z, 2)*pow(Y_H, 2)+17.777597*pow(z, 2)*pow(Y_H, 3)-96.471672*pow(z, 3)*pow(Y_H, 0)+178.633133*pow(z, 3)*pow(Y_H, 1)-38.253731*pow(z, 3)*pow(Y_H, 2)+124.799034*pow(z, 4)*pow(Y_H, 0)-186.443405*pow(z, 4)*pow(Y_H, 1)-12.444495*pow(z, 5)*pow(Y_H, 0)"; addcol kcorr_J "0-0.286459*pow(z, 1)*pow(J_K, 0)+2.827541*pow(z, 1)*pow(J_K, 1)-0.445788*pow(z, 1)*pow(J_K, 2)+0.342958*pow(z, 1)*pow(J_K, 3)-7.434520*pow(z, 2)*pow(J_K, 0)-43.319092*pow(z, 2)*pow(J_K, 1)-15.636494*pow(z, 2)*pow(J_K, 2)-2.740727*pow(z, 2)*pow(J_K, 3)+72.026680*pow(z, 3)*pow(J_K, 0)+262.249664*pow(z, 3)*pow(J_K, 1)+27.496885*pow(z, 3)*pow(J_K, 2)-323.571259*pow(z, 4)*pow(J_K, 0)-340.672150*pow(z, 4)*pow(J_K, 1)+395.254425*pow(z, 5)*pow(J_K, 0)"; addcol kcorr_H "0+0.359437*pow(z, 1)*pow(H_K, 0)+5.440102*pow(z, 1)*pow(H_K, 1)-3.286579*pow(z, 1)*pow(H_K, 2)-18.617926*pow(z, 1)*pow(H_K, 3)-6.899263*pow(z, 2)*pow(H_K, 0)-158.156555*pow(z, 2)*pow(H_K, 1)-24.999006*pow(z, 2)*pow(H_K, 2)+359.200012*pow(z, 2)*pow(H_K, 3)+53.301483*pow(z, 3)*pow(H_K, 0)+1283.672363*pow(z, 3)*pow(H_K, 1)-159.837448*pow(z, 3)*pow(H_K, 2)-1539.324707*pow(z, 3)*pow(H_K, 3)-236.131577*pow(z, 4)*pow(H_K, 0)-4242.764160*pow(z, 4)*pow(H_K, 1)+1433.016479*pow(z, 4)*pow(H_K, 2)+1853.301758*pow(z, 4)*pow(H_K, 3)+502.851013*pow(z, 5)*pow(H_K, 0)+6217.591309*pow(z, 5)*pow(H_K, 1)-2105.408447*pow(z, 5)*pow(H_K, 2)-527.163269*pow(z, 6)*pow(H_K, 0)-3288.293945*pow(z, 6)*pow(H_K, 1)+241.101410*pow(z, 7)*pow(H_K, 0)"; addcol kcorr_K "0-2.732966*pow(z, 1)*pow(J_K, 0)+1.711451*pow(z, 1)*pow(J_K, 1)-3.227258*pow(z, 1)*pow(J_K, 2)+1.290244*pow(z, 1)*pow(J_K, 3)+11.497313*pow(z, 2)*pow(J_K, 0)-16.103697*pow(z, 2)*pow(J_K, 1)+3.783451*pow(z, 2)*pow(J_K, 2)-0.315761*pow(z, 2)*pow(J_K, 3)-59.601894*pow(z, 3)*pow(J_K, 0)+91.185181*pow(z, 3)*pow(J_K, 1)-9.114658*pow(z, 3)*pow(J_K, 2)+97.592857*pow(z, 4)*pow(J_K, 0)-104.714813*pow(z, 4)*pow(J_K, 1)-36.927883*pow(z, 5)*pow(J_K, 0)";' > galaxies_11bands_0007_003_kcorr.fits