** 10 June 2002 on ev6 ls -ld /*/*/*F22*I* /*/*/*/*F22*I* | grep -v .fits> fred rm fred2 awk '{print "du -s" ,$9,">>fred2"}' fred > fredx source fredx awk '{print $1}' fred2 > fred3 paste fred3 fred | awk '$1>1000' >relfiles clean (edit by hand) (cutting the Decoupe stuff out: it is all stacked images) 995786 drwxr-xr-x 2 mellier crap 1536 Apr 28 18:14 /evsix2/mellier/DecoupeMarch2001/F22P5-I/ 995786 drwxr-xr-x 2 mellier crap 1536 Apr 28 18:14 /i2/mellier/DecoupeMarch2001/F22P5-I/ 1014370 drwxr-xr-x 2 mellier crap 1536 May 1 11:42 /evsix2/mellier/DecoupeMarch2001/F22P7-I/ 1014370 drwxr-xr-x 2 mellier core 1536 May 1 11:42 /i2/mellier/DecoupeMarch2001/F22P7-I/ 1013546 drwxr-xr-x 2 mellier crap 1536 May 1 12:06 /evsix2/mellier/DecoupeMarch2001/F22P8-I/ 1013546 drwxr-xr-x 2 mellier crap 1536 May 1 12:06 /i2/mellier/DecoupeMarch2001/F22P8-I/ 1012330 drwxr-xr-x 2 mellier crap 1536 May 1 12:31 /evsix2/mellier/DecoupeMarch2001/F22P9-I/ 1012330 drwxr-xr-x 2 mellier crap 1536 May 1 12:31 /i2/mellier/DecoupeMarch2001/F22P9-I/ 1019898 drwxr-xr-x 2 mellier crap 1536 May 2 09:06 /evsix2/mellier/DecoupeMarch2001/F22P10-I/ 1019898 drwxr-xr-x 2 mellier crap 1536 May 2 09:06 /i2/mellier/DecoupeMarch2001/F22P10-I/ 1011306 drwxr-xr-x 2 mellier core 1536 May 1 10:52 /i2/mellier/DecoupeMarch2001/F22P11-I/ 1011306 drwxr-xr-x 2 mellier core 1536 May 1 10:52 /evsix2/mellier/DecoupeMarch2001/F22P11-I/ 1012010 drwxr-xr-x 2 mellier crap 1536 May 1 11:17 /i2/mellier/DecoupeMarch2001/F22P12-I/ 1012010 drwxr-xr-x 2 mellier crap 1536 May 1 11:17 /evsix2/mellier/DecoupeMarch2001/F22P12-I/ The good stuff: 5331340 drwxr-xr-x 3 mellier core 512 Jan 16 2000 /samplix3/mellier/F22P3-I-May99/ 5360250 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/ 5343722 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/ 6008898 drwxr-xr-x 4 mellier core 512 Jan 14 2000 /samplix3/mellier/F22P4-I-Nov99/ 5350106 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/ 7519582 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/ 3936970 drwxr-xr-x 3 mellier system 512 Jun 21 2000 /samplix3/mellier/F22P6b-I-Apr00-2/ 7509296 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/ 4724364 drwxr-xr-x 3 mellier system 512 Jun 21 2000 /samplix3/mellier/F22P11b-I-Apr00-2/ Note: mellier and dantel files *not* necessarily the same. P1 Lost due to moon P2 Lost due to moon P3 /samplix3/mellier/F22P3-I-May99/notcoadded/*MDFG.fits.gz (5*12) /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded/*MOBDFG.fits.gz (5*12) P4 /samplix3/mellier/F22P4-I-Nov99/notcoadded/*MDFG.fits.gz (5*12) /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded/*MOBDFG.fits.gz (5*12) P5 /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded/*MOBDFG.fits.gz (5*12) P6 /samplix3/mellier/F22P6b-I-Apr00-2/notcoadded/flatnotcoadded/*MODFOfl.fits (5*12) /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded/*MOBDFG.fits.gz (7*12) P7 /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded/*MOBDFG.fits.gz (7*12) P8 nada P9 nada P10 nada P11 /samplix3/mellier/F22P11b-I-Apr00-2/notcoadded/flatnotcoadded/*MODFOfl.fits (6*12) P12 nada New note added: 22 June 2002: Somwhere in here you should run /samplix1/gwyn/findim.pl In theory: P1: 1.485930/0720/0.9/May99(01)[Lost: Moon] 2.485931/0720/0.9/May99(01)[Lost: Moon] 3.485932/0720/1.0/May99(01)[Lost: Moon] 4.485933/0720/0.9/May99(01)[Lost: Moon] 5.485934/0720/0.9/May99(01)[Lost: Moon] 6.486004/0720/1.2/May99(01)[Lost: Moon] 7.486005/0720/1.2/May99(01)[Lost: Moon] P2: 1.486006/0720/1.0/May99(02)[Lost: Moon] 2.486007/0720/0.9/May99(02)[Lost: Moon] 3.486008/0720/0.9/May99(02)[Lost: Moon] 4.486009/0720/0.9/May99(02)[Lost: Moon] 5.486090/0720/0.8/May99(02)[Lost: Moon] 6.486091/0720/0.7/May99(02)[Lost: Moon] 7.486092/0720/0.8/May99(02)[Lost: Moon] P3: 1.557282/0720/0.80/Nov00(03b)[ignore] 2.550361/0720/0.60/Sept00(03b) 3.550362/0720/0.60/Sept00(03b) 4.550363/0720/0.60/Sept00(03b) 5.550364/0720/0.75/Sept00(03b)[Lost: elongated] 6.550366/0720/0.65/Sept00(03b) 7.486093/0720/0.7/May99(03)[Lost: Moon] 8.486094/0720/0.6/May99(03)[Lost: Moon] 9.486370/0720/1.1/May99(03) 10.486371/0720/1.1/May99(03) 11.486372/0720/1.1/May99(03) 12.486373/0720/1.1/May99(03) 13.486374/0720/1.1/May99(03) P4: 1.550367/0720/0.75/Sept00(04b) 2.550368/0720/0.75/Sept00(04b) 3.550369/0720/0.80/Sept00(04b) 4.550373/0720/0.60/Sept00(04b) 5.550374/0720/0.60/Sept00(04b) 6.508265/0720/0.7/Nov99(04) 7.508266/0720/0.8/Nov99(04) 8.508267/0720/0.8/Nov99(04) 9.508268/0720/0.8/Nov99(04) 10.508269/0720/0.8/Nov99(04) P5: 1.550375/0720/0.55/Sept00(05b) 2.550376/0720/0.55/Sept00(05b) 3.550377/0720/0.60/Sept00(05b) 4.550378/0720/0.60/Sept00(05b) 5.550379/0720/0.65/Sept00(05b) P6: 1.531288/0720/0.8/Apr00-2(06b) 2.531289/0720/0.7/Apr00-2(06b) 3.531290/0720/0.8/Apr00-2(06b) 4.531182/0720/1.0/Apr00-2(06b) 5.531183/0720/1.0/Apr00-2(06b) 6.551371/0720/0.7/Oct00-1(06b) 7.551372/0720/0.6/Oct00-1(06b) 8.551373/0720/0.6/Oct00-1(06b) 9.551374/0720/0.6/Oct00-1(06b) 10.551375/0720/0.7/Oct00-1(06b) 11.551376/0720/0.7/Oct00-1(06b) 12.551377/0720/0.8/Oct00-1(06b) P7: 1.551379/0720/0.6/Oct00-1(07b) 2.551380/0720/0.6/Oct00-1(07b) 3.551381/0720/0.65/Oct00-1(07b) 4.551382/0720/0.65/Oct00-1(07b) 5.551383/0720/0.65/Oct00-1(07b) 6.551384/0720/0.6/Oct00-1(07b) 7.551385/0720/0.8/Oct00-1(07b)[Double images] P8: 1.552797/0720/0.6/Oct00-3(8b)[Cirrus] 2.552798/0720/0.6/Oct00-3(8b)[Cirrus] 3.552799/0720/0.5/Oct00-3(8b)[Cirrus] 4.552800/0720/0.5/Oct00-3(8b)[Cirrus] 5.552801/0720/0.5/Oct00-3(8b)[Cirrus] 6.552855/0720/0.8/Oct00-3(8b)[Photometric] P9: 1.552802/0720/0.5/Oct00-3(09b) 2.552803/0720/0.5/Oct00-3(09b) 3.552804/0720/0.5/Oct00-3(09b) 4.552805/0720/0.5/Oct00-3(09b) 5.552806/0720/0.5/Oct00-3(09b) 6.552856/0720/1.0/Oct00-3(09b)[Photometric] P10: 1.552850/0720/0.8/Oct00-3(10b) 2.552851/0720/0.6/Oct00-3(10b) 3.552852/0720/0.6/Oct00-3(10b) 4.552853/0720/0.7/Oct00-3(10b) 5.552854/0720/0.8/Oct00-3(10b) P11: 1.557283/0720/0.9/Nov00(11b)[ignore] 2.552953/0720/0.6/Oct00-3(11b) 3.552954/0720/0.6/Oct00-3(11b) 4.552858/0720/0.8/Oct00-3(11b)[Photometric] 5.552859/0720/0.8/Oct00-3(11b)[Photometric] 6.552860/0720/0.8/Oct00-3(11b)[Photometric] 7.531630/0720/1.0/Apr00-2(11b) 8.531631/0720/0.9/Apr00-2(11b) 9.531632/0720/0.8/Apr00-2(11b) 10.531633/0720/0.8/Apr00-2(11b) 11.531477/0720/1.3/Apr00-2(11b)[Lost: Clouds] 12.531478/0720/1.5/Apr00-2(11b)[Lost: Clouds] P12: 1.552955/0720/0.6/Oct00-3(12b) 2.552956/0720/0.7/Oct00-3(12b) 3.552957/0720/0.7/Oct00-3(12b) 4.552958/0720/0.8/Oct00-3(12b) 5.552959/0720/1.0/Oct00-3(12b) /samplix1/dantel/DATAprocessed-avril2000-2.DEF/astrometry/F22P11b-I-Apr00-2: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/astrometry/F22P6b-I-Apr00-2: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded/flatnotcoadded: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/cut: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/flatnotcoadded: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF: /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded/flatnotcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded/flatnotcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded/flatnotcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded/flatnotcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded/flatnotcoadded/bad: /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded/flatnotcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/astrometry/oct2000/F22P6b-I-Oct00: /samplix1/dantel/cfh12k_sept_oct_2000/astrometry/oct2000/F22P7b-I-Oct00: /samplix1/dantel/cfh12k_sept_oct_2000/astrometry/sept2000/F22P3b-I-Sep00: /samplix1/dantel/cfh12k_sept_oct_2000/astrometry/sept2000/F22P4b-I-Sep00: /samplix1/dantel/cfh12k_sept_oct_2000/astrometry/sept2000/F22P5b-I-Sep00: /samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I: /samplix3/dantel/cfh12k_oct2000-MOBDF: /samplix3/dantel/cfh12k_sep2000-MOBDF: /samplix3/dantel/coadd-oct2000_2_3/F22/F22P10b-I-Oct00_3: /samplix3/dantel/coadd-oct2000_2_3/F22/F22P11b-I-Oct00_3: /samplix3/dantel/coadd-oct2000_2_3/F22/F22P12b-I-Oct00_3: /samplix3/dantel/coadd-oct2000_2_3/F22/F22P8b-I-Oct00_3: /samplix3/dantel/coadd-oct2000_2_3/F22/F22P9b-I-Oct00_3: /samplix3/mellier/F22P11b-I-Apr00-2/notcoadded/flatnotcoadded: /samplix3/mellier/F22P3-I-May99/notcoadded/flatnotcoadded: /samplix3/mellier/F22P3-I-May99/notcoadded: /samplix3/mellier/F22P4-I-Nov99/notcoadded/flatnotcoadded: /samplix3/mellier/F22P4-I-Nov99/notcoadded: /samplix3/mellier/F22P6b-I-Apr00-2/notcoadded/flatnotcoadded: /u2/morin/terapix/dbterapix/v2/src/dbtest/tpx/data: (less the *fl.fits) /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/cut: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded: /samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF: /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded: /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded: /samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I: /samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I: /samplix3/dantel/cfh12k_oct2000-MOBDF: /samplix3/dantel/cfh12k_sep2000-MOBDF: /samplix3/mellier/F22P3-I-May99/notcoadded: /samplix3/mellier/F22P4-I-Nov99/notcoadded: /u2/morin/terapix/dbterapix/v2/src/dbtest/tpx/data: ** 11 June 2002 Let's go with /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/cut( = cropped) /samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF (includes other images) /samplix3/dantel/cfh12k_oct2000-MOBDF (other images) /samplix3/dantel/cfh12k_sep2000-MOBDF (other images) /u2/morin/terapix/dbterapix/v2/src/dbtest/tpx/data (No good: just headers) /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded /samplix3/mellier/F22P3-I-May99/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded /samplix3/mellier/F22P4-I-Nov99/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded /samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded /samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I (Fav's) /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded All available in /samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF According to Seb, don't use the Fav's... But according to Mirelle Dantel, do use them... Final list: cat >F22I.dirlist /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded /samplix3/mellier/F22P3-I-May99/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded /samplix3/mellier/F22P4-I-Nov99/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded /samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded /samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded /samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF run: /samplix/gwyn/track/linkall.pl Put the links in F22I.links 552959o11oct003p12.fits ^ ^^ ^ ^ ^^^ | || | | |||____ pointing number: 2 digits 03 to 12 (should be 01 to 16) | || | | ||_____ the letter 'p' | || | | |______ date: 1 digit for sub-run number | || | |________ date: 2 digits for year | || |___________ date: | ||_____________ CCD number: 2 digits 00 to 11 | |______________ the letter 'o' |____________________ the exposure number: 6 digits 3 letters for month, 2 digits for year 1 digit switch to pix cd /raid mkdir gwyn cd gwyn mkdir F22 cd F22 mkdir F22I cd F22I mkdir data cd data ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22I get * bye gzip -d *gz ** 13 June 2002 Redid most the above with the Fav images. in iraf: cd F22I/data !ls *.fits | perl -nle 'print $_."\[1:2048,1:4096\]"' >listin !ls *.fits | perl -nle 's/\.fits/C.fits/ ; print $_' >listout imcopy @listin @listout ** 14 June 2002 instead of all that: rm imlist* stack.pl stackall ** 17 June 2002 ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' mv mask* .. mv imlist* .. cd .. sex.pl ** 19 June 2002 Bug in linkall.pl means I have to do a redo from start. On ev6: cd /samplix1/gwyn/track linkall.pl On pix1: cd /raid/gwyn/F22/F22I/data ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22I get * bye gzip -d *gz mv * back ln -s back/*.fits . delhead CONSWV *.fits rm imlist* imlist.pl stackall 2551.540u 197.060s 53:57.22 84.9% 0+0k 0+0io 245pf+0w rm *p??.fits ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd .. cp data/mask*.box . imlist.pl sex.pl >sex.log (both imlist.pl and sex.pl are now pix1-based, not lasa28-based) ** 20 June 2002 maskch.pl on lasa28: scat -r -4000 -n 1000000 -d -c ua2 334.5 .33333 | awk '{print $2,$3,$4,$5}' >USNO also using getout.f to get the outlines this is checked with sm.1 on Morgan: /home/gwyn/red/F22/F22I/inithead Check all the masks with the help of Wayne & Schuster (1963) 00 images: most have a big box on the upper left that should not be included in the final images, but could be used in the astrometry part. Also there are a lot very faint baddish columns 531182o04apr002p06.fits shows a flatfield gradient 0 486370o04may990p03.fits (which looks like the normal CCD00) has the faint bad cols and the big bad square. CCD05 problems with tvm... 486370o06may990p03.fits streaky stars 486370o07may990p03.fits strange effects: examine image without boxes 486370o10may990p03.fits lots of extensions CCD00 leave as is: DONE the blurry area in left side. might vary from exposure to exposure CCD01 leave as is: DONE gaps in the bad columns in one exposures CCD02 No problems DONE CCD03 No problems DONE CCD04 Leave as is DONE Blurry area in May99 CCD05 homogenize all but 990 DONE may990 completely different CCD06 All but may990 : homogenize: DONE may990 fine as is. CCD07 homogenize all by cating boxes DONE CCD08 homogenize all: DONE (could have been ignored. The correction is the strip along the bottom. CCD09 homogenize all: DONE CCD10 No problems DONE CCD11 may990: added marginal bad column: keep it DONE ** June 21 2002 Script for dealing with CCD09 disp 486370o09may990p03.fits[700:850,3900:4096] 1 erase+ xcen=.25 ycen=.3333 tvmark 1 '' commands=add09may990 col=205 tvmark 1 '' commands=add09may990 col=204 disp 531182o09apr002p06.fits[700:850,3900:4096] 1 erase- xcen=.50 ycen=.3333 tvmark 1 '' commands=add09may990 col=205 tvmark 1 '' commands=add09apr002 col=204 disp 551371o09oct001p06.fits[700:850,3900:4096] 1 erase- xcen=.75 ycen=.3333 tvmark 1 '' commands=add09may990 col=205 tvmark 1 '' commands=add09oct001 col=204 disp 508265o09nov990p04.fits[700:850,3900:4096] 1 erase- xcen=.25 ycen=.6666 tvmark 1 '' commands=add09may990 col=205 tvmark 1 '' commands=add09nov990 col=204 disp 550361o09sep000p03.fits[700:850,3900:4096] 1 erase- xcen=.50 ycen=.6666 tvmark 1 '' commands=add09may990 col=205 tvmark 1 '' commands=add09sep000 col=204 disp 552797o09oct003p08.fits[700:850,3900:4096] 1 erase- xcen=.75 ycen=.6666 tvmark 1 '' commands=add09may990 col=205 tvmark 1 '' commands=add09oct003 col=204 Script for dealing with CCD08 disp 486370o08may990p03.fits[1:100,1:4096] 1 erase+ xcen=.15 ycen=.5 ymag=.2 tvmark 1 '' commands=mask08may990.tvm col=204 tvmark 1 '' commands=add08may990 col=205 disp 531182o08apr002p06.fits[1:100,1:4096] 1 erase- xcen=.30 ycen=.5 ymag=.2 tvmark 1 '' commands=mask08may990.tvm col=204 tvmark 1 '' commands=add08apr002 col=205 disp 551371o08oct001p06.fits[1:100,1:4096] 1 erase- xcen=.45 ycen=.5 ymag=.2 tvmark 1 '' commands=mask08may990.tvm col=204 tvmark 1 '' commands=add08oct001 col=205 disp 508265o08nov990p04.fits[1:100,1:4096] 1 erase- xcen=.60 ycen=.5 ymag=.2 tvmark 1 '' commands=mask08may990.tvm col=204 tvmark 1 '' commands=add08nov990 col=205 disp 550361o08sep000p03.fits[1:100,1:4096] 1 erase- xcen=.75 ycen=.5 ymag=.2 tvmark 1 '' commands=mask08may990.tvm col=204 tvmark 1 '' commands=add08sep000 col=205 disp 552797o08oct003p08.fits[1:100,1:4096] 1 erase- xcen=.90 ycen=.5 ymag=.2 tvmark 1 '' commands=mask08may990.tvm col=204 tvmark 1 '' commands=add08oct003 col=205 script for dealing with CCD07 disp 486370o07may990p03.fits[1550:1730,100:200] 1 erase+ xcen=.25 ycen=.3333 tvmark 1 '' commands=fred col=205 tvmark 1 '' commands=add07may990 col=204 disp 531182o07apr002p06.fits[1550:1730,100:200] 1 erase- xcen=.50 ycen=.3333 tvmark 1 '' commands=fred col=205 tvmark 1 '' commands=add07apr002 col=204 disp 551371o07oct001p06.fits[1550:1730,100:200] 1 erase- xcen=.75 ycen=.3333 tvmark 1 '' commands=fred col=205 tvmark 1 '' commands=add07oct001 col=204 disp 508265o07nov990p04.fits[1550:1730,100:200] 1 erase- xcen=.25 ycen=.6666 tvmark 1 '' commands=fred col=205 tvmark 1 '' commands=add07nov990 col=204 disp 550361o07sep000p03.fits[1550:1730,100:200] 1 erase- xcen=.50 ycen=.6666 tvmark 1 '' commands=fred col=205 tvmark 1 '' commands=add07sep000 col=204 disp 552797o07oct003p08.fits[1550:1730,100:200] 1 erase- xcen=.75 ycen=.6666 tvmark 1 '' commands=fred col=205 tvmark 1 '' commands=add07oct003 col=204 script for dealing with CCD01 disp 486370o00may990p03.fits[140:575,1400:4096] 1 erase+ xcen=.15 ycen=.5 tvmark 1 '' commands=mask00may990.tvm col=204 tvmark 1 '' commands=add00may990 col=205 disp 531182o00apr002p06.fits[140:575,1400:4096] 1 erase- xcen=.30 ycen=.5 tvmark 1 '' commands=mask00apr002.tvm col=204 tvmark 1 '' commands=add00apr002 col=205 disp 551371o00oct001p06.fits[140:575,1400:4096] 1 erase- xcen=.45 ycen=.5 tvmark 1 '' commands=mask00oct001.tvm col=204 tvmark 1 '' commands=add00oct001 col=205 disp 508265o00nov990p04.fits[140:575,1400:4096] 1 erase- xcen=.60 ycen=.5 tvmark 1 '' commands=mask00nov990.tvm col=204 tvmark 1 '' commands=add00nov990 col=205 disp 550361o00sep000p03.fits[140:575,1400:4096] 1 erase- xcen=.75 ycen=.5 tvmark 1 '' commands=mask00sep000.tvm col=204 tvmark 1 '' commands=add00sep000 col=205 disp 552797o00oct003p08.fits[140:575,1400:4096] 1 erase- xcen=.90 ycen=.5 tvmark 1 '' commands=mask00oct003.tvm col=204 tvmark 1 '' commands=add00oct003 col=205 look at the .gif's of the masks: maskapr002.gif maskmay990.gif masknov990.gif maskoct001.gif maskoct003.gif masksep000.gif maskpic: takes the mask*tvm files (as specified by imlistall) and makes cute .pgm files of the masks. mog.pl: takes these pgm and assembles them into a big gif, one per epoch. back to pix1 cd /raid/gwyn/F22/F22I/data ncftp -u gwyn lasa28.astrsp-mrs.fr cd /home/gwyn/red/F22/F22I/maskcheck/ get add* bye imlist.pl stackall >stacklog ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd .. cp data/mask*.box . imlist.pl sex.pl >sex.log dofit1.pl dofit1.log ** 22 June 2002 531478o??apr002p11.fits is completely blank it shall be discarded But it should never have been read: Ditto with the following: There was a problem with linkall.pl It has been fixed. sep000 03 5.550364/0720/0.75/Sept00(03b)[Lost: elongated] oct001 07 7.551385/0720/0.8/Oct00-1(07b)[Double images] oct003 08 1.552797/0720/0.6/Oct00-3(8b)[Cirrus] oct003 08 2.552798/0720/0.6/Oct00-3(8b)[Cirrus] oct003 08 3.552799/0720/0.5/Oct00-3(8b)[Cirrus] oct003 08 4.552800/0720/0.5/Oct00-3(8b)[Cirrus] oct003 08 5.552801/0720/0.5/Oct00-3(8b)[Cirrus] apr002 11 11.531477/0720/1.3/Apr00-2(11b)[Lost: Clouds] apr002 11 12.531478/0720/1.5/Apr00-2(11b)[Lost: Clouds] These images have been removed... The process has been restarted though: backup of linkall.pl is on Morgan: /home/gwyn/red/scripts/linkall.pl run: rm 550364* rm 551385* rm 531477* rm 531478* Keeping the cirrus images (reluctantly) on ev6 cd /samplix1/gwyn/track linkall.pl On pix1: cd /raid/gwyn/F22/F22I/data ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22I get * bye gzip -d *gz mv * back ln -s back/*.fits . delhead CONSWV *.fits rm imlist* imlist.pl stackall ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd .. ** 23 June 2002 cd /raid/gwyn/F22/F22I cp data/mask*box imlist.pl sex.pl >sex.log dofit1.pl >dofit1.log check results... checkdo.pl dofit1.log worst centering errors: 486371o06may990p03.fits 75 -1 65 53 1.255 486373o06may990p03.fits 76 -1 64 48 1.241 486374o06may990p03.fits 72 -1 60 49 1.233 486372o06may990p03.fits 76 -1 65 52 1.217 486370o06may990p03.fits 75 -1 64 51 1.216 552954o06oct003p11.fits 75 -1 66 67 1.044 fewest initial matches: 550376o05sep000p05.fits 50 -1 42 43 0.645 550375o05sep000p05.fits 51 -1 43 45 0.717 550377o05sep000p05.fits 52 -1 43 45 0.694 550361o09sep000p03.fits 52 -1 30 30 0.570 550362o09sep000p03.fits 53 -1 32 32 0.557 550362o05sep000p03.fits 53 -1 50 50 0.781 fewest final matches: 550361o09sep000p03.fits 52 -1 30 30 0.570 550362o09sep000p03.fits 53 -1 32 32 0.557 552955o02oct003p12.fits 54 -1 33 33 0.674 550377o09sep000p05.fits 55 -1 33 34 0.725 550376o09sep000p05.fits 55 -1 33 34 0.654 552956o09oct003p12.fits 57 -1 35 35 0.529 Not too shabby. Results of inspecting images: 486370o00may990p03.fits multiple dots 486371o00may990p03.fits multiple dots 486372o00may990p03.fits multiple dots 486373o00may990p03.fits multiple dots 486374o00may990p03.fits multiple dots generally, the swath of crap is of variable width 551372o00oct001p06.fits streak 551383o00oct001p06.fits streak 552854o00oct003p10.fits streak 552958o00oct003p12.fits faint streak 551372o01oct001p06.fits streak 551383o01oct001p07.fits streak 552854o01oct003p10.fits streak 551372o02oct001p06.fits streak 551383o02oct001p06.fits streak 552854o02oct003p10.fits streak 552956o02oct003p12.fits faint streak 531288o01apr002p06.fits faint streak 551372o03oct001p06.fits streak 531288o04apr002p06.fits bad cosmic ray 550361o04sep000p03.fits streak 550362o04sep000p03.fits thread? 550363o04sep000p03.fits thread? 550366o04sep000p03.fits bad cosmic ray 550378o04sep000p05.fits scattered light 552806o04oct003p09.fits bad cosmic ray 552850o04oct003p10.fits tumbling satellite 552851o04oct003p10.fits bumpy cosmic rays 552856o04oct003p09.fits bad cosmic ray 552957o04oct003p12.fits big streak 486370o05may990p03.fits thread? 550379o05sep000p05.fits internal reflection 551372o05oct001p06.fits streak 551383o05oct001p07.fits streak 552801o05oct003p08.fits streak 552853o05oct003p10.fits streak 552954o05oct003p11.fits streak 486370o06may990p03.fits stars streaked 486371o06may990p03.fits stars streaked 486372o06may990p03.fits stars streaked 486373o06may990p03.fits stars streaked 486374o06may990p03.fits stars streaked 531183o06apr002p06.fits scattered light 531290o06apr002p06.fits scattered light 552802o06oct003p09.fits bright star just off frame 552803o06oct003p09.fits bright star just off frame 552804o06oct003p09.fits bright star just off frame 552805o06oct003p09.fits bright star just off frame 552806o06oct003p09.fits bright star just off frame 552855o06oct003p08.fits bright star just off frame 552856o06oct003p09.fits bright star just off frame 486371o07may990p03.fits slightly trailed images 486373o07may990p03.fits slightly trailed images some other very slightly trailed image in may... 531183o07apr002p06.fits two "asteroids" 552802o07oct003p09.fits very bright star 552803o07oct003p09.fits very bright star 552804o07oct003p09.fits very bright star 552805o07oct003p09.fits very bright star 552806o07oct003p09.fits very bright star 552856o07oct003p09.fits very bright star 486370o08may990p03.fits very slightly trailed images 486371o08may990p03.fits very slightly trailed images 486372o08may990p03.fits slightly trailed images 486373o08may990p03.fits slightly trailed images 486374o08may990p03.fits slightly trailed images 550368o08sep000p04.fits bright star off to one side. 552806o08oct003p09.fits two horizontal lines 486370o09may990p03.fits streak on lower right corner 486371o09may990p03.fits very slightly trailed images 486373o09may990p03.fits slightly trailed images 531630o09apr002p11.fits bright star 531631o09apr002p11.fits bright star 531632o09apr002p11.fits bright star 531633o09apr002p11.fits bright star 550367o09sep000p04.fits bright star off edge 550368o09sep000p04.fits bright star off edge 550369o09sep000p04.fits bright star off edge 550373o09sep000p04.fits bright star off edge 550374o09sep000p04.fits bright star off edge 552956o09oct003p12.fits faint streak 486370o10may990o03.fits bright streak 550367o10sep000p04.fits scattered light at bottom 550368o10sep000p04.fits scattered light at bottom 550369o10sep000p04.fits scattered light at bottom 550373o10sep000p04.fits scattered light at bottom 550374o10sep000p04.fits scattered light at bottom 552957o10oct003p12.fits bright streak may990o11 some more bits might need masking 486370o11may990o03.fits bright streak *o11apr002p11.fits don't look good 531632o11apr002p11.fits doubled bright star 550ls 361o11sep000p03.fits faint streaks 550369o11sep000p04.fits ghostly doubled stars 550375o11sep000p05.fits ghostly residual stars 551372o11oct001p06.fits streak 551373o11oct001p06.fits bright streak 552954o11oct003p11.fits streak And looking at the results of the first fit: 508269o03nov990p04.fits corner star unmatched 551371o06oct001p06.fits missing stars in bottom left 551372o06oct001p06.fits a few stars in bottom left 551373o06oct001p06.fits missing stars in bottom left 551374o06oct001p06.fits missing stars in bottom left 551375o06oct001p06.fits missing stars in bottom left 551376o06oct001p06.fits missing stars in bottom left 551377o06oct001p06.fits missing stars in bottom left ** 24 June 2002 Photometry: I am assuming: Y-date Header-date Table-date Zero-point may990 1999-05-13 May 99 25.802 nov990 1999-11-09 Nov 99 25.894 apr002 2000-04-30 Apr00-2 26.3386 2000-05-01 2000-05-03 sep000 2000-09-22 Sept+Oct00 26.2906 oct001 2000-10-03 Sept+Oct00 26.2906 oct003 2000-10-23 Oct00-3 26.0078 2000-10-24 2000-10-25 dofit2.pl is being maintained on pix1 with this version of getzp sub getzp { my $head=shift @_; my $expkey ='EXPTIME'; my $airkey ='AIRMASS'; my $datekey ='DATE-OBS'; my $colour='I'; my %fact= qw (B 0.20 V 0.12 R 0.10 I 0.08); my %zpoint = qw { I1999-05-13 25.802 I1999-11-09 25.894 I2000-04-30 26.3386 I2000-05-01 26.3386 I2000-05-03 26.3386 I2000-09-22 26.2906 I2000-10-03 26.2906 I2000-10-23 26.0078 I2000-10-24 26.0078 I2000-10-25 26.0078 }; open HEAD,$head; while () { chomp; ($key,$value)= split /=/; ($value) = split /\//,$value; $key=~s/\s//g; $value=~s/\s//g; $value=~s/\'//g; if ($key eq $expkey ) {$exptime = $value} if ($key eq $airkey ) {$airmass = $value} if ($key eq $datekey) {$date = $value} } close HEAD; $fact=$fact{$colour}; unless (defined($zpoint{$colour.$date})) { print "ERROR: couldn't find zero-point\n"; print "$head $colour $date\n"; die; } $zpoint =$zpoint{$colour.$date} + 2.5*&log10($exptime) - $airmass*$fact; } ** 25 June 2002 Seb says: Oct 00-2 No calibrations because of Supernovae dudes Apr-00-2 first reduction run of Seb: Probably not good. ** 28 June 2002 checked gradients accross images. Get the minibackground images from SExtractor. use /home/virmos/gwyn/scripts/back.pl on pix1 Run backch.f on lasa28.astrsp-mrs.fr on all the images. backch.f prints out image_name median peak-to-peak horizontal_gradient vertical_gradient diagonal_gradient (sw-ne) diagonal_gradient (se-nw) The gradients should not be more than 2% The peak to peak is actually affected by bright stars more than anything else. All images have less than 2% gradients. Some have high medians: 531633o??apr002p11.fits 40000+ It shall be tossed... ** 1 July 2002 Checking the gain: run /home/virmos/gwyn/scripts/noise.pl (it runs SExtractor, throws away the catlogs and outputs a object-less images with a mb.fits suffix) run /home/virmos/gwyn/pg/stddev.f ls *mb.fits | perl -nle 'system "stddev $_"' | grep -v read > stddev.log This gives you the std. dev. by examining 10x10 pixel sub-rasters, tossing those contaminated by objects or bad columns, and averaging all the std-devs. Using a small sub-raster means even if there is fringing, you are probably either on a fringe or completely off it. The average will change slightly, but the std-dev won't. Copy back to lasa28 ~/red/F22/F22I/backcheck run sm.1 and sm.2 to get: I seem to have deleted all the add-files I don't think they got included in the last stackall. Redo from start: cp stackall >stackall.log ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd .. cd /raid/gwyn/F22/F22I cp data/mask* . imlist.pl sex.pl >sex.log Problems with one image. Turns out sex.pl got cut off in mid stream. Not enough stars --> not a good fit. dofit1.pl dofit1.log mv *06may* toss dofit2.pl dofit2nomay06.log crossref dofit3.pl dofit3nomay06.log ** 9 July 2002 I'm supposed to be on vacation set PHOT_APERTURES 20 # MAG_APER aperture diameter(s) in pixels sex.pl dofit.pl USNO 1 dofit1.log checkdo.pl dofit1.log ls *.fits | perl -nle 'system "uphead $_ fead"' dofit.pl USNO 2 dofit2.log checkdo.pl dofit2.log 550376o03sep000p05.fits 74 2 52 65 0.727 508268o00nov990p04.fits 72 2 45 54 0.709 550377o03sep000p05.fits 75 2 49 63 0.701 508268o05nov990p04.fits 70 2 55 62 0.698 550376o04sep000p05.fits 57 2 42 52 0.690 486373o10may990p03.fits 61 2 46 56 0.689 fewest initial matches: 550366o05sep000p03.fits 49 2 42 48 0.574 550377o05sep000p05.fits 50 2 38 45 0.508 550376o05sep000p05.fits 50 2 38 42 0.498 550375o05sep000p05.fits 51 2 38 45 0.566 550362o09sep000p03.fits 51 2 30 32 0.390 550362o05sep000p03.fits 51 2 43 48 0.502 fewest final matches: 550361o09sep000p03.fits 52 2 30 32 0.393 552955o02oct003p12.fits 52 2 30 32 0.506 550376o09sep000p05.fits 55 2 28 32 0.415 550362o09sep000p03.fits 51 2 30 32 0.390 550377o09sep000p05.fits 55 2 29 33 0.409 552957o02oct003p12.fits 53 2 31 33 0.440 No real worries. crossref dofit.pl master 2 500 dofit3.log cd astrocheck ln -s ../*fead . ln -s ../master ln -s ../rdall.used \ls *.fead | perl -nle 'system "getout $_"' ls *out >outlist cat *out | sort -n | head > outrel cat *out | sort -nr | head >> outrel cat *out | sort -nr +1 | head >> outrel cat *out | sort -n +1 | head >> outrel ** 10 July 2002 Added a feature to stackall.f Every pixel less than the median (less 2000) is replaced with the saturation value: 65535 This will hopefully avoid the problem of solarized saturated stars being identified as two stars. stackall ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cp mask*box .. cd .. mv *o06may*.fits toss mv 531633o??apr002p11.fits toss imlist.pl sex.pl>sex.log ** 12 July 2002 It looks like the Oct00-3 I zero point is off sub getzp { my $head=shift @_; my $expkey ='EXPTIME'; my $airkey ='AIRMASS'; my $datekey ='DATE-OBS'; my $colour='I'; my %fact= qw (B 0.20 V 0.12 R 0.10 I 0.08); my %zpoint = qw { I1999-05-13 25.802 I1999-11-09 25.894 I2000-04-30 26.3386 I2000-05-01 26.3386 I2000-05-03 26.3386 I2000-09-22 26.2906 I2000-10-03 26.2906 I2000-10-23 26.260 I2000-10-24 26.260 I2000-10-25 26.260 }; open HEAD,$head; while () { chomp; ($key,$value)= split /=/; ($value) = split /\//,$value; $key=~s/\s//g; $value=~s/\s//g; $value=~s/\'//g; if ($key eq $expkey ) {$exptime = $value} if ($key eq $airkey ) {$airmass = $value} if ($key eq $datekey) {$date = $value} } close HEAD; $fact=$fact{$colour}; unless (defined($zpoint{$colour.$date})) { print "ERROR: couldn't find zero-point\n"; print "$head $colour $date\n"; die; } $zpoint =$zpoint{$colour.$date} + 2.5*&log10($exptime) - $airmass*$fact; } Back to astrometry. sex.pl >sex.log (did this before) dofit.pl USNO 1 100 dofit1.log checkdo.pl dofit1.log dofit.pl USNO 2 100 dofit2.log checkdo.pl dofit2.log ls *.fits | perl -nle 'system "uphead $_ fead"' dofit.pl USNO 2 100 dofit2.log /raid/gwyn/F22/F22I #52 > checkdo.pl dofit2.log worst centering errors: 551373o06oct001p06.fits 74 2 21 57 1.322 550366o10sep000p03.fits 67 2 50 64 0.730 550376o03sep000p05.fits 73 2 51 64 0.727 508267o09nov990p04.fits 73 2 52 65 0.723 550379o04sep000p05.fits 59 2 46 55 0.709 508268o00nov990p04.fits 72 2 45 54 0.709 fewest initial matches: 550366o05sep000p03.fits 49 2 42 48 0.574 550377o05sep000p05.fits 50 2 38 45 0.508 550376o05sep000p05.fits 50 2 38 42 0.498 550362o05sep000p03.fits 51 2 43 48 0.502 550375o05sep000p05.fits 51 2 38 45 0.566 550361o05sep000p03.fits 52 2 41 49 0.590 fewest final matches: 550376o01sep000p05.fits 58 2 32 34 0.568 550377o01sep000p05.fits 60 2 34 36 0.569 550375o01sep000p05.fits 60 2 36 38 0.536 552957o03oct003p12.fits 52 2 36 39 0.541 550376o05sep000p05.fits 50 2 38 42 0.498 550361o01sep000p03.fits 59 2 40 42 0.528 crossref rdall 2 8 .5 master loaded 576169 objects from 859 files sorting them... done cross correlating, iteration: 1 found 449487 new matches tossed 2350 bad objects 93916 reference objects cross correlating, iteration: 2 found 1 new matches tossed 10 bad objects 93906 reference objects cross correlating, iteration: 3 found 0 new matches tossed 0 bad objects 93906 reference objects 93906 dofit.pl master 2 500 dofit3.log checkdo.pl dofit3.log worst centering errors: 531632o02apr002p11.fits 433 2 321 414 0.163 531631o02apr002p11.fits 430 2 324 413 0.161 531630o02apr002p11.fits 381 2 298 367 0.159 531290o00apr002p06.fits 373 2 270 337 0.155 531288o00apr002p06.fits 364 2 257 323 0.151 531182o00apr002p06.fits 377 2 257 333 0.147 fewest initial matches: 531632o05apr002p11.fits 281 2 204 268 0.081 531632o04apr002p11.fits 282 2 201 265 0.048 531631o04apr002p11.fits 284 2 204 269 0.051 531631o05apr002p11.fits 285 2 207 270 0.080 531630o05apr002p11.fits 285 2 208 271 0.084 531630o04apr002p11.fits 288 2 210 275 0.055 fewest final matches: 531632o04apr002p11.fits 282 2 201 265 0.048 531289o10apr002p06.fits 296 2 230 267 0.068 531632o05apr002p11.fits 281 2 204 268 0.081 531631o04apr002p11.fits 284 2 204 269 0.051 531631o05apr002p11.fits 285 2 207 270 0.080 531630o05apr002p11.fits 285 2 208 271 0.084 crossref rdall 2 8 .5 master2 loaded 576169 objects from 859 files sorting them... done cross correlating, iteration: 1 found 449670 new matches tossed 2300 bad objects 93932 reference objects cross correlating, iteration: 2 found 2 new matches tossed 4 bad objects 93931 reference objects cross correlating, iteration: 3 found 0 new matches tossed 0 bad objects 93931 reference objects 93931 dofit.pl master2 2 500 dofit4.log checkdo.pl dofit4.log worst centering errors: 531631o02apr002p11.fits 430 2 323 414 0.122 531630o02apr002p11.fits 381 2 300 367 0.119 531290o07apr002p06.fits 454 2 300 409 0.114 531288o00apr002p06.fits 368 2 258 324 0.113 531632o02apr002p11.fits 433 2 320 415 0.113 531182o05apr002p06.fits 399 2 275 348 0.112 fewest initial matches: 531632o05apr002p11.fits 281 2 204 267 0.070 531631o05apr002p11.fits 284 2 206 270 0.066 531632o04apr002p11.fits 284 2 201 267 0.047 531630o05apr002p11.fits 285 2 209 271 0.070 531631o04apr002p11.fits 286 2 204 269 0.047 531630o04apr002p11.fits 290 2 210 275 0.050 fewest final matches: 531632o05apr002p11.fits 281 2 204 267 0.070 531632o04apr002p11.fits 284 2 201 267 0.047 531289o10apr002p06.fits 297 2 231 268 0.060 531631o04apr002p11.fits 286 2 204 269 0.047 531631o05apr002p11.fits 284 2 206 270 0.066 531630o05apr002p11.fits 285 2 209 271 0.070 tar -cvf phot.tar *rd *fead (at least for the nonce) crossref rdall 2 8 .5 master3 loaded 576169 objects from 859 files sorting them... done cross correlating, iteration: 1 found 449681 new matches tossed 2294 bad objects 93932 reference objects cross correlating, iteration: 2 found 4 new matches tossed 2 bad objects 93930 reference objects cross correlating, iteration: 3 found 0 new matches tossed 0 bad objects 93930 reference objects 93930 dofit.pl master2 2 500 dofit5.log worst centering errors: 531631o02apr002p11.fits 430 2 323 414 0.122 531630o02apr002p11.fits 381 2 300 367 0.119 531290o07apr002p06.fits 454 2 300 409 0.114 531288o00apr002p06.fits 368 2 258 324 0.113 531632o02apr002p11.fits 433 2 320 415 0.113 531182o05apr002p06.fits 399 2 275 348 0.112 fewest initial matches: 531632o05apr002p11.fits 281 2 204 267 0.070 531631o05apr002p11.fits 284 2 206 270 0.066 531632o04apr002p11.fits 284 2 201 267 0.047 531630o05apr002p11.fits 285 2 209 271 0.070 531631o04apr002p11.fits 286 2 204 269 0.047 531630o04apr002p11.fits 290 2 210 275 0.050 fewest final matches: 531632o05apr002p11.fits 281 2 204 267 0.070 531632o04apr002p11.fits 284 2 201 267 0.047 531289o10apr002p06.fits 297 2 231 268 0.060 531631o04apr002p11.fits 286 2 204 269 0.047 531631o05apr002p11.fits 284 2 206 270 0.066 531630o05apr002p11.fits 285 2 209 271 0.070 ** 13 July 2002 The zero-point for Oct 00-3 is 26.260+/-0.006 bring *.fead and *.rd over to lasa28 /home/gwyn/red/F22/F22I/phot Run erd.pl erdall Exposure.erd zero-point UTC-date Hawaii calculated from (DATE-OBS) Standard the .fead file Time 486370may990p03.erd 32.8002912379574 1999-05-13 03:42 486371may990p03.erd 32.8133312379574 1999-05-13 03:56 486372may990p03.erd 32.8225312379574 1999-05-13 04:09 486373may990p03.erd 32.8272512379574 1999-05-13 04:23 486374may990p03.erd 32.8334112379574 1999-05-13 04:36 508265nov990p04.erd 32.9521312379574 1999-11-09 19:10 etc. run ~/pg/magmatch.f (which now uses the median, instead of the average) run ~/pg/connph.pl which calculates zero-points for every image based on every photometric image. It outputs the average. The graphical results are in: photcalF22I.gif photcalF22I.ps which were outputed by sm.5 Lower left: layout of the 22hr pointings, outlined in red labelled in black. Note the CCD06 of May 1999 has been excluded. Upper left: layout of the 22hr pointings. "Photometric" exposures are in heavy black and labelled "0". Exposures that "touch" (more than 150 common objects) photometric exposures are in red and labelled "1". Exposures that "touch" level 1 exposures are outlined in blue and labelled "2". Big right figure: For every exposure, I calculated the zero-point you would get if you assumed it was photometric, using the airmass and exposure time from the header and the zero-point for the run from Yannick's web page. (Except for Oct-00-3 where I used my own zero-point) I selected some exposures as really being photometric. These ones are shown in black. For each "photometric" exposure, I derived zero-points for all the other exposures, by comparing magnitudes in the overlap regions. The difference between these magnitudes is plotted. The exposures are labelled on the left and the UTC date and HST time are shown on the right. The colour-coding is the same as the upper left panel. sm.6 produces magmatchdiff.ps which shows the difference between Average (3-sigma clipped) Average (2-sigma clipped) Median Only the worst case scenarios (for Noverlaps>150) are shown Not much of a difference really. I checked a few bad cases. In every case, the median is better. They are often cases where CCD's with internal zero-point offsets are being compared. The zero-point shifts internal to each pointing are shown in interncheck.ps generated with sm.2 sm.1 is for comparing the results of matchex... and has the macro comp in it. There are night to night variations in the zero-points. Look at ~/red/stan/*oct003*d.ps SA92 2000-10-23 26.259 SA92 2000-10-24 26.238 SA95 2000-10-24 26.254 ->26.257 SA101 2000-10-24 26.280 SA101 2000-10-25 26.309 ** 14 July 2002 photometry finalized final astrometry checked using lasa28:/home/gwyn/ag/rdmatch.f which matches all the .rd files together matching everything with r<2" and deltamag=.5 It runs on rdall (image zero-point) mkrdall.pl makes rdall, working from the .fead files. Examine the worst case scenarios: lasa28:/home/gwyn/red/F22/F22I/dofit5check/sm.1 which plots up various parameter from rdall.sum name name #overlap #bad .95% fraction_bad 552856o05oct003p09.rd 486370o02may990p03.rd 303 10 0.232 0.033 486372o10may990p03.rd 486370o10may990p03.rd 251 8 0.225 0.032 486370o10may990p03.rd 486372o10may990p03.rd 250 8 0.236 0.032 551376o06oct001p06.rd 531182o00apr002p06.rd 373 12 0.255 0.032 551377o06oct001p06.rd 531182o00apr002p06.rd 362 17 0.285 0.047 550369o05sep000p04.rd 550373o05sep000p04.rd 384 13 0.228 0.034 550378o05sep000p05.rd 550375o05sep000p05.rd 473 15 0.181 0.032 531182o00apr002p06.rd 551376o06oct001p06.rd 380 12 0.255 0.032 bad here means worse than 0.3" (or 1.5 pixels) and lasa28:/home/gwyn/red/F22/F22I/dofit5check/sm.2 which shows the bad region green <0.3" blue <0.6" red >0.6" lasa28:/home/gwyn/ag/applyzp.pl corrects to AB magnitude, then puts the zeropoint into the headers. copy the headers back to pix1 (after doing some checking on the astrometry) cd /raid/gwyn/F22/F22I/ mkdir swarp mkdir stacks cd swarp mkdir vmem mkdir resamp doswarp.pl ** 15 July 2002 For future reference, swarp resamples all the images it can lay it's hands on, regardless if they are relevant. The resampled images come out 1x1 pixel... The speed is not greatly affected. SWarp fails (on pix1) if the center is specified and there are more than 5 images. ** 16 July 2002 Successfully SWarped *oct003p08 nov990p04 oct003p09 F22t2.. sex -c bigimage.sex F22It2.fits -CATALOG_NAME F22It2.cat -WEIGHT_IMAGE F22It2.weight.fits swarping *sep000* Saturation value is about 1000, not 40000 [1] doswarp.pl F22I.list & /raid/gwyn/F22/F22I/swarp #39 > > *Error*: cannot open for reading /raid/gwyn/F22/F22I/swarp/resamp//551382o09oct001p07.resamp.weight.fits ** 17 July 2002 F22It2.fits 0.830 square degrees F22It3.fits 0.965 square degrees sex -c bigimage.sex F22It3.fits -CATALOG_NAME F22It3.cat -WEIGHT_IMAGE F22It3.weight.fits Redo SExtractor because I used the wrong threshold. More data for P1, P2 taken on 20-21 September 2000 setting gain to 0.04 (that is the average of the FSCALES time the average gain of the original images) 0.033 * 1.2 SATUR_LEVEL 1200. # level (in ADUs) at which arises saturation bigimage.sex looks like this: # Default configuration file for SExtractor V1.2b14 - > 2.0 # EB 23/07/98 # (*) indicates parameters which can be omitted from this config file. #-------------------------------- Catalog ----------------------- CATALOG_NAME F22It3.cat # name of the output catalog CATALOG_TYPE ASCII_HEAD # "NONE","ASCII_HEAD","ASCII","FITS_1.0" # or "FITS_LDAC" PARAMETERS_NAME header.param # name of the file # containing catalog contents #------------------------------- Extraction ------------------------- DETECT_TYPE CCD # "CCD" or "PHOTO" (*) DETECT_MINAREA 3 # minimum number of pixels above threshold DETECT_THRESH 1 # or , in mag.arcsec-2 ANALYSIS_THRESH 1 # or , in mag.arcsec-2 THRESH_TYPE RELATIVE FILTER Y FILTER_NAME gauss_3.0_7x7.conv # name of the file containing the filter DEBLEND_NTHRESH 32 # Number of deblending sub-thresholds DEBLEND_MINCONT 0.002 # Minimum contrast parameter for deblending CLEAN Y CLEAN_PARAM 1.0 # Cleaning efficiency MASK_TYPE CORRECT # type of detection MASKing: can be one of # "NONE", "BLANK" or "CORRECT" #------------------------------ Photometry ----------------------------------- PHOT_APERTURES 15.,25. # MAG_APER aperture diameter(s) in pixels PHOT_FLUXFRAC 0.2, 0.5, 0.8 # Fraction of FLUX_AUTO defining FLUX_RADIUS PHOT_AUTOPARAMS 2.5, 3.5 # MAG_AUTO parameters: , SATUR_LEVEL 1000. # level (in ADUs) at which arises saturation PHOT_AUTOAPERS 15.0, 15.0 # MAG_AUTO minimum apertures: estimation,photometry MAG_ZEROPOINT 30.00 # magnitude zero-point MAG_GAMMA 4.0 # gamma of emulsion (for photographic scans) GAIN 0.04 # detector gain in e-/ADU. PIXEL_SCALE 0. # size of pixel in arcsec (0=use FITS WCS info). #------------------------- Star/Galaxy Separation ------------------- SEEING_FWHM 0.9 # stellar FWHM in arcsec STARNNW_NAME default.nnw # Neural-Network_Weight table filename #------------------------------ Background --------------------------- BACK_SIZE 512 # Background mesh: or , BACK_FILTERSIZE 9 # Background filter: or , BACKPHOTO_TYPE LOCAL # can be "GLOBAL" or "LOCAL" (*) BACKPHOTO_THICK 30 # thickness of the background LOCAL annulus (*) BACK_TYPE AUTO BACK_VALUE 0 #------------------------------ Check Image ----------------------- CHECKIMAGE_TYPE NONE # can be one of "NONE", "IDENTICAL", # "BACKGROUND", "-BACKGROUND","BACKGROUND_RMS", # "MINIBACKGROUND", "MINIBACK_RMS", # "FILTERED", "OBJECTS", "-OBJECTS", # "SEGMENTATION", or "APERTURES", CHECKIMAGE_NAME check.fits # Filename for the check-image (*) #--------------------- Memory (change with caution!) -------------------- MEMORY_OBJSTACK 2000 # number of objects in stack MEMORY_PIXSTACK 5000000 # number of pixels in stack MEMORY_BUFSIZE 512 # number of lines in buffer #----------------------------- Miscellaneous -------------------------- VERBOSE_TYPE NORMAL # can be "QUIET", "NORMAL" or "FULL" (*) #------------------------------- New Stuff ---------------------------- WEIGHT_TYPE MAP_WEIGHT #WEIGHT_TYPE BACKGROUND While header.param looks like this: NUMBER X_IMAGE Y_IMAGE ALPHA_J2000 DELTA_J2000 A_IMAGE ERRA_IMAGE B_IMAGE ERRB_IMAGE THETA_IMAGE ERRTHETA_IMAGE X2_IMAGE ERRX2_IMAGE Y2_IMAGE ERRY2_IMAGE XY_IMAGE ERRXY_IMAGE A_WORLD ERRA_WORLD B_WORLD ERRB_WORLD THETA_J2000 ERRTHETA_J2000 X2_WORLD ERRX2_WORLD Y2_WORLD ERRY2_WORLD XY_WORLD ERRXY_WORLD ISOAREA_IMAGE THRESHOLD MU_THRESHOLD BACKGROUND FLUX_MAX MU_MAX KRON_RADIUS FLUX_RADIUS(3) MAG_ISO MAGERR_ISO FLUX_ISO FLUXERR_ISO MAG_ISOCOR MAGERR_ISOCOR FLUX_ISOCOR FLUXERR_ISOCOR MAG_AUTO MAGERR_AUTO FLUX_AUTO FLUXERR_AUTO MAG_BEST MAGERR_BEST FLUX_BEST FLUXERR_BEST MAG_APER(2) MAGERR_APER(2) FLUX_APER(2) FLUXERR_APER(2) FLAGS CLASS_STAR Exciting stuff, eh? The mapix job finished. FTP'ing to lasa28 cut at: RA=334.531 and Dec= 0.409371 into four images see: /home/gwyn/red/F22/F22Icut4.ps generated with /home/gwyn/red/F22/F22I/swarpcheck/sm.1 This works out to: 20006.784 3381.134 334.531 -.4 20006.713 12119.821 334.531 0.1 20006.584 27848.442 334.531 1 for x and 1988.075 17527.351 333.5 0.409371 19464.935 17526.346 334.5 0.409371 32572.268 17527.022 335.25 0.409371 for y call it (20000,17500) NAXIS1 = 37573 / Number of pixels along this axis NAXIS2 = 34187 / Number of pixels along this axis what we want is imcopy F22I[1:20000,1:17500] imcopy F22I[1:20000,17501:34187] imcopy F22I[20001:37573,1:17500] imcopy F22I[20001:37573,17501:34187] From Seb: voila comment on utilise toFITS2d: /u3/hjmcc/bin/toFITS2d -x 9990 -X 15000 -y 3709 -Y 10000 F10I-center+P11_12_new2.fits >! F10I-center+P11_12_new2_sub.fits ceci va produire une sous-image de 15000*10000 dont le coin bas a gauche a pour coordonnees (9990,3709) dans l'ancienne image ... cheers seb /u3/hjmcc/bin/toFITS2d -x 9990 -X 15000 -y 3709 -Y 10000 F10I-center+P11_12_new2.fits >! F10I-center+P11_12_new2_sub.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.fits >! F22Isw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.fits >! F22Ise.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.fits >! F22Inw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.fits >! F22Ine.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.weight.fits >! F22Isw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.weight.fits >! F22Ise.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.weight.fits >! F22Inw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.weight.fits >! F22Ine.weight.fits The north fiel The north fields don't work. There is a problem somewhere 72% into the image. Actually, some of the image is missing. ** 19 July 2002 Working on the quadrants that work. sex -c bigimage.sex F22Ise.fits -CATALOG_NAME F22Ise.cat -WEIGHT_IMAGE F22Ise.weight.fits sex -c bigimage.sex F22Isw.fits -CATALOG_NAME F22Isw.cat -WEIGHT_IMAGE F22Isw.weight.fits exsex F22Ise.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ise.xyrdm exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdm swarp (F22Itry2) ^[[1A^[[1M> Co-adding line: 29360 / 34187 ^[[1A^[[1M> Co-adding line: 29376 / 34187 ^[[1A^[[1M> Co-adding line: 29392 / 34187 ^[[1A^[[1M> Co-adding line: 29408 / 34187 ^[[1Afloating exception Sigh. The part that was written is the same as before imstat F22Itry1/stacks/F22Isw.fits # IMAGE NPIX MEAN MODE MIDPT MIN MAX F22Itry1/stacks/F22Isw.fits 350000000 3.394 -1.397 -0.9151 -350.8 3163. cl> e imstat F22Itry2/stacks/F22Isw.fits # IMAGE NPIX MEAN MODE MIDPT MIN MAX F22Itry2/stacks/F22Isw.fits 350000000 3.394 -1.397 -0.9151 -350.8 3163. cl> Recompiling with -ieee at Henry's suggestion. I don't think it will work. Regions in DS9: "phyiscal" does the correct thing on binned images ( image;box(2807,3512.25,3,10.5,0) # color=red ^ ^ ^ ^ ^ | | | | |____ God knows | | | |_________ length of the box sides | | |___________ length of the box top and bottom | |________________ center of box |_____________________ center of box ** 20 July 2002 Looking for new images in the central part of the F22... /evsix /evsix1/cfh12k_mai_1999/485930o: /sampl /samplix1/DATA/cfh12k_oct2_2000/552568o: 485930 485930o/ bad 485931 485931o.phu bad 485932 485932o1.pl bad 485933 485933o/ bad 485934 485934o/ bad 486004 486004o.phu bad 486005 486005o.phu bad 486006 486006o.phu bad 486007 486007o.phu bad 486008 486008o.phu bad 486009 486009o.phu bad 486090 486090o.phu bad 486091 486091o.phu bad 486092 486092o.phu bad 486093 486093o.phu bad 486176 486176o.phu check it out no fits files 486177 486177o.phu check it out no fits files 486370 486370o10MDFGOfl.fits used 486372 486372o/ used 486373 486373o/ used 486374 486374o10MDFGOfl.fits used 486663 486663o/ B 486664 486664o/ B 486665 486665o/ B 486830 486830o/ B 486831 486831o/ B 486832 486832o/ B 486833 486833o/ B 508265 508265o00MDFG.fits.gz used 508267 508267o00MDFG.fits.gz used 508268 508268o00MDFG.fits.gz used 508269 508269o00MDFG.fits.gz used 508364 508364o10MDFOfl.fits B 508365 508365o00MDF.fits.gz B 508496 508496o00MDF.fits.gz B 508497 508497o10MDFOfl.fits.gz B 508498 508498o00MDF.fits.gz B 508594 508594o10MDFOfl.fits V 508595 508595o00MDF.fits.gz V 508596 508596o00MDF.fits.gz V 508597 508597o00MDF.fits.gz V 508670 508670o00MDF.fits.gz V 552568 552568o/ unreduced 552569 552569o09_hdr.fits unreduced 552570 552570o00.fits unreduced 552571 552571o.phu unreduced 552572 552572o.phu unreduced 552573 552573o00_hdr.fits unreduced 552574 552574o00.fits unreduced 552576 552576o.phu unreduced 552577 552577o/ unreduced 552578 552578o00_hdr.fits unreduced 552579 552579o00.fits unreduced 552632 552632o/ unreduced 552633 552633o/ unreduced 552634 552634o00_hdr.fits unreduced 552635 552635o00.fits unreduced 552637 552637o.phu unreduced 552638 552638o/ unreduced 552639 552639o00_hdr.fits unreduced 552640 552640o.phu unreduced 552641 552641o.phu unreduced 552642 552642o/ unreduced 552643 552643o00_hdr.fits unreduced unreduced 552644 552644o00.fits unreduced 552960 552960o11_hdr.fits unreduced 552962 552962o/ R 552963 552963o09_hdr.fits R 552964 552964o00.fits R 553047 553047f.phu unreduced 553049 553049o/ R 553050 553050d00_hdr.fits unreduced 553051 553051d00.fits unreduced 553052 553052d.phu unreduced 553053 553053d/ unreduced 553054 553054d00_hdr.fits unreduced 553055 553055d.phu unreduced 587773 587773~2.jpg R 594985 594985~1.jpg no fits files 603105 603105~1.hdr V ** 21 July 2002 The latest stack worked! Maybe there is something to this -ieee stuff after all. I doubt it. /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.fits >! F22Isw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.fits >! F22Ise.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.fits >! F22Inw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.fits >! F22Ine.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.weight.fits >! F22Isw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.weight.fits >! F22Ise.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.weight.fits >! F22Inw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.weight.fits >! F22Ine.weight.fits boxes: F22Iseall.reg all the boxes, plus doubles, plus targets F22Iseb4.reg don't know F22Isebox.reg just the boxes, saturated stars and bad regions F22Isestars.reg early star boxes fred.reg USNO.reg USNO catalog tvm.pl does tvmark like stuff. Running imstat to see if everything is hunky-dory IMAGE NPIX MEAN MODE MIDPT MIN MAX F22Isw.fits 350000000 3.394 -1.397 -0.9151 -350.8 3163. It is. Compare with above. SExtraction is happening in pix1:/raid/gwyn/F22/F22Imapix The catalogs from try1 (which failed 72% through the final write) are in pix1:/raid/gwyn/F22/F22Imapix/try1 The good stuff is in pix1:/raid/gwyn/F22/F22Imapix sex -c bigimage.sex F22Ise.fits -CATALOG_NAME F22Ise.cat -WEIGHT_IMAGE F22Ise.weight.fits sex -c bigimage.sex F22Ine.fits -CATALOG_NAME F22Ine.cat -WEIGHT_IMAGE F22Ine.weight.fits sex -c bigimage.sex F22Inw.fits -CATALOG_NAME F22Inw.cat -WEIGHT_IMAGE F22Inw.weight.fits sex -c bigimage.sex F22Isw.fits -CATALOG_NAME F22Isw.cat -WEIGHT_IMAGE F22Isw.weight.fits exsex F22Ise.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ise.xyrdm exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdm exsex F22Ine.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ine.xyrdm exsex F22Inw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Inw.xyrdm Can't display F22I.fits on mapix ** 3 August 2002 Updated F22I.imlist in /samplix1/gwyn/track to include the P13 images. ** 9 August 2002 Rexamine areas: 22:21:31 -00:32:00 ** 13 August 2002 (working from Victoria) updated F22I.dirlist: first line has P13... /samplix1/foucaud/Oct2000_2/Results/f22/ /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded /samplix3/mellier/F22P3-I-May99/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded /samplix3/mellier/F22P4-I-Nov99/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded /samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I /samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded /samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded /samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded Tried to run linkall.pl didn't work because samplix1 isn't online modified linkall.pl to do just P13 (I have now modified it back) On pix1: moved /raid/gwyn/F22/F22I to /raid/gwyn/F22/F22Iv1 working in F22I mkdir data cd data cp /raid/gwyn/F22/F22Iv1/data/*.fits . ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22IP13 get *.fits bye cp /raid/gwyn/F22/F22Iv1/data/add* . cp add00oct001 add00oct002 cp add01oct001 add01oct002 cp add02oct001 add02oct002 cp add03oct001 add03oct002 cp add04oct001 add04oct002 cp add05oct001 add05oct002 cp add06oct001 add06oct002 cp add07oct001 add07oct002 cp add08oct001 add08oct002 cp add09oct001 add09oct002 cp add10oct001 add10oct002 cp add11oct001 add11oct002 imlist.pl stackall >stackall.log ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cp mask*box .. cd .. mkdir toss mv *o06may*.fits toss mv 531633o??apr002p11.fits toss imlist.pl sex.pl>sex.log (get the USNO catalog...) ** 13 August 2002 Using 26.26 for the oct002 zero-point: my %zpoint = qw { I1999-05-13 25.802 I1999-11-09 25.894 I2000-04-30 26.3386 I2000-05-01 26.3386 I2000-05-03 26.3386 I2000-09-22 26.2906 I2000-10-03 26.2906 I2000-10-20 26.26 I2000-10-21 26.26 I2000-10-23 26.260 I2000-10-24 26.260 I2000-10-25 26.260 }; dofit.pl USNO 1 100 dofit1.log checkdo.pl dofit1.log dofit.pl USNO 2 100 dofit2.log problems with some images: Flatfielding problems: look at: 552632o10oct002p13.fits (for example) ** 22 August 2002 On pix1: cd /raid/gwyn/F22/F22I/data/ imcopy the [1:2048,1:4096] section of all the images to backoct002/ (use a ls *.fits | awk '{print imcopy $1....}' >fred.cl style of thing in Linux: cd backoct002 back.pl tar -cvf mb.tar *mb.fits switch to lasa28: cd /home/gwyn/red/F22/F22I/backcheck/mbdir ncftp -u virmos pix1.iap.fr cd /raid/gwyn/F22/F22I/data/backoct002 get mb.tar quit tar -xvf mb.tar ls *mb.fits | perl -nle 'system "backch $_"' > backch.res results: Yup. They're wrong. examine /home/gwyn/red/F22/F22I/backcheck/mbdir/backch.res (morgan or lasa28) Things like: 552642o01oct002p13mb.fits 23635. 15.26 0.53 13.67 11.56 12.23 bad 552641o01oct002p13mb.fits 23103. 15.46 0.27 13.75 11.70 12.27 bad 552643o01oct002p13mb.fits 24005. 15.60 0.53 13.75 11.54 12.49 bad 552644o01oct002p13mb.fits 24335. 16.62 0.88 14.39 11.64 13.21 bad 552641o06oct002p13mb.fits 24123. 21.48 9.25 1.09 2.86 8.36 bad 552642o06oct002p13mb.fits 24665. 20.81 9.25 1.25 2.70 8.44 bad 552644o06oct002p13mb.fits 25684. 21.51 9.29 1.66 2.15 9.15 bad = no good at all ** 23 August 2002 /home/gwyn/mask (on Morgan) will contain a bunch of programs for dealing with the masks. Or at least one: dedub.pl will take F22Iseall.reg and F22Ise.xyrdm and make fred.reg which has all the boxes of F22Iseall.reg (no doubles) and all the usuable stars of F22Ise.xyrdm (flag<4 and mag<22.5) as red boxes and blue circular points F22Isefin.reg is the definitive version. ** 24 August 2002 F22I??.xyrdm (x,y,RA,Dec,Mag_auto,Magerr_auto,flux_radius) fluxradius in pixels... The general masking procedure - make some boxes with ds9 regions - also make them using ~/mask/mksat.pl - save them to F22Iswall.reg - merge them with the relevant object list fro F22Isw.xyrdm using ~/mask/dedub.pl. Output is fred.reg - plot them again with ds9 - add more boxes to fred.reg using a distinctive colour - add these new boxes (by hand) to F22Iswall.reg ** 25 August 2002 Rexamine areas: 22:21:31 -00:32:00 22:15:37 +00:11:53 (sw) Big circle-ish area... exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_APER MAGERR_APER 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdma exsex F22Isw.cat ALPHA_j2000 DELTA_j2000 MAG_APER MAG_AUTO >F22Isw.xymm There are a number of problems.... Which will be full resolved at a later date. ~/mask/dedub.f takes over from ~/mask/dedub.pl It will also take the "lines" which it blankets with boxes grep cyan fred.reg > fred ; cat F22Ineall.reg fred >fred2 ; \mv fred2 F22Ineall.reg ; dedub ** 26 August 2002 awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Isw.xyrdm >fredsw awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Ise.xyrdm >fredse awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Inw.xyrdm >frednw awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Ine.xyrdm >fredne cat fredsw fredse frednw fredne > F22I.rdm conmask.pl takes the F22I??all.reg files and converts them from DS9 format physical coordinates to pairs of RA,Dec (2000). It uses the F22I??.head and xy2skypv to do the transforming. Using the F22I??b4.fits headers (and multiplying the coordinates by 4) doesn't work as well. Compare using ds9 to do the transformation (x,y)->(RA,Dec) and chmask.pl to make convert from DS9 format to corners appmask takes up to three aruguments: box list (assumes ra1 dec1 ra2 dec2 in decimal degrees) inital object list (ra,dec in columns 3,4) output object list (same format as inital, but only the objects outside the boxes) appmask F22Ise.box F22Ise.xyrdm appmask F22Isw.box F22Isw.xyrdm appmask F22Ine.box F22Ine.xyrdm appmask F22Inw.box F22Inw.xyrdm appmask F22I.box F22I.xyrdm (Should be obvious at this point that the .box files contain the masks...) There are patches of badness, but what can you do? New data from Seb /samplix1/foucaud/Oct2000_2/Results/f22_2 Now linked correctly in: /samplix1/gwyn/F22/F22IP13 on pix1: cd /raid/gwyn/F22/F22I/data/ mkdir backoct0002bis cd backoct0002bis ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22IP13 get * bye ls *.fits | awk '{print "imcopy " $1 "[1:2048,1:4096] " $1}' >fred.cl (in IRAF) cl backch.res ** 27 August 2002 Decision: include 552573o??oct002p13 552574o??oct002p13 552576o??oct002p13 552577o??oct002p13 552578o??oct002p13 552579o??oct002p13 except for 00 Others have large flat field gradients... ncftp -u gwyn samplix cd /samplix1/gwyn/F22/F22IP13 get 552573*.fits get 552574*.fits get 552576*.fits get 552577*.fits get 552578*.fits get 552579*.fits bye imlist.pl stackall >stackall.log ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cp mask*box .. cd .. mkdir toss imlist.pl sex.pl>sex.log mv *o06may*.* toss mv 531633o??apr002p11.* toss mv 55257?o00oct002p13.* toss ** 28 August 2002 dofit.pl USNO 1 100 dofit1.log checkdo.pl dofit1.log dofit.pl USNO 2 100 dofit2.log checkdo.pl dofit2.log Looks good... ls *.fits | perl -nle 'system "uphead $_ fead"' dofit.pl USNO 2 100 dofit2.log checkdo.pl dofit2.log crossref (checking... megamatch rdall get rdall.sum plot number of overlapping galaxies vs. relative centering error do this on Morgan: /home/gwyn/red/F22/F22I/dofit6check sm.1 ) dofit.pl master 2 500 dofit3.log checkdo.pl dofit3.log Woah! checkdo.pl dofit3.log worst centering errors: 552576o10oct002p13.fits 421 2 56 158 1.446 531632o02apr002p11.fits 434 2 322 415 0.166 531631o02apr002p11.fits 431 2 325 414 0.164 531630o02apr002p11.fits 382 2 299 368 0.160 531182o05apr002p06.fits 399 2 275 349 0.144 531183o07apr002p06.fits 456 2 311 413 0.143 fewest initial matches: 531632o04apr002p11.fits 285 2 203 269 0.059 531632o05apr002p11.fits 287 2 208 273 0.090 531631o04apr002p11.fits 290 2 211 275 0.073 531631o05apr002p11.fits 294 2 214 277 0.093 531630o04apr002p11.fits 294 2 216 280 0.072 531630o05apr002p11.fits 297 2 218 281 0.094 fewest final matches: 552576o10oct002p13.fits 421 2 56 158 1.446 531632o04apr002p11.fits 285 2 203 269 0.059 531632o05apr002p11.fits 287 2 208 273 0.090 550361o05sep000p03.fits 316 2 203 274 0.080 531631o04apr002p11.fits 290 2 211 275 0.073 To fix: Examine image: looks OK plot detections: also looks OK dofit 552576o10oct002p13.fits .cex master fead 2 y 500 WCS: CRVALs 334.4591667 0.4119167 CRPIXs -2161.23 4218.18 CD: 5.667E-05 -7.500E-07 1.177E-06 5.676E-05 order: 2 PV1_: -1.15E-03 1.01E+00 -4.45E-03 0.00E+00 -2.18E-02 2.16E-02 -1.26E-03 PV2_: 1.16E-03 1.01E+00 -4.27E-03 0.00E+00 3.58E-02 -2.03E-02 3.23E-03 read 500 lines from 552576o10oct002p13.cex read 96044 lines from master retreived 497 stars from memory In tol: 4.0 sec Out tol: 0.0 sec found 421 matches from 500 to 497 Doing linear part In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 497 retreived 495 stars from memory In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 495 In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 495 In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 495 Doing 2nd order In tol: 4.0 sec Out tol: 8.0 sec found 316 matches from 500 to 495 In tol: 1.0 sec Out tol: 8.0 sec found 313 matches from 500 to 495 In tol: 1.0 sec Out tol: 8.0 sec found 312 matches from 500 to 495 In tol: 1.0 sec Out tol: 8.0 sec found 56 matches from 500 to 495 In tol: 2.0 sec Out tol: 4.0 sec found 158 matches from 500 to 495 552576o10oct002p13.fits 421 2 56 158 1.446 WCS: CRVALs 334.4591667 0.4119167 CRPIXs -2161.06 4215.46 CD: 5.666E-05 -7.827E-07 1.150E-06 5.677E-05 order: 2 PV1_: -2.91E-03 1.01E+00 -3.83E-02 0.00E+00 -2.78E-02 1.65E-02 -1.38E-01 PV2_: -6.78E-04 9.79E-01 -2.75E-03 0.00E+00 -9.89E-02 -2.37E-02 -2.59E-04 Probably too many .cex detections Fix it with: dofit 552576o10oct002p13.fits .cex master fead 2 n 400 552576o10oct002p13.fits 332 2 257 318 0.060 We now resume our regualar programming: crossref dofit.pl master 2 500 dofit4.log checkdo.pl dofit4.log megamatch (check...) crossref dofit.pl master 2 500 dofit5.log checkdo.pl dofit5.log megamatch (checking... megamatch rdall get rdall.sum plot number of overlapping galaxies vs. relative centering error do this on Morgan: /home/gwyn/red/F22/F22I/dofit6check sm.1 ) ** 30 August 2002 Briefly re-ran sex.pl dofit.pl master 2 500 dofit6.log with PHOT_APER 15 Then I switched it back again. Having problems matching results from before. Magnitudes have larger errors it seems. Switching the magnitude threshold lower... DETECT_MINAREA 5 # minimum number of pixels above threshold DETECT_THRESH 5 # or , in mag.arcsec-2 ANALYSIS_THRESH 5 # or , in mag.arcsec-2 becomes DETECT_MINAREA 3 # minimum number of pixels above threshold DETECT_THRESH 3 # or , in mag.arcsec-2 ANALYSIS_THRESH 3 # or , in mag.arcsec-2 ** 3 September 2002 To do: From Henry: - Figure out flags - done - square edges - yes - polygon code - done From Vincent: - delivery mechanism: Put it somewhere and tell him - images: F22I.fits and the subrasters - catalogs: with the # thingies - Henry flags: Two columns: ID, flag - Stephen flags - putting stuff on vcom: I have an account and /vcom/gwyn From Stephen: - mask edges: DONE - document everything: - star/galaxy separation - bright star masking - surface area - new flags - check catalog parameters parameters: Done Bottom: -00:31:45 =~ 1122 -0.529166666 Top: +01:18:00 =~ 33095 1.3 West: 22:13:48 =~ 1115 333.45 East: 22:21:52 =~ 36360 335.46666666 Which becomes: # Region file format: DS9 version 3.0 global color=red font="helvetica 10 normal" select=1 edit=1 move=1 delete=1 include=1 fixed=0 source fk5;box(334.4583333,0.38541666,2.0166666,1.8291666666,-0) cat>corners 333.45 -0.529166666 335.46666666 -0.529166666 335.46666666 1.3 333.45 1.3 /home/gwyn/red/F22/F22I/swarpcheck #8 > sky2xypv F22I.head corners 1112.108 1121.920 333.45 -0.529166666 36361.445 1121.937 335.46666666 -0.529166666 36357.678 33095.594 335.46666666 1.3 1115.896 33095.634 333.45 1.3 seeing.f ** 4 September 2002 Dicking with the boxes, just to be sure... on Morgan /home/gwyn/red/F22/F22I/swarpcheck #17 > conmask.pl se sw ne nw cat F22I??.box > F22I.box Had another look at F22Iseb4.fits New and improved version of F22Iseall.fits /home/gwyn/red/F22/F22I/swarpcheck/sm.o produces F22Iboxes.gif black boxes should be retained. blue boxes are outside the boundaries and should be ditched. Two more boxes to mask the seams... RA : 334:31:40 +/- 5" Dec: 0:24:28 +/- 5" 334.5291666666 334:31:45 22:18:07.00 334.5347222222 333:31:55 22:17:08.33 0:24:23.0 0.4063888888 0:24:33.0 0.4091666666 The boxes are: F22Icr.box 333.45 0.4063888888 335.46666666 0.4091666666 334.5291666666 -0.529166666 334.5347222222 1.3 ** OK, we were extracting the wrong -ing half light radius with exsex: This has been fixed on Morgan exsex F22Ise.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ise.xyrdm exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdm exsex F22Ine.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ine.xyrdm exsex F22Inw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Inw.xyrdm cat F22I??.xyrdm >F22I.xyrdm appmask F22Ise.box F22Ise.xyrdm appmask F22Isw.box F22Isw.xyrdm appmask F22Ine.box F22Ine.xyrdm appmask F22Inw.box F22Inw.xyrdm appmask F22I.box F22I.xyrdm The area of the F22I is: 3.349333 (dx=.01 degrees) 3.248901 (dx=.001 degrees) ** 9 September 2002 Belated documenation of star-galaxy separation /home/gwyn/mask/seeing (on Morgan) /home/gwyn/red/F22/F22I/swarpcheck #20 > seeing read 107540 lines fromF22I.xyrdmcl 333.440521 335.479919 -0.543314993 1.31409574 nailing done the last squares 1.70000005 3.20000005 0 1.70000005 1.88750005 1 1.88750005 2.07500005 2 2.07500005 2.26250005 3 2.26250005 2.45000005 4 2.45000005 2.63750005 5 2.63750005 2.82500005 6 2.82500005 3.01250005 7 3.01250005 3.20000005 8 3.20000005 3.38750005 9 3.38750005 3.57500005 ** 10 September 2002 Using lasa28:/home/gwyn/mask/area2.f We get: 3.3426 sq deg at a resolution of 20000 3.3382 sq deg at a resolution of 10000 3.3293 sq deg at a resolution of 5000 3.3110 sq deg at a resolution of 2500 2500 3.3110 5000 3.3293 10000 3.3382 15000 3.3408 20000 3.3426 Running maskfits.f (vcom:/vcom/gwyn/mask or lasa28:/home/gwyn/mask) F22Ine.fits WCS: CRVALs 334.4612061 0.3846313 CRPIXs -1213.00 -406.00 CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05 reading F22Ine.fits done oneing array done reading F22I.box done read 773 boxes -18887.892-16378.080 16357.678 15595.594 doing right doing top writing F22IneMask.fits done area (sq deg): 0.7365 F22Ise.fits WCS: CRVALs 334.4612061 0.3846313 CRPIXs -1213.00 17093.00 CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05 reading F22Ise.fits done oneing array done reading F22I.box done read 468 boxes -18887.892 1120.920 16357.678 33094.594 doing right doing bottom writing F22IseMask.fits done area (sq deg): 0.7834 F22Isw.fits WCS: CRVALs 334.4612061 0.3846313 CRPIXs 18786.00 17093.00 CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05 reading F22Isw.fits done oneing array done reading F22I.box done read 1205 boxes 1111.108 1120.920 36356.678 33094.594 doing left doing bottom writing F22IswMask.fits done area (sq deg): 0.9182 F22Inw.fits WCS: CRVALs 334.4612061 0.3846313 CRPIXs 18786.00 -406.00 CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05 reading F22Inw.fits done oneing array done reading F22I.box done read 1031 boxes 1111.108-16378.080 36356.678 15595.594 doing left doing top writing F22InwMask.fits done area (sq deg): 0.9071 F22Ine area (sq deg): 0.7365 F22Ise area (sq deg): 0.7834 F22Isw area (sq deg): 0.9182 F22Inw area (sq deg): 0.9071 3.3452 sq deg total (adopted) ** 11 September 2002 Belated description of software... area.f calculates the area by making small steps in RA,Dec and counting the ones that that are outside all the masks. It goes through every point and checks to see if it is in every box. SLOW ares2.f (faster) sets up a big grid, sets them all to 1 goes through every box and sets all points inside the box to 0. maskfits.f reads the 4 quartered .fits images, and makes a mask image which is 0 inside a mask (or outside the borders) and 1 outside a mask and inside the borders. It also calculates the surface area. The princples of star galaxy separation are documented at http://www.astrsp-mrs.fr/gwyn/cat/starsep.html seeing.f:(morgan:/home/gwyn/mask/seeing.f) Given a file with ra,dec,magnitude,half-light-radius Figures out typical stellar half-light-radius by finding the peak of the half-light-radius distribution in a given magnitude range. Figures out the width of the locus by following it down until falls below 20% of the peak of the histogram. Run it with ns=40 (40x40 boxes) to get a high resolution seeing map of the F22I field. Run it with ns=20 to get good star-galaxy separation You need more counts in your histogram to get the drop off than you do to get the peak reads: F22I.xyrdmcl writes: radecseeing (Ra,Dec,mode of half-light radius) seeing.fort (binary file for sm.j) histout (histogram (x,y,seeingblock), for reading by sm.n) maskout (mag, half-light, useable, seeing blck, for reading by sm.n) F22I.zones (ra1,ra2,dec1,dec2,hmax,hon,hoff) objects with positions in the box (ra1:ra2,dec1:dec2) with magntiudes in the range 19erdmall #magmatch erdmall erd.pl magmatch erdall ls *o09*fead | perl -nle 'system "getoutex $_"' ls *exp >expall sort -n *exp | head -1 >fred1 sort -n *exp | tail -1 >fred2 sort -n +1 *exp | head -1 >fred3 sort -n +1 *exp | head -1 >fred4 cat fred1 fred2 fred3 fred4 > exprel connph.pl (works now) sm.5 (works) /raid/gwyn/F22/F22I #21 > dofit.pl master 2 500 dofit7.log update dofit on pix1 so it puts in the AstroGwyn header cards. and on pix1:/raid/gwyn/F22/F22I dofit.pl master 2 500 dofit7.log tar -cvf feadrd.tar *.fead *.rd gzip feadrd.tar on lasa28:/home/gwyn/pg/F22I/ ncftp -u gwyn pix1.iap.fr cd /raid/gwyn/F22/F22I get feadrd.tar.gz bye gzip -dc feadrd.tar.gz | tar -xvf - erd.pl magmatch erdall (mag range 1850 Connectivity level Photometric image 1 2 3 4 5 6 7 8 9 531182apr002p06 29 19 11 3 5 10 531183apr002p06 30 18 11 3 5 10 552850oct003p10 9 8 7 16 16 15 6 552851oct003p10 9 8 7 16 16 15 6 552852oct003p10 9 8 7 16 16 15 6 552853oct003p10 9 8 7 16 16 15 6 552854oct003p10 9 8 7 16 16 15 6 552855oct003p08 22 24 21 10 552856oct003p09 19 27 25 6 552858oct003p11 17 7 16 16 15 6 552859oct003p11 16 8 16 16 15 6 552860oct003p11 17 7 16 16 15 6 18150 Doesn't work 1860 Connectivity level Photometric image 1 2 3 4 5 6 7 8 9 531182apr002p06 26 18 15 3 5 10 531183apr002p06 21 23 15 3 5 10 552850oct003p10 9 8 5 12 14 20 9 552851oct003p10 9 8 5 12 14 20 9 552852oct003p10 9 8 5 12 14 20 9 552853oct003p10 9 8 5 12 14 20 9 552854oct003p10 9 8 5 12 14 20 9 552855oct003p08 21 24 15 7 10 552856oct003p09 16 20 20 21 552858oct003p11 17 5 12 14 20 9 552859oct003p11 16 6 12 14 20 9 552860oct003p11 17 5 12 14 20 9 1850 Connectivity level Photometric image 1 2 3 4 5 6 7 8 9 531182apr002p06 29 19 11 3 5 10 531183apr002p06 30 18 11 3 5 10 552850oct003p10 9 8 7 16 16 15 6 552851oct003p10 9 8 7 16 16 15 6 552852oct003p10 9 8 7 16 16 15 6 552853oct003p10 9 8 7 16 16 15 6 552854oct003p10 9 8 7 16 16 15 6 552855oct003p08 22 24 21 10 552856oct003p09 19 27 25 6 552858oct003p11 17 7 16 16 15 6 552859oct003p11 16 8 16 16 15 6 552860oct003p11 17 7 16 16 15 6 ** 2 October 2002 Forget the above. In fact, there was a bug in magmatch.f (which has been fixed) connph.pl sm.5 (just to check) examine ~/log/F22Ip13zp.gif applyzp.pl ls *.fead | perl -nle '($a=$_)=~s/fead/head/; system "cp $_ $a"' tar -cvf head.tar *head move to pix1 cd /raid/gwyn/F22/F22I tar -xvf head.tar final astrometry checked using lasa28:/home/gwyn/ag/rdmatch.f which matches all the .rd files together matching everything with r<2" and deltamag=.5 It runs on rdall (image zero-point) mkrdall.pl makes rdall, working from the .fead files. ** 3 October 2002 rerun crossref and dofit one more time... on pix1:/raid/gwyn/F22/F22I crossref dofit.pl master 2 500 dofitX.log mkdir rdcheck cd rdcheck/ cp ../*.fead . cp ../*.cex . ls *.cex | perl -nle '($r=$_)=~s/cex/rd/; ($f=$_)=~s/cex/fead/; system "xy2skypv $f $_ $r"' cp ../rdall . rdmatch rdall awk '$3>200 && $6>0.03' rdall.sum /raid/gwyn/F22/F22I #46 > awk '$3>200 && $6>0.03' rdall.sum 552856o05oct003p09.rd 486371o02may990p03.rd 354 11 0.229 0.031 552574o03oct002p13.rd 552573o03oct002p13.rd 413 16 0.224 0.039 552576o03oct002p13.rd 552573o03oct002p13.rd 410 14 0.229 0.034 552577o03oct002p13.rd 552573o03oct002p13.rd 379 18 0.249 0.047 552578o03oct002p13.rd 552573o03oct002p13.rd 391 14 0.228 0.036 552579o03oct002p13.rd 552573o03oct002p13.rd 414 17 0.269 0.041 551384o05oct001p07.rd 551379o05oct001p07.rd 317 10 0.136 0.032 552576o05oct002p13.rd 552573o05oct002p13.rd 247 10 0.226 0.040 552577o05oct002p13.rd 552573o05oct002p13.rd 230 8 0.225 0.035 552578o05oct002p13.rd 552573o05oct002p13.rd 265 12 0.249 0.045 552579o05oct002p13.rd 552573o05oct002p13.rd 257 11 0.238 0.043 552576o05oct002p13.rd 552574o05oct002p13.rd 257 8 0.224 0.031 552573o05oct002p13.rd 552576o05oct002p13.rd 247 10 0.226 0.040 552574o05oct002p13.rd 552576o05oct002p13.rd 255 8 0.223 0.031 552578o05oct002p13.rd 552576o05oct002p13.rd 267 10 0.239 0.037 552578o05oct002p13.rd 552577o05oct002p13.rd 259 8 0.226 0.031 552573o05oct002p13.rd 552578o05oct002p13.rd 265 12 0.249 0.045 552576o05oct002p13.rd 552578o05oct002p13.rd 269 9 0.237 0.033 552577o05oct002p13.rd 552578o05oct002p13.rd 261 9 0.227 0.034 552573o05oct002p13.rd 552579o05oct002p13.rd 254 9 0.229 0.035 486370o02may990p03.rd 552805o05oct003p09.rd 276 13 0.250 0.047 486371o02may990p03.rd 552805o05oct003p09.rd 321 14 0.259 0.044 486373o02may990p03.rd 552805o05oct003p09.rd 323 10 0.221 0.031 486374o02may990p03.rd 552805o05oct003p09.rd 333 12 0.227 0.036 552803o05oct003p09.rd 552805o05oct003p09.rd 406 14 0.160 0.034 552856o05oct003p09.rd 552805o05oct003p09.rd 388 17 0.263 0.044 486370o02may990p03.rd 552856o05oct003p09.rd 314 10 0.232 0.032 552956o05oct003p12.rd 552959o05oct003p12.rd 401 13 0.220 0.032 550368o05sep000p04.rd 550373o05sep000p04.rd 399 14 0.224 0.035 550369o05sep000p04.rd 550373o05sep000p04.rd 388 15 0.236 0.039 552579o06oct002p13.rd 552574o06oct002p13.rd 286 9 0.226 0.031 552574o06oct002p13.rd 552579o06oct002p13.rd 286 9 0.226 0.031 552855o06oct003p08.rd 552798o06oct003p08.rd 477 19 0.273 0.040 552798o06oct003p08.rd 552855o06oct003p08.rd 467 18 0.273 0.039 486371o10may990p03.rd 486370o10may990p03.rd 259 9 0.220 0.035 486372o10may990p03.rd 486370o10may990p03.rd 252 8 0.222 0.032 486370o10may990p03.rd 486371o10may990p03.rd 259 9 0.220 0.035 486370o10may990p03.rd 486372o10may990p03.rd 251 8 0.222 0.032 looks good. let update the zero-points # move to lasa28 cd /raid/gwyn/F22/F22I tar -cvf - *.fead | gzip > fead.tar.gz on lasa28: cd /home/gwyn/pg/F22I ncftp -u virmos pix1.iap.fr cd /raid/gwyn/F22/F22I get fead.tar.gz bye gzip -dc fead.tar.gz | tar -xvf - # apply the zero-points. applyzp.pl ls *.fead | perl -nle '($a=$_)=~s/fead/head/; system "cp $_ $a"' tar -cvf - *.head| gzip > head.tar.gz #move to pix1 ncftp -u virmos pix1.iap.fr cd /raid/gwyn/F22/F22I put head.tar.gz gzip -dc head.tar.gz | tar -xvf - # get ready to swarp telnet samplix.iap.fr cd /samplix1/gwyn/F22 mkdir swarpF22I cd swarpF22I ncftp -u virmos pix1.iap.fr cd /raid/gwyn/F22/F22I get *.head get *.fits bye mkdir masks cd masks ncftp -u virmos pix1.iap.fr cd /raid/gwyn/F22/F22I/data/ get mask*.fits bye cd .. ls *.fits >F22I.swarplist mkdir stacks mkdir resamp mkdir vmem get F22I.swarp it looks like this: (note this may not be the same as was used for the first image! mapix died, taking the original with it) # Default configuration file for SWarp 1.33# EB 30/04/2002 # #----------------------------------- Output ----------------------------------- IMAGEOUT_NAME F22Isw.fits # Output filename WEIGHTOUT_NAME F22Isw.weight.fits # Output weight-map filename HEADER_ONLY N # Only a header as an output file (Y/N)? HEADER_SUFFIX .head # Filename extension for additional headers #------------------------------- Input Weights -------------------------------- WEIGHT_TYPE MAP_WEIGHT # "BACKGROUND","MAP_RMS","MAP_VARIANCE" # or "MAP_WEIGHT" WEIGHT_SUFFIX .weight.fits # Suffix to use for weight-maps WEIGHT_IMAGE # Weightmap filename if suffix not used # (all or for each weight-map) #------------------------------- Co-addition ---------------------------------- COMBINE Y # Combine resampled images (Y/N)? COMBINE_TYPE MEDIAN # "MEDIAN","AVERAGE","MIN","MAX" # or "WEIGHTED" #-------------------------------- Astrometry ---------------------------------- CELESTIAL_TYPE NATIVE # "NATIVE", "EQUATORIAL", "GALACTIC", # or "ECLIPTIC" PROJECTION_TYPE TAN # Any WCS projection code or "NONE" CENTER_TYPE ALL # "MANUAL", "ALL" or "MOST" CENTER 00:00:00.0, +00:00:00.0 # Coordinates of the image center #CENTER_TYPE MANUAL #CENTER 22:20:00.0 +00:00:00.0 PIXELSCALE_TYPE MANUAL # "MANUAL", "FIT", "MIN", "MAX" # or "MEDIAN" PIXEL_SCALE 0.206 # Pixel scale IMAGE_SIZE 0 # Image size (0 = AUTOMATIC) #IMAGE_SIZE 1000,1000 #-------------------------------- Resampling ---------------------------------- RESAMPLE Y # Resample input images (Y/N)? RESAMPLE_DIR /raid/gwyn/F22/F22I/swarp/resamp/ # Directory path for resampled images RESAMPLE_SUFFIX .resamp.fits # filename extension for resampled images RESAMPLING_TYPE LANCZOS3 # "NEAREST","BILINEAR","LANCZOS2", # "LANCZOS3" or "LANCZOS4" (1 per axis) OVERSAMPLING 0 # Oversampling in each dimension # (0 = automatic) INTERPOLATE N # Interpolate bad input pixels (Y/N)? # (all or for each image) FSCALASTRO_TYPE FIXED # "NONE" or "FIXED" FSCALE_KEYWORD FSCALE # FITS keyword for the multiplicative # factor to apply to each input image FSCALE_DEFAULT 1.0 # Default FSCALE value if not in header GAIN_KEYWORD GAIN # FITS keyword for eff. gain (e-/ADU) GAIN_DEFAULT 0.0 # Default gain if no FITS keyword found # 0 = infty (all or for each image) #--------------------------- Background subtraction --------------------------- SUBTRACT_BACK Y # Subtraction sky background (Y/N)? # (all or for each image) BACK_TYPE AUTO # "AUTO" or "MANUAL" # (all or for each image) BACK_DEFAULT 0.0 # Default background value in MANUAL # (all or for each image) BACK_SIZE 128 # Background mesh size (pixels) # (all or for each image) BACK_FILTERSIZE 3 # Background map filter range (meshes) # (all or for each image) #------------------------- Virtual memory management -------------------------- VMEM_DIR /raid/gwyn/F22/F22I/swarp/vmem/ # Directory path for swap files VMEM_MAX 2047 # Maximum amount of virtual memory (MB) MEM_MAX 512 # Maximum amount of usable RAM (MB) #------------------------------ Miscellaneous --------------------------------- DELETE_TMPFILES N # Delete temporary resampled FITS files # (Y/N)? WRITE_FILEINFO N # Write information about each input # file in the output image header? VERBOSE_TYPE FULL # "QUIET","NORMAL" or "FULL" NTHREADS 1 # Number of simultaneous threads for # the SMP version of SWarp The actual command to be run is: doswarp.pl F22I.swarplist F22Ihole.fits which generates a swarp command swarp 486370o00may990p03.fits 486370o01may990p03.fits 486370o02may990p03.fits 486370o03may990p03.fits 486370o04may990p03.fits 486370o05may990p03.fits 486370o07may990p03.fits 486370o08may990p03.fits 486370o09may990p03.fits 486370o10may990p03.fits 486370o11may990p03.fits 486371o00may990p03.fits 486371o01may990p03.fits 486371o02may990p03.fits 486371o03may990p03.fits 486371o04may990p03.fits 486371o05may990p03.fits 486371o07may990p03.fits 486371o08may990p03.fits 486371o09may990p03.fits 486371o10may990p03.fits 486371o11may990p03.fits 486372o00may990p03.fits 486372o01may990p03.fits 486372o02may990p03.fits 486372o03may990p03.fits 486372o04may990p03.fits 486372o05may990p03.fits 486372o07may990p03.fits 486372o08may990p03.fits 486372o09may990p03.fits 486372o10may990p03.fits 486372o11may990p03.fits 486373o00may990p03.fits 486373o01may990p03.fits 486373o02may990p03.fits 486373o03may990p03.fits 486373o04may990p03.fits 486373o05may990p03.fits 486373o07may990p03.fits 486373o08may990p03.fits 486373o09may990p03.fits 486373o10may990p03.fits 486373o11may990p03.fits 486374o00may990p03.fits 486374o01may990p03.fits 486374o02may990p03.fits 486374o03may990p03.fits 486374o04may990p03.fits 486374o05may990p03.fits 486374o07may990p03.fits 486374o08may990p03.fits 486374o09may990p03.fits 486374o10may990p03.fits 486374o11may990p03.fits 508265o00nov990p04.fits 508265o01nov990p04.fits 508265o02nov990p04.fits 508265o03nov990p04.fits 508265o04nov990p04.fits 508265o05nov990p04.fits 508265o06nov990p04.fits 508265o07nov990p04.fits 508265o08nov990p04.fits 508265o09nov990p04.fits 508265o10nov990p04.fits 508265o11nov990p04.fits 508266o00nov990p04.fits 508266o01nov990p04.fits 508266o02nov990p04.fits 508266o03nov990p04.fits 508266o04nov990p04.fits 508266o05nov990p04.fits 508266o06nov990p04.fits 508266o07nov990p04.fits 508266o08nov990p04.fits 508266o09nov990p04.fits 508266o10nov990p04.fits 508266o11nov990p04.fits 508267o00nov990p04.fits 508267o01nov990p04.fits 508267o02nov990p04.fits 508267o03nov990p04.fits 508267o04nov990p04.fits 508267o05nov990p04.fits 508267o06nov990p04.fits 508267o07nov990p04.fits 508267o08nov990p04.fits 508267o09nov990p04.fits 508267o10nov990p04.fits 508267o11nov990p04.fits 508268o00nov990p04.fits 508268o01nov990p04.fits 508268o02nov990p04.fits 508268o03nov990p04.fits 508268o04nov990p04.fits 508268o05nov990p04.fits 508268o06nov990p04.fits 508268o07nov990p04.fits 508268o08nov990p04.fits 508268o09nov990p04.fits 508268o10nov990p04.fits 508268o11nov990p04.fits 508269o00nov990p04.fits 508269o01nov990p04.fits 508269o02nov990p04.fits 508269o03nov990p04.fits 508269o04nov990p04.fits 508269o05nov990p04.fits 508269o06nov990p04.fits 508269o07nov990p04.fits 508269o08nov990p04.fits 508269o09nov990p04.fits 508269o10nov990p04.fits 508269o11nov990p04.fits 531182o00apr002p06.fits 531182o01apr002p06.fits 531182o02apr002p06.fits 531182o03apr002p06.fits 531182o04apr002p06.fits 531182o05apr002p06.fits 531182o06apr002p06.fits 531182o07apr002p06.fits 531182o08apr002p06.fits 531182o09apr002p06.fits 531182o10apr002p06.fits 531182o11apr002p06.fits 531183o00apr002p06.fits 531183o01apr002p06.fits 531183o02apr002p06.fits 531183o03apr002p06.fits 531183o04apr002p06.fits 531183o05apr002p06.fits 531183o06apr002p06.fits 531183o07apr002p06.fits 531183o08apr002p06.fits 531183o09apr002p06.fits 531183o10apr002p06.fits 531183o11apr002p06.fits 531288o00apr002p06.fits 531288o01apr002p06.fits 531288o02apr002p06.fits 531288o03apr002p06.fits 531288o04apr002p06.fits 531288o05apr002p06.fits 531288o06apr002p06.fits 531288o07apr002p06.fits 531288o08apr002p06.fits 531288o09apr002p06.fits 531288o10apr002p06.fits 531288o11apr002p06.fits 531289o00apr002p06.fits 531289o01apr002p06.fits 531289o02apr002p06.fits 531289o03apr002p06.fits 531289o04apr002p06.fits 531289o05apr002p06.fits 531289o06apr002p06.fits 531289o07apr002p06.fits 531289o08apr002p06.fits 531289o09apr002p06.fits 531289o10apr002p06.fits 531289o11apr002p06.fits 531290o00apr002p06.fits 531290o01apr002p06.fits 531290o02apr002p06.fits 531290o03apr002p06.fits 531290o04apr002p06.fits 531290o05apr002p06.fits 531290o06apr002p06.fits 531290o07apr002p06.fits 531290o08apr002p06.fits 531290o09apr002p06.fits 531290o10apr002p06.fits 531290o11apr002p06.fits 531630o00apr002p11.fits 531630o01apr002p11.fits 531630o02apr002p11.fits 531630o03apr002p11.fits 531630o04apr002p11.fits 531630o05apr002p11.fits 531630o06apr002p11.fits 531630o07apr002p11.fits 531630o08apr002p11.fits 531630o09apr002p11.fits 531630o10apr002p11.fits 531630o11apr002p11.fits 531631o00apr002p11.fits 531631o01apr002p11.fits 531631o02apr002p11.fits 531631o03apr002p11.fits 531631o04apr002p11.fits 531631o05apr002p11.fits 531631o06apr002p11.fits 531631o07apr002p11.fits 531631o08apr002p11.fits 531631o09apr002p11.fits 531631o10apr002p11.fits 531631o11apr002p11.fits 531632o00apr002p11.fits 531632o01apr002p11.fits 531632o02apr002p11.fits 531632o03apr002p11.fits 531632o04apr002p11.fits 531632o05apr002p11.fits 531632o06apr002p11.fits 531632o07apr002p11.fits 531632o08apr002p11.fits 531632o09apr002p11.fits 531632o10apr002p11.fits 531632o11apr002p11.fits 550361o00sep000p03.fits 550361o01sep000p03.fits 550361o02sep000p03.fits 550361o03sep000p03.fits 550361o04sep000p03.fits 550361o05sep000p03.fits 550361o06sep000p03.fits 550361o07sep000p03.fits 550361o08sep000p03.fits 550361o09sep000p03.fits 550361o10sep000p03.fits 550361o11sep000p03.fits 550362o00sep000p03.fits 550362o01sep000p03.fits 550362o02sep000p03.fits 550362o03sep000p03.fits 550362o04sep000p03.fits 550362o05sep000p03.fits 550362o06sep000p03.fits 550362o07sep000p03.fits 550362o08sep000p03.fits 550362o09sep000p03.fits 550362o10sep000p03.fits 550362o11sep000p03.fits 550363o00sep000p03.fits 550363o01sep000p03.fits 550363o02sep000p03.fits 550363o03sep000p03.fits 550363o04sep000p03.fits 550363o05sep000p03.fits 550363o06sep000p03.fits 550363o07sep000p03.fits 550363o08sep000p03.fits 550363o09sep000p03.fits 550363o10sep000p03.fits 550363o11sep000p03.fits 550366o00sep000p03.fits 550366o01sep000p03.fits 550366o02sep000p03.fits 550366o03sep000p03.fits 550366o04sep000p03.fits 550366o05sep000p03.fits 550366o06sep000p03.fits 550366o07sep000p03.fits 550366o08sep000p03.fits 550366o09sep000p03.fits 550366o10sep000p03.fits 550366o11sep000p03.fits 550367o00sep000p04.fits 550367o01sep000p04.fits 550367o02sep000p04.fits 550367o03sep000p04.fits 550367o04sep000p04.fits 550367o05sep000p04.fits 550367o06sep000p04.fits 550367o07sep000p04.fits 550367o08sep000p04.fits 550367o09sep000p04.fits 550367o10sep000p04.fits 550367o11sep000p04.fits 550368o00sep000p04.fits 550368o01sep000p04.fits 550368o02sep000p04.fits 550368o03sep000p04.fits 550368o04sep000p04.fits 550368o05sep000p04.fits 550368o06sep000p04.fits 550368o07sep000p04.fits 550368o08sep000p04.fits 550368o09sep000p04.fits 550368o10sep000p04.fits 550368o11sep000p04.fits 550369o00sep000p04.fits 550369o01sep000p04.fits 550369o02sep000p04.fits 550369o03sep000p04.fits 550369o04sep000p04.fits 550369o05sep000p04.fits 550369o06sep000p04.fits 550369o07sep000p04.fits 550369o08sep000p04.fits 550369o09sep000p04.fits 550369o10sep000p04.fits 550369o11sep000p04.fits 550373o00sep000p04.fits 550373o01sep000p04.fits 550373o02sep000p04.fits 550373o03sep000p04.fits 550373o04sep000p04.fits 550373o05sep000p04.fits 550373o06sep000p04.fits 550373o07sep000p04.fits 550373o08sep000p04.fits 550373o09sep000p04.fits 550373o10sep000p04.fits 550373o11sep000p04.fits 550374o00sep000p04.fits 550374o01sep000p04.fits 550374o02sep000p04.fits 550374o03sep000p04.fits 550374o04sep000p04.fits 550374o05sep000p04.fits 550374o06sep000p04.fits 550374o07sep000p04.fits 550374o08sep000p04.fits 550374o09sep000p04.fits 550374o10sep000p04.fits 550374o11sep000p04.fits 550375o00sep000p05.fits 550375o01sep000p05.fits 550375o02sep000p05.fits 550375o03sep000p05.fits 550375o04sep000p05.fits 550375o05sep000p05.fits 550375o06sep000p05.fits 550375o07sep000p05.fits 550375o08sep000p05.fits 550375o09sep000p05.fits 550375o10sep000p05.fits 550375o11sep000p05.fits 550376o00sep000p05.fits 550376o01sep000p05.fits 550376o02sep000p05.fits 550376o03sep000p05.fits 550376o04sep000p05.fits 550376o05sep000p05.fits 550376o06sep000p05.fits 550376o07sep000p05.fits 550376o08sep000p05.fits 550376o09sep000p05.fits 550376o10sep000p05.fits 550376o11sep000p05.fits 550377o00sep000p05.fits 550377o01sep000p05.fits 550377o02sep000p05.fits 550377o03sep000p05.fits 550377o04sep000p05.fits 550377o05sep000p05.fits 550377o06sep000p05.fits 550377o07sep000p05.fits 550377o08sep000p05.fits 550377o09sep000p05.fits 550377o10sep000p05.fits 550377o11sep000p05.fits 550378o00sep000p05.fits 550378o01sep000p05.fits 550378o02sep000p05.fits 550378o03sep000p05.fits 550378o04sep000p05.fits 550378o05sep000p05.fits 550378o06sep000p05.fits 550378o07sep000p05.fits 550378o08sep000p05.fits 550378o09sep000p05.fits 550378o10sep000p05.fits 550378o11sep000p05.fits 550379o00sep000p05.fits 550379o01sep000p05.fits 550379o02sep000p05.fits 550379o03sep000p05.fits 550379o04sep000p05.fits 550379o05sep000p05.fits 550379o06sep000p05.fits 550379o07sep000p05.fits 550379o08sep000p05.fits 550379o09sep000p05.fits 550379o10sep000p05.fits 550379o11sep000p05.fits 551371o00oct001p06.fits 551371o01oct001p06.fits 551371o02oct001p06.fits 551371o03oct001p06.fits 551371o04oct001p06.fits 551371o05oct001p06.fits 551371o06oct001p06.fits 551371o07oct001p06.fits 551371o08oct001p06.fits 551371o09oct001p06.fits 551371o10oct001p06.fits 551371o11oct001p06.fits 551372o00oct001p06.fits 551372o01oct001p06.fits 551372o02oct001p06.fits 551372o03oct001p06.fits 551372o04oct001p06.fits 551372o05oct001p06.fits 551372o06oct001p06.fits 551372o07oct001p06.fits 551372o08oct001p06.fits 551372o09oct001p06.fits 551372o10oct001p06.fits 551372o11oct001p06.fits 551373o00oct001p06.fits 551373o01oct001p06.fits 551373o02oct001p06.fits 551373o03oct001p06.fits 551373o04oct001p06.fits 551373o05oct001p06.fits 551373o06oct001p06.fits 551373o07oct001p06.fits 551373o08oct001p06.fits 551373o09oct001p06.fits 551373o10oct001p06.fits 551373o11oct001p06.fits 551374o00oct001p06.fits 551374o01oct001p06.fits 551374o02oct001p06.fits 551374o03oct001p06.fits 551374o04oct001p06.fits 551374o05oct001p06.fits 551374o06oct001p06.fits 551374o07oct001p06.fits 551374o08oct001p06.fits 551374o09oct001p06.fits 551374o10oct001p06.fits 551374o11oct001p06.fits 551375o00oct001p06.fits 551375o01oct001p06.fits 551375o02oct001p06.fits 551375o03oct001p06.fits 551375o04oct001p06.fits 551375o05oct001p06.fits 551375o06oct001p06.fits 551375o07oct001p06.fits 551375o08oct001p06.fits 551375o09oct001p06.fits 551375o10oct001p06.fits 551375o11oct001p06.fits 551376o00oct001p06.fits 551376o01oct001p06.fits 551376o02oct001p06.fits 551376o03oct001p06.fits 551376o04oct001p06.fits 551376o05oct001p06.fits 551376o06oct001p06.fits 551376o07oct001p06.fits 551376o08oct001p06.fits 551376o09oct001p06.fits 551376o10oct001p06.fits 551376o11oct001p06.fits 551377o00oct001p06.fits 551377o01oct001p06.fits 551377o02oct001p06.fits 551377o03oct001p06.fits 551377o04oct001p06.fits 551377o05oct001p06.fits 551377o06oct001p06.fits 551377o07oct001p06.fits 551377o08oct001p06.fits 551377o09oct001p06.fits 551377o10oct001p06.fits 551377o11oct001p06.fits 551379o00oct001p07.fits 551379o01oct001p07.fits 551379o02oct001p07.fits 551379o03oct001p07.fits 551379o04oct001p07.fits 551379o05oct001p07.fits 551379o06oct001p07.fits 551379o07oct001p07.fits 551379o08oct001p07.fits 551379o09oct001p07.fits 551379o10oct001p07.fits 551379o11oct001p07.fits 551380o00oct001p07.fits 551380o01oct001p07.fits 551380o02oct001p07.fits 551380o03oct001p07.fits 551380o04oct001p07.fits 551380o05oct001p07.fits 551380o06oct001p07.fits 551380o07oct001p07.fits 551380o08oct001p07.fits 551380o09oct001p07.fits 551380o10oct001p07.fits 551380o11oct001p07.fits 551381o00oct001p07.fits 551381o01oct001p07.fits 551381o02oct001p07.fits 551381o03oct001p07.fits 551381o04oct001p07.fits 551381o05oct001p07.fits 551381o06oct001p07.fits 551381o07oct001p07.fits 551381o08oct001p07.fits 551381o09oct001p07.fits 551381o10oct001p07.fits 551381o11oct001p07.fits 551382o00oct001p07.fits 551382o01oct001p07.fits 551382o02oct001p07.fits 551382o03oct001p07.fits 551382o04oct001p07.fits 551382o05oct001p07.fits 551382o06oct001p07.fits 551382o07oct001p07.fits 551382o08oct001p07.fits 551382o09oct001p07.fits 551382o10oct001p07.fits 551382o11oct001p07.fits 551383o00oct001p07.fits 551383o01oct001p07.fits 551383o02oct001p07.fits 551383o03oct001p07.fits 551383o04oct001p07.fits 551383o05oct001p07.fits 551383o06oct001p07.fits 551383o07oct001p07.fits 551383o08oct001p07.fits 551383o09oct001p07.fits 551383o10oct001p07.fits 551383o11oct001p07.fits 551384o00oct001p07.fits 551384o01oct001p07.fits 551384o02oct001p07.fits 551384o03oct001p07.fits 551384o04oct001p07.fits 551384o05oct001p07.fits 551384o06oct001p07.fits 551384o07oct001p07.fits 551384o08oct001p07.fits 551384o09oct001p07.fits 551384o10oct001p07.fits 551384o11oct001p07.fits 552573o01oct002p13.fits 552573o02oct002p13.fits 552573o03oct002p13.fits 552573o04oct002p13.fits 552573o05oct002p13.fits 552573o06oct002p13.fits 552573o07oct002p13.fits 552573o08oct002p13.fits 552573o09oct002p13.fits 552573o10oct002p13.fits 552573o11oct002p13.fits 552574o01oct002p13.fits 552574o02oct002p13.fits 552574o03oct002p13.fits 552574o04oct002p13.fits 552574o05oct002p13.fits 552574o06oct002p13.fits 552574o07oct002p13.fits 552574o08oct002p13.fits 552574o09oct002p13.fits 552574o10oct002p13.fits 552574o11oct002p13.fits 552576o01oct002p13.fits 552576o02oct002p13.fits 552576o03oct002p13.fits 552576o04oct002p13.fits 552576o05oct002p13.fits 552576o06oct002p13.fits 552576o07oct002p13.fits 552576o08oct002p13.fits 552576o09oct002p13.fits 552576o10oct002p13.fits 552576o11oct002p13.fits 552577o01oct002p13.fits 552577o02oct002p13.fits 552577o03oct002p13.fits 552577o04oct002p13.fits 552577o05oct002p13.fits 552577o06oct002p13.fits 552577o07oct002p13.fits 552577o08oct002p13.fits 552577o09oct002p13.fits 552577o10oct002p13.fits 552577o11oct002p13.fits 552578o01oct002p13.fits 552578o02oct002p13.fits 552578o03oct002p13.fits 552578o04oct002p13.fits 552578o05oct002p13.fits 552578o06oct002p13.fits 552578o07oct002p13.fits 552578o08oct002p13.fits 552578o09oct002p13.fits 552578o10oct002p13.fits 552578o11oct002p13.fits 552579o01oct002p13.fits 552579o02oct002p13.fits 552579o03oct002p13.fits 552579o04oct002p13.fits 552579o05oct002p13.fits 552579o06oct002p13.fits 552579o07oct002p13.fits 552579o08oct002p13.fits 552579o09oct002p13.fits 552579o10oct002p13.fits 552579o11oct002p13.fits 552797o00oct003p08.fits 552797o01oct003p08.fits 552797o02oct003p08.fits 552797o03oct003p08.fits 552797o04oct003p08.fits 552797o05oct003p08.fits 552797o06oct003p08.fits 552797o07oct003p08.fits 552797o08oct003p08.fits 552797o09oct003p08.fits 552797o10oct003p08.fits 552797o11oct003p08.fits 552798o00oct003p08.fits 552798o01oct003p08.fits 552798o02oct003p08.fits 552798o03oct003p08.fits 552798o04oct003p08.fits 552798o05oct003p08.fits 552798o06oct003p08.fits 552798o07oct003p08.fits 552798o08oct003p08.fits 552798o09oct003p08.fits 552798o10oct003p08.fits 552798o11oct003p08.fits 552799o00oct003p08.fits 552799o01oct003p08.fits 552799o02oct003p08.fits 552799o03oct003p08.fits 552799o04oct003p08.fits 552799o05oct003p08.fits 552799o06oct003p08.fits 552799o07oct003p08.fits 552799o08oct003p08.fits 552799o09oct003p08.fits 552799o10oct003p08.fits 552799o11oct003p08.fits 552800o00oct003p08.fits 552800o01oct003p08.fits 552800o02oct003p08.fits 552800o03oct003p08.fits 552800o04oct003p08.fits 552800o05oct003p08.fits 552800o06oct003p08.fits 552800o07oct003p08.fits 552800o08oct003p08.fits 552800o09oct003p08.fits 552800o10oct003p08.fits 552800o11oct003p08.fits 552801o00oct003p08.fits 552801o01oct003p08.fits 552801o02oct003p08.fits 552801o03oct003p08.fits 552801o04oct003p08.fits 552801o05oct003p08.fits 552801o06oct003p08.fits 552801o07oct003p08.fits 552801o08oct003p08.fits 552801o09oct003p08.fits 552801o10oct003p08.fits 552801o11oct003p08.fits 552802o00oct003p09.fits 552802o01oct003p09.fits 552802o02oct003p09.fits 552802o03oct003p09.fits 552802o04oct003p09.fits 552802o05oct003p09.fits 552802o06oct003p09.fits 552802o07oct003p09.fits 552802o08oct003p09.fits 552802o09oct003p09.fits 552802o10oct003p09.fits 552802o11oct003p09.fits 552803o00oct003p09.fits 552803o01oct003p09.fits 552803o02oct003p09.fits 552803o03oct003p09.fits 552803o04oct003p09.fits 552803o05oct003p09.fits 552803o06oct003p09.fits 552803o07oct003p09.fits 552803o08oct003p09.fits 552803o09oct003p09.fits 552803o10oct003p09.fits 552803o11oct003p09.fits 552804o00oct003p09.fits 552804o01oct003p09.fits 552804o02oct003p09.fits 552804o03oct003p09.fits 552804o04oct003p09.fits 552804o05oct003p09.fits 552804o06oct003p09.fits 552804o07oct003p09.fits 552804o08oct003p09.fits 552804o09oct003p09.fits 552804o10oct003p09.fits 552804o11oct003p09.fits 552805o00oct003p09.fits 552805o01oct003p09.fits 552805o02oct003p09.fits 552805o03oct003p09.fits 552805o04oct003p09.fits 552805o05oct003p09.fits 552805o06oct003p09.fits 552805o07oct003p09.fits 552805o08oct003p09.fits 552805o09oct003p09.fits 552805o10oct003p09.fits 552805o11oct003p09.fits 552806o00oct003p09.fits 552806o01oct003p09.fits 552806o02oct003p09.fits 552806o03oct003p09.fits 552806o04oct003p09.fits 552806o05oct003p09.fits 552806o06oct003p09.fits 552806o07oct003p09.fits 552806o08oct003p09.fits 552806o09oct003p09.fits 552806o10oct003p09.fits 552806o11oct003p09.fits 552850o00oct003p10.fits 552850o01oct003p10.fits 552850o02oct003p10.fits 552850o03oct003p10.fits 552850o04oct003p10.fits 552850o05oct003p10.fits 552850o06oct003p10.fits 552850o07oct003p10.fits 552850o08oct003p10.fits 552850o09oct003p10.fits 552850o10oct003p10.fits 552850o11oct003p10.fits 552851o00oct003p10.fits 552851o01oct003p10.fits 552851o02oct003p10.fits 552851o03oct003p10.fits 552851o04oct003p10.fits 552851o05oct003p10.fits 552851o06oct003p10.fits 552851o07oct003p10.fits 552851o08oct003p10.fits 552851o09oct003p10.fits 552851o10oct003p10.fits 552851o11oct003p10.fits 552852o00oct003p10.fits 552852o01oct003p10.fits 552852o02oct003p10.fits 552852o03oct003p10.fits 552852o04oct003p10.fits 552852o05oct003p10.fits 552852o06oct003p10.fits 552852o07oct003p10.fits 552852o08oct003p10.fits 552852o09oct003p10.fits 552852o10oct003p10.fits 552852o11oct003p10.fits 552853o00oct003p10.fits 552853o01oct003p10.fits 552853o02oct003p10.fits 552853o03oct003p10.fits 552853o04oct003p10.fits 552853o05oct003p10.fits 552853o06oct003p10.fits 552853o07oct003p10.fits 552853o08oct003p10.fits 552853o09oct003p10.fits 552853o10oct003p10.fits 552853o11oct003p10.fits 552854o00oct003p10.fits 552854o01oct003p10.fits 552854o02oct003p10.fits 552854o03oct003p10.fits 552854o04oct003p10.fits 552854o05oct003p10.fits 552854o06oct003p10.fits 552854o07oct003p10.fits 552854o08oct003p10.fits 552854o09oct003p10.fits 552854o10oct003p10.fits 552854o11oct003p10.fits 552855o00oct003p08.fits 552855o01oct003p08.fits 552855o02oct003p08.fits 552855o03oct003p08.fits 552855o04oct003p08.fits 552855o05oct003p08.fits 552855o06oct003p08.fits 552855o07oct003p08.fits 552855o08oct003p08.fits 552855o09oct003p08.fits 552855o10oct003p08.fits 552855o11oct003p08.fits 552856o00oct003p09.fits 552856o01oct003p09.fits 552856o02oct003p09.fits 552856o03oct003p09.fits 552856o04oct003p09.fits 552856o05oct003p09.fits 552856o06oct003p09.fits 552856o07oct003p09.fits 552856o08oct003p09.fits 552856o09oct003p09.fits 552856o10oct003p09.fits 552856o11oct003p09.fits 552858o00oct003p11.fits 552858o01oct003p11.fits 552858o02oct003p11.fits 552858o03oct003p11.fits 552858o04oct003p11.fits 552858o05oct003p11.fits 552858o06oct003p11.fits 552858o07oct003p11.fits 552858o08oct003p11.fits 552858o09oct003p11.fits 552858o10oct003p11.fits 552858o11oct003p11.fits 552859o00oct003p11.fits 552859o01oct003p11.fits 552859o02oct003p11.fits 552859o03oct003p11.fits 552859o04oct003p11.fits 552859o05oct003p11.fits 552859o06oct003p11.fits 552859o07oct003p11.fits 552859o08oct003p11.fits 552859o09oct003p11.fits 552859o10oct003p11.fits 552859o11oct003p11.fits 552860o00oct003p11.fits 552860o01oct003p11.fits 552860o02oct003p11.fits 552860o03oct003p11.fits 552860o04oct003p11.fits 552860o05oct003p11.fits 552860o06oct003p11.fits 552860o07oct003p11.fits 552860o08oct003p11.fits 552860o09oct003p11.fits 552860o10oct003p11.fits 552860o11oct003p11.fits 552953o00oct003p11.fits 552953o01oct003p11.fits 552953o02oct003p11.fits 552953o03oct003p11.fits 552953o04oct003p11.fits 552953o05oct003p11.fits 552953o06oct003p11.fits 552953o07oct003p11.fits 552953o08oct003p11.fits 552953o09oct003p11.fits 552953o10oct003p11.fits 552953o11oct003p11.fits 552954o00oct003p11.fits 552954o01oct003p11.fits 552954o02oct003p11.fits 552954o03oct003p11.fits 552954o04oct003p11.fits 552954o05oct003p11.fits 552954o06oct003p11.fits 552954o07oct003p11.fits 552954o08oct003p11.fits 552954o09oct003p11.fits 552954o10oct003p11.fits 552954o11oct003p11.fits 552955o00oct003p12.fits 552955o01oct003p12.fits 552955o02oct003p12.fits 552955o03oct003p12.fits 552955o04oct003p12.fits 552955o05oct003p12.fits 552955o06oct003p12.fits 552955o07oct003p12.fits 552955o08oct003p12.fits 552955o09oct003p12.fits 552955o10oct003p12.fits 552955o11oct003p12.fits 552956o00oct003p12.fits 552956o01oct003p12.fits 552956o02oct003p12.fits 552956o03oct003p12.fits 552956o04oct003p12.fits 552956o05oct003p12.fits 552956o06oct003p12.fits 552956o07oct003p12.fits 552956o08oct003p12.fits 552956o09oct003p12.fits 552956o10oct003p12.fits 552956o11oct003p12.fits 552957o00oct003p12.fits 552957o01oct003p12.fits 552957o02oct003p12.fits 552957o03oct003p12.fits 552957o04oct003p12.fits 552957o05oct003p12.fits 552957o06oct003p12.fits 552957o07oct003p12.fits 552957o08oct003p12.fits 552957o09oct003p12.fits 552957o10oct003p12.fits 552957o11oct003p12.fits 552958o00oct003p12.fits 552958o01oct003p12.fits 552958o02oct003p12.fits 552958o03oct003p12.fits 552958o04oct003p12.fits 552958o05oct003p12.fits 552958o06oct003p12.fits 552958o07oct003p12.fits 552958o08oct003p12.fits 552958o09oct003p12.fits 552958o10oct003p12.fits 552958o11oct003p12.fits 552959o00oct003p12.fits 552959o01oct003p12.fits 552959o02oct003p12.fits 552959o03oct003p12.fits 552959o04oct003p12.fits 552959o05oct003p12.fits 552959o06oct003p12.fits 552959o07oct003p12.fits 552959o08oct003p12.fits 552959o09oct003p12.fits 552959o10oct003p12.fits 552959o11oct003p12.fits -c /samplix1/gwyn/F22/swarpF22I/F22I.swarp -VMEM_DIR /samplix1/gwyn/F22/swarpF22I/vmem/ -RESAMPLE_DIR /samplix1/gwyn/F22/swarpF22I/resamp/ -VERBOSE_TYPE NORMAL -IMAGEOUT_NAME /samplix1/gwyn/F22/swarpF22I/stacks/F22Ihole.fits -WEIGHTOUT_NAME /samplix1/gwyn/F22/swarpF22I/stacks/F22Ihole.weight.fits -WEIGHT_TYPE MAP_WEIGHT Coadding 925 CCDs... here goes! shall I proceed (y/n) ** 7 October 2002 The Swarp failed due to lack of disk space. Start again on mapix cd /mapxi2/gwyn mkdir F22 ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22 get -R swarpF22I bye cd /mapix4/gwyn/ mkdir F22 cd F22 mkdir swarpF22I cd swarpF22I mkdir stacks mkdir vmem mkdir resamp cd ln -s /mapix4/gwyn/F22/swarpF22I/vmem/ . ln -s /mapix4/gwyn/F22/swarpF22I/stacks ln -s /mapix4/gwyn/F22/swarpF22I/resamp doswarp.pl F22I.swarplist F22Ihole.fits >swarp.log "no ident" unweighted no ext. header 37573x34187 32 bits (floats) Center: 22:17:50.69 +00:23:04.7 129'x117' Scale: 0.206 ''/pixel Gain: 0 e-/ADU Flux scaling (astrom/photom): 1 X / 1 X > Framing and scaling images... ** 9 October 2002 -------------- File 552959o11oct003p12.fits: "no ident" WEIGHTED EXT. HEADER 2048x4096 16 bits (integers) Center: 22:21:40.66 -00:25:36.8 6.92'x13.8' Scale: 0.2027 ''/pixel Gain: 1.45 e-/ADU Flux scaling (astrom/photom): 1.032 X / 0.03112 X Background: 23245.43 RMS: 168.4411 Weight scale: 168.4411 > Resampling line: 4076 / 4078 > Not enough memory for multibuf (nlines*multiwidth elements) ! FUCK! Henry claims the virtual memory limits need to be set larger... Gilles reboots, and we submit the job again... ** 10 October 2002 Swarp fails again: "no ident" WEIGHTED EXT. HEADER 2048x4096 16 bits (integers) Center: 22:21:40.66 -00:25:36.8 6.92'x13.8' Scale: 0.2027 ''/pixel Gain: 1.5 e-/ADU Flux scaling (astrom/photom): 1 X / 0.03112 X Background: 0 RMS: 0 Weight scale: 1 > Co-adding line: 944 / 34187 /mapix2/gwyn/F22/swarpF22I #13 > /mapix2/gwyn/F22/swarpF22I #13 > > Co-adding line: 4096 / 34187 > *Error* while reading 486370o00may990p03.fits /u3/hjmcc/bin/toFITS2d -x 19000 -X 1000 -y 19000 -Y 1000 F22Ihole1.fits >! fred.fits ** 14 October 2002 The last swarp didn't work once again.... doswarp.pl F22I.swarplist F22Ihole.fits > & swarp.err ** 15 October 2002 resampled images and stacks should be on different disks. ** 17 October 2002 /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Ihole.fits >! F22Iholesw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Ihole.fits >! F22Iholese.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Ihole.fits >! F22Iholenw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Ihole.fits >! F22Iholene.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Ihole.weight.fits >! F22Iholesw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Ihole.weight.fits >! F22Iholese.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Ihole.weight.fits >! F22Iholenw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Ihole.weight.fits >! F22Iholene.weight.fits sex -c bigimage.sex F22Iholese.fits -CATALOG_NAME F22Iholese.cat -WEIGHT_IMAGE F22Iholese.weight.fits sex -c bigimage.sex F22Iholene.fits -CATALOG_NAME F22Iholene.cat -WEIGHT_IMAGE F22Iholene.weight.fits sex -c bigimage.sex F22Iholenw.fits -CATALOG_NAME F22Iholenw.cat -WEIGHT_IMAGE F22Iholenw.weight.fits sex -c bigimage.sex F22Iholesw.fits -CATALOG_NAME F22Iholesw.cat -WEIGHT_IMAGE F22Iholesw.weight.fits ** 18 October 2002 on lasa28: working in /home/gwyn/red/F22/F22I/hole get the catalogs... exsex F22Iholene.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholene.rdm exsex F22Iholenw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholenw.rdm exsex F22Iholese.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholese.rdm exsex F22Iholesw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholesw.rdm ln -s /home/gwyn/red/F22/F22I/swarpcheck/F22I??.cat exsex F22Ine.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Ine.rdm exsex F22Inw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Inw.rdm exsex F22Ise.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Ise.rdm exsex F22Isw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Isw.rdm 142430 9254635 105928561 F22Ine.cat 181433 11789830 134946785 F22Inw.cat 123094 7997795 91542551 F22Ise.cat 157748 10250305 117325141 F22Isw.cat 154528 10041005 114929450 F22Iholene.cat 181551 11797500 135034562 F22Iholenw.cat 132558 8612955 98583752 F22Iholese.cat 158213 10280530 117671077 F22Iholesw.cat matchex F22Iholene.rdm F22Ine.rdm F22Ine.mat 1 4 matchex F22Iholenw.rdm F22Inw.rdm F22Inw.mat 1 4 matchex F22Iholese.rdm F22Ise.rdm F22Ise.mat 1 4 matchex F22Iholesw.rdm F22Isw.rdm F22Isw.mat 1 4 about .07" and .01 magntitudes... marginal check out F22Idswarp.gif generated with /home/gwyn/red/F22/F22I/hole/sm.2 exsex F22Iholene.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholene.rdm exsex F22Iholenw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholenw.rdm exsex F22Iholese.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholese.rdm exsex F22Iholesw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholesw.rdm exsex F22Ine.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Ine.rdm exsex F22Inw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Inw.rdm exsex F22Ise.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Ise.rdm exsex F22Isw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Isw.rdm 30708 122832 1197612 F22Iholene.rdm 33344 133376 1300416 F22Iholenw.rdm 28899 115596 1127061 F22Iholese.rdm 30364 121456 1184196 F22Iholesw.rdm 28009 112036 1092351 F22Ine.rdm 33044 132176 1288716 F22Inw.rdm 26250 105000 1023750 F22Ise.rdm 29657 118628 1156623 F22Isw.rdm matchex F22Iholene.rdm F22Ine.rdm F22Ine.mat 1 4 matchex F22Iholenw.rdm F22Inw.rdm F22Inw.mat 1 4 matchex F22Iholese.rdm F22Ise.rdm F22Ise.mat 1 4 matchex F22Iholesw.rdm F22Isw.rdm F22Isw.mat 1 4 ave_sep RMS 68% F22Iholene.rdm F22Ine.rdm 23800 0.044 0.075 0.047 F22Iholenw.rdm F22Inw.rdm 28202 0.066 0.100 0.110 F22Iholese.rdm F22Ise.rdm 22681 0.046 0.081 0.096 F22Iholesw.rdm F22Isw.rdm 25630 0.070 0.109 0.110 doswarp.pl F22I.swarplist F22Iholem.fits > & swarpmatch.err using a a configuration file F22I.matchswarp which looks like this: # Default configuration file for SWarp 1.33# EB 30/04/2002 # #----------------------------------- Output ----------------------------------- IMAGEOUT_NAME F22Isw.fits # Output filename WEIGHTOUT_NAME F22Isw.weight.fits # Output weight-map filename HEADER_ONLY N # Only a header as an output file (Y/N)? HEADER_SUFFIX .head # Filename extension for additional headers #------------------------------- Input Weights -------------------------------- WEIGHT_TYPE MAP_WEIGHT # "BACKGROUND","MAP_RMS","MAP_VARIANCE" # or "MAP_WEIGHT" WEIGHT_SUFFIX .weight.fits # Suffix to use for weight-maps WEIGHT_IMAGE # Weightmap filename if suffix not used # (all or for each weight-map) #------------------------------- Co-addition ---------------------------------- COMBINE Y # Combine resampled images (Y/N)? COMBINE_TYPE MEDIAN # "MEDIAN","AVERAGE","MIN","MAX" # or "WEIGHTED" #-------------------------------- Astrometry ---------------------------------- CELESTIAL_TYPE NATIVE # "NATIVE", "EQUATORIAL", "GALACTIC", # or "ECLIPTIC" PROJECTION_TYPE TAN # Any WCS projection code or "NONE" CENTER_TYPE MANUAL # "MANUAL", "ALL" or "MOST" CENTER 22:17:50.6894,0:23:04.67271 # Coordinates of the image center #CENTER_TYPE MANUAL #CENTER 22:20:00.0 +00:00:00.0 PIXELSCALE_TYPE MANUAL # "MANUAL", "FIT", "MIN", "MAX" # or "MEDIAN" PIXEL_SCALE 0.206 # Pixel scale IMAGE_SIZE 37573,34187 # Image size (0 = AUTOMATIC) #IMAGE_SIZE 1000,1000 #-------------------------------- Resampling ---------------------------------- RESAMPLE Y # Resample input images (Y/N)? RESAMPLE_DIR /raid/gwyn/F22/F22I/swarp/resamp/ # Directory path for resampled images RESAMPLE_SUFFIX .resamp.fits # filename extension for resampled images RESAMPLING_TYPE LANCZOS3 # "NEAREST","BILINEAR","LANCZOS2", # "LANCZOS3" or "LANCZOS4" (1 per axis) OVERSAMPLING 0 # Oversampling in each dimension # (0 = automatic) INTERPOLATE N # Interpolate bad input pixels (Y/N)? # (all or for each image) FSCALASTRO_TYPE FIXED # "NONE" or "FIXED" FSCALE_KEYWORD FSCALE # FITS keyword for the multiplicative # factor to apply to each input image FSCALE_DEFAULT 1.0 # Default FSCALE value if not in header GAIN_KEYWORD GAIN # FITS keyword for eff. gain (e-/ADU) GAIN_DEFAULT 0.0 # Default gain if no FITS keyword found # 0 = infty (all or for each image) #--------------------------- Background subtraction --------------------------- SUBTRACT_BACK Y # Subtraction sky background (Y/N)? # (all or for each image) BACK_TYPE AUTO # "AUTO" or "MANUAL" # (all or for each image) BACK_DEFAULT 0.0 # Default background value in MANUAL # (all or for each image) BACK_SIZE 128 # Background mesh size (pixels) # (all or for each image) BACK_FILTERSIZE 3 # Background map filter range (meshes) # (all or for each image) #------------------------- Virtual memory management -------------------------- VMEM_DIR /raid/gwyn/F22/F22I/swarp/vmem/ # Directory path for swap files VMEM_MAX 2047 # Maximum amount of virtual memory (MB) MEM_MAX 512 # Maximum amount of usable RAM (MB) #------------------------------ Miscellaneous --------------------------------- DELETE_TMPFILES N # Delete temporary resampled FITS files # (Y/N)? WRITE_FILEINFO N # Write information about each input # file in the output image header? VERBOSE_TYPE FULL # "QUIET","NORMAL" or "FULL" NTHREADS 1 # Number of simultaneous threads for # the SMP version of SWarp ** 21 October 2002 The result of the last swarp was: -rw-r--r-- 1 gwyn system 5138040960 Oct 21 08:10 F22Iholem.weight.fits -rw-r--r-- 1 gwyn system 5138040960 Oct 21 08:10 F22Iholem.fits gethead CRVAL1 CRVAL2 F22Iholem.fits 3.344612058e+02 3.846313083e-01 CENTER 22:17:50.6894,0:23:04.67271 # Coordinates of the image center Try yet again with: CENTER 3.344612061e+02,3.846313087e-01 # Coordinates of the image center doswarp.pl F22I.swarplist F22Iholem2.fits >& swarpmatch2.err meanwhile... /mapix2/gwyn/F22/swarpF22I/stacks mkdir try2 mv F22Iholem.weight.fits try2 mv F22Iholem.fits try2 cd try2 /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem.fits >! F22Iholemse.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem.fits >! F22Iholemsw.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem.fits >! F22Iholemnw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem.fits >! F22Iholemne.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem.weight.fits >! F22Iholemse.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem.weight.fits >! F22Iholemsw.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem.weight.fits >! F22Iholemnw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem.weight.fits >! F22Iholemne.weight.fits sex -c bigimage.sex F22Iholemse.fits -CATALOG_NAME F22Iholemse.cat -WEIGHT_IMAGE F22Iholemse.weight.fits sex -c bigimage.sex F22Iholemne.fits -CATALOG_NAME F22Iholemne.cat -WEIGHT_IMAGE F22Iholemne.weight.fits sex -c bigimage.sex F22Iholemnw.fits -CATALOG_NAME F22Iholemnw.cat -WEIGHT_IMAGE F22Iholemnw.weight.fits sex -c bigimage.sex F22Iholemsw.fits -CATALOG_NAME F22Iholemsw.cat -WEIGHT_IMAGE F22Iholemsw.weight.fits exsex F22Iholemne.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemne.rdm exsex F22Iholemnw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemnw.rdm exsex F22Iholemse.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemse.rdm exsex F22Iholemsw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemsw.rdm matchex F22Iholemne.rdm F22Ine.rdm F22Imne.mat 1 4 matchex F22Iholemnw.rdm F22Inw.rdm F22Imnw.mat 1 4 matchex F22Iholemse.rdm F22Ise.rdm F22Imse.mat 1 4 matchex F22Iholemsw.rdm F22Isw.rdm F22Imsw.mat 1 4 ** 22 October 2002 /mapix2/gwyn/F22/swarpF22I/stacks #3 > gethead CRVAL1 CRVAL2 F22Iholem2.fits 3.344612061e+02 3.846313087e-01 /raid/gwyn/F22/F22Imapix #6 > gethead CRVAL1 CRVAL2 F22I.fits 3.344612061e+02 3.846313087e-01 this means that: CENTER 3.344612061e+02,3.846313087e-01 # Coordinates of the image center works.... /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem2.fits >! F22Iholem2se.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem2.fits >! F22Iholem2sw.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem2.fits >! F22Iholem2nw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem2.fits >! F22Iholem2ne.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem2.weight.fits >! F22Iholem2se.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem2.weight.fits >! F22Iholem2sw.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem2.weight.fits >! F22Iholem2nw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem2.weight.fits >! F22Iholem2ne.weight.fits sex -c bigimage.sex F22Iholem2.fits -CATALOG_NAME F22Iholem2.cat -WEIGHT_IMAGE F22Iholem2.weight.fits sex -c bigimage.sex F22Iholem2se.fits -CATALOG_NAME F22Iholem2se.cat -WEIGHT_IMAGE F22Iholem2se.weight.fits sex -c bigimage.sex F22Iholem2ne.fits -CATALOG_NAME F22Iholem2ne.cat -WEIGHT_IMAGE F22Iholem2ne.weight.fits sex -c bigimage.sex F22Iholem2nw.fits -CATALOG_NAME F22Iholem2nw.cat -WEIGHT_IMAGE F22Iholem2nw.weight.fits sex -c bigimage.sex F22Iholem2sw.fits -CATALOG_NAME F22Iholem2sw.cat -WEIGHT_IMAGE F22Iholem2sw.weight.fits sex -c bigimage.sex F22Iholem2ce.fits -CATALOG_NAME F22Iholem2ce.cat -WEIGHT_IMAGE F22Iholem2ce.weight.fits exsex F22Iholem2ne.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2ne.rdm exsex F22Iholem2nw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2nw.rdm exsex F22Iholem2se.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2se.rdm exsex F22Iholem2sw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2sw.rdm area near the hole: 334.506 0.156675 334.794 0.665195 22:18:01.44 22.19:10.56 00:09:24.03 00:39:54.702 This includes all the borders image section --> [19569:24603,13110:21998] or /u3/hjmcc/bin/toFITS2d -x 19569 -X 5034 -y 13110 -Y F22Iholem2.fits >! F22Iholem2ce.fits ** 29 October 2002 sex -c bigimage.sex F22Iholem2ce.fits -CATALOG_NAME F22Iholem2ce.cat -WEIGHT_IMAGE F22Iholem2ce.weight.fits move to Morgan:/home/gwyn/red/F22/F22I/hole exsex F22Iholem2ce.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Iholem2ce.xyrdm ** 30 October 2002 Made boxes for F22Iholem2ce.fits Make a big box: cat > F22Iceout.reg # Region file format: DS9 version 3.0 global color=green font="helvetica 10 normal" select=1 edit=1 move=1 delete=1 include=1 fixed=0 source fk5;box(334.65,0.410935,0.288,0.50852,0) # color=red (Actually, this is off by one pixel) Everything inside this box should be deleted. The objects to be deleted are: Morgan:/home/gwyn/red/F22/F22I/hole/vdb.rdxym Morgan:/home/gwyn/mask/seehole.f Calculates seeing for the hole alone. The fairly ad hoc perl program Morgan:/home/gwyn/mask/patchzone.pl takes F22I.zones.old (the old seeing list) and replaces the hole parts with the new values from the new F22I.zones. patchzone.pl >F22Ihole.zones Hmm. Apparently /home/gwyn/mask/flag.f wasn't documented. lasa28:/home/gwyn/mask/flag.f merges all the catalogs and makes nice flags started on Morgan, edited on lasa28: this version is the definitive with-hole catalog generator edited again on Morgan: becomes: flaghole.f: This version will be the definitive hole-filling versin ** 31 October 2002 Morgan:/home/gwyn/mask/dedub.f now also outputs x1,y1,x2,y2 of the masks. It remains to convert them to RA,Dec awk '{print $1,$2}' F22Ice.boxy >fred xy2skypv F22Iholem2ce.head fred | awk '{print $1,$2}' >fred1 awk '{print $3,$4}' F22Ice.boxy >fred xy2skypv F22Iholem2ce.head fred | awk '{print $1,$2}' >fred2 paste -d ' ' fred1 fred2 >F22Ice.box awk '$1<334.506 || $2<0.156675 || $3>334.794 || $4>0.665195' F22I.box >F22Ib.box mv F22I.box F22Iold.box cat F22Ib.box F22Ice.box >F22I.box ** 1 November 2002 exsex F22Ice.cat ALPHA_J2000 DELTA_J2000 MAG_AUTO MAGERR_AUTO 'MAG_APER(1)' 'MAGERR_APER(2)' FLAGS > F22Ice.rdmmf exsex F22Iholem2.cat ALPHA_J2000 DELTA_J2000 MAG_AUTO MAGERR_AUTO 'MAG_APER(1)' 'MAGERR_APER(2)' FLAGS > F22I.rdmmf ** 4 November 2002 working in virmos:/vcom/gwyn/F22Iastro ln -s ../F22I/F22I.cat ln -s ../F22Ipatch/F22Ice.cat awk '$4<334.506 || $4>334.794 || $5<0.156675 || $5>0.665195' F22I.cat > F22Iex.cat awk '($4<334.506 || $4>334.794 || $5<0.156675 || $5>0.665195) && $56<22.5 ' F22I.cat >F22Iexm.cat awk '$56<22.5' F22Ice.cat > F22Icem.cat grep -v # F22Icem.cat >fred cat F22Iexm.cat fred >F22Im.cat grep -v # F22Ice.cat > fred2 cat F22Iex.cat fred2 > F22Ipa.cat awk '$56<22.5 && $48>$56-1 && $64==0' F22Ipa.cat > F22Ia.cat (that is to say, MAG_APER(1) <22.5 MAG_AUTO > MAG_APER(1) -1 FLAGS=0 ) alternatively, one could say that it has reasonably compact fairly bright uncrowded, unsaturated sources. awk '$0!~/#/ {print $4,$5}' /vcom/gwyn/F22Iastro/F22Ia.cat > F22Ia.rd awk '$0!~/#/ {print $4,$5,$56}' /vcom/gwyn/F22Iastro/F22Ia.cat > F22Ia.rdm