2000.11
  重新制订了PIP2流程, 使更多的星用PSF求出,而不是用j_phot求出.
  生成一个星表文件, 例: T359i_15.mag  (大写字母打头), 放在同一地方, 星表很大.
包括了制作过程中的该删但没删的东西.
  星表格式同于以前, 但在第二行"测光误差"的前面加了一个序号, 使用者不必理会, 它
是一开始,统一找到星的编号.
  具体做法:
  1. 对各颜色做j_find,j_phot 得到各颜色的待求星. 转成赤经赤纬
  2. 合并各颜色的待求星, 凡出现3次, 即被选中. 做成一个统一的待求星表.
  3. 待求星表转成各颜色的XY坐标, 同时生成"边缘星编号", (靠边20pixel 为边缘星)
  4. 生成统一的"边缘星编号", 待后面总星表形成后, 在总星表上贴上"e"
  5. 用pip21_8处理各颜色, 免去j_find步骤, 从j_phot开始做.
  6. 根据各颜色的结果*.als的赤经赤纬, 合并成"总星表", 原则也是出现3次.
  7. 转成定标星等, 贴上"e"
  8. 用mo4目视检查所有星象, 贴上"bdgmnst"
  9. 用不贴标记的星, 对各颜色分别做一次j_phot, 求得5象元的孔径测光星等及误差,
     和已知的总星表上的对应星等做较差统计, 得"中值", 修改孔径测光星等, 贴在
     总星表上缺数据的地方. 其测光误差前冠以"X", 以此区分何种测光结果.
     虽然补充了不少数据, 但仍有空缺, 原因是孔径测光也测不到结果.
说明:
在2角秒内找到星, 认为是同一天体.
做了pip2后, 星的位置有变化, 因此合并总星表时, 不用早期的星编号.
尤其在一些"不深"的颜色上, 位置变了许多, 因此该颜色星等被剔除.
使用者检查星表图像程序 mo4   (新老计算机版本都有)
本计算过程, 全部在新机器上完成.
pipeline 2 的程序未作改动.

---------------------------------
2000.7
t329j_16.* 中标记说明

  b  亮星(饱和), 坐标不准
  d  受污染(鬼像,CCD缺陷等)
  e  靠边, 星像出框架
  f  星太暗, 但可见
  g  星系
  m  混淆 (邻近星的亮度差不多)
  n  星看不见
  s  受亮星污染(在晕中,或spark)
  t  受星系污染
  x  错误的星像, (亮星系引起)
  X  在PIP2中,在本颜色上没找到此星, 利用已知赤经赤纬, 反推出X,Y座标
     人为生成一个*.coo文件, 用j_phot孔径测光, 注意, 有的星误差很大.
     在星表上, 仍有一部分星没有星等, (j_phot测不出来), 属于两种情况,
     A. 太暗;  B. 出框架, (合并图像的框架<2048)
--------------
t329j.mag       所有星, 没修订 (星等均已定标, 除非有说明readme, 见T534)
t329j_16.mag    有8个以上颜色,
t329j_16.mag1   有8个以上颜色, 目视检查过. 做了以上类型标记.
t329j_16.mag2   在t329j_16.mag基础上, 反推X,Y座标, 加上孔径测光, 剔除"坏"者.

检查过, 在"基准"颜色中没有的星, 而在其它颜色中出现了2次以上的天体, 均为"坏"者.

文件名中字母含义:
  j   以该颜色框架坐标为基准, (挑星多, 合并后图像,质优者, 有时挑 i )
  16  共有多少颜色.
  建议使用 t329j_16.mag2

--------------
星表读出的FORTRAN例子:
        character*133 d0,d1,d2
        real mag(16),err(16)
        character*1 mark(16)
        open(1,file='/base/photo/T329/t329j_16.mag2',status='old')
        read(1,'(a)')d0
10      read(1,'(a)',end=19)d1
        read(1,'(a)')d2
        do i=1,16
        read(d1(i*6+33:),'(i5)')k
        mag(i)=k*0.001
        read(d2(i*6+33:),'(a1,i4)'),mark(i),k
        err(i)=k*0.001
        enddo
        goto 10
19      close(1)
        end