### -- extract RA,DEC,flux,redshift from IIFSCz catalog
gawk '{print $2,$3,$5,$6,$7,$8,$37}' ~/data/IIFSCz/IIFSCz.v2.dat > IIFSCz.cat


### -- make IIFSCz mask
make_geometry make_geometry.param


### -- copy Akari data 
#    - objects detected with all bands
cp ~/data/Akari_FIS_BSC/public1.0/AKARI-FIS_BSC_V1.txt2 .
#    - all objects
cp ~/data/Akari_FIS_BSC/public1.0/AKARI-FIS_BSC_V1.txt .



### -- Make flux limited sample with S60
#    - make positive data to avoid divergence when logalithmic binning
gawk '$4>0 {print}' AKARI-FIS_BSC_V1.txt > temp

#    - monitor the logN-logS relation
count temp 4 50 log > S60.cnt2
rm temp -f

#    - cut flux at S60>0.6 Jy
gawk '$4>= 0.6 {print}' AKARI-FIS_BSC_V1.txt > AKARI-FIS_BSC_V1.S60flim.txt

#    - make flux_completeness map to find the quadrupole pattern
flux_completeness  flux_completeness.param



### -- apply IIFSCz mask to Akari catalog
apply_mask apply_mask.param




### -- matching Akari vs IIFSCz
matching AKARI-FIS_BSC_V1.S60flim.IIFSCzMask.txt 1 2 IIFSCz.cat 1 2 30 v | tee matching_Akari_IIFSCz_30sec.log

#    - extract catalog which is matched
match2ext AKARI-FIS_BSC_V1.S60flim.IIFSCzMask.txt IIFSCz.cat matching_Akari_IIFSCz_30sec.log 
mv -f fort.100 matching_Akari_IIFSCz_30sec.Akari.cat
mv -f fort.101 matching_Akari_IIFSCz_30sec.IIFSCz.cat

#    - extract catalog which is not matched
match2ext AKARI-FIS_BSC_V1.S60flim.IIFSCzMask.txt IIFSCz.cat matching_Akari_IIFSCz_30sec.log inv
mv -f fort.100 matching_Akari_IIFSCz_30sec.invAkari.cat
mv -f fort.101 matching_Akari_IIFSCz_30sec.invIIFSCz.cat



### -- matching with SDSS for objects not mached with IIFSCz
#    - copy SDSS DR7 spec-z sample RA.DEC.redshift and removing NF!=3 line
sed -e 's/,/ /g' ~/data/SDSS/DR7/result.csv | gawk '{print $2,$3,$15}' | gawk 'NF==3 {print}'> sdss.radec.z.txt

#    - matching within 30 arcsec
matching matching_Akari_IIFSCz_30sec.invAkari.cat 1 2 sdss.radec.z.txt 1 2 30 v | tee matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.log

#    - extract catalog which is matched
match2ext matching_Akari_IIFSCz_30sec.invAkari.cat sdss.radec.z.txt matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.log
mv -f fort.100 matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.Akari.cat
mv -f fort.101 matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.sdss.cat

#    - extract catalog which is not matched
match2ext matching_Akari_IIFSCz_30sec.invAkari.cat sdss.radec.z.txt matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.log inv
mv -f fort.100 matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.invAkari.cat
mv -f fort.101 matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.invsdss.cat



### -- find unmatched objects in the SIMBAD
URL:     http://simbad.u-strasbg.fr/simbad/sim-fcoo
click the output options and check the list display column as
   output format     : Ascii
   Max. # of objects : maximum 200000
   Identifier        : off
   Object type       : on/ standard description
   Coordinates1      : on/FK5 epoch2000
   Coordinates2-4    : off
   Proper motions    : off
   Parallax          : off
   Radial velocity   : on/ check only redshift
   Flux/Magnitudes   : on/ check J,H,K
   Spectral type     : on
   Morphological type: on
   Angular size      : on
   Hierarchy links   : off
   Bibliography      : off
   Notes             : on
   Measurements      : off

And click "save" button. Then jump to "Coordinate query" and upload file named 

  matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.invAkari.SIMBAD.entry

coordinate system is FK5, serch radius is 30-40 arcsec.
output on the web browser screen is saved as

  matching_Akari_IIFSCz_30sec.invAkari.matching_Akari_sdss_30sec.invAkari.SIMBAD.result

The result file seems somewhat awkward and can be well arranged by SIMBAD.arrange.f90.

  SIMBAD.arrange matching_AKA...SIMBAD.result > output
  


# -- star/galaxy separation for object not mathed with IIFSCz