2000.11,  2002.5
  重新制订了PIP3流程, 使更多的星用PSF求出. (pip3.step 描述做星表的流程)
  最终结果:
  例: T517
  在目录 /base/photo/T517下,  (在vegao下有相同的版本)
  T517i_15.mag  星表, 整数, 老格式. 千分之一星等. 有X,Y.
  T517i_15.mag2 星表, 实数, 百分之一星等, 有序号,USNOSA20对应星(2角秒)
  T517i_15.mag3 星表, 同于*.mag2, 增加了孔径测光值,
  con           颜色定标文件.
  cob           有各个颜色的合并图像的制作过程说明.
  cob_fit       存放主要的一幅合并图像, 在此例中为i颜色,
文件名中字母含义:
  i   以该颜色框架坐标为基准, (挑星多, 合并后图像,质优者, 有时挑 j )
  15  共有多少颜色.

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

---------------------------------
T517i_15.* 中标记说明

  b  亮星(饱和), 坐标不准
  d  受污染(鬼像,CCD缺陷等)
  e  靠边, 星像出框架
  f  星太暗, 但可见
  g  星系
  m  混淆 (邻近星的亮度差不多)
  n  星看不见
  s  受亮星污染(在晕中,或spark)
  t  受星系污染
  x  错误的星像, (亮星系引起)
  X  在PIP2中,在本颜色上没找到此星, 利用已知赤经赤纬, 反推出X,Y座标
     人为生成一个*.coo文件, 用j_phot孔径测光, 注意, 有的星误差很大.
     在星表上, 仍有一部分星没有星等, (j_phot测不出来), 属于两种情况,
     A. 太暗;  B. 出框架, (合并图像的框架<2048)

--------------
星表读出的FORTRAN例子1:  (整型)
        character*133 d0,d1,d2
        real mag(16),err(16)
        character*1 mark(16)
        open(1,file='/base/photo/T517/T517i_15.mag',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
----------------------------
星表读出的FORTRAN例子2:  (整型)
        character*126 d0,d1,d2
        real mag(16),err(16)
        character*1 mark(16)
        open(1,file='/base/photo/T517/T517i_15.mag2',status='old')
        read(1,'(a)')d0
        nc=
10      read(1,'(a)',end=19)d1
        read(1,'(a)')d2
        do i=1,nc
        read(d1(i*6+26:),'(f5.2)')mag(i)
        read(d2(i*6+26:),'(a1,f4.2)'),mark(i),err(i)
        enddo
        goto 10
19      close(1)
        end
----------------------------
作颜色星等"较差"图程序:  在/vega2/rhbin下:
1. dim     处理*.mag
2. dim2    处理*.mag2


----------------------
             深浅暴光图像之间的星等常数确定
                                        2000,4
对于合并的DEEP_BATC天区图像, 请用PIP21处理.
(分9块做PSF, 最后得*.als--BATC星等单色星表, *.lst--全部PSF星, deep_psf.lst)

发表星等(i)= BATC星等(i) + CONST
    星等误差(i)= F ( BATC星等误差(i) , CONST 误差)

下面叙述如何得到 CONST 及 CONST 误差

    CONST = CONST1 + CONST2
    CONST1: 用测光夜的测光标准星 定出当天 浅暴光_BATC天区的值, (周旭做)
CONST2:
1. 用j_find, j_phot 处理浅暴光_BATC天区图像, 得到*.ap
2. 用pip6处理*.ap, 分三步:
  A. 类似j_auto自动找PSF星的方法, 找到一批"孤立星", 星等范围5个星等.
     按照cob50的做法合并. 生成的图像是,中央星的信噪比极高, 背景极平的图像.
     按DAO_PHOT的做法, 对该星做不同半径的孔径测光(3--20 PIXEL)
     求出5pixel与15pixel的流量比.
     取这批"孤立星" 5pixel的星等值, 根据流量比, 归算到15pixel 的星等
  B. 深暴光的PSF星的星等, 可以在*.als中找到.
  C. 取两个框架上坐标对应的30颗亮星, 去3sigma点, 求出平均差值及误差

注: 在/base/photo 的星表中列出的星等误差, 仅为pip2给出的误差
    在已发表的文章中, 有叙述, 此误差为主要误差, 其它误差可忽略.
    定标误差在文件 av.dat中.