pro astDiam, id



  triad = pds_TRIAD_getCatalog(id)
  iras  = pds_IRAS_getCatalog(id)
  msx   = pds_MSX_getCatalog(id)
  occ   = pds_occ_getCatalog(id,/DIAM)
  rw    = irasRW_getCatalog(id)
  akari = Akari_getCatalog(id)
  wise  = WISE_getCatalog(id)

;  help, iras, msx, rw, akari, wise
print, ''


  dimT = size(triad)
  if dimT(dimT(0)+1) eq 8 then begin
    print, ' '+string(triad.num, format='(I6)')+' '+$
           string(triad.des, format='(A-17)')+' '+$
           string(triad.d.val, format='(F7.2)')+'   '+$
           string(triad.d.unc, format='(F7.2)')+$
           '     STM      1    PDSSBN-TRIAD'
  endif




  dimI = size(iras)
  if dimI(dimI(0)+1) eq 8 then begin
    print, ' '+string(iras.num, format='(I6)')+' '+$
           string(iras.des, format='(A-17)')+' '+$
           string(iras.d.val, format='(F7.2)')+'   '+$
           string(iras.d.unc, format='(F7.2)')+$
           '     STM      1    PDSSBN-IRAS'
  endif




  dimM = size(msx)
  if dimM(dimM(0)+1) eq 8 then begin
    print, ' '+string(msx.num, format='(I6)')+' '+$
           string(msx.des, format='(A-17)')+' '+$
           string(msx.d.val, format='(F7.2)')+'   '+$
           string(msx.d.unc, format='(F7.2)')+$
           '     STM      1    PDSSBN-MSX'
  endif


  dimO = size(occ)
  if dimO(dimO(0)+1) eq 8 then begin
    circ= where(occ.b.val eq -9999.9)
    if circ(0) ne -1 then occ.b(circ).val=occ.a(circ).val
    diam = mean(sqrt( occ.a.val^2 + occ.b.val^2 ))
    v=where( occ.a.unc ne -99.9 and occ.b.unc ne -99.9)
    if v(0) ne -1 then $
      unc = diam*mean( sqrt( (occ(v).a.unc/occ(v).a.val)^2 + (occ(v).b.unc/occ(v).b.val)^2 ) ) $
    else $
      unc = 0.10*diam
    print, ' '+string(occ(0).num, format='(I6)')+' '+$
           string(occ(0).des, format='(A-17)')+' '+$
           string(diam, format='(F7.2)')+'   '+$
           string(unc, format='(F7.2)')+$
           '     OCC      1    PDSSBN-OCC'
  endif




  dimR = size(rw)
  if dimR(dimR(0)+1) eq 8 then begin
    nRyan = dimR(1)
    for k=0, nRyan-1 do begin

      if rw(k).stm.d.val ne -1000. then begin
        print, ' '+string(rw(k).num, format='(I6)')+' '+$
               string(rw(k).des, format='(A-17)')+' '+$
               string(rw(k).stm.d.val, format='(F7.2)')+'   '+$
               string(rw(k).stm.d.unc, format='(F7.2)')+$
               '     STM      1    2010-AJ-140-Ryan'
      endif
      if rw(k).neatm.d.val ne -1000. then begin
        print, ' '+string(rw(k).num, format='(I6)')+' '+$
               string(rw(k).des, format='(A-17)')+' '+$
               string(rw(k).neatm.d.val, format='(F7.2)')+'   '+$
               string(rw(k).neatm.d.unc, format='(F7.2)')+$
               '     NEATM    1    2010-AJ-140-Ryan'
      endif
    endfor
  endif




  dimA = size(akari)
  if dimA(dimA(0)+1) eq 8 then begin
    print, ' '+string(akari.num, format='(I6)')+' '+$
           string(akari.des, format='(A-17)')+' '+$
           string(akari.d.val, format='(F7.2)')+'   '+$
           string(akari.d.unc, format='(F7.2)')+$
           '     STM      1    2011-PASJ-63-Usui'
  endif




  dimW = size(wise)
  if dimW(dimW(0)+1) eq 8 then begin
    nWise = dimW(1)
    for k=0, nWise-1 do begin
      split = strsplit(wise(k).src,'.',/extract)
      print, ' '+string(wise(k).num, format='(I6)')+' '+$
             string(wise(k).des, format='(A-17)')+' '+$
             string(wise(k).d.val, format='(F7.2)')+'   '+$
             string(wise(k).d.unc, format='(F7.2)')+$
             '     STM      1    '+split(0)
    endfor
  endif







end