** 12 September 2002 Not a good start. As usual, the terapix disks are in a shambles. copy over morgan:/home/gwyn/red/scripts/findim.pl working in ev6:/u2/gwyn/find ls -R / > files findim.pl F22B.imlist results in: /samplix3/mellier/F22P1-B-Mai99/notcoadded /samplix3/mellier/F22P1-B-Nov99/notcoadded /samplix3/mellier/F22P2-B-Nov99/notcoadded /samplix5/dantel/process-nov2000/images-B/F22P2b-B-Nov00/notcoadded (the last one you can only find on samplix...) F22B.imlist looks like this: P1: 1.508437/1440/0.8/Nov99(01) gone 2.508438/1440/0.8/Nov99(01) gone 3.508439/1440/0.9/Nov99(01) gone 4.508364/1440/1.0/Nov99(01) found 5.508365/1440/1.0/Nov99(01) found 6.486663/1080/1.0/May99(01) found 7.486664/1080/1.0/May99(01) found 8.486665/1080/1.0/May99(01) found 9.486830/1080/0.9/May99(01) found 10.486831/1080/0.9/May99(01) found 11.486832/1080/1.1/May99(01) found 12.486833/1080/1.1/May99(01) found P2: 1.557365/1440/0.75/Nov00(02) found 2.557366/1440/0.75/Nov00(02) found 3.508496/1440/1.0/Nov99(02) found 4.508497/1440/1.4/Nov99(02) found 5.508498/1440/1.2/Nov99(02) found P3: 1.592159/1380/0.97/Jun01(03) 2.592160/1380/0.91/Jun01(03) 3.592161/1380/0.94/Jun01(03) 4.592162/1380/0.93/Jun01(03) 5.592163/1380/0.98/Jun01(03) P4: 1.592464/1380/1.05/Jun01(04) 2.592465/1380/1.03/Jun01(04) 3.592466/1380/1.07/Jun01(04) 4.592467/1380/1.15/Jun01(04) 5.592468/1380/1.26/Jun01(04) A bunch of headers in: /evsix5/mellier/O1AF51-O.leFevre/header P3 was taken on DATE-OBS= '2001-06-19' / UTC Start of observation P4 was taken on DATE-OBS= '2001-06-21' / UTC Start of observation ** 23 September 2002 Found some data on tape. on samplix:/u3/gwyn/find (using files.samplix) findim.pl F22B.imlist gets you: /samplix1/gwyn/raw /samplix3/mellier/F22P1-B-Mai99/notcoadded /samplix3/mellier/F22P1-B-Nov99/notcoadded /samplix3/mellier/F22P2-B-Nov99/notcoadded /samplix5/dantel/process-nov2000/images-B/F22P2b-B-Nov00/notcoadded cat >F22B.dirlist /samplix1/gwyn/raw /samplix3/mellier/F22P1-B-Mai99/notcoadded /samplix3/mellier/F22P1-B-Nov99/notcoadded /samplix3/mellier/F22P2-B-Nov99/notcoadded /samplix5/dantel/process-nov2000/images-B/F22P2b-B-Nov00/notcoadded and giving up on: 1.508437/1440/0.8/Nov99(01) gone 2.508438/1440/0.8/Nov99(01) gone 3.508439/1440/0.9/Nov99(01) gone linkall.pl F22B.imlist F22B.dirlist /samplix1/gwyn/F22/F22B on pix1: cd /raid/gwyn/F22/F22B/data ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22B bye gzip -d *gz ** 7 November 2002 imlist.pl ln -s /raid/gwyn/masks/mask*.fits . maskall problems: FITSIO Error Status = 412: Data type conversion overflow WHY ME OH LORD! ** 12 November 2002 the Nov99 data has a poorly set BZERO keyword gethead bzero *.fits ... 486833o09may990p01.fits 32768.0 486833o10may990p01.fits 32768.0 486833o11may990p01.fits 32768.0 508364o00nov990p01.fits 32767.0 508364o01nov990p01.fits 32767.0 508364o02nov990p01.fits 32767.0 508364o03nov990p01.fits 32767.0 508364o04nov990p01.fits 32767.0 508364o05nov990p01.fits 32767.0 ... Causing IRAF to think they are real (not ushort) and causing maskall.f to print rude messages. This problem doesn't seem to be present in the F22I data. boxall.f does the same thing as maskall.f, but reads box files rather than mask??????.fits ls -1 5*nov99*fits > fredin ls -1 5*nov99*fits > fredout imstat @fredin >fredbefore imarith @fredin + 0 @fredout pixtype=ushort imstat @fredout >fredafter gethead bzero *.fits | awk '$2!=32768' (no answers) diff fredbefore fredafter32c32 < 508496o06nov990p02.fits 8586240 2497. 2524. 2526. 1296. -1. 65534. --- > 508496o06nov990p02.fits 8586240 2497. 2525. 2527. 1296. 0. 65535. /raid/gwyn/F22/F22B/data #74 > head fredbefore rm imlist* imlist.pl maskall > maskall.log check maskall.log and imlistall.log for peculilarities.... Nope. Everything looks OK. Fair range in medians, but then there is a fair range in exposure times... ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd /raid/gwyn/F22/F22B imlist.pl ln -s /raid/gwyn/masks/sexbox/*box . sex.pl >sex.log ln -s /raid/gwyn/USNO/USNOF22 USNO Add some lines to dofit.pl B1999-05-14 25.548 Measured by me ~/log/STAN.log B1999-05-15 25.636 Measured by me ~/log/STAN.log B1999-11-10 26.034 From Yannick's page B1999-11-12 26.034 From Yannick's page B2000-11-27 26.2 Pulled out of my ass (well, a rough estimate from the run) B2001-06-19 25.504 PHOT_C from the header B2001-06-21 26.061 PHOT_C from the header 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"' ln -s /raid/gwyn/USNO/F22Ia.rdm dofit.pl F22Ia.rdm 2 300 dofit3.log ** 13 November 2002 Looking at the data... 486663o00may990p01.fits brickwall 486664o00may990p01.fits brickwall streak 486665o00may990p01.fits brickwall 486830o00may990p01.fits brickwall 486831o00may990p01.fits brickwall 486832o00may990p01.fits brickwall 486833o00may990p01.fits brickwall 486663o01may990p01.fits bright star 486664o01may990p01.fits bright star 486665o01may990p01.fits bright star faint streak 486830o01may990p01.fits bright star 486831o01may990p01.fits bright star 486832o01may990p01.fits bright star 486833o01may990p01.fits bright star 508498o01nov990p02.fits two faint streaks 557365o01nov000p02.fits bright short streak 557366o01nov000p02.fits bright short streak 592163o01jun010p03.fits bright long streak 486663o03may990p01.fits mild brickwall 486664o03may990p01.fits mild brickwall 486665o03may990p01.fits mild brickwall 486830o03may990p01.fits very mild brickwall 486831o03may990p01.fits very mild brickwall 486832o03may990p01.fits very mild brickwall 486833o03may990p01.fits mild brickwall 508496o03nov990p02.fits bright streak 592464o03jun010p04.fits faint streak 486830o04may990p01.fits short bright streak 486831o04may990p01.fits short bright streak 486832o04may990p01.fits short bright streak 508497o04nov990p02.fits two faint streaks 592161o04jun010p03.fits streak 592464o04jun010p04.fits streak 508496o05nov990p02.fits streak 508497o05nov990p02.fits streak 592159o05jun010p03.fits very bright star 592160o05jun010p03.fits very bright star 592161o05jun010p03.fits very bright star 592162o05jun010p03.fits very bright star 592163o05jun010p03.fits very bright star *06may* crap 508497o06nov990p02.fits two streaks 592163o06jun010p03.fits bright streak 508364o07nov990p01.fits short bright streak 508365o07nov990p01.fits short bright streak 508497o07nov990p02.fits two streaks 486663o08may990p01.fits short bright streak 486664o08may990p01.fits short bright streak 486665o08may990p01.fits short bright streak 508497o08nov990p02.fits bright streak 486832o09may990p01.fits scattered light 508496o09nov990p02.fits streak 508497o09nov990p02.fits streaks 592159o09jun010p03.fits bright star 592160o09jun010p03.fits bright star 592161o09jun010p03.fits bright star streak 592162o09jun010p03.fits bright star 592163o09jun010p03.fits bright star 486663o10may990p01.fits scattered light 486664o10may990p01.fits scattered light 486665o10may990p01.fits scattered light 486830o10may990p01.fits scattered light 486831o10may990p01.fits scattered light 486832o10may990p01.fits scattered light 486833o10may990p01.fits scattered light 508496o10nov990p02.fits bright streak 508497o10nov990p02.fits streak 486663o11may990p01.fits scattered light 486664o11may990p01.fits scattered light 486665o11may990p01.fits scattered light 486830o11may990p01.fits scattered light 486831o11may990p01.fits scattered light 486832o11may990p01.fits scattered light 486833o11may990p01.fits scattered light 508364o11nov990p01.fits ripples 508365o11nov990p01.fits ripples 508496o11nov990p02.fits streak 508497o11nov990p02.fits ripples 508498o11nov990p02.fits ripples 557365o11nov000p02.fits ripples 557366o11nov000p02.fits ripples 592159o11jun010p03.fits ripples 592160o11jun010p03.fits ripples 592161o11jun010p03.fits ripples 592162o11jun010p03.fits ripples 592163o11jun010p03.fits ripples 592464o11jun010p04.fits ripples 592465o11jun010p04.fits ripples 592466o11jun010p04.fits ripples 592467o11jun010p04.fits ripples 592468o11jun010p04.fits ripples o03 extra bad columns after may99 FIXED for jun010 o04may extra bad columns FIXED o05 extra bad columns after may99 FIXED o06may complete crap (no shit) o10may extra bad columns FIXED o11may extra bad columns FIXED jun gradients ** 15 November 2002 fixed all the masks for F22B cd /raid/gwyn/F22/F22B/data maskall ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd .. rm *ex sex.pl >sex.log 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 F22Ia.rdm 2 500 dofit3.log ** 11 December 2002 checkdo.pl dofit3.log Hmm. Not many hits on the 11may990 Redo masks cd /raid/gwyn/F22/F22B/data maskall ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd .. rm *ex sex.pl >sex.log 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 F22Ia.rdm 2 500 dofit3.log checkdo.pl dofit3.log mkdir toss *mv 486833* toss *mv 486663* toss mv *may* toss rm imlist* imlist.pl dofit.pl F22Ia.rdm 2 500 dofit3.log checkdo.pl dofit3.log 486832o10may990p01.fits 91 2 75 70 0.079 486664o10may990p01.fits 88 2 74 73 0.106 486665o10may990p01.fits 89 2 73 75 0.130 486832o06may990p01.fits 87 2 77 77 0.134 592159o05jun010p03.fits 98 2 76 77 0.149 486830o10may990p01.fits 106 2 78 78 0.141 ** 12 December Looks good: check out: Morgan:/home/gwyn/red/F22/F22B/astrocheck/sm.1 which produces F22Bastrores.ps going back to replace the tossed images: except mv *o06may* toss Note: I've lowered the thesholds in xymf.sex: DETECT_MINAREA 5 # minimum number of pixels above threshold DETECT_THRESH 3 # or , in mag.arcsec-2 ANALYSIS_THRESH 3 # or , in mag.arcsec-2 I've added a parameter to xymf.param X_IMAGE Y_IMAGE MAG_APER MAGERR_APER FLAGS MAG_AUTO cleansex.f now tosses objects whose MAG_AUTO and MAG_APER don't match: if (am-amauto.gt.0.2) goto 100 rm imlist* imlist.pl sex.pl dofit.pl F22Ia.rdm 2 500 dofit3.log checkdo.pl dofit3.log A quick look at the photometry: examine: F22Bprephot.ps generated with Morgran:/home/gwyn/red/F22/F22B/photcheck/sm.1 Mag cut off is 18 exprel cat *exp | sort -n | head -1 >> exprel cat *exp | sort -n +1 | tail -1 >> exprel cat *exp | sort -n +1 | head -1 >> exprel ls *exp >expall ** 16 December 2002 check over magmatch.f rdm.pl (takes rdall and the .rd files, and makes .rdm and rdmall files) magmatch rdmall checkphot.pl | grep z z jun010p04may990p01 z may990p01nov000p02 z may990p01nov990p02 z nov000p02nov990p02 The bad flatfield problem is still present... magmatch erdall sm.5 (which runs connph.pl) hem, haw hrmm... Take all the jun01 exposures as photometric: F22B.exlist: #486663may990p01.erd 32.7665593854043 1999-05-14 03:46 #486664may990p01.erd 32.8073593854043 1999-05-14 04:06 #486665may990p01.erd 32.8377593854043 1999-05-14 04:28 #486830may990p01.erd 32.8263593854043 1999-05-15 03:30 #486831may990p01.erd 32.8937593854043 1999-05-15 04:02 #486832may990p01.erd 32.9255593854043 1999-05-15 04:22 #486833may990p01.erd 32.9443593854043 1999-05-15 04:41 #508364nov990p01.erd 33.9355062267886 1999-11-10 19:16 #508365nov990p01.erd 33.9371062267886 1999-11-10 19:41 #508496nov990p02.erd 33.9345062267886 1999-11-12 19:04 #508497nov990p02.erd 33.9365062267886 1999-11-12 19:30 #508498nov990p02.erd 33.9357062267886 1999-11-12 19:55 #557365nov000p02.erd 33.8805062267886 2000-11-27 19:16 #557366nov000p02.erd 33.8753062267886 2000-11-27 19:41 592159jun010p03.erd 33.0382977125737 2001-06-19 02:02 592160jun010p03.erd 33.0696977125737 2001-06-19 02:26 592161jun010p03.erd 33.0926977125737 2001-06-19 02:50 592162jun010p03.erd 33.1098977125737 2001-06-19 03:15 592163jun010p03.erd 33.1224977125737 2001-06-19 03:39 592464jun010p04.erd 33.6276977125737 2001-06-21 02:20 592465jun010p04.erd 33.6502977125737 2001-06-21 02:44 592466jun010p04.erd 33.6670977125737 2001-06-21 03:09 592467jun010p04.erd 33.6794977125737 2001-06-21 03:33 592468jun010p04.erd 33.6884977125737 2001-06-21 03:57 And done! F22B.zpoint #### the following exposures were assumed photometric #### # 592159jun010p03 33.6563127125737 # 592160jun010p03 33.6809127125737 # 592161jun010p03 33.6930127125737 # 592162jun010p03 33.7081127125737 # 592163jun010p03 33.7211127125737 # 592464jun010p04 33.7165127125737 # 592465jun010p04 33.7391127125737 # 592466jun010p04 33.7437127125737 # 592467jun010p04 33.7418127125737 # 592468jun010p04 33.7347127125737 486663may990p01 32.792 7.6390e-02 486664may990p01 32.820 7.4486e-02 486665may990p01 32.857 7.2010e-02 486830may990p01 32.774 7.7731e-02 486831may990p01 32.768 7.8126e-02 486832may990p01 32.802 7.5738e-02 486833may990p01 32.844 7.2837e-02 508364nov990p01 33.908 2.7345e-02 508365nov990p01 33.908 2.7342e-02 508496nov990p02 33.930 2.6789e-02 508497nov990p02 33.831 2.9338e-02 508498nov990p02 33.905 2.7426e-02 557365nov000p02 33.809 2.9950e-02 557366nov000p02 33.773 3.0954e-02 592159jun010p03 33.656 3.4473e-02 592160jun010p03 33.681 3.3700e-02 592161jun010p03 33.693 3.3327e-02 592162jun010p03 33.708 3.2867e-02 592163jun010p03 33.721 3.2475e-02 592464jun010p04 33.717 3.2613e-02 592465jun010p04 33.739 3.1941e-02 592466jun010p04 33.744 3.1806e-02 592467jun010p04 33.742 3.1862e-02 592468jun010p04 33.735 3.2071e-02 ** 18 December 2002 ls *fead | perl -nle 's/.fead//; rename "$_.fead","$_.head"' ** 10 January 2003 on mapix... cd /mapix2/gwyn/F22 mkdir swarpF22B cd swarpF22B mkdir /mapix4/gwyn/F22/swarpF22B mkdir /mapix4/gwyn/F22/swarpF22B/vmem mkdir /mapix4/gwyn/F22/swarpF22B/resamp ln -s /mapix4/gwyn/F22/swarpF22B/vmem ln -s /mapix4/gwyn/F22/swarpF22B/resamp mkdir stacks mkdir masks ncftp -u virmos pix1 cd /raid/gwyn/F22/F22B get *.fits get head.tar ## which comes from lasa28:/home/gwyn/red/F22/F22B/photcheck bye rm *o00may990* ## cuz of brickwall... cd masks ncftp -u virmos pix1 cd /raid/gwyn/masks get mask*jun010*.fits get mask*may990*.fits get mask*nov990*.fits get mask*nov000*.fits bye cd .. ls *.fits >F22B.swarplist cp /mapix2/gwyn/F22/swarpF22I/F22I.matchswarp F22B.matchswarp imhead /mapix2/gwyn/F22/swarpF22I/stacks/F22Iholem2.fits > F22B.head doswarp.pl F22B.swarplist F22B.fits >&F22B.swarplog ** 12 January 2003 /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22B.fits >! F22Bsw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22B.fits >! F22Bse.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22B.fits >! F22Bnw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22B.fits >! F22Bne.fits /u3/hjmcc/bin/toFITS2d -x 19569 -X 5034 -y 13110 -Y 8888 F22B.fits >! F22Bce.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22B.weight.fits >! F22Bsw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22B.weight.fits >! F22Bse.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22B.weight.fits >! F22Bnw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22B.weight.fits >! F22Bne.weight.fits /u3/hjmcc/bin/toFITS2d -x 19569 -X 5034 -y 13110 -Y 8888 F22B.weight.fits >! F22Bce.weight.fits ** 13 January 2003 sex -c bigimage.sex F22Bse.fits -CATALOG_NAME F22Bse.cat -WEIGHT_IMAGE F22Bse.weight.fits sex -c bigimage.sex F22Bne.fits -CATALOG_NAME F22Bne.cat -WEIGHT_IMAGE F22Bne.weight.fits sex -c bigimage.sex F22Bnw.fits -CATALOG_NAME F22Bnw.cat -WEIGHT_IMAGE F22Bnw.weight.fits sex -c bigimage.sex F22Bsw.fits -CATALOG_NAME F22Bsw.cat -WEIGHT_IMAGE F22Bsw.weight.fits Mow the field into chunks, to compare the F22B and F22I astrometry... F22Bsw[12000:20000,9000:17500] F22Bnw[12000:20000, 1:9000] F22Bne[1:17573,1:9000] F22Bse[1:17573,9000:17500] lasa28:/home/gwyn/ag/mow.pl cuts all quartered images into 1200x850 chunks. The size is determined by the size of lasa28's screen. This way we can blink all the images section by section... cd sub ls *B*.fits | perl -nle '($a=$_)=~s/B/I/; system "ds9 -geometry 1300x1200+0+0 -nopanner -nomagnifier -noinfo -histequ $a -histequ -zscale $_ -single"' exsex F22Bse.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bse.rdm exsex F22Bsw.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bsw.rdm exsex F22Bne.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bne.rdm exsex F22Bnw.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bnw.rdm /home/gwyn/red/F22/F22B/swarpcheck #22 > match F22Bnw.rdm F22Ia.rdm F22Bnw.rdm F22Ia.rdm 2097 0.145 0.225 0.143 1426 /home/gwyn/red/F22/F22B/swarpcheck #23 > match F22Bne.rdm F22Ia.rdm F22Bne.rdm F22Ia.rdm 5844 0.126 0.206 0.124 3974 /home/gwyn/red/F22/F22B/swarpcheck #24 > match F22Bse.rdm F22Ia.rdm F22Bse.rdm F22Ia.rdm 5852 0.136 0.202 0.134 3979 /home/gwyn/red/F22/F22B/swarpcheck #25 > match F22Bsw.rdm F22Ia.rdm F22Bsw.rdm F22Ia.rdm 2192 0.131 0.189 0.138 1491 matchex F22Bsw.rdm F22Ia.rdm matchsw2 2 4 matchex F22Bse.rdm F22Ia.rdm matchse2 2 4 matchex F22Bnw.rdm F22Ia.rdm matchnw2 2 4 matchex F22Bne.rdm F22Ia.rdm matchne2 2 4 F22Bsw.rdm F22Ia.rdm 2192 0.131 0.189 0.138 F22Bse.rdm F22Ia.rdm 5852 0.136 0.202 0.134 F22Bnw.rdm F22Ia.rdm 2097 0.145 0.225 0.143 F22Bne.rdm F22Ia.rdm 5844 0.126 0.206 0.124 matchex F22Bsw.rdm F22Ia.rdm matchsw1 1 2 matchex F22Bse.rdm F22Ia.rdm matchse1 1 2 matchex F22Bnw.rdm F22Ia.rdm matchnw1 1 2 matchex F22Bne.rdm F22Ia.rdm matchne1 1 2 F22Bsw.rdm F22Ia.rdm 2323 0.128 0.166 0.138 F22Bse.rdm F22Ia.rdm 6593 0.132 0.180 0.135 F22Bnw.rdm F22Ia.rdm 2737 0.135 0.178 0.143 F22Bne.rdm F22Ia.rdm 7016 0.117 0.164 0.124 cat match??2 >match2 cat match??1 >match1 examined patches se021 vertical shifts se025 no problems se107 horizontal shifts ne148 vertical shifts ne104 barely perceptible diagonal shifts ** 14 January 2003 # clean the files... ls *rd | perl -nle '($a=$_)=~s/rd/rdm/; system "awk \x27 \$6 < .5 && \$8 > \$5-1 && \$7==0 \{print \$1,\$2,\$5\}\x27 $_ > $a"' #match to the standard... ls *rdm | perl -nle '($a=$_)=~s/rdm/rdma/; system "matchex $_ F22Ia.rdm $a"' > matchlog awk '{print $0,FILENAME}' *.rdma > allrdma matchex F22Ia.rdm USNOF22 F22IaUSNO awk '{print $1,$2,$3}' master >master.rdm matchex F22Ia.rdm master.rdm F22Iamaster ** 15 January 2003 dofit 486663o08may990p01.fits .cex F22Ia.rdm fead 2 n 500 0 0 cp 486663o08may990p01.mat 486663o08may990p01.mat500 dofit 486663o08may990p01.fits .cex F22Ia.rdm fead 2 n 100 0 0 cp 486663o08may990p01.mat 486663o08may990p01.mat100 dofit 486663o08may990p01.fits .cex F22Ia.rdm fead 2 n 150 0 0 cp 486663o08may990p01.mat 486663o08may990p01.mat150 dofit 486663o08may990p01.fits .cex F22Ia.rdm fead 2 n 200 0 0 cp 486663o08may990p01.mat 486663o08may990p01.mat200 dofit 486663o08may990p01.fits .cex F22Ia.rdm fead 3 n 500 0 0 cp 486663o08may990p01.mat 486663o08may990p01.mat5003 ** 16 January 2003 Well, astrometry problems exist. Solution for the F22 data: mask them and the end Now for extracting catalogs... on mapix: cd /mapix2/gwyn/F22/swarpF22B/stacks # get the old quadrent images... ncftp -u virmos pix1 cd /raid/gwyn/F22/F22Imapix get F22I??.fits get F22I??.weight.fits bye # get the centre, hole plugging image... ln -s ../../swarpF22I/stacks/F22Iholem2ce*fits . cp ~/config/header.param . cp ~/config/bigimage.sex . cp ~/config/default.nnw . cp ~/config/gauss_3.0_7x7.conv . sex -c bigimage.sex F22Ine.fits,F22Bne.fits -WEIGHT_IMAGE F22Ine.weight.fits,F22Bne.weight.fits -CATALOG_NAME F22Bne.cat sex -c bigimage.sex F22Ise.fits,F22Bse.fits -WEIGHT_IMAGE F22Ise.weight.fits,F22Bse.weight.fits -CATALOG_NAME F22Bse.cat sex -c bigimage.sex F22Isw.fits,F22Bsw.fits -WEIGHT_IMAGE F22Isw.weight.fits,F22Bsw.weight.fits -CATALOG_NAME F22Bsw.cat sex -c bigimage.sex F22Inw.fits,F22Bnw.fits -WEIGHT_IMAGE F22Inw.weight.fits,F22Bnw.weight.fits -CATALOG_NAME F22Bnw.cat sex -c bigimage.sex F22Iholem2ce.fits,F22Bce.fits -WEIGHT_IMAGE F22Iholem2ce.weight.fits,F22Bce.weight.fits -CATALOG_NAME F22Bce.cat /mapix2/gwyn/F22/swarpF22B/stacks #41 > wc F22Ine.cat 142419 9253920 105945832 F22Ine.cat /mapix2/gwyn/F22/swarpF22B/stacks #42 > wc F22Bne.cat 142419 9253920 107818405 F22Bne.cat sex -c bigimage.sex F22Ine.fits,F22Ine.fits -WEIGHT_IMAGE F22Ine.weight.fits,F22Ine.weight.fits -CATALOG_NAME F22Ine.cat sex -c bigimage.sex F22Ine.fits -CATALOG_NAME F22Ine.cat -WEIGHT_IMAGE F22Ine.weight.fits /mapix2/gwyn/F22/swarpF22B/stacks #44 > wc F22Ine.cat 142419 9253920 105920349 F22Ine.cat sex -c bigimage.sex F22Ine.fits,F22Ine.fits -WEIGHT_IMAGE F22Ine.weight.fits,F22Ine.weight.fits -CATALOG_NAME F22Ine.cat sex -c bigimage.sex F22Ise.fits,F22Ise.fits -WEIGHT_IMAGE F22Ise.weight.fits,F22Ise.weight.fits -CATALOG_NAME F22Ise.cat sex -c bigimage.sex F22Isw.fits,F22Isw.fits -WEIGHT_IMAGE F22Isw.weight.fits,F22Isw.weight.fits -CATALOG_NAME F22Isw.cat sex -c bigimage.sex F22Inw.fits,F22Inw.fits -WEIGHT_IMAGE F22Inw.weight.fits,F22Inw.weight.fits -CATALOG_NAME F22Inw.cat sex -c bigimage.sex F22Iholem2ce.fits,F22Iholem2ce.fits -WEIGHT_IMAGE F22Iholem2ce.weight.fits,F22Iholem2ce.weight.fits -CATALOG_NAME F22Ice.cat ** 20 January 2003 on mapix: Extract some matching parameters: exsex F22Ice.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Ice.rdxyim exsex F22Ise.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Ise.rdxyim exsex F22Ine.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Ine.rdxyim exsex F22Isw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Isw.rdxyim exsex F22Inw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Inw.rdxyim exsex F22Bce.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bce.rdxyim exsex F22Bse.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bse.rdxyim exsex F22Bne.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bne.rdxyim exsex F22Bsw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bsw.rdxyim exsex F22Bnw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bnw.rdxyim on lasa28: get the files from mapix: ncftp -u gwyn mapix.iap.fr cd /mapix2/gwyn/F22/swarpF22B/stacks get *rdxyim get the virmos database off the web: awk '$0!~/#/ {print $2,$3,$4,$5,$1,$6}' Jan20-2003_14_24_04.virmos | sort -n +4 >vdb awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Ise.rdxyim > fred \mv fred F22Ise.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Isw.rdxyim > fred \mv fred F22Isw.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Ine.rdxyim > fred \mv fred F22Ine.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Inw.rdxyim > fred \mv fred F22Inw.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Bse.rdxyim > fred \mv fred F22Bse.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Bsw.rdxyim > fred \mv fred F22Bsw.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Bne.rdxyim > fred \mv fred F22Bne.rdxyim awk '$1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195' F22Bnw.rdxyim > fred \mv fred F22Bnw.rdxyim check to see if <= vs < matters. awk '$1==334.506 || $1==334.794 || $2==0.156675 || $2==0.665195' vdb awk '$1==334.506 || $1==334.794 || $2==0.156675 || $2==0.665195' *rdxyim It does not. The big long thin boxes from F22I.boxes: 333.45 0.4063888888 335.46666666 0.4091666666 334.5291666666 -0.529166666 334.5347222222 1.3 awk '$1<334.5291666666 && $2>.4091666666 && ($1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195)' vdb > vdbInw awk '$1>334.5347222222 && $2>.4091666666 && ($1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195)' vdb > vdbIne awk '$1<334.5291666666 && $2<.4063888888 && ($1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195)' vdb > vdbIsw awk '$1>334.5347222222 && $2<.4063888888 && ($1<334.506 || $1> 334.794 || $2<0.156675 || $2>0.665195)' vdb > vdbIse matchall F22Inw.rdxyim vdbInw vdbInwm .01 0 matchall F22Ine.rdxyim vdbIne vdbInem .01 0 matchall F22Isw.rdxyim vdbIsw vdbIswm .01 0 matchall F22Ise.rdxyim vdbIse vdbIsem .01 0 lasa28:/home/gwyn/ag/matchall.f is similar to matchex.f but if it can't find a match, prints the line with "no match" instead skipping. wc F22Inw.rdxyim vdbInw vdbInwm 180591 1083546 10329758 F22Inw.rdxyim 180050 1080300 10057378 vdbInw 178543 2142516 20544350 vdbInwm awk '$0!~no {print $6,$12}' vdb*m >fred ** 21 January 2003 awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIse >vdbIsef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Ise.rdxyim >F22Isef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIne >vdbInef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Ine.rdxyim >F22Inef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIsw >vdbIswf awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Isw.rdxyim >F22Iswf.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbInw >vdbInwf awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Inw.rdxyim >F22Inwf.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIce >vdbIcef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Ice.rdxyim >F22Icef.rdxyim matchall F22Isef.rdxyim vdbIsef vdbIsefm .01 matchall F22Inef.rdxyim vdbInef vdbInefm .01 matchall F22Iswf.rdxyim vdbIswf vdbIswfm .01 matchall F22Inwf.rdxyim vdbInwf vdbInwfm .01 F22Isef.rdxyim vdbIsef 120379 280 0.000 0.001 0.000 F22Inef.rdxyim vdbInef 139653 288 0.000 0.001 0.000 F22Iswf.rdxyim vdbIswf 155890 455 0.000 0.001 0.000 F22Inwf.rdxyim vdbInwf 179677 373 0.000 0.001 0.000 grep no vdbIsefm > vdbIsefmno grep no vdbIswfm > vdbIswfmno grep no vdbInefm > vdbInefmno grep no vdbInwfm > vdbInwfmno rd2reg.pl vdbIsefmno rd2reg.pl vdbIswfmno rd2reg.pl vdbInefmno rd2reg.pl vdbInwfmno ** 22 January 2003 The cuts are actually at: RA : 334:31:40 +/- 5" Dec: 0:24:28 +/- 5" And of course, they're curved in RA,Dec... awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdb >vdbf awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIse >vdbIsef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIne >vdbInef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIsw >vdbIswf awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbInw >vdbInwf awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIce >vdbIcef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Ise.rdxyim >F22Isef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Ine.rdxyim >F22Inef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Isw.rdxyim >F22Iswf.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Inw.rdxyim >F22Inwf.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Ice.rdxyim >F22Icef.rdxyim awk '$1==334.506 || $1==334.794 || $2==0.156675 || $2==0.665195' vdb awk '$1==334.506 || $1==334.794 || $2==0.156675 || $2==0.665195' F22I*rdxyim no results: good ** 23 January 2003 awk '334.506 < $1 && $1 < 334.794 && 0.15667 < $2 && $2 < 0.66519' vdb> vdbIce awk '334.506 >=$1 || $1 >=334.794 || 0.15667 >=$2 || $2 >=0.66519' vdb> vdbInc awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbIce >vdbIcef awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' vdbInc >vdbIncf cat F22Isef.rdxyim F22Iswf.rdxyim F22Inef.rdxyim F22Inwf.rdxyim > F22Incf.rdxyim matchf22: from: 16103 F22Icef.rdxyim to: 16103 vdbIcef matches: 16103 cefound missed to: 0 celost missed from: 0 cerej Centring statistics: Average: 0.00035736 RMS: 0.00114148 68 %-tile: 0.00000000 from: 599355 F22Incf.rdxyim to: 599340 vdbIncf matches: 598920 ncfound missed to: 420 nclost missed from: 435 ncrej Centring statistics: Average: 0.00033884 RMS: 0.00111225 68 %-tile: 0.00000000 cat F22Bse.rdxyim F22Bsw.rdxyim F22Bne.rdxyim F22Bnw.rdxyim > F22Bnc.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Bce.rdxyim >F22Bcef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Bnc.rdxyim >F22Bncf.rdxyim from: 16103 F22Bcef.rdxyim to: 16103 vdbIcef matches: 16103 cefound missed to: 0 celost missed from: 0 cerej Centring statistics: Average: 0.00035736 RMS: 0.00114148 68 %-tile: 0.00000000 from: 599355 F22Bncf.rdxyim to: 599340 vdbIncf matches: 598920 ncfound missed to: 420 nclost missed from: 435 ncrej Centring statistics: Average: 0.00033884 RMS: 0.00111225 68 %-tile: 0.00000000 ** 24 January 2003 blkavg F22Bce.fits F22Bceb4.fits 4 4 blkavg F22Bse.fits F22Bseb4.fits 4 4 blkavg F22Bne.fits F22Bneb4.fits 4 4 blkavg F22Bnw.fits F22Bnwb4.fits 4 4 blkavg F22Bsw.fits F22Bswb4.fits 4 4 334.112500 nw 0:52:09 22:16:27 ne 0:51:00 22:21:52 se -0:03:31 22:21:52 sw -0:03:20 22:16:28 --> B-area: 22:16:28 < R.A. < 22:21:52 -0:03:20 < Dec. < 0:51:00 334.11666666 < R.A. < 335.4666666 -0.05555555 < Dec. < 0.85 F22Ilim.reg # Region file format: DS9 version 3.0 global color=blue font="helvetica 10 normal" select=1 edit=1 move=1 delete=1 include=1 fixed=0 source fk5;line(334.116666,-0.055555,335.466666,-0.055555) # line=0 0 width=4 fk5;line(335.466666,-0.055555,335.466666,00.850000) # line=0 0 width=4 fk5;line(335.466666,00.850000,334.116666,00.850000) # line=0 0 width=4 fk5;line(334.116666,00.850000,334.116666,-0.055555) # line=0 0 width=4 save the new boxes to fred.reg ... using emacs, edit from: # 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.71702,0.38007835,44.908",49.028",0) fk5;box(334.67547,0.21699777,49.028",49.852",0) fk5;box(335.41801,0.17051158,381.1",26.78",0) fk5;box(335.09087,-0.044496036,44.084",40.788",0) to: 334.71702 0.38007835 44.908 49.028 334.67547 0.21699777 49.028 49.852 335.41801 0.17051158 381.1 26.78 335.09087 -0.044496036 44.084 40.788 awk '{dx=$3/2/3600 ; dy =$4/2/3600 ; print $1-dx,$2-dy,$1+dx,$2+dy}' fred.reg > F22Badd.box cat F22I.box F22Badd.box >F22B.box box2reg.pl F22Badd.box red F22Badd.box.reg box2reg.pl F22B.box red F22B.box.reg ** 27 January 2003 Masking the bad astrometry regions: lasa28://home/gwyn/red/F22/F22B/swarpcheck/sub making blink.reg ds9 -geometry 1200x81000+0+0 -histequ *se021*.fits -single ds9 -geometry 1200x81000+0+0 -histequ *se107*.fits -single ds9 -geometry 1200x81000+0+0 -histequ *se108*.fits -single ds9 -geometry 1200x81000+0+0 -histequ *ne148*.fits -single ds9 -geometry 1200x81000+0+0 -histequ *ne104*.fits -single # Region file format: DS9 version 3.0 global color=green font="helvetica 10 normal" select=1 edit=1 move=1 delete=1 in clude=1 fixed=0 source fk5;box(335.17996,0.62625866,161.607",48.41",-0) fk5;box(335.45752,0.62248013,177.057",23.587",-0) fk5;box(334.70099,-0.052938054,196.318",119.068",-0) fk5;box(335.24791,0.24984967,28.432971",135.54799",-0) fk5;box(335.16531,0.63894848,42.539",61.079",-0) which (after clipping for the actual borders of the F22B) becomes. 335.158 0.619535 335.202 0.632982 335.433 0.619204 335.4666666 0.625756 334.674 -0.05555555 334.728 -0.0364008 335.244 0.231024 335.252 0.268676 335.159 0.630465 335.171 0.647432 add these lines to F22Badd.box and F22B.box awk '{a+=(($3-$1)*($4-$2))} ; END {print a}' 335.158 0.619535 335.202 0.632982 335.433 0.619204 335.4666666 0.625756 334.674 -0.05555555 334.728 -0.0364008 335.244 0.231024 335.252 0.268676 335.159 0.630465 335.171 0.647432 0.00235143 awk '{print(($3-$1)*($4-$2))}' 334.11666666 -0.05555555 335.4666666 0.85 rough area: 1.2225 total affected area: 8.5 sq arcmin. or a percentage of: 0.19% awk '{print $1,$2}' F22Bcef.rdxyim > F22Bcef.rd awk '{print $1,$2}' F22Bncf.rdxyim > F22Bncf.rd grep -v # F22Bce.cat > F22Bce.fat awk '$1!~/#/ && (334.506 >=$4 || $4 >=334.794 || 0.15667 >=$5 || $5 >=0.66519) ' F22Bse.cat F22Bsw.cat F22Bne.cat F22Bnw.cat > F22Bnc.fat wc *fat *f.rd 16103 1046695 11981680 F22Bce.fat 599355 38958075 456151111 F22Bnc.fat 16103 32206 322060 F22Bcef.rd 599355 1198710 12133452 F22Bncf.rd paste F22Bcef.rd F22Bce.fat > F22Bce.rdfat paste F22Bncf.rd F22Bnc.fat > F22Bnc.rdfat ** 28 January 2003 508437 508438 508439 the raw frames look like crap. based on the FLIPS info in the other headers and searches in mapix:/u3/gwyn/find files* ncftp -u gwyn ev6 cd /evsix4/jcc/CFH12K-F156/SkyFlat/B get *.fits cd /evsix5/dantel/VirmosMasterFiles/Nov99-01/Mask get *.fits (actually gzip them into /u3/gwyn/q first...) cd cd /evsix4/jcc/CFH12K-F156/Dark/1080 get *.fits HISTORY = imred: Input image: 508364o00.fits HISTORY = imred: Output image: 508364o00MDF.fits HISTORY = imred: Mask bad pixels: * Mask-CFH12K-99II.00.fits HISTORY = imred: Dark subtraction: - Dark18mn_121199.00.fits * 1.0 HISTORY = imred: Flat-field: * 7101.00 / SkyFlatB_1199.00.fits HISTORY = imred: Flat-field scale factor = max[mode flats] = 7101 from CCD 4 COMMENT = imred: F22B-nov99 @ IAP - Mireille so... operand1= 508439o11.fits Operand image or numerical constant op = * Operator operand2= Mask-CFH12K-99II.11.fits Operand image or numerical constant result = fred.fits Resultant image (title = ) Title for resultant image (divzero= 0.) Replacement value for division by zero (hparams= ) List of header parameters (pixtype= double) Pixel type for resultant image (calctyp= double) Calculation data type (verbose= yes) Print operations? (noact = no) Print operations without performing them? (mode = ql) imarith 508437o00.fits * Mask-CFH12K-99II.00.fits 508437o00M.fits imarith 508437o01.fits * Mask-CFH12K-99II.01.fits 508437o01M.fits imarith 508437o02.fits * Mask-CFH12K-99II.02.fits 508437o02M.fits imarith 508437o03.fits * Mask-CFH12K-99II.03.fits 508437o03M.fits imarith 508437o04.fits * Mask-CFH12K-99II.04.fits 508437o04M.fits imarith 508437o05.fits * Mask-CFH12K-99II.05.fits 508437o05M.fits imarith 508437o06.fits * Mask-CFH12K-99II.06.fits 508437o06M.fits imarith 508437o07.fits * Mask-CFH12K-99II.07.fits 508437o07M.fits imarith 508437o08.fits * Mask-CFH12K-99II.08.fits 508437o08M.fits imarith 508437o09.fits * Mask-CFH12K-99II.09.fits 508437o09M.fits imarith 508437o10.fits * Mask-CFH12K-99II.10.fits 508437o10M.fits imarith 508437o11.fits * Mask-CFH12K-99II.11.fits 508437o11M.fits imarith 508438o00.fits * Mask-CFH12K-99II.00.fits 508438o00M.fits imarith 508438o01.fits * Mask-CFH12K-99II.01.fits 508438o01M.fits imarith 508438o02.fits * Mask-CFH12K-99II.02.fits 508438o02M.fits imarith 508438o03.fits * Mask-CFH12K-99II.03.fits 508438o03M.fits imarith 508438o04.fits * Mask-CFH12K-99II.04.fits 508438o04M.fits imarith 508438o05.fits * Mask-CFH12K-99II.05.fits 508438o05M.fits imarith 508438o06.fits * Mask-CFH12K-99II.06.fits 508438o06M.fits imarith 508438o07.fits * Mask-CFH12K-99II.07.fits 508438o07M.fits imarith 508438o08.fits * Mask-CFH12K-99II.08.fits 508438o08M.fits imarith 508438o09.fits * Mask-CFH12K-99II.09.fits 508438o09M.fits imarith 508438o10.fits * Mask-CFH12K-99II.10.fits 508438o10M.fits imarith 508438o11.fits * Mask-CFH12K-99II.11.fits 508438o11M.fits imarith 508439o00.fits * Mask-CFH12K-99II.00.fits 508439o00M.fits imarith 508439o01.fits * Mask-CFH12K-99II.01.fits 508439o01M.fits imarith 508439o02.fits * Mask-CFH12K-99II.02.fits 508439o02M.fits imarith 508439o03.fits * Mask-CFH12K-99II.03.fits 508439o03M.fits imarith 508439o04.fits * Mask-CFH12K-99II.04.fits 508439o04M.fits imarith 508439o05.fits * Mask-CFH12K-99II.05.fits 508439o05M.fits imarith 508439o06.fits * Mask-CFH12K-99II.06.fits 508439o06M.fits imarith 508439o07.fits * Mask-CFH12K-99II.07.fits 508439o07M.fits imarith 508439o08.fits * Mask-CFH12K-99II.08.fits 508439o08M.fits imarith 508439o09.fits * Mask-CFH12K-99II.09.fits 508439o09M.fits imarith 508439o10.fits * Mask-CFH12K-99II.10.fits 508439o10M.fits imarith 508439o11.fits * Mask-CFH12K-99II.11.fits 508439o11M.fits NOTE WELL!! No overscan substraction!! NOTE WELL!! No bias substraction!! NOTE WELL!! The dark has the wrong length of time, in keeping with what people have done before!! The modes of the overscan sectors of the images and the darks more or less the same. imarith 508437o00M.fits - Dark18mn_121199.00.fits 508437o00MD.fits imarith 508437o01M.fits - Dark18mn_121199.01.fits 508437o01MD.fits imarith 508437o02M.fits - Dark18mn_121199.02.fits 508437o02MD.fits imarith 508437o03M.fits - Dark18mn_121199.03.fits 508437o03MD.fits imarith 508437o04M.fits - Dark18mn_121199.04.fits 508437o04MD.fits imarith 508437o05M.fits - Dark18mn_121199.05.fits 508437o05MD.fits imarith 508437o06M.fits - Dark18mn_121199.06.fits 508437o06MD.fits imarith 508437o07M.fits - Dark18mn_121199.07.fits 508437o07MD.fits imarith 508437o08M.fits - Dark18mn_121199.08.fits 508437o08MD.fits imarith 508437o09M.fits - Dark18mn_121199.09.fits 508437o09MD.fits imarith 508437o10M.fits - Dark18mn_121199.10.fits 508437o10MD.fits imarith 508437o11M.fits - Dark18mn_121199.11.fits 508437o11MD.fits imarith 508438o00M.fits - Dark18mn_121199.00.fits 508438o00MD.fits imarith 508438o01M.fits - Dark18mn_121199.01.fits 508438o01MD.fits imarith 508438o02M.fits - Dark18mn_121199.02.fits 508438o02MD.fits imarith 508438o03M.fits - Dark18mn_121199.03.fits 508438o03MD.fits imarith 508438o04M.fits - Dark18mn_121199.04.fits 508438o04MD.fits imarith 508438o05M.fits - Dark18mn_121199.05.fits 508438o05MD.fits imarith 508438o06M.fits - Dark18mn_121199.06.fits 508438o06MD.fits imarith 508438o07M.fits - Dark18mn_121199.07.fits 508438o07MD.fits imarith 508438o08M.fits - Dark18mn_121199.08.fits 508438o08MD.fits imarith 508438o09M.fits - Dark18mn_121199.09.fits 508438o09MD.fits imarith 508438o10M.fits - Dark18mn_121199.10.fits 508438o10MD.fits imarith 508438o11M.fits - Dark18mn_121199.11.fits 508438o11MD.fits imarith 508439o00M.fits - Dark18mn_121199.00.fits 508439o00MD.fits imarith 508439o01M.fits - Dark18mn_121199.01.fits 508439o01MD.fits imarith 508439o02M.fits - Dark18mn_121199.02.fits 508439o02MD.fits imarith 508439o03M.fits - Dark18mn_121199.03.fits 508439o03MD.fits imarith 508439o04M.fits - Dark18mn_121199.04.fits 508439o04MD.fits imarith 508439o05M.fits - Dark18mn_121199.05.fits 508439o05MD.fits imarith 508439o06M.fits - Dark18mn_121199.06.fits 508439o06MD.fits imarith 508439o07M.fits - Dark18mn_121199.07.fits 508439o07MD.fits imarith 508439o08M.fits - Dark18mn_121199.08.fits 508439o08MD.fits imarith 508439o09M.fits - Dark18mn_121199.09.fits 508439o09MD.fits imarith 508439o10M.fits - Dark18mn_121199.10.fits 508439o10MD.fits imarith 508439o11M.fits - Dark18mn_121199.11.fits 508439o11MD.fits # IMAGE NPIX MEAN MIDPT MODE STDDEV MIN MAX SkyFlatB_1199.00.fits 8586240 4386. 4481. 4596. 1157. 0. 65534. SkyFlatB_1199.01.fits 8586240 5179. 5339. 5344. 1078. 0. 65534. SkyFlatB_1199.02.fits 8586240 4268. 4347. 4347. 653.3 0. 16714. SkyFlatB_1199.03.fits 8586240 5411. 5498. 5495. 1000. 0. 65534. SkyFlatB_1199.04.fits 8586240 6948. 7082. 7081. 1127. 0. 65534. SkyFlatB_1199.05.fits 8586240 3601. 3660. 3723. 1120. 0. 65534. SkyFlatB_1199.06.fits 8586240 4410. 4531. 4610. 712. 0. 9504. SkyFlatB_1199.07.fits 8586240 5121. 5225. 5226. 784.9 0. 44155. SkyFlatB_1199.08.fits 8586240 6405. 6536. 6541. 1050. 0. 65534. SkyFlatB_1199.09.fits 8586240 5922. 6046. 6049. 897.3 0. 30894. SkyFlatB_1199.10.fits 8586240 4669. 4769. 4777. 709.2 0. 20447. SkyFlatB_1199.11.fits 8586240 4061. 4177. 4253. 660.5 0. 53026. imarith 508437o00MD.fits * 7101 508437o00MDF.fits imarith 508437o01MD.fits * 7101 508437o01MDF.fits imarith 508437o02MD.fits * 7101 508437o02MDF.fits imarith 508437o03MD.fits * 7101 508437o03MDF.fits imarith 508437o04MD.fits * 7101 508437o04MDF.fits imarith 508437o05MD.fits * 7101 508437o05MDF.fits imarith 508437o06MD.fits * 7101 508437o06MDF.fits imarith 508437o07MD.fits * 7101 508437o07MDF.fits imarith 508437o08MD.fits * 7101 508437o08MDF.fits imarith 508437o09MD.fits * 7101 508437o09MDF.fits imarith 508437o10MD.fits * 7101 508437o10MDF.fits imarith 508437o11MD.fits * 7101 508437o11MDF.fits imarith 508438o00MD.fits * 7101 508438o00MDF.fits imarith 508438o01MD.fits * 7101 508438o01MDF.fits imarith 508438o02MD.fits * 7101 508438o02MDF.fits imarith 508438o03MD.fits * 7101 508438o03MDF.fits imarith 508438o04MD.fits * 7101 508438o04MDF.fits imarith 508438o05MD.fits * 7101 508438o05MDF.fits imarith 508438o06MD.fits * 7101 508438o06MDF.fits imarith 508438o07MD.fits * 7101 508438o07MDF.fits imarith 508438o08MD.fits * 7101 508438o08MDF.fits imarith 508438o09MD.fits * 7101 508438o09MDF.fits imarith 508438o10MD.fits * 7101 508438o10MDF.fits imarith 508438o11MD.fits * 7101 508438o11MDF.fits imarith 508439o00MD.fits * 7101 508439o00MDF.fits imarith 508439o01MD.fits * 7101 508439o01MDF.fits imarith 508439o02MD.fits * 7101 508439o02MDF.fits imarith 508439o03MD.fits * 7101 508439o03MDF.fits imarith 508439o04MD.fits * 7101 508439o04MDF.fits imarith 508439o05MD.fits * 7101 508439o05MDF.fits imarith 508439o06MD.fits * 7101 508439o06MDF.fits imarith 508439o07MD.fits * 7101 508439o07MDF.fits imarith 508439o08MD.fits * 7101 508439o08MDF.fits imarith 508439o09MD.fits * 7101 508439o09MDF.fits imarith 508439o10MD.fits * 7101 508439o10MDF.fits imarith 508439o11MD.fits * 7101 508439o11MDF.fits imarith 508437o00MDF.fits / SkyFlatB_1199.00.fits 508437o00MDF.fits imarith 508437o01MDF.fits / SkyFlatB_1199.01.fits 508437o01MDF.fits imarith 508437o02MDF.fits / SkyFlatB_1199.02.fits 508437o02MDF.fits imarith 508437o03MDF.fits / SkyFlatB_1199.03.fits 508437o03MDF.fits imarith 508437o04MDF.fits / SkyFlatB_1199.04.fits 508437o04MDF.fits imarith 508437o05MDF.fits / SkyFlatB_1199.05.fits 508437o05MDF.fits imarith 508437o06MDF.fits / SkyFlatB_1199.06.fits 508437o06MDF.fits imarith 508437o07MDF.fits / SkyFlatB_1199.07.fits 508437o07MDF.fits imarith 508437o08MDF.fits / SkyFlatB_1199.08.fits 508437o08MDF.fits imarith 508437o09MDF.fits / SkyFlatB_1199.09.fits 508437o09MDF.fits imarith 508437o10MDF.fits / SkyFlatB_1199.10.fits 508437o10MDF.fits imarith 508437o11MDF.fits / SkyFlatB_1199.11.fits 508437o11MDF.fits imarith 508438o00MDF.fits / SkyFlatB_1199.00.fits 508438o00MDF.fits imarith 508438o01MDF.fits / SkyFlatB_1199.01.fits 508438o01MDF.fits imarith 508438o02MDF.fits / SkyFlatB_1199.02.fits 508438o02MDF.fits imarith 508438o03MDF.fits / SkyFlatB_1199.03.fits 508438o03MDF.fits imarith 508438o04MDF.fits / SkyFlatB_1199.04.fits 508438o04MDF.fits imarith 508438o05MDF.fits / SkyFlatB_1199.05.fits 508438o05MDF.fits imarith 508438o06MDF.fits / SkyFlatB_1199.06.fits 508438o06MDF.fits imarith 508438o07MDF.fits / SkyFlatB_1199.07.fits 508438o07MDF.fits imarith 508438o08MDF.fits / SkyFlatB_1199.08.fits 508438o08MDF.fits imarith 508438o09MDF.fits / SkyFlatB_1199.09.fits 508438o09MDF.fits imarith 508438o10MDF.fits / SkyFlatB_1199.10.fits 508438o10MDF.fits imarith 508438o11MDF.fits / SkyFlatB_1199.11.fits 508438o11MDF.fits imarith 508439o00MDF.fits / SkyFlatB_1199.00.fits 508439o00MDF.fits imarith 508439o01MDF.fits / SkyFlatB_1199.01.fits 508439o01MDF.fits imarith 508439o02MDF.fits / SkyFlatB_1199.02.fits 508439o02MDF.fits imarith 508439o03MDF.fits / SkyFlatB_1199.03.fits 508439o03MDF.fits imarith 508439o04MDF.fits / SkyFlatB_1199.04.fits 508439o04MDF.fits imarith 508439o05MDF.fits / SkyFlatB_1199.05.fits 508439o05MDF.fits imarith 508439o06MDF.fits / SkyFlatB_1199.06.fits 508439o06MDF.fits imarith 508439o07MDF.fits / SkyFlatB_1199.07.fits 508439o07MDF.fits imarith 508439o08MDF.fits / SkyFlatB_1199.08.fits 508439o08MDF.fits imarith 508439o09MDF.fits / SkyFlatB_1199.09.fits 508439o09MDF.fits imarith 508439o10MDF.fits / SkyFlatB_1199.10.fits 508439o10MDF.fits imarith 508439o11MDF.fits / SkyFlatB_1199.11.fits 508439o11MDF.fits rm *M.fits rm *MD.fits mkdir back cp *MDF.fits back PACKAGE = imutil TASK = imarith operand1= 508439o01MDF.fits Operand image or numerical constant op = / Operator operand2= SkyFlatB_1199.01.fits Operand image or numerical constant result = 508439o01MDF.fits Resultant image (title = ) Title for resultant image (divzero= 0.) Replacement value for division by zero (hparams= ) List of header parameters (pixtype= ushort) Pixel type for resultant image (calctyp= double) Calculation data type (verbose= yes) Print operations? (noact = no) Print operations without performing them? (mode = ql) imarith 508437o00MDF.fits * 1 508437o00MDF.fits imarith 508437o01MDF.fits * 1 508437o01MDF.fits imarith 508437o02MDF.fits * 1 508437o02MDF.fits imarith 508437o03MDF.fits * 1 508437o03MDF.fits imarith 508437o04MDF.fits * 1 508437o04MDF.fits imarith 508437o05MDF.fits * 1 508437o05MDF.fits imarith 508437o06MDF.fits * 1 508437o06MDF.fits imarith 508437o07MDF.fits * 1 508437o07MDF.fits imarith 508437o08MDF.fits * 1 508437o08MDF.fits imarith 508437o09MDF.fits * 1 508437o09MDF.fits imarith 508437o10MDF.fits * 1 508437o10MDF.fits imarith 508437o11MDF.fits * 1 508437o11MDF.fits imarith 508438o00MDF.fits * 1 508438o00MDF.fits imarith 508438o01MDF.fits * 1 508438o01MDF.fits imarith 508438o02MDF.fits * 1 508438o02MDF.fits imarith 508438o03MDF.fits * 1 508438o03MDF.fits imarith 508438o04MDF.fits * 1 508438o04MDF.fits imarith 508438o05MDF.fits * 1 508438o05MDF.fits imarith 508438o06MDF.fits * 1 508438o06MDF.fits imarith 508438o07MDF.fits * 1 508438o07MDF.fits imarith 508438o08MDF.fits * 1 508438o08MDF.fits imarith 508438o09MDF.fits * 1 508438o09MDF.fits imarith 508438o10MDF.fits * 1 508438o10MDF.fits imarith 508438o11MDF.fits * 1 508438o11MDF.fits imarith 508439o00MDF.fits * 1 508439o00MDF.fits imarith 508439o01MDF.fits * 1 508439o01MDF.fits imarith 508439o02MDF.fits * 1 508439o02MDF.fits imarith 508439o03MDF.fits * 1 508439o03MDF.fits imarith 508439o04MDF.fits * 1 508439o04MDF.fits imarith 508439o05MDF.fits * 1 508439o05MDF.fits imarith 508439o06MDF.fits * 1 508439o06MDF.fits imarith 508439o07MDF.fits * 1 508439o07MDF.fits imarith 508439o08MDF.fits * 1 508439o08MDF.fits imarith 508439o09MDF.fits * 1 508439o09MDF.fits imarith 508439o10MDF.fits * 1 508439o10MDF.fits imarith 508439o11MDF.fits * 1 508439o11MDF.fits ls *MDF.fits | perl -nle 'system "ds9 -geometry 1200x1000+0+0 -histequ $_"' 508438o02 bright streak on samplix cd /samplix5/gwyn mkdir Nov99B-reproc back on lasa28: ncftp -u gwyn samplix.iap.fr cd /samplix5/gwyn/Nov99B-reproc put *MDF.fits F22B.imlist now looks like: P1: 1.486663/1080/1.0/May99(01b) 2.486664/1080/1.0/May99(01b) 3.486665/1080/1.0/May99(01b) 4.486830/1080/0.9/May99(01b) 5.486831/1080/0.9/May99(01b) 6.486832/1080/1.1/May99(01b) 7.486833/1080/1.1/May99(01b) 8.557365/1440/0.75/Nov00(01) 9.557366/1440/0.75/Nov00(01) 10.508496/1440/1.0/Nov99(01) 11.508497/1440/1.4/Nov99(01) 12.508498/1440/1.2/Nov99(01) P2: 1.508437/1440/0.8/Nov99(02) 2.508438/1440/0.8/Nov99(02) 3.508439/1440/0.9/Nov99(02) 4.508364/1440/1.0/Nov99(02) 5.508365/1440/1.0/Nov99(02) P3: 1.592159/1380/0.97/Jun01(03) 2.592160/1380/0.91/Jun01(03) 3.592161/1380/0.94/Jun01(03) 4.592162/1380/0.93/Jun01(03) 5.592163/1380/0.98/Jun01(03) P4: 1.592464/1380/1.05/Jun01(04) 2.592465/1380/1.03/Jun01(04) 3.592466/1380/1.07/Jun01(04) 4.592467/1380/1.15/Jun01(04) 5.592468/1380/1.26/Jun01(04) ls -R / > files findim.pl F22B.imlist results in: /samplix1/gwyn/raw /samplix3/mellier/F22P1-B-Mai99/notcoadded /samplix3/mellier/F22P1-B-Nov99/notcoadded /samplix3/mellier/F22P2-B-Nov99/notcoadded /samplix5/dantel/process-nov2000/images-B/F22P2b-B-Nov00/notcoadded /samplix5/gwyn/Nov99B-reproc linkall.pl F22B.imlist F22B.dirlist /samplix1/gwyn/F22/F22B ncftp -u gwyn samplix.iap.fr cd /samplix1/gwyn/F22/F22B get * bye gzip -d *gz ls -1 5*nov99*fits > fredin ls -1 5*nov99*fits > fredout imstat @fredin >fredbefore imarith @fredin + 0 @fredout pixtype=ushort imstat @fredout >fredafter rm imlist* imlist.pl ln -s /raid/gwyn/masks/mask*.fits . maskall > maskall.log check maskall.log and imlistall.log for peculilarities.... . ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"' cd /raid/gwyn/F22/F22B imlist.pl ln -s /raid/gwyn/masks/sexbox/*box . Also, retreive the masked images and have a look on lasa28: cd /north/gwyn/q ncftp -u virmos pix1.iap.fr cd /raid/gwyn/F22/F22B/ get 50843*.fits bye ls 508437*nov990p02.fits | perl -nle 'system "ds9 -geometry 1200x1000+0+0 -histequ $_"' ls 508438*nov990p02.fits | perl -nle 'system "ds9 -geometry 1200x1000+0+0 -histequ $_"' ls 508439*nov990p02.fits | perl -nle 'system "ds9 -geometry 1200x1000+0+0 -histequ $_"' Meanwhile back in pix1:/raid/gwyn/F22/F22B sex.pl >sex.log cp /raid/gwyn/USNO/USNOF22 USNO 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 F22Ia.rdm 2 500 dofit3.log checkdo.pl dofit3.log mkdir toss mv *o06may* toss mv *o00may* toss (charge transfer and brickwall) rm imlist* imlist.pl dofit.pl F22Ia.rdm 2 500 dofit3.log checkdo.pl dofit3.log 486663o08may990p01.fits 238 2 226 230 0.324 508497o10nov990p01.fits 297 2 279 280 0.311 486664o08may990p01.fits 248 2 226 234 0.290 486832o08may990p01.fits 243 2 225 225 0.285 486665o08may990p01.fits 251 2 237 234 0.278 508496o10nov990p01.fits 300 2 285 284 0.264 fewest initial matches: 486663o11may990p01.fits 71 2 62 63 0.092 486833o00may990p01.fits 77 2 62 63 0.167 486832o11may990p01.fits 78 2 70 69 0.096 486833o10may990p01.fits 84 2 70 73 0.139 486831o11may990p01.fits 84 2 69 70 0.150 486663o10may990p01.fits 84 2 72 70 0.093 fewest final matches: 486833o00may990p01.fits 77 2 62 63 0.167 486663o11may990p01.fits 71 2 62 63 0.092 486664o11may990p01.fits 85 2 67 68 0.109 486832o11may990p01.fits 78 2 70 69 0.096 486663o10may990p01.fits 84 2 72 70 0.093 486831o11may990p01.fits 84 2 69 70 0.150 looks OK except for May99... rm imlist* imlist.pl dofit.pl F22Ia.rdm 2 500 dofit4.log checkdo.pl dofit4.log look at the data... reinstate the May99 except CCD00 and CCD06 rdmatch rdall send rdall.sum to Morgan:/home/gwyn/red/F22/F22B/astrocheck (mv previous version of this directory to astrocheckv1) Morgan:/home/gwyn/red/F22/F22B/astrocheck/sm.1 which produces F22Bastrores.ps looks acceptable on to photometry still working in pix1:/raid/gwyn/F22/F22B rdm.pl (takes rdall and the .rd files, and makes .rdm and rdmall files) Use the version in pix1:/home/virmos/gwyn/pg (which is the same as lasa28:/home/gwyn/pg except for the perl path) erd.pl magmatch erdmall ls *o02*fead | perl -nle '`getoutex $_`' tar -cvf phot.tar erdall.sum erdall *exp *fead in Morgan:/home/gwyn/red/F22/F22B/photocheck fp cd /raid/gwyn/F22/F22B get phot.tar bye tar -xvf phot.tar connph.pl ls *exp >expall cat *exp | sort -n | head -1 > exprel cat *exp | sort -n | tail -1 >> exprel cat *exp | sort -n +1 | head -1 >> exprel cat *exp | sort -n +1 | tail -1 >> exprel sm.5 things look nominal (except for the May99 outlines) (edit applyzp.pl) applyzp.pl ls *fead | perl -nle 's/.fead//; rename "$_.fead","$_.head"' tar -cvf head.tar *.head fp cd /raid/gwyn/F22/F22B put head.tar bye back to pix1:/raid/gwyn/F22/F22B tar -xvf head.tar on mapix cd /mapix1/gwyn/F22 mkdir swarpF22B cd swarpF22B mkdir /mapix4/gwyn/F22/swarpF22B mkdir /mapix4/gwyn/F22/swarpF22B/vmem mkdir /mapix4/gwyn/F22/swarpF22B/resamp ln -s /mapix4/gwyn/F22/swarpF22B/vmem ln -s /mapix4/gwyn/F22/swarpF22B/resamp mkdir stacks mkdir masks ncftp -u virmos pix1 cd /raid/gwyn/F22/F22B get *.fits get head.tar ## which comes from Morgan:/home/gwyn/red/F22/F22B/photocheck bye cd masks ncftp -u virmos pix1 cd /raid/gwyn/masks get mask*jun010*.fits get mask*may990*.fits get mask*nov990*.fits get mask*nov000*.fits bye cd .. ls *.fits >F22B.swarplist tar -xvf head.tar cp /mapix2/gwyn/F22/swarpF22I/F22I.matchswarp F22B.matchswarp imhead /mapix2/gwyn/F22/swarpF22I/stacks/F22Iholem2.fits > F22B.head doswarp.pl F22B.swarplist F22B.fits >&F22B.swarplog ** 31 January 2003 /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22B.fits >! F22Bsw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22B.fits >! F22Bse.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22B.fits >! F22Bnw.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22B.fits >! F22Bne.fits /u3/hjmcc/bin/toFITS2d -x 19569 -X 5034 -y 13110 -Y 8888 F22B.fits >! F22Bce.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22B.weight.fits >! F22Bsw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22B.weight.fits >! F22Bse.weight.fits /u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22B.weight.fits >! F22Bnw.weight.fits /u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22B.weight.fits >! F22Bne.weight.fits /u3/hjmcc/bin/toFITS2d -x 19569 -X 5034 -y 13110 -Y 8888 F22B.weight.fits >! F22Bce.weight.fits cp ~/config/header.param . cp ~/config/bigimage.sex . cp ~/config/default.nnw . cp ~/config/gauss_3.0_7x7.conv . ln -s /mapix2/gwyn/F22/swarpF22B/stacks/F22I*.fits . sex -c bigimage.sex F22Ine.fits,F22Bne.fits -WEIGHT_IMAGE F22Ine.weight.fits,F22Bne.weight.fits -CATALOG_NAME F22Bne.cat sex -c bigimage.sex F22Ise.fits,F22Bse.fits -WEIGHT_IMAGE F22Ise.weight.fits,F22Bse.weight.fits -CATALOG_NAME F22Bse.cat sex -c bigimage.sex F22Isw.fits,F22Bsw.fits -WEIGHT_IMAGE F22Isw.weight.fits,F22Bsw.weight.fits -CATALOG_NAME F22Bsw.cat sex -c bigimage.sex F22Inw.fits,F22Bnw.fits -WEIGHT_IMAGE F22Inw.weight.fits,F22Bnw.weight.fits -CATALOG_NAME F22Bnw.cat sex -c bigimage.sex F22Iholem2ce.fits,F22Bce.fits -WEIGHT_IMAGE F22Iholem2ce.weight.fits,F22Bce.weight.fits -CATALOG_NAME F22Bce.cat sex -c bigimage.sex F22Bse.fits -CATALOG_NAME F22Bose.cat -WEIGHT_IMAGE F22Bse.weight.fits sex -c bigimage.sex F22Bne.fits -CATALOG_NAME F22Bone.cat -WEIGHT_IMAGE F22Bne.weight.fits sex -c bigimage.sex F22Bnw.fits -CATALOG_NAME F22Bonw.cat -WEIGHT_IMAGE F22Bnw.weight.fits sex -c bigimage.sex F22Bsw.fits -CATALOG_NAME F22Bosw.cat -WEIGHT_IMAGE F22Bsw.weight.fits -rw-r--r-- 1 gwyn system 5138040960 Oct 24 04:39 F22Iholem2.fits -rw-r--r-- 1 gwyn system 5138040960 Oct 24 04:39 F22Iholem2.weight.fits -rw-r--r-- 1 gwyn system 178974720 Oct 24 14:46 F22Iholem2ce.fits -rw-r--r-- 1 gwyn system 178974720 Oct 24 14:49 F22Iholem2ce.weight.fits -rw-r--r-- 1 gwyn system 1172900160 Oct 24 12:14 F22Iholem2ne.fits -rw-r--r-- 1 gwyn system 1172900160 Oct 24 12:40 F22Iholem2ne.weight.fits -rw-r--r-- 1 gwyn system 1334885760 Oct 24 12:03 F22Iholem2nw.fits -rw-r--r-- 1 gwyn system 1334885760 Oct 24 12:34 F22Iholem2nw.weight.fits -rw-r--r-- 1 gwyn system 1230117120 Oct 24 11:51 F22Iholem2se.fits -rw-r--r-- 1 gwyn system 1230117120 Oct 24 12:21 F22Iholem2se.weight.fits -rw-r--r-- 1 gwyn system 1400008320 Oct 24 11:57 F22Iholem2sw.fits -rw-r--r-- 1 gwyn system 1400008320 Oct 24 12:28 F22Iholem2sw.weight.fits blkavg F22Bnw.fits F22Bnwb4.fits 4 4 blkavg F22Bce.fits F22Bceb4.fits 4 4 blkavg F22Bse.fits F22Bseb4.fits 4 4 blkavg F22Bne.fits F22Bneb4.fits 4 4 blkavg F22Bsw.fits F22Bswb4.fits 4 4 (ship to Morgan) exsex F22Ice.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Ice.rdxyim exsex F22Ise.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Ise.rdxyim exsex F22Ine.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Ine.rdxyim exsex F22Isw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Isw.rdxyim exsex F22Inw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Inw.rdxyim exsex F22Bce.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bce.rdxyim exsex F22Bse.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bse.rdxyim exsex F22Bne.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bne.rdxyim exsex F22Bsw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bsw.rdxyim exsex F22Bnw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bnw.rdxyim ** 3 Febuary 2003 exsex F22Bose.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bose.rdm exsex F22Bosw.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bosw.rdm exsex F22Bone.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bone.rdm exsex F22Bonw.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6==0 {print $1,$2,$3}' > F22Bonw.rdm match F22Bone.rdm F22Ia.rdm match F22Bonw.rdm F22Ia.rdm match F22Bose.rdm F22Ia.rdm match F22Bosw.rdm F22Ia.rdm matchex F22Bosw.rdm F22Ia.rdm matchsw1 1 2 matchex F22Bose.rdm F22Ia.rdm matchse1 1 2 matchex F22Bonw.rdm F22Ia.rdm matchnw1 1 2 matchex F22Bone.rdm F22Ia.rdm matchne1 1 2 cat match??1 >match1 matchex F22Bosw.rdm F22Ia.rdm matchsw2 2 4 matchex F22Bose.rdm F22Ia.rdm matchse2 2 4 matchex F22Bonw.rdm F22Ia.rdm matchnw2 2 4 matchex F22Bone.rdm F22Ia.rdm matchne2 2 4 cat match??2 >match2 working in:/home/gwyn/F22/F22B/swarpcheck mkdir sub mow.pl cd sub ls *.fits | perl -nle 'system "getout $_"' ls *out >outall One more shifty region: # 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(335.17996,0.62625866,161.607",48.41",-0) fk5;box(335.45752,0.62248013,177.057",23.587",-0) fk5;box(334.70099,-0.052938054,196.318",119.068",-0) fk5;box(335.24791,0.24984967,28.432971",135.54799",-0) fk5;box(335.16531,0.63894848,42.539",61.079",-0) fk5;box(334.62467,0.40302315,202.086",20.394",0) # color=red edit with awk 335.17996 0.62625866 161.607 48.41 335.45752 0.62248013 177.057 23.587 334.70099 -0.052938054 196.318 119.068 335.24791 0.24984967 28.432971 135.54799 335.16531 0.63894848 42.539 61.079 334.62467 0.40302315 202.086 20.394 awk '{dx=$3/2/3600 ; dy =$4/2/3600 ; print $1-dx,$2-dy,$1+dx,$2+dy}' > F22Badd.box edit for the limits and you get: and you get: 335.158 0.619535 335.202 0.632982 335.433 0.619204 335.4666666 0.625756 334.674 -0.0555555 334.728 -0.0364008 335.244 0.231024 335.252 0.268676 335.159 0.630465 335.171 0.647432 334.597 0.400191 334.653 0.405856 edit for the limits and you get: (the limits are: 22:16:28 < RA < 22:21:52 -00:03:20 F22B.box OK this is very disorganized here is F22Badd.box it includes: - bad imaging masking from version 1 - bad astrometry from version 1 - bad astrometry from version 2 334.711 0.373269 334.723 0.386888 334.669 0.210074 334.682 0.223922 335.365 0.166792 335.471 0.174231 335.085 -0.050161 335.097 -0.038831 335.075 -0.0449545 335.092 -0.0327662 335.065 -0.0382604 335.078 -0.025557 335.057 -0.0364298 335.07 -0.0191486 335.045 -0.0286483 335.061 -0.0127405 335.032 -0.0208669 335.048 -0.00587466 334.645 0.430893 334.657 0.432152 334.774 0.62338 334.784 0.633108 334.879 0.623602 334.883 0.627607 334.722 0.411436 334.726 0.415212 335.367 0.616673 335.386 0.623196 335.436 0.615056 335.483 0.621808 334.324 0.268899 334.345 0.272447 334.283 0.291673 334.292 0.293618 334.31 0.234222 334.328 0.238915 334.172 0.829829 334.179 0.833377 334.243 0.433868 334.25 0.440162 334.506 0.75225 334.506 0.759002 334.518 0.652 334.526 0.655776 334.375 0.696632 334.382 0.703384 334.37 0.67237 334.378 0.679122 334.373 0.654517 334.383 0.662414 334.392 0.642615 334.399 0.648452 334.514 0.684501 334.521 0.690338 334.519 0.708763 334.525 0.713684 334.291 0.685643 334.298 0.691937 334.307 0.66802 334.312 0.674314 334.401 0.676033 334.407 0.680267 335.158 0.619535 335.202 0.632982 335.433 0.619204 335.4666666 0.625756 334.674 -0.05555555 334.728 -0.0364008 335.244 0.231024 335.252 0.268676 335.159 0.630465 335.171 0.647432 334.597 0.400191 334.653 0.405856 ds9 -geometry 1200x800+0+0 -histequ F22Bceb4.fits -region F22Badd.reg -region F22I.box.reg -region F22Blim.reg ds9 -geometry 1200x800+0+0 -histequ F22Bneb4.fits -region fred.reg ds9 -geometry 1200x800+0+0 -histequ F22Bnwb4.fits -region fred.reg ds9 -geometry 1200x800+0+0 -histequ F22Bseb4.fits -region fred.reg ds9 -geometry 1200x800+0+0 -histequ F22Bswb4.fits -region fred.reg grep cyan fred.reg (and edit...) 334.70742 0.76649108 39.964 725.532 335.37139 0.73516334 39.964 848.308 335.24746 0.031199798 18.952 646.84 335.24982 0.19557986 20.188 240.196 335.24776 0.13081147 10.3 80.34 335.24995 0.3000571 14.42 241.02 335.25035 0.37037645 16.48 269.448 335.37248 0.21273797 38.316 291.284 335.3734 0.27624615 31.724 203.116 335.37318 0.32968484 35.02 192.404 335.37296 0.37316822 36.668 133.076 awk '{dx=$3/2/3600 ; dy =$4/2/3600 ; print $1-dx,$2-dy,$1+dx,$2+dy}' >fred.box 334.702 0.665723 334.713 0.867259 335.366 0.617343 335.377 0.852984 335.245 -0.0586391 335.25 0.121039 335.247 0.162219 335.253 0.22894 335.246 0.119653 335.249 0.14197 335.248 0.266582 335.252 0.333532 335.248 0.332953 335.253 0.4078 335.367 0.172282 335.378 0.253194 335.369 0.248036 335.378 0.304457 335.368 0.302962 335.378 0.356408 335.368 0.354685 335.378 0.391651 cat F22I.box F22Badd.box > F22BI.box ~/red/F22/F22B/swarpcheck/trim.pl removes all the boxes and parts of boxes outside of the F22B bounding box. It reads F22BI.box and write F22B.box. goto lasa28:/home/gwyn/red/F22/F22B/catcheck cp ../catcheckv1/* . ncftp -u gwyn mapix.iap.fr cd /mapix1/gwyn/F22/swarpF22B/stacks/ get F22B??.cat wc F22B*cat 16162 1047215 11985169 F22Bce.cat 142419 9253920 108057147 F22Bne.cat 181446 11790675 137897750 F22Bnw.cat 123103 7998380 93642166 F22Bse.cat 157749 10250370 121230159 F22Bsw.cat grep -v # F22Bce.cat > F22Bce.fat awk '$1!~/#/ && (334.506 >=$4 || $4 >=334.794 || 0.15667 >=$5 || $5 >=0.66519) ' F22Bse.cat F22Bsw.cat F22Bne.cat F22Bnw.cat > F22Bnc.fat exsex F22Bce.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bce.rdxyim exsex F22Bse.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bse.rdxyim exsex F22Bne.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bne.rdxyim exsex F22Bsw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bsw.rdxyim exsex F22Bnw.cat ALPHA_J2000 DELTA_J2000 X_IMAGE Y_IMAGE NUMBER MAG_AUTO > F22Bnw.rdxyim cat F22Bse.rdxyim F22Bsw.rdxyim F22Bne.rdxyim F22Bnw.rdxyim | awk '334.506 >=$1 || $1 >=334.794 || 0.15667 >=$2 || $2 >=0.66519' > F22Bnc.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Bce.rdxyim >F22Bcef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Bnc.rdxyim >F22Bncf.rdxyim awk '{print $1,$2}' F22Bcef.rdxyim > F22Bcef.rd awk '{print $1,$2}' F22Bncf.rdxyim > F22Bncf.rd paste F22Bcef.rd F22Bce.fat > F22Bce.rdfat paste F22Bncf.rd F22Bnc.fat > F22Bnc.rdfat wc *.fat *f.rd Actually, probably best to redo all that tomorrow... from: 16103 F22Bcef.rdxyim to: 16103 vdbIcef matches: 16103 cefound missed to: 0 celost missed from: 0 cerej Centring statistics: Average: 0.00035736 RMS: 0.00114148 68 %-tile: 0.00000000 from: 599355 F22Bncf.rdxyim to: 599340 vdbIncf matches: 598920 ncfound missed to: 420 nclost missed from: 435 ncrej Centring statistics: Average: 0.00033884 RMS: 0.00111225 68 %-tile: 0.00000000 ** 4 Febuary 2003 exsex F22Bose.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6<4 && $1>334.11666666 &&$1<335.4666666 && $2>-0.05555555 && $2<0.85 {print $1,$2,$3}' > F22Bose.rdm exsex F22Bosw.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6<4 && $1>334.11666666 &&$1<335.4666666 && $2>-0.05555555 && $2<0.85 {print $1,$2,$3}' > F22Bosw.rdm exsex F22Bone.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6<4 && $1>334.11666666 &&$1<335.4666666 && $2>-0.05555555 && $2<0.85 {print $1,$2,$3}' > F22Bone.rdm exsex F22Bonw.cat alpha_j2000 delta_j2000 mag_auto mag_aper magerr_auto flags | awk '$3> $4-1 && $5<.5 && $6<4 && $1>334.11666666 &&$1<335.4666666 && $2>-0.05555555 && $2<0.85 {print $1,$2,$3}' > F22Bonw.rdm cat F22Bo??.rdm > F22Bo.rdm Working out the area with /vcom/gwyn/mask/fitsmask.f ln -s /vcom/gwyn/F22B/F22Bne.fits . ln -s /vcom/gwyn/F22B/F22Bnw.fits . ln -s /vcom/gwyn/F22B/F22Bse.fits . ln -s /vcom/gwyn/F22B/F22Bsw.fits . F22Bnw.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 image dimensions: 20000 16686 array dimensions: 20000 16686 reading F22Bnw.fits done oneing array done reading F22B.box done read 1028 boxes 12764.675 -8098.755 36357.605 7729.123 doing left doing top writing F22BnwMask.fits done area (sq deg): 0.1706 F22Bsw.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 image dimensions: 20000 17500 array dimensions: 20000 17500 reading F22Bsw.fits done oneing array done reading F22B.box done read 1199 boxes 12764.675 9400.245 36357.605 25228.123 doing left doing bottom writing F22BswMask.fits done area (sq deg): 0.1774 F22Bse.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 image dimensions: 17573 17500 array dimensions: 17573 17500 reading F22Bse.fits done oneing array done reading F22B.box done read 534 boxes -7234.325 9400.245 16358.605 25228.123 doing right doing bottom writing F22BseMask.fits done area (sq deg): 0.3989 F22Bne.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 image dimensions: 17573 16686 array dimensions: 17573 16686 reading F22Bne.fits done oneing array done reading F22B.box done read 873 boxes -7234.325 -8098.755 16358.605 7729.123 doing right doing top writing F22BneMask.fits done area (sq deg): 0.3898 In summary: area (sq deg): 0.1706 area (sq deg): 0.1774 area (sq deg): 0.3989 area (sq deg): 0.3898 Total area (adopted) 1.1367 preliminary number counts: /home/gwyn/red/F22/F22B/catcheck/sm.7 awk '334.506 >=$1 || $1 >=334.794 || 0.15667 >=$2 || $2 >=0.66519 {$5+=100000; print $0}' F22Bsw.rdxyim > fredsw awk '334.506 >=$1 || $1 >=334.794 || 0.15667 >=$2 || $2 >=0.66519 {$5+=300000; print $0}' F22Bse.rdxyim > fredse awk '334.506 >=$1 || $1 >=334.794 || 0.15667 >=$2 || $2 >=0.66519 {$5+=500000; print $0}' F22Bnw.rdxyim > frednw awk '334.506 >=$1 || $1 >=334.794 || 0.15667 >=$2 || $2 >=0.66519 {$5+=700000; print $0}' F22Bne.rdxyim > fredne cat fredsw fredse frednw fredne >fred awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' F22Bce.rdxyim >F22Bcef.rdxyim awk '{printf "%10.6f %10.6f %10.3f %10.3f %12d %8.3f\n",$1,$2,$3,$4,$5,$6}' fred >F22Bncf.rdxyim grep -v # Feb04-2003_14_39_09.virmos >vdb awk '334.506 < $2 && $2 < 334.794 && 0.15667 < $3 && $3 < 0.66519' vdb> vdbIce awk '334.506 >=$2 || $2 >=334.794 || 0.15667 >=$3 || $3 >=0.66519' vdb> vdbInc 220067845 334.180778 .742550 13886.597 23349.127 10110 21.824 2.30500E+00 .000 220067847 333.444565 .741192 1019.738 23327.279 00110 20.419 4.44300E+00 .000 220067853 333.642723 .741214 4483.362 23326.937 10110 21.199 2.80200E+00 .000 220067857 334.048210 .742432 11569.933 23347.251 10110 22.348 3.02200E+00 .000 220067860 333.866536 .741670 8395.006 23334.293 11010 20.152 2.17900E+00 3.000 220067620 334.305610 .741455 16068.018 23329.889 10110 22.681 2.51000E+00 .000 ra dec x y ma hlr sex_flag db_flag ID awk '{printf "%10.6f %10.6f %10.3f %10.3f %8.3f %8.4f %4d %6d %12d\n",$2,$3,$4,$5,$7,$8,$9,$6,$1}' vdbIce >vdbIcef awk '{printf "%10.6f %10.6f %10.3f %10.3f %8.3f %8.4f %4d %6d %12d\n",$2,$3,$4,$5,$7,$8,$9,$6,$1}' vdbInc >vdbIncf making sense of the VDB flags: stars resolved/faint edge masked saturated 10110 ^^^^^ |||||__ saturated (1 = saturated) ||||___ masked (1=not masked) |||____ galaxy? (1=too fuzzy, too bright, too faint) ||_____ star? (1=is definitely a star) |______ edge (1=not on edge) matchf22.f from: 16103 F22Bcef.rdxyim to: 16103 vdbIcef matches: 16103 cefound missed to: 0 celost missed from: 0 cerej Centring statistics: Average: 0.00035736 RMS: 0.00114148 68 %-tile: 0.00000000 from: 599355 F22Bncf.rdxyim to: 599340 vdbIncf matches: 598920 ncfound missed to: 420 nclost missed from: 435 ncrej Centring statistics: Average: 0.00033884 RMS: 0.00111225 68 %-tile: 0.00000000 cat cefound celost ncfound nclost | sort -n +8 >allmatch cat cefound celost ncfound nclost | awk '{print $9}' >fred wc fred sort -nu fred |wc ## make files that have all the data and start with RA,Dec awk '{print $1,$2}' F22Bcef.rdxyim > F22Bcef.rd awk '{print $1,$2}' F22Bncf.rdxyim > F22Bncf.rd ## modify matchf22 so it takes the longer lines awk '{print length}' *.rdfat | sort -nu | tail filefrom='F22B'//orient//'f.rdfat' character*830 linefrom(nmax) character*830 crr(n),tempc,c matchf22 (with any luck, for the last time) reading FROM data pixels or arcseconds? 16103 arcseconds sorting FROM data matching TO data from: 16103 F22Bce.rdfat to: 16103 vdbIcef matches: 16103 cefound missed to: 0 celost missed from: 0 cerej Centring statistics: Average: 0.00035736 RMS: 0.00114148 68 %-tile: 0.00000000 reading FROM data pixels or arcseconds? 599355 arcseconds sorting FROM data matching TO data from: 599355 F22Bnc.rdfat to: 599340 vdbIncf matches: 598920 ncfound missed to: 420 nclost missed from: 435 ncrej Centring statistics: Average: 0.00033884 RMS: 0.00111225 68 %-tile: 0.00000000 wc ncfound nclost cefound 598920 45517920 521494319 ncfound 420 31920 92820 nclost 16103 1223828 13719804 cefound 615443 46773668 535306943 total wc vdb 615443 5538986 89239235 vdb cat ncfound nclost cefound >fred sort -T./ -n +8 -o F22B.allmatch fred (do this on vcom) ** 5 Febuary 2003 wc F22B.allmatch 615443 46773668 535306943 F22B.allmatch awk '{print length}' F22B.allmatch > fred sort -nu fred | tail 901 902 903 904 906 907 909 911 912 915 awk '{print NF}' allmatch | sort -nu 76 awk '{print length}' 334.306245 0.407949 16078.997 17501.514 25.751 1.7430 24 10100 220000009 86 334.243270 0.407949 14978.488 17501.542 25.872 1.6950 24 10100 220000008 Got that? 615443 lines 76 fields per line, up to 915 characters per line. the VDB ID number is (77:86) the original part of the line (108:950) the stuff we aren't taking from the VDB (i.e. everything but ID,X,Y,RA,DEC) is: (165:950) in 60 fields MADE a mistake somewhere.... grep -v # F22Bce.cat > F22Bce.fat awk '$1!~/#/ && (334.506 < $4 && $4 < 334.794 && 0.15667 < $5 && $5 < 0.66519)' F22Bce.cat > fred diff fred F22Bce.fat (just checking) awk '$1!~/#/ && (334.506 >=$4 || $4 >=334.794 || 0.15667 >=$5 || $5 >=0.66519) ' F22Bse.cat F22Bsw.cat F22Bne.cat F22Bnw.cat > F22Bnc.fat awk '{printf "%10.6f %10.6f\n",$4,$5}' F22Bce.fat >F22Bcef.rd awk '{printf "%10.6f %10.6f\n",$4,$5}' F22Bnc.fat >F22Bncf.rd wc *.fat *f.rd 16103 1046695 11980679 F22Bce.fat 599355 38958075 456995489 F22Bnc.fat 16103 32206 354266 F22Bcef.rd 599355 1198710 13185810 F22Bncf.rd paste -d' ' F22Bcef.rd F22Bce.fat > F22Bce.rdfat paste -d' ' F22Bncf.rd F22Bnc.fat > F22Bnc.rdfat awk '{print length}' *.rdfat | sort -nu | tail cp -f *.rdfat /vcom/gwyn/sort goto virmos:/vcom/gwyn/sort matchf22 reading FROM data pixels or arcseconds? 16103 arcseconds sorting FROM data matching TO data from: 16103 F22Bce.rdfat to: 16103 vdbIcef matches: 16103 cefound missed to: 0 celost missed from: 0 cerej Centring statistics: Average: 0.00035736 RMS: 0.00114148 68 %-tile: 0.00000000 reading FROM data pixels or arcseconds? 599355 arcseconds sorting FROM data matching TO data cat ncfound nclost cefound >fred sort -T./ -n +8 -o F22B.allmatch fred from: 599355 F22Bnc.rdfat to: 599340 vdbIncf matches: 598920 ncfound missed to: 420 nclost missed from: 435 ncrej Centring statistics: Average: 0.00033884 RMS: 0.00111225 68 %-tile: 0.00000000 cat ncfound nclost cefound >fred sort -T./ -n +8 -o F22B.allmatch fred wc F22B.allmatch 615443 46773668 536392142 F22B.allmatch wk '{print length}' F22B.allmatch > fred sort -nu fred | tail 900 901 902 904 905 908 911 913 914 917 awk '{print NF}' allmatch | sort -nu 76 ** 6 Febuary 2002 awk '{print $15,$16,$59,$5}' F22B.allmatch > F22.rdbi awk '$10!=0 {print $1-333,$2,$10-333,$11,$15-333,$16}' F22B.allmatch >F22.rdrdrd (the 333 is to deal with the round off error in sm) Checking B-I over the field lasa28:/home/gwyn/red/F22/F22B/catcheck/sm.9 F22BIv.gif median B-I (suitable clipped) as a function of sky poistion lasa28:/home/gwyn/red/F22/F22B/catcheck/sm.d F22BIvhist.gif Histograms for the same 25 RA/Dec bins.. dashed line is median solid line is 1.6 Is there a problem?? Can't really say. in lasa28:/home/gwyn/red/F22/F22B/swarpcheck exsex F22Bone.cat alpha_j2000 delta_j2000 mag_auto magerr_auto mag_aper flags >F22Bone.rdmf exsex F22Bonw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto mag_aper flags >F22Bonw.rdmf exsex F22Bose.cat alpha_j2000 delta_j2000 mag_auto magerr_auto mag_aper flags >F22Bose.rdmf exsex F22Bosw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto mag_aper flags >F22Bosw.rdmf awk '$6==0 && $4<.5 && $3>$5-1' *.rdmf >F22Bo.rdm cd ../swarpcheck get all the .head files (with the zero points) get all the .rd files (no correction for zero points) lasa28:/home/gwyn/pg/retrophot.pl lasa28:/home/gwyn/red/F22/F22B/swarpcheck/sm.6 produces F22Bswarprawerr.gif lasa28:/home/gwyn/red/F22/F22B/swarpcheck/sm.6 produces F22Bswarpraw.gif (we are plotting differences in mag_auto) There are some offsets between individual images and the swarped image. This offset is mostly constant (within the errorbars) at about 0.05 The offset is correlates with chip number in some epochs. (things only get worse if you use an average instead of median BTW) ** 7 Febuary 2003 no way to avoid star-galaxy separation... lasa28:/home/gwyn/F22/F22B/catcheck/sm.g does it using Besancon counts from http://www.obs-besancon.fr/www/modele/ exsex F22Bone.cat alpha_j2000 delta_j2000 mag_auto 'flux_radius(2)' flags >F22Bone.rdmh exsex F22Bonw.cat alpha_j2000 delta_j2000 mag_auto 'flux_radius(2)' flags >F22Bonw.rdmh exsex F22Bose.cat alpha_j2000 delta_j2000 mag_auto 'flux_radius(2)' flags >F22Bose.rdmh exsex F22Bosw.cat alpha_j2000 delta_j2000 mag_auto 'flux_radius(2)' flags >F22Bosw.rdmh awk '$1>334.11666666 && $1<335.4666666 && $2>-0.05555555 && $2<0.85 && $3<25' *.rdmh >F22Bo.rdmh ship to Morgan Morgan:/home/gwyn/mask/seeingB.f plus Morgan:/home/gwyn/red/F22/F22B/catcheck/sm.j produce F22Bseeing.gif we now have F22B.zones... flagF22Bo.f works like flag.f for the F22B. It pretends there is no I image.. it produces F22Bo.mag and nothing else really, although it could be kludged... F22Bo.mag has one column with all the magnitude of all the galaxies+we're not sures detected in B in the F22.